虚拟变量要怎样加入程序中呢？
比如我想将虚拟变量加入到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)
}
这个程序要怎么改呢？
到目前为止我还没有找到有关虚拟变量的程序例子，如果有的话，求推荐！
小妹分少，难以悬赏，求高手帮助啊！