R/logMCMC.post.R

.logMCMC.post <- function (x,n,theta){

 if(length(x)!=length(n) | length(x)!=length(theta)){
	print("data lengths do not match!")
	stop()
 }
 sum(dbinom(x,n,prob=theta,log=T)/length(x))

}

Try the IFP package in your browser

Any scripts or data that you put into this service are public.

IFP documentation built on Nov. 26, 2020, 1:07 a.m.