[疑难杂症]这个Winbugs程序为什么不能编译

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这里就有错误

试运行了一下，出错信息是：variable x is not defined.
我猜是因为你的x数据最后5个数据是NA，这样L[k]等等就无法运算。可能要把那5个NA改成数据才行。

TimeT 发表于 2013-12-21 12:21
试运行了一下，出错信息是：variable x is not defined.
我猜是因为你的x数据最后5个数据是NA，这样L[k]等 ...

我是想试一下，这几个数据是缺失的情况能不能运行。

你这想x[i]用错了 哦