全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 winbugs及其他软件专版
2507 3
2013-12-20
悬赏 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这里就有错误

二维码

扫码加我 拉你入群

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

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

全部回复
2013-12-21 12:21:07
试运行了一下,出错信息是:variable x is not defined.
我猜是因为你的x数据最后5个数据是NA,这样L[k]等等就无法运算。可能要把那5个NA改成数据才行。
二维码

扫码加我 拉你入群

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

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

2013-12-21 18:28:24
TimeT 发表于 2013-12-21 12:21
试运行了一下,出错信息是:variable x is not defined.
我猜是因为你的x数据最后5个数据是NA,这样L[k]等 ...
我是想试一下,这几个数据是缺失的情况能不能运行。
二维码

扫码加我 拉你入群

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

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

2014-6-7 22:40:00
你这想x[i]用错了 哦
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群