I am Jupita, undergraduate student department of statistic Sepuluh Nopember Institute of Technology, Surabaya, Indonesia. I am using Bayesian Hierarchy for my final project. When using openbugs, I have obstacles when using this program. There is error message "initial value loaded but chain uninitialized variables" when I load inits, I continue the step whit gen inits, but there is error message "error for node y[1,1] of type graphGPD.StdNode first aqrgument must be data". So, I want to ask you,:
1) What this error message mean?
2) How can I solve this problem?
This is my model, initial chain and data:
model{
for( j in 1 : N ) {
for( i in 1 : T ) {
y[i , j] ~ dgpar(mu[j], tau[j], teta[j])
}
tau[j] ~ dgamma(0.001, 0.001)
teta[j] ~ dgamma(0.001, 0.001)
mu[j] <- beta[1] + beta[2] * x1[j] + beta[3] * x2[j] + beta[4] * x3[j] + beta[5] * x4[j]
sigma[j] <- 1 / sqrt(tau[j])
}
beta[1] ~ dnorm(0.0, 1.0E-6)
beta[2] ~ dnorm(0.0, 1.0E-6)
beta[3] ~ dnorm(0.0, 1.0E-6)
beta[4] ~ dnorm(0.0, 1.0E-6)
beta[5] ~ dnorm(0.0, 1.0E-6)
tau.r ~ dgamma(0.001, 0.001)
sigma.r <- 1 / sqrt(tau.r)
}
INISIALISASI
list(beta=c(1, 1, 1, 1, 1), teta=c(0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5), tau=c(1, 1, 1, 1, 1, 1, 1, 1, 1), tau.r=0.8)
DATA
list(x1=c(27.1, 26.8, 26.3, 28.75, 28.35, 27.4, 27.05, 25.7, 27.35), x2=c(14, 100, 46, 2, 4, 0, 8, 207, 26), x3=c(106.38, 106.11, 106.57, 106.88, 106.84, 106.88, 106.65, 106.75, 106.75), x4=c(-6.10, -6.11, -6.29, -6.11, -6.16, -6.26, -6.12, -6.50, -6.26), N=9, T=985,
y=structure(
.Data=c(16,20,94,54.5,19.1,49,29,40,45
,
...
,NA,NA,NA,NA,NA,NA,NA,76.7,NA),
.Dim=c(985,9)))
I really need your help. Thanks before
-------------------------------------------------------------------
This list is for discussion of modelling issues and the BUGS software.
For help with crashes and error messages, first mail [log in to unmask]
To mail the BUGS list, mail to [log in to unmask]
Before mailing, please check the archive at www.jiscmail.ac.uk/lists/bugs.html
Please do not mail attachments to the list.
To leave the BUGS list, send LEAVE BUGS to [log in to unmask]
If this fails, mail [log in to unmask], NOT the whole list
|