时空模型code

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

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:NL]~spatial.exp(mu[],x[],y[],tau.u,phi,1)
for(i in 1:NL){
mu<-0
}
#AR(1)prior distribution for temporal random effects：
v[1]~dnorm(0.0,tau2)
for(i in 2:NY){
vmean[i-1]<-ro*v[i-1];
v~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 prediction
for(j in l:m){
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:m1){
u.p[j]~spatial.unipred(mu.p[j],x.p[j],y.p[j],u[])
mu.p[j]<-0
}
#priors
Tau.u~dgamma(0.1,0.1)
#wior range for correlation at min distance(O．6kin)is O．01 to O．99
#prior range for correlation at max distance(49kin)is 0 to 0．2
for(i in 1:10){
beta~dnorm(0,0.1)
}
}

这个是二项分布的时空模型，想做poisson怎么修改呀，请各位大神赐教啊

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：code COD ODE model logit 模型

自己写的代码吗？？

请问会写成泊松分布了吗？我也想做泊松分布的，求助