Hi,
We've been working with some model code which is getting very long, but
it is pretty modular. I can see a clutzy way of managing this by
keeping separate text files and pasting them together line-by-line with
an R package where you could select components, but that won't keep the
problem at bay for long.
For example we might have a model where:
z[i,t] <- dbern(z[i,t-1] * phi[i,t])
phi[i,t] is then handled separately depending on what type of model is
intended. For example it might just be a single term:
phi[i,t] <- interceptPhi
Or it might be a linear model, spline, etc...
My question is whether anyone has experience managing complexity with
BUGS/JAGS code the way you would with other programing languages. I
know there are some references in JAGS to user-created functions but I
believe this means implemented in C++ (?).
My apologies if this is a little off-topic.
Krzysztof
-----------------------------------------------
Krzysztof Sakrejda-Leavitt
Organismic and Evolutionary Biology
University of Massachusetts, Amherst
319 Morrill Science Center South
611 N. Pleasant Street
Amherst, MA 01003
work #: 413-325-6555
email: [log in to unmask]
-----------------------------------------------
-------------------------------------------------------------------
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
|