空间相关性

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

model
{
for(i in 1:m)
{
y~dpois(mu)
log(mu)<-alpha+v+u+beta1*quark+beta2*bed
theta<-exp(alpha+v+u+beta1*quark+beta2*bed)
res_theta<-exp(alpha+v+u)
pp<-step(res_theta+eps)
v~dnorm(0,tau.v)
RR_exp<-exp(beta1*quark+beta2*bed)
RR_het<-exp(v)
RR_clust<-exp(u)
}
eps<-1.0E-6
u[1:m]~car.normal(adj[],weights[],num[],tau.u)
for(k in 1:SumNumNeigh)
{
weights[k]<-1
}
alpha~dflat()
mean<-exp(alpha)
beta1~dnorm(0.0,1.0E-5)
beta2~dnorm(0.0,1.0E-5)
tau.u~dgamma(0.5,0.0005)
tau.v~dgamma(0.5,0.0005)
}



list(m=44,y=c(449,437,418,428,304,439,412,295,372,298,955,997,647,542,430,849,865,588,574,458,689,793,624,578,440,615,571,475,343,297,282,376,261,234,195,555,643,521,402,369,489,367,262,283,198),
bed=c(7.36,7.76,7.71,6.81,7.41,11.3,12.24,12.11,12.78,12.78,7.68,8.08,7.97,7.63,7.84,11.13,10.66,12.25,13.23,11.82,5.06,4.8,4.72,4.71,4.71,2.57,2.74,2.74,4.08,3.03,2.95,3.12,3.14,3.71,3.23,3.52,3.37,3.62,3.92,3.61,5.98,5.88,5.91,6.48,6.06),
quark=c(1,1,2,1,1,1,1,2,1,1,1,1,2,1,1,1,1,2,1,1,1,1,2,1,1,1,1,2,1,1,1,1,2,1,1,1,1,2,1,1,1,1,2,1,1),
u=c(30.6,30.6,30.6,30.6,30.6,30.68,30.68,30.68,30.68,30.68,30.73,30.73,30.73,30.73,30.73,30.61,30.61,30.61,30.61,30.61,30.69,30.69,30.69,30.69,30.69,30.6,30.6,30.6,30.6,30.6,30.79,30.79,30.79,30.79,30.79,30.84,30.84,30.84,30.84,30.84,30.73,30.73,30.73,30.73,30.73),
v=c(0,1,2,3,4,0,1,2,3,4,0,1,2,3,4,0,1,2,3,4,0,1,2,3,4,0,1,2,3,4,0,1,2,3,4,0,1,2,3,4,0,1,2,3,4),
adj=c(0,1,1,1,1,1,1,1,1),
num=c(15,8.7, 13.6,14.5, 17.5,11.6,11.7,31.1),
SumNumNeigh=8
)

请教各位高手此代码哪儿错了？

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：空间相关性 空间相关 相关性 model mode 相关性 空间

adj=c(1,2,2,2,2,2,2,2,1),
num=c(30.68,30.6,30.73,30.68,30.61,30.73,30.69,30.61,30.60,30.69,30.79,30.60,30.84,30.79,30.73,30.73),
SumNumNeigh=16

学习了。。。

o牧孑 发表于 2014-8-21 16:56
学习了。。。

程序编译出现错误，不知道哪儿错了

这是什么软件啊大神