全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 winbugs及其他软件专版
1474 3
2016-01-04
参数x服从beta分布,这里拟合beta分布惨回溯,在Update时出错了
model
{
   for(i in 1:N)
           {

              x[i]~dbeta(a,b)
           }
                      a~dunif(0,2)
                      b~dunif(0,2)
}
# Data
        list(N=54,x=c(0.239,0.111,0.865,0.481,0.1488,0.378,0.3,0.3595,0.3776,0.454,
0.36,0.2553,0.6126,0.796,0.454,1,0.343,0.426,0.631,0.556,0.426,0.372,0.387,0.248,0.372,0.2505,0.5367,0.3736,0.672,0.267,0.579,0.582,0.387,0.3396,0.7098,0.6288,0.6094,0.5816,0.5522,0.2068,
0.528,0.316,0.1189,0.357,0.4187,0.491,0.371,0.498,0.102,0.2068,0,0.283,0.641,0.387))
#initial data
        list(a=2,b=2);

二维码

扫码加我 拉你入群

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

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

全部回复
2016-1-4 20:57:54
the upper limit of beta distribution should be less than 1,and the lower limit of beta distribution should be large than 0!so,the data you estimated the model may be wrong.Replacing 1 and 0 in the data you used by others in (0,1)
二维码

扫码加我 拉你入群

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

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

2016-1-5 22:19:10
同意楼上,你可以把数据中的那个1改成0.99999999,把其中那个0改成0.00000001,就能正常运行(我试过了)。
二维码

扫码加我 拉你入群

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

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

2016-1-5 22:19:10
同意楼上,你可以把数据中的那个1改成0.99999999,把其中那个0改成0.00000001,就能正常运行(我试过了)。
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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