Hi all,
I'm using winbugs for the first time and I'm doing some examples just to
learn more about how winbugs works.
I processed the model below without problema.
model;
{
for( i in 1 : N ) {
Y[i] ~ dt(mu[i],1,20)
mu[i] <- beta * x[i]
}
beta ~ dflat()
}
list(x = c(1, 2, 3, 4, 10), Y = c(1, 1, 4, 4, 1), N=5) # Data
list(beta = 0) # Inital Values
Now, I would like to process the model varying the value of one
obersations, namely, Y[5], i.e., I would like to consider the data
"list(x = c(1, 2, 3, 4, 10), Y = c(1, 1, 4, 4, k), N=5)", where k is the
value I would like to vary, say k=1,2,...,50.
How is the most efficient way to do it? is it necessary to insert K as a
node?
Thank you.
Ailton.
-------------------------------------------------------------------
To mail the BUGS list, mail to [log in to unmask]
You can search old messages at www.jiscmail.ac.uk/lists/bugs.html
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
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]
|