z**********i 发帖数: 12276 | 1 最近用zero inflated beta binomial来做我的DATA.
但,不确定我的Likelihood function对不对.
r 是numerator, n 是denominator.因为,r有很多0,所以想用ZIBB来做.
if r=0 then ll=log(p_0+(1-p_0)*beta(A,B+n)/beta(A,B));
else ll = log(1-p_0)+(lgamma(n + 1)-lgamma(r + 1)-lgamma(n-r+1))+log(beta(A+
r,B+n-r))-log(beta(A,B));
多谢指教! |
|