<James.Dell <at> csiro.au> writes:
...
> 1) I can't seem to send variables classed as factors (Month), is there a
> way do this?
You can not use factors per se in BUGS. You have to convert them to numeric
(integer) variables before.
> 2) Checking the Log in WinBUGS I can see that the model is Syntactically
> correct, but Bugs is not able to recognise the the initial values for
> each of the chains.
...
> model
> {
> #Centre variables
> mSeaW <- mean(SeaW[])
> s_dSeaW <- sd(SeaW[])
>
> #normalise Variables
> nSeaWiFS <- mSeaW/s_dSeaW
>
> for(i in 1:N) {
> log(lambda[i]) <- delta0 + alpha1 * Month[i] + alpha2 * Lat[i] +
> beta1 * (SeaW[i] - nSeaW)
I guess you want to use mSeaW here instead of nSeaW!
^ ^
> # recalculate the original intercept term
> Intercept <- delta0 - beta1 * nSeaW
I guess you want to use nSeaWiFS here instead of nSeaW!
...
>
> ##And the Error Message from WinBUGS
>
> model is syntactically correct
>
> data(C:/Program Files/R/R-2.6.0/data.txt)
>
> data loaded
>
> compile(3)
>
> made use of undefined node nSeaW
This warned you that the model is not OK!
Gregor
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
