全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
7985 17
2010-06-22
那位大侠行行好,帮小妹一下?
想做一个sv模型的分析,基本的几个程序已经有了,可是winbug一些使用过程还是看不太懂。。。
不知道如何操作?
谁能帮帮我啊?
不胜感激啊 ~~T T
二维码

扫码加我 拉你入群

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

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

全部回复
2010-6-22 14:02:10
model
{                                                               
for (i in 1:n) {
Ymean<-rho/tau*exp(0.5*theta)*(theta[i+1]-mu-phi*(theta-mu));
Yisigma2 <- 1/(exp(theta)*(1-rho*rho));
y~ dnorm(Ymean,Yisigma2);   
    }     
isigma2<-(1-phi*phi)*itau2;                                                        
theta0 ~ dnorm(mu,isigma2);  
thetamean[1] <- mu + phi*(theta0-mu);   
theta[1] ~ dnorm(thetamean[1],itau2);                                 
for (i in 2:(n+1)) {                                                                    
thetamean <- mu + phi*(theta[i-1]-mu);                                   
theta~dnorm(thetamean,itau2);}                                 
phi1 ~ dbeta(20,1.5);     
phi <- 2*phi1-1;     
mu ~ dnorm(0,0.04);                                                              
itau2 ~ dgamma(2.5,0.025);                                                      
tau <- sqrt(1/itau2) ;
rho ~ dunif(-1,1)                                                        
}
list(Y=c(0.13,0.23,0.4,0.53,0.21,0.23,0.43,0.21,0.4,0.21,0.22,0.22,0.11,0.22,0.33,0.23)
n=16)
list(mu=0, phi=0, tau=0)


这是我的程序 上方黑色的部分已经检验没有问题了 可是我后边的list部分却检验部通过
是哪里错了 啊  winbugs的教程我也看的一知半解的。。。
二维码

扫码加我 拉你入群

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

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

2010-6-22 16:52:17
model
{                                                               
for (i in 1:n) {
Ymean<-rho/tau*exp(0.5*theta)*(theta[i+1]-mu-phi*(theta-mu));
Yisigma2 <- 1/(exp(theta)*(1-rho*rho));
Y~ dnorm(Ymean,Yisigma2);   
    }     
isigma2<-(1-phi*phi)*itau2;                                                        
theta0 ~ dnorm(mu,isigma2);  
thetamean[1] <- mu + phi*(theta0-mu);   
theta[1] ~ dnorm(thetamean[1],itau2);                                 
for (i in 2:(n+1)) {                                                                    
thetamean <- mu + phi*(theta[i-1]-mu);                                   
theta~dnorm(thetamean,itau2)}                                 
phi1 ~ dbeta(20,1.5);     
phi <- 2*phi1-1;     
mu ~ dnorm(0,0.04);                                                              
itau2 ~ dgamma(2.5,0.025);                                                      
tau <- sqrt(1/itau2) ;
rho ~ dunif(-1,1)                                                        
}
#data
list(Y=c(0.13,0.23,0.4,0.53,0.21,0.23,0.43,0.21,0.4,0.21,0.22,0.22,0.11,0.22,0.33,0.23),
n=16)
#inits
list(mu=0, phi1=0.975, itau2=50)

#load inits
#gen inits
二维码

扫码加我 拉你入群

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

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

2010-6-22 17:00:02
真的是真的是太感谢了。。。
终于可以走一步了 T T
二维码

扫码加我 拉你入群

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

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

2010-6-24 03:41:00
上面的程序为什么我check的时候总是提示,expected right parenthesis?
二维码

扫码加我 拉你入群

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

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

2010-7-7 20:10:25
5# autozhao
我也遇到同样问题,楼上你现在解决了吗,麻烦告知一下,谢谢了
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

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