我正在尝试在OpenBUGS中运行模型。但我收到错误"变量TT未定义"。这是我的BUGS代码:
model{
for(t in 2:TT){
for(i in 1:N){
y[t,i] ~ dmnorm(y1[t,i], tau[i])
}
y1[t,1] <- (p[1,1] * y[t-1,1]) + (p[1,2] *y[t-1,2]) + (p[1,3] *y[t-1,3])
y1[t,2] <- (p[2,1] * y[t-1,1]) + (p[2,2] *y[t-1,2]) + (p[2,3] *y[t-1,3])
y1[t,3] <- (p[3,1] * y[t-1,1]) + (p[3,2] *y[t-1,2]) + (p[3,3] *y[t-1,3])
}
for (s1 in 1:N) {
for (s2 in 1:N) {
p[s1, s2]~ dnorm(1, 1)
}
}
for(i in 1:N){
tau[i] ~ dgamma(0.001, 0.001)
}
}
这是运行BUGS代码的数据和R代码:
Data=structure(.Data=c(0.1, 0.2, 0.5, 0.2, 0.2, 0.8, 0.2, 0.4, 0.3, 0.2, 0.1, 0.6, 0.1, 0.5, 0.6), .Dim=c(5,3))
TT=nrow(Data)
N=ncol(Data)
y=as.matrix(Data)
sp.data = list(y, N, TT)
Sp.sim0<- bugs (sp.data, inits=NULL,
parameters=c('p'), model.file= "model2.txt",n.chains=2, DIC=TRUE,n.iter=1000, codaPkg=FALSE,debug=TRUE)
有任何请帮助解决这个问题吗?提前谢谢。