悬赏 10 个论坛币 未解决
model{
d<-100000
for( k in 1 : N ) {
zeros[k]<-0
phi[k]<- d-log(L[k])
zeros[k]~dpois(phi[k])
L[k]<- beta*theta/pow(1+beta*x[k],theta+1)
}
theta~dgamma(0.5,5.0E-04)
beta~dgamma(0.5,5.0E-04)
}
list(N=30,x=c(0.3702472,0.3947406,0.4149750,0.4706422,0.4805167,0.4825024
,0.4847931,0.4872666,0.4952600,0.5018638,0.5189860,0.5247481,
0.5341415,0.5368265,0.5380381,0.5435172,0.5580318,0.5615500,0.5645411,0.5650451,0.5811612,0.5826908,0.5829479 ,0.5829557,0.5838025,NA,NA,NA,NA,NA))
list(theta=0.5,beta=0.5)
每次到了compile这里就有错误