winbugs中的一个错误

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

model
{
for(i in 1:M)
{
for(j in 1:T)
{
y[i,j]~dbin(p[i,j],50)
logit(p[i,j])<-inprod(x[i,],beta[])+lambda[j]*kethe[i]
}
kethe[i]<-kethe[i-1]*phi+delta[i]
}
for(i in 1:M)   {    x[i,1:T]~dmnorm(mu[], Omega[ , ])    }
for(j in 1:T)   {    mu[j]<-0.0     }
Omega[1:T,1:T]~ dwish(R[ , ], 3)
Sigma[1:T,1:T] <- inverse(Omega[ , ])
lambda[1:T]~dmnorm(mu[], Omega[ , ])
for(j in 1:T)   {    mu[j]<-0.0     }
Omega[1:T,1:T]~ dwish(R[ , ], 3)
Sigma[1:T,1:T] <- inverse(Omega[ , ])

for(i in 1:M)   {   
delta[i]~dnorm(mean[i],prec[i])
mean[i]<-mu[i]*(1-2*p)/(p*(1-p))
prec[i]<-2*mu[i]/(tau*p*(1-p))
mu[i]~dexp(tau)     }
phi~dunif(-1,1)
tau~dgamma(alpha0,gamma0)
kethe[0]~dnorm(0,1.0E-4)
}




winbugs始终在mean[i]<-mu[i]*(1-2*p)/(p*(1-p))
prec[i]<-2*mu[i]/(tau*p*(1-p))
报错，哪位好心人帮我找找错~？
万分感激啊~

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS bug Win winbugs

p change to P

for(i in 1:M)   {   
delta[i]~dnorm(mean[i],prec[i])
mean[i]<-mu[i]*(1-(2*P))/(P*(1-P))
prec[i]<-2*mu[i]/(tau*P*(1-P))
mu[i]~dexp(tau)}

2# epoh

万分感谢~
现在没有错误了~
但是为什么换成大写字母就可以了呢？~