betabinexch0=function(theta , data)
{
al=theta[,1]
be=theta[,2]
y=data[,1]
n=data[,2]
N=length(y)
val=0*be;
for (i in 1:N){
val=val+lbeta(al+y[i],be+n[i]-y[i])
}
val=val-sum(lbeta(al,be))
val=val-5/2*log(al+be)
return(val)
}
mycontour(betabinexch0,c(100,200,100,200),data)


雷达卡




京公网安备 11010802022788号







