RT。其实就是WinBUGS,我用来做SV-T模型,但是却报错:
model
{
for( i in 1 : n )
{
y~dt(0,p,omega)
p<-exp(-theta)
}
theta[1]~dnorm(mu,tau)
for(j in 2:n)
{
theta[j]~dnorm(theta2[j],tau)
theta2[j]<-mu+phi*(theta[j-1]-mu)
}
phi<-2*phi1-1
mu~dnorm(0,0.01)
tau~dgamma(2.5,0.025)
phi1~dbeta(20,1.5)
omega~dchisqr(8)
}
list(phi1= 0.5,tau =0.5,mu=0)
我数据是从文件中载入的。
在gen 初值时报错:unable to generate initial values for node[02318080H] of type GrahpT.Mixing.
请问什么问题?谢谢!