Hi All,
I am trying run the model below for an AR(1) process
but I get the “undefined real result” error message.
Does anyone have a suggestion on how to solve this
problem? I have already tried different initial values
but this does not seem to make any difference.
Any suggestion will be dearly appreciated.
Best regards,
Enrique Cuervo Guzmán
[log in to unmask]
[log in to unmask]
model
{
for(i in 2:N)
{
y[i] ~ dnorm(mu[i],tau)
mu[i] <- alpha+rho*y[i-1]
}
rho <- B + (1 - B)*U
B ~ dbern(0.5)
U ~ dunif(0,1)
tau ~ dgamma(0.1, 0.1)
alpha ~ dnorm(0, tau)
}
Initial values:
list(tau=1500, alpha=0, B=0,U=1)
Enrique Cuervo Guzmán.
[log in to unmask]
5
__________________________________________________
Do You Yahoo!?
Tired of spam? Yahoo! Mail has the best spam protection around
http://mail.yahoo.com
-------------------------------------------------------------------
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
|