When I try to use the following expression in my BUGS program, I get a trap
message.
mu[i,j] <-
log(1000000*(((k*f*dose[i,j])*((km*exp(-k*t[i,j]))*(exp(k)-exp(-k))-
(k*exp(-km*t[i,j]))*(exp(km)-exp(-km)))) / (k*(km-k)*(0.06*2))))
The program runs fine if I do not use the log of the formula. The formula is
never zero or negative unless dose is zero and I have excluded this
possibility by constraining dose. If anyone has any ideas, please let me
know.
Bonita Singal
University of Cincinnati
-------------------------------------------------------------------
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
|