WinBUGS时空模型代码

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

小弟初学Bayesian，完全不懂，下面是一个二项分布时空模型，在compile阶段出错：Variable is not difined。请帮忙看看哪儿出错了呀，谢谢！model{for(j in l:N){z[j]~dbin(p[j],n[j])logit(p[j])<-beta[1]+beta[2]*lst[j]+beta[3]*ndvi[j]+beta[4]*dist[j]+beta[5]*pipe[j]+u[loc[j]]+v[year[j]]}#Spatial random effects?u[1:9]~spatial.exp(mu[],x[],y[],tau.u,phi,1)for(i in 1:NL){mu[i]<-0}#AR(1)prior distribution for temporal random effects?v[1]~dbin(0.0,tau2)for(i in 2:5){vmean[i-1]<-ro*v[i-1];v[i]~dnorm(vmean[i-1],tau.v)}tau2<-1/((1-pow(ro,2))*pow(sigma.v,2))ro~dunif(-1,1)tau.v~dgamma(0.1,0.1)sigma.v<-1/sqrt(tau.v)#validation or predictionfor(j in l:44){z.p[j]~dbin(p.p[j],n.p[j])logit(p[j])<-beta[1]+beta[2]*lst.p[j]+beta[3]*ndvi.p[j]+beta[4]*dist.p[j]+beta[5]*pipe.p[j]+ u.p[loc.p[j]]+v[year.p[j]]}for(j in 1:9){u.p[j]~spatial.unipred(mu.p[j],x.p[j],y.p[j],u[])mu.p[j]<-0}#priorsTau.u~dgamma(0.1,0.1)for(i in 1:10){beta[i]~dnorm(0,0.1)}}

list(z=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),n=c(441968,438851,436556,436575,438763,548258,544790,541934,542007,544168,884737,817890,870929,868917,870117,778772,671350,764461,761944,762196,725694,720874,715930,715299,717247,485215,481381,477763,477449,478438,400178,397388,394597,394916,396658,610793,606985,603531,603831,606470,322986,321379,319846,320180,321656),lst=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),ndvi=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),dist=c(29.17,29.17,29.17,29.17,29.17,30.54,30.54,30.54,30.54,30.54,33.94,33.94,33.94,33.94,33.94,36.43,36.43,36.43,36.43,36.43,29.97,29.97,29.97,29.97,29.97,26.48,26.48,26.48,26.48,26.48,22.57,22.57,22.57,22.57,22.57,23.39,23.39,23.39,23.39,23.39,25.52,25.52,25.52,25.52),pipe=c(0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,1,1,1,1,1,0,0,0,0,0,1,1,1,1,1,0,1,1,1,1,1,1,1,1,1,1,1,1,1),loc=c(1,1,1,1,1,2,2,2,2,2,3,3,3,3,3,4,4,4,4,4,5,5,5,5,5,6,6,6,6,6,7,7,7,7,7,8,8,8,8,8,9,9,9,9),year=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),N=44 )

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS bug Win 二项分布 模型

求指教winbugs贝叶斯model

同求指教贝叶斯的model编写

你好，我是学医学的，有贝叶斯统计建模问题（主要是WINBUGS）软件请教你可以么？qq 1400236514，我可以交学费。