虚拟变量要怎样加入程序中呢?
比如我想将虚拟变量加入到SV-T模型中,使它变成yt=exp(theta(t)/2)*epslion(t)+a*D(i),epslion(t)服从t(0,1,w)分布,theta(t)=mu+phi(theta(t-1)-mu)+b*D(i)+eta(t),eta(t)服从N(0,tau^2)分布。
不加虚拟变量的SV-T模型是
model{
for(i in 1:n)
{y~dt(0,p,omega)
p<-exp(-theta)
}
theta[1]~dnorm(mu,itau2)
for(j in 2:n)
{theta[j]~dnorm(theta2[j],itau2)
theta2[j]<-mu+phi*(theta[j-1]-mu)}
phi<-2*phi1-1
tau<-sqrt(1/itau2)
mu~dnorm(0,0.01)
itau2~dgamma(2.5,0.025)
phi1~dbeta(20,1.5)
omega~dchisqr(8)
}
这个程序要怎么改呢?
到目前为止我还没有找到有关虚拟变量的程序例子,如果有的话,求推荐!
小妹分少,难以悬赏,求高手帮助啊!


雷达卡


京公网安备 11010802022788号







