全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
2128 0
2011-10-11
complie 的时候出现 array index is not an integer?怎么解决,我的索引值,明显是整数的,请高手赐教


#model
{
  
C <- 10
      for (i in 1:N) {
          ones[i] <- 1
          ones[i] ~ dbern( ones.p[i])
          ones.p[i] <- exp( l[i] - C )
          # write here the expression of the log-likelihood for i observation
          l[i] <-   0.5*( log( lambda ) - log(2*3.14) - 3*log(y[i]) ) - 0.5* lambda * pow( (y[i]-mu[i])/mu[i], 2 )/y[i]
           lambda<-2#
          mu[i]<-4###加入了一个逆高斯分布,还不完全肯定是否正确
                 
    delta2<-y[i]#逆高斯值
    rdaoshu<-1/r[i]
     tsf<-1/365
     sqts<-sqrt(tsf)
     sqeb<-sqrt(eb)
     powrg<-pow(r[i],gamma)
    ds<-delta*powrg*sqeb*sqts

logit[r[i+1] -r[i]]<-(a0+a1*r[i]+a2*r[i]*r[i]+a3*rdaoshu)*tsf+ds

  }
a0~dnorm(0,1)
a1~dnorm(0,1)
a2~dnorm(0,1)
a3~dnorm(0,1)
eb~dnorm(0,1)
delta2<-delta*delta
gamma~dnorm(0,0.5)
}
#data
list(  N=165,
r=c(0.126530567,0.123957562,0.124254475,0.123264735,0.114353403,0.121482998,0.12059203,0.120394028,0.112471275,
0.113560964,0.114749602,0.113164725,0.114254351,0.11395719,0.113660021,0.112768473,0.110192515,0.110390685,
0.108408836,0.106426659,0.110886096,0.104444154,0.094625931,0.087877445,0.08757963,0.085196845,0.0785424,
0.069697203,0.065819163,0.064327275,0.048800636,0.056068808,0.054675253,0.049099415,0.048601446,0.042125972,
0.041827019,0.042723854,0.052783741,0.042125972,0.049796536,0.042026322,0.033951924,0.028266682,0.028665734,
0.027867616,0.026869895,0.02696967,0.024774409,0.025572732,0.024874202,0.023976032,0.023776429,0.023676627,
0.023676627,0.023576824,0.023377215,0.023377215,0.023377215,0.024774409,0.026371003,0.02617144,0.025672519,
0.025273367,0.025073786,0.024475024,0.024475024,0.024375227,0.023976032,0.02347702,0.023277409,0.023776429,
0.023576824,0.022977986,0.022379119,0.021580583,0.021181295,0.020282849,0.019883518,0.020083185,0.020981646,
0.020881821,0.020981646,0.022379119,0.021580583,0.021281119,0.020582339,0.019683847,0.020183017,0.020981646,
0.021480763,0.021979858,0.026770118,0.028266682,0.025672519,0.022379119,0.023576824,0.022678557,0.020781994,
0.023077795,0.022379119,0.024075832,0.023277409,0.023377215,0.022578745,0.022079675,0.022079675,0.020881821,
0.020582339,0.022578745,0.019783683,0.017886663,0.016388804,0.015689739,0.015489998,0.015589869,0.015489998,
0.014790881,0.015989343,0.017287542,0.015390127,0.013792071,0.014890757,0.016388804,0.016788251,0.019584011,
0.023177602,0.024175631,0.021380941,0.023576824,0.032057143,0.021580583,0.020282849,0.030261811,0.019883518,
0.033253882,0.023576824,0.030860284,0.031059769,0.023676627,0.047406237,0.038338691,0.032755255,0.033054433,
0.034450501,0.032755255,0.025273367,0.031957409,0.033652768,0.034350787,0.030261811,0.032555798,0.031558466,
0.030760541,0.026770118,0.015589869,0.010295582,0.010095752,0.009895918,0.009995836,0.009995836,0.010595321,
0.016189075,0.014890757,0.015689739,0.015689739   ),
y = c(41.273, 72.628, 3.895, 56.98, 5.552, 5.744, 6.399, 10.324, 19.688, 14.923, 70.551, 36.876, 73.614, 0.63, 6.016, 14.83, 5.941, 8.226, 32.997, 5.993, 11.129, 6.265, 1.09, 34.514, 22.988, 13.241, 1.984, 1.587, 3.534, 125.5, 85.96, 15.998, 305.689,  68.492, 72.521, 7.445, 2.273, 29.444, 11.045, 3.131, 5.312, 33.617, 23.875, 41.723, 52.23, 5.472, 1.327, 3.533, 2.698, 34.683, 3.868, 20.195, 10.076, 8.767, 9.262, 28.838, 8.238, 8.719, 22.366, 15.731, 6.479, 5.372, 1150.242, 2.373, 5.696, 27.043, 10.092, 9.669, 21.733, 1.151, 11.928, 1.313, 3.558, 0.867, 8.658, 23.43, 0.361, 37.741, 6.074, 10.722, 5.246, 21.654, 4.344, 2.405, 12.9, 4.277,  4.003, 0.833, 3.549, 0.091, 12.15, 9.041, 2.137, 2.016, 8.747, 1.256, 25.758, 16.548, 9.812, 3.415,6.016, 14.83, 5.941, 8.226, 32.997, 5.993, 11.129, 6.265, 1.09, 34.514, 22.988, 13.241, 1.984, 1.587, 3.534, 125.5, 85.96, 15.998, 305.689,  68.492, 72.521, 7.445, 2.273, 29.444, 11.045, 3.131, 5.312, 33.617, 23.875,20.195,3.558, 0.867, 8.658, 23.43, 0.361, 37.741, 6.074, 10.722, 5.246, 21.654, 4.344, 2.405, 12.9, 4.277,  4.003,0.833, 3.549, 0.091, 12.15, 9.041, 2.137, 2.016, 8.747, 1.256, 25.758, 16.548, 9.812, 3.415,6.016, 14.83,7.445, 2.273, 29.444, 11.045, 3.131)
)

#initial values
list(a0=0,a1=0,a2=0,a3=0,eb=0,delta=1,gamma=0)

二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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