Dear BUGS-users,
I am using a similar model explained in the Kidney example. I need to
estimate an aggregated frailty (for a categorical variable with 7 categories
called country) and I would like to prove another distribution instead of a
normal for this random effect.
It is possible to fit the same model with a gamma frailty?, e.g:
log(mu[i,j]) <- alpha + beta.age*age[i,j] + b.frail[country[i]];
}
# Random effects:
b.frail[i] ~ dgamma(0.0001, 0.0001)
}
I have compiled this model and I have obtained an error message:
"Invalid value for node expected positive value".
Note: When I have use
# Random effects:
b.frail[i] ~ dnormal(0, tau.p)
}
with
tau.p~dgamma(0.001,0.001)
it works fine
Any suggestion?
Thanks
Juan R Gonzalez
-------------------------------------------------------------------
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
|