Dear Bugsies,
Is anybody else concerned by the cyclic construction of this graph?
Should we assume that bugs still works when the models are not DAGs or
Chain graphs? Does the compiler check for cyclic structures? Enquiring
minds want to know.
Finn
[snip]
Solution by Andrew Millard
--------------------------
I should have seen this first time - your censoring was not symmetrical.
This updates successfully:
model{
for (i in 1:N) {
a[i] ~ dweib(gammaa, lambdaa)I(leftem[i], aupper[i])
aupper[i] <- min(rightem[i], c[i]) ## ADD THIS
LINE
blower[i] <- max(a[i], leftcar[i])
c[i] ~ dweib(gammab, lambdab)I(blower[i], rightcar[i])
}
You had forced new cvalues of c to be more than a, but new values of a
did
not have to be less than c.
The Jama Valley example on my website is rather more complicated, but
has
the same multiple censoring of unknowns by one another.
Andrew
-------------------------------------------------------------------
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
|