全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 winbugs及其他软件专版
916 0
2018-02-23
这是winbugs的model,但是运行第一步出现错误,显示“expected right parenthesis”,大家有可以找出问题的吗?

model
  {
     for(i in 1:N){
                Y[i]~dpois(lambda[i])
                log(lambda[i])<-beta0
                        +X1[i]*beta
                        +X2[i]*beta
                        +theta[i]
                        +phi[i]
                        theta[i]~dnorm(0,1/tau)
        }
        beta0~dnorm(0.0,1.0E-5)
        beta~dnorm(0,1000)
        tau~dgamma(0.001,0.001)
        for(j in 1:29){
                phi[j]~dnorm(Q[j],1/tau[j])
        }
        for(j in 1:N){
                for(k in 1:N){
                        A=A+w[j,k]
                        B=B+phi[k]*w[j,k]
                }
                        Q[j,k]<-B/A
                        tau[j]<-tau[c]/A
                }
        }
        tau[c]~dgamma(0.1,0.1)
}

Data list(N=29,Y=c(7674394,985200,2692194,2417510,399400,750000,4055704,645400,22804610,2110700,680300,
                                2463754,581520,9368938,465352,1697500,863450,667900,1671300,3936600,1510300,3525301,337500,1776700,
                                3836500,4829825,838600,622500,4000)
                                X1=c(404.655,547.202,476.307,455.824,614.270,253.589,249.091,409.959,4222.248,930.228,387.508,270.839,
                                                258.339,3557.005,383.271,371.705,332.901,639.661,2009.795,501.011,537.530,498.393,152.509,366.851,
                                                361.744,114.875,104.050,82.930,115.226)
                                X2=c(42.890,57.999,50.485,34.273,65.108,19.067,26.402,43.452,157.749,76.878,32.025,22.383,21.350,150.171,
                                                25.311,24.547,21.984,42.243,100.771,33.086,22.650,32.913,10.072,24.226,27.301,8.670,4.769,3.801,5.281)
                                w=structure(
                                        .Data=c(0,1,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
                                                                1,0,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
                                                                0,1,0,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
                                                                1,1,1,0,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
                                                                0,0,1,1,0,0,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
                                                                0,0,0,1,0,0,1,0,0,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
                                                                0,0,0,1,1,1,0,0,1,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
                                                                0,0,0,0,1,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
                                                                0,0,0,0,0,0,1,1,0,0,1,1,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
                                                                0,0,0,0,0,1,0,0,0,0,1,0,1,0,1,1,0,1,0,0,0,0,0,0,0,0,0,0,0,
                                                                0,0,0,0,0,0,1,0,1,1,0,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
                                                                0,0,0,0,0,0,0,0,1,0,1,0,1,1,1,0,0,1,1,0,0,0,0,0,0,0,0,0,0,
                                                                0,0,0,0,0,0,0,0,0,1,1,1,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
                                                                0,0,0,0,0,0,0,0,1,0,0,1,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,
                                                                0,0,0,0,0,0,0,0,0,1,0,1,1,0,0,0,0,1,1,0,0,0,0,0,0,0,0,0,0,
                                                                0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,1,1,0,0,0,0,0,0,0,0,0,0,0,
                                                                0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,1,1,1,0,0,0,0,0,0,0,0,0,
                                                                0,0,0,0,0,0,0,0,0,1,0,0,0,0,1,1,1,0,1,0,0,0,0,0,0,0,0,0,0,
                                                                0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,0,1,1,0,1,1,0,0,0,0,0,0,0,0,
                                                                0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,0,1,0,1,1,0,1,1,0,0,0,0,
                                                                0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,0,0,0,0,1,0,0,0,0,
                                                                0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,1,1,0,0,0,0,0,
                                                                0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,1,0,1,0,0,0,
                                                                0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,1,1,0,1,1,0,0,0,
                                                                0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,0,0,1,0,1,1,0,0,
                                                                0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,0,0,0,0,
                                                                0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,0,0,1,0,0,0,0,
                                                                0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,1,
                                                                0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0)))
       
Inits  list(beta0 = 0.5, beta= 0, A=0,B=0,tau = 1)


二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

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