急急急急急！！关于winbugs做sv模型

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

那位大侠行行好，帮小妹一下？
想做一个sv模型的分析，基本的几个程序已经有了，可是winbug一些使用过程还是看不太懂。。。
不知道如何操作？
谁能帮帮我啊？
不胜感激啊 ~~T T

扫码加我 拉你入群

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

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

分享0 收藏3 回帖

关键词：winbugs WINBUG BUGS SV模型 bug 模型 winbugs

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的教程我也看的一知半解的。。。

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

真的是真的是太感谢了。。。
终于可以走一步了 T T

上面的程序为什么我check的时候总是提示，expected right parenthesis？

5# autozhao
我也遇到同样问题，楼上你现在解决了吗，麻烦告知一下，谢谢了

楼主的程序写的有问题
楼上的问题是因为参数index的值设定不正确
问题段程序修改后如下
model
{                                                              
        for (i in 1:n){
                Ymean[i] <- (rho/tau)*exp(0.5*theta[i])*(theta[i+1]-mu-phi*(theta[i]-mu))
                Yisigma2[i] <- 1/(exp(theta[i])*(1-rho*rho))
                Y[i]~ dnorm(Ymean[i], Yisigma2[i])
            }     
        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[i] <- mu + phi*(theta[i-1]-mu)                                 
                theta[i]~dnorm(thetamean[i], 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)                                                        
}

楼上的 你好 本人刚开始学SV的WINbugs 不知道能不能帮个忙给下你这个程序的原模型形式和各变量都代表的那个参数啊 重谢 万分感激！！

1# 雪雪冰
你好 我也正开始学用WINbugs解SV 听说你那几个程序都有了 能分享下吗 当万分感激 ~~QQ517818168

7# lhdavidmice 你好 ，弱弱问下按你这个程序在输入初始值后，winbugs怎么说 this chain contais unintialzed varibles    ???谢了