大家帮帮忙看看，面板ordered probit （winbugs ）代码

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

大家好，我刚开始用winbugs做贝叶斯，很多地方都不懂，下面我写了一个面板数据的ordered probit模型，大家帮我看看哪里不对，老是不对报错。
model{
for (i in 1:N){for (t in 1:T){
fVie[i,t]~dcat(p[i,t,])
p[i,t,1]<-  phi(r-mu[i,t])
p[i,t,2]<-  phi(r+d-mu[i,t])        -phi(r-mu[i,t])
p[i,t,3]<-  phi(r+d+e-mu[i,t])      -phi(r+d-mu[i,t])
p[i,t,4]<-  phi(r+d+e+f-mu[i,t])    -phi(r+d+e-mu[i,t])
p[i,t,5]<-  phi(r+d+e+f+g-mu[i,t])  -phi(r+d+e+f-mu[i,t])
p[i,t,6]<-  phi(r+d+e+f+g+h-mu[i,t])-phi(r+d+e+f+g-mu[i,t])
p[i,t,7]<-1-phi(r+d+e+f+g+h-mu[i,t])
mu[i,t]<-alpha+b1*marriage[i,t]+b2*ad[i,t]+b3*kid[i,t]+b4*heal[i,t]+b5*dlRev[i,t]+b6*my[i,t]
}}


for (i in 1:N){alpha~dnorm(alphabar,phi)}

alphabar~dnorm(0,0.0001)
b1~dnorm(0,0.0001)
b2~dnorm(0,0.0001)
b3~dnorm(0,0.0001)
b4~dnorm(0,0.0001)
b5~dnorm(0,0.0001)
b6~dnorm(0,0.0001)
phi~dgamma(0.1,0.1)
r~dnorm(0,0.001)
d~dgamma(0.1,0.1)
e~dgamma(0.1,0.1)
f~dgamma(0.1,0.1)
g~dgamma(0.1,0.1)
h~dgamma(0.1,0.1)
}

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs ordered WINBUG Probit Order 模型

面板也可以做有序 probit回归？