Hi to all!
I intend to use a OpenBUGS as a realization of Metropolis-Hasting algorithm for my projects in bioinformatics. I have a problem formulation as a task of finding maximum of join probability function, what has been decomposed into product of independent of each other distributions.
Example: P(Z,B,A,G|W)=P(A|B,W)P(G|A,B)P(Z|B)P(B|W)
max_{B,A,G} P(Z,B,A,G|W)-?
The variables have a complex nature and could be represented as objects and not as simple data types (integers, real etc.). Some PDFs from I need to sample, are also complex (where isn’t a analytical expression for them). For example, I need to sample a one variable "B" (the string of states) from the - a Hidden Markov Model. I need a way to define my own PDFs and sampling algorithms in BUGS. I would like to use a mechanism of writing BUGS Extensions as mentioned in Developer Manual. Can I define on Component Pascal abstract classes derived from classes in Graph and Updater class hierarchies for my variables and functions and write realizations of them on C++? May realization be compiled in DLL and linked with CP codes? I don’t wand to write all project on CP wholly because no one from my team so I am knows it. May any one know a better way to solve such problems? I will be very appreciated if you give me some examples of simple projects with user-defined stochastic functions and datatypes for BUGS. Is there a more detailed manual how to use Graph and Updater class hierarchies?
Sincerely, Alexey Neverov.
-------------------------------------------------------------------
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
|