Hi All,
I have written a c-code for for generating samples from joint posterior density using Gibbs sampling technique. I have used Adaptive Rejection sampling (Wild & Gilks (1993)) to generate
samples from the univariate conditional density. However, the code is working fine for some cases,
for some cases it is not working. Is anybody aware of the source code for ARS? I can share my code if you can help me to find out the bug in the code.
Regards,
swarup.
Post the relevent code sections... also whatever warnings and errors your compiler is coming up with...
Is this a "produces the wrong answer" kind of doesn't work, or "crashes with segfault" kind of doesn't work?
Plus, looking at someone else's code (assuming there is any to find) won't tell you much about how you messed up your code.
If you dance barefoot on the broken glass of undefined behaviour, you've got to expect the occasional cut.
If at first you don't succeed, try writing your phone number on the exam paper.
I love these kinds of thread starters. Analogous perhaps to:
"Doctor, I have a pain - somewhere, sometimes. Can you tell me what's wrong with me, over the phone?"
Sure thing, let me just check with the Gypsy lady down the block here. . . < lol >
Originally Posted by Adak
LOL...
A guy walks into a doctor's office, hooks his thumb over the edge of the desk and whacks it with a hammer...
Guy: "Doc, it hurts like he** every time I do this, can you help me?"
Doc: "Stop doing that."
You make it sound like there's only one implementation on the planet.
Perhaps you mean the algorithm?
We need a few things to get anywhere at all:
1. A description of what "not working" means. We have to know what on earth we are looking for.
2. Some source code containing the bug.
2. A nice-to-have: List of the compiler warnings that the code produces. Whilst we'll look over the code, most of us won't bother to actually compile it.
My homepage
Advice: Take only as directed - If symptoms persist, please see your debugger
Linus Torvalds: "But it clearly is the only right way. The fact that everybody else does it some other way only means that they are wrong"
3 Minutes Googling got me this link
Adaptive Rejection Sampling
Tim S.
Googling, what's that?
Warning: Some or all of my posted code may be non-standard and as such should not be used and in no case looked at.
