As explicitely stated in the manual, logical nodes cannot be associated
to data values. As a consequence the following example will fail at
compilation time.
---------------------------
# model
model {
a ~ dnorm(10,1)
b ~ dnorm(10,1)
c <- a * b
}
# data
list(c=75)
---------------------------
However it is can be useful to solve such models. Is the trick of a very
concentrated variate around the logical node valuable? (see an example
below), or do exist other possibilities?
---------------------------
# model
model {
a ~ dnorm(10,1)
b ~ dnorm(10,1)
c <- a * b
d ~ dnorm(c,100000)
}
# data
list(d=75)
---------------------------
Thanks for your help.
Jean-Baptiste DENIS
mél : [log in to unmask] unité de recherche BIA
tél : ##.33.(0)1.34.65.22.23 INRA
fax : ##.33.(0)1.34.65.22.17 F-78352 Jouy-en-Josas
------------< http://www.inra.fr/bia/J/ >--------------
-------------------------------------------------------------------
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
|