Print

Print


Just to follow up, Martyn helped answer my question. This works:

model {
r ~ dunif(0, 5)
F(x) <- r*exp(-r*x)
z <- integral(F(x), 0, 1, 0.001)
}

Apparently the function must be named F.

Thanks Martyn,
Richard




On Fri, Oct 1, 2010 at 11:20 AM, Martyn Plummer <[log in to unmask]> wrote:
> Sorry, there are no functionals in JAGS.
>
> Richard, I think your problem is that the first argument needs to be a
> logical node.
>
> Martyn
>
> On Fri, 2010-10-01 at 10:09 -0400, Krzysztof Sakrejda-Leavitt wrote:
>> I am also interested in whether this can be done in JAGS, if anyone has
>> insight on that.
>>
>> Krzysztof
>>
>> On 10/01/2010 09:35 AM, Richard Chandler wrote:
>> > Hi,
>> >
>> > Is it possible to integrate over intervals of pdfs using the new integral function in OpenBUGS? I thought I would be able to do something like:
>> >
>> > model {
>> >
>> > rate ~ dunif(0, 1)
>> >
>> > z <- integral(dexp(rate), 0, 1, 0.001)
>> >
>> > }
>> >
>> > But I get the error message: "Sorry something went wrong in procedure Node.Parents in module GraphFunctional"
>> >
>> > Any help would be greatly appreciated.
>> >
>> > I'm using OpenBUGS 3.1.1 rev 524 on Windows.
>> >
>> >
>> > Thanks,
>> > Richard
>> >
>> > -------------------------------------------------------------------
>> > 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
>> >
>>
>
>
> -----------------------------------------------------------------------
> This message and its attachments are strictly confidential. If you are
> not the intended recipient of this message, please immediately notify
> the sender and delete it. Since its integrity cannot be guaranteed,
> its content cannot involve the sender's responsibility. Any misuse,
> any disclosure or publication of its content, either whole or partial,
> is prohibited, exception made of formally approved use
> -----------------------------------------------------------------------
>

-------------------------------------------------------------------
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