model{
for (i in 1:n){
y~dnorm(mean.y,v)
mu~dnorm(mean.mu,w1)
beta~dnorm(mean.beta,w2)
mean.y<-mean.mu[i-1]+mean.beta[i-1]
mean.beta<-mean.beta[i-1]
}
v~dnorm(mu1,tau1)
w1~dnorm(mu2,tau2)
w2~dnorm(mu3,tau3)
for ( k in 1:3){
mu[k]~dnorm(0,1.0E-6)
tau[k]~dgamma(0.001,0.0001)}
}
#data
list(n=35,y=c(0.00,-0.70,-1.85,-3.00,-4.60,-5.55,-6.55,-13.90,-16.90,-19.25,-21.60,-23.30,-24.85,-28.80,-33.35,-35.40,-37.10,-36.90,
-36.80,-37.40,-37.70,-37.60,-37.60,-38.40,-38.50,-38.00,-38.45,-39.05,-39.65,-40.70,-40.45,-40.65,-40.85,-40.70,-40.75))
#inite
list(mu=1,mean.beta=c(1),mean.mu=c(1),v=0.1,w1=0.1,w2=0.1)
出现错误mean.mu没有定义或在数据中没有设置
本文来自: 人大经济论坛 S-Plus&R专版 版,详细出处参考:http://www.pinggu.org/bbs/viewthread.php?tid=928006&page=1&from^^uid=1382521