model
{
a ~ dunif (0, 1000) # uninformative prior,constrained to be positive
b ~ dnorm (0, 1.0E-6) # uninformative prior
c ~ dunif(0, 10) # uninformative prior, constrained to be positive
prec ~ dgamma(0.001, 0.001) # uninformative prior
mTree <- mean(TreeDens[]) # mean of the explanatory variable
for (i in 1:16) # for each lake sampled
{
pred[i] <- a/(1 + exp(-(b+c*(TreeDens[i]- mTree)))) # predicted CWD
logpred[i] <- log(pred[i]) # take the logarithm of the prediction
CWD[i] ~ dlnorm(logpred [i], prec) # CWD drawn from lognormal
}
for (j in 8:21) # make predictions fro densities 800- 2100
{
pred[j] <- a/(1 + exp(-(b+c*(j*100-mTree))))
}
}
list(TreeDens = c(1270, 1210, 1800, 1875, 1300, 2150, 1330, 964, 961, 1400, 1280, 976, 771, 833, 883, 956),
CWD = c(121, 41, 183, 130, 127, 134, 65, 52, 12, 46, 54, 97, 1, 4, 1, 4))
list(a=0,b=0,c=0,prec=100)
在compile时,出现“multiple definitions of node pred[8]”错误信息。该如何修改呢?谢谢了!


雷达卡



京公网安备 11010802022788号







