全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 winbugs及其他软件专版
1632 1
2014-05-24
悬赏 100 个论坛币 未解决

模型语法没错,数据加载没出问题。就是编译时说变量change没定义,求大神指点!!!!!!

Model{

   for(i in 1:N-2)

  {

     punif<-1/(N-2)

  }

  for(j in 1:6)

  {

     tao[j]~dgamma(1.001,0.001)

  }

K[1]~dcat(punif[])

K[2]~dcat(punif[])

K[3]~dcat(punif[])

K[4]~dcat(punif[])

K[5]~dcat(punif[])

  for( j in 1:5)

  {

     x.change[j]<-ranked(k[],j)+1

   }

  change[1]<-0

  change[7]<-300

for(j in 1:6)

{

mu[j]~dnorm(0,mtao[j])

mtao[j]<-mean(y[change[j]:change[j+1]])*mean(y[change[j]:change[j+1]])

}

for(j in 1:6)

{

   for(i in change[j]+1:change[j+1])

    {

       y~dnorm(mu[j],tao[j])

    }

}

}

list(y=c(0.007346,-0.01004,0.00428,0.006131,0.004033,-0.00043,-0.00048,-0.00855,-0.02353,0.007209,-0.00601,-0.01132,0.02467,0.017108,-0.00123,-0.00984,-0.01302,-0.01002,-0.01817,-0.00243,0.052088,-0.00817,0.021262,0.007752,-0.00215,-0.01354,0.020061,0.019121,0.000151,0.000362,0.001904,-0.00492,-0.00199,-0.01232,0.005337,-0.00513,0.004668,0.018785,-0.00112,-0.02348,0.002879,0.002755,-0.01901,-0.01074,-0.00969,-0.01301,0.002964,-0.01979,-0.00438,0.00543,-0.00212,-0.00459,-0.00425,-0.01896,0.00736,0.035208,-0.01698,-0.00937,0.019217,0.008703,0.018621,-0.00606,-0.01744,0.023255,-0.01117,-0.01417,-0.01593,0.001808,-0.01318,-0.00988,-0.00288,-0.00964,0.006339,-0.01433,0.017569,-0.00857,-0.02473,0.004171,-0.00999,-0.0182,0.003286,-0.01104,0.003801,0.003019,0.000462,-0.00382,-0.0262,0.0028,-0.00159,-0.01245,-0.00656,0.00778,0.000543,-0.02053,-0.0225,-0.00239,0.020255,-0.00333,0.078904,0.013772,-0.02033,-0.0018,-0.01166,-0.01905,0.012205,-0.00037,0.027268,-0.01275,0.000496,-0.00759,0.007451,0.020443,-0.01438,-0.00325,-0.02201,-0.02373,-0.00791,-0.00789,-0.00531,0.004876,-0.02005,-0.00639,0.033694,-0.00842,0.002909,-0.01384,-0.01376,0.002229,0.006581,-0.00041,0.024858,-0.00088,0.025879,0.015816,-0.00309,-0.00324,0.005447,0.013762,0.003311,-0.00462,0.023444,0.008942,0.011916,0.010786,0.015793,-0.01332,0.016214,-0.00823,0.013127,-0.03828,0.00187,0.007288,-0.00749,0.014829,0.004563,-0.00052,-0.01499,-0.00699,0.014216,0.018857,0.003299,0.006185,-0.01933,0.017729,0.0000854,-0.00402,-0.00119,0.015824,0.008329,0.0000164,-0.00357,0.006312,-0.00658,-0.02052,-0.02388,-0.0069,0.002655,-0.01906,-0.00129,0.020734,0.000108,-0.01452,0.01589,0.00402,-0.00799,-0.0114,-0.0072,0.008648,-0.00579,0.001181,0.004715,-0.00013,-0.01702,-0.02231,0.000566,-0.01553,0.010995,-0.00263,0.013525,-0.00865,0.004347,0.00055,0.008592,-0.0018,-0.01808,0.001739,-0.00142,-0.00105,0.007678,-0.00053,0.019608,0.002626,-0.01918,0.006434,0.006875,0.001384,-0.00368,-0.01253,0.002072,-0.00047,-0.00407,-0.01389,0.007935,0.010808,-0.00118,0.013703,0.002592,0.001324,0.006723,-0.00163,0.003512,0.003756,0.004048,-0.00493,0.003953,0.00845,0.010385,-0.00219,0.002376,0.011026,-0.00755,0.016999,0.013713,0.010099,0.005151,0.004064,-0.00787,0.012841,-0.00428,-0.01534,0.004619,0.020405,0.01467,0.002977,0.000371,-0.00296,0.004768,0.02324,-0.0043,0.006189,-0.0161,0.010346,-0.00235,0.005212,0.009932,-0.02221,-0.00254,0.000106,0.016336,-0.00326,0.003592,0.006203,0.000247,0.001421,0.0058,-0.01614,0.005918,-0.00336,-0.02279,-0.00759,-0.00418,0.000391,0.011193,-0.0005,0.012446,-0.00049,-0.00372,0.014822,0.001533,0.004879,0.004425,-0.00597,-0.00349,0.006644,0.005204,-0.00834,0.002759,0.016179,0.007796,0.007763),N=300)

二维码

扫码加我 拉你入群

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

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

全部回复
2014-5-30 23:10:23
很多问题,不可能check model时没有错
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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