1374 0

关于openbugs数据输入的问题 [推广有奖]

  • 0关注
  • 0粉丝

硕士生

89%

还不是VIP/贵宾

-

威望
0
论坛币
1243 个
通用积分
7.4500
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
48331 点
帖子
48
精华
0
在线时间
379 小时
注册时间
2016-9-4
最后登录
2024-4-26

10论坛币
我的模型是时间序列ARMA(1,1)模型,用openbugs的check model 和load data 都通过了,但是compile的时候显示array index is greater than array upper bound than Y,请教是我输入的数据有问题还是代码有问题。附代码:model
{
for(t in 1:N+1){
    Y[1]~dnorm(mu[T],tau);}
    mu[1]<-a[1]+M
for(t in 2:N){
     mu[t]<-a[t]+theta1*Y[t-1]+phi1*a[t-1];
                   }

M~dnorm(0,1.0E-2)
tau~dgamma(0.1,0.1);

for(t in 1:N+1){
       a[t]~dnorm(0,tau)
                   }
beta1~dbeta(1,1);
theta1<-2*beta1-1
gamma1~dbeta(1,1);
phi1<-2*gamma1-1

}

输入数据的格式就是list(N=2276,
mu=c(0.0·······))


想问一下哪里需要修改



关键词:时间序列 问题。
您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注jltj
拉您入交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-27 17:30