 |
 |
You most likely did two different things without noticing. You must supply the full code used in both cases. Bendix Carstensen > -----Original Message----- > From: (The BUGS software mailing list) [mailto:[log in to unmask]] On > Behalf Of Everton B. da Rocha > Sent: 4. oktober 2013 20:06 > To: [log in to unmask] > Subject: [BUGS] Poisson-Normal regression model > > Dear all, > > I'm fitting a regression model poisson-normal. I made the adjustment using > the OpenBUGS then made the adjustment by R using the interface between R > and OpenBUGS through the library BRugs by the function BRugsFit. In both > cases I used the same model specification, i.e., initial values, thin, > burn-in, number of chains, prior distributions, however the results are > different. > > The output of the two adjustments are as follows: > > #------------------------------------------------------------------------- > - > OpenBUGS output > > mean sd MC_error val2.5pc median val97.5pc > start sample > D 1.211 0.1996 0.00252 0.8799 1.19 1.656 5001 > 10000 > int0 0.835 0.1532 0.009751 0.529 0.8325 1.134 5001 > 10000 > int1 0.6538 0.1659 0.01105 0.3185 0.6535 > 0.9736 5001 10000 > slope0 -0.01431 0.004407 8.038E-5 -0.02305 -0.01433 > -0.005729 5001 10000 > slope1 -0.0121 0.004291 8.08E-5 -0.0205 -0.01212 > -0.003624 5001 10000 > tau 0.8479 0.1366 0.001715 0.6038 0.8405 > 1.137 5001 10000 > > #------------------------------------------------------------------------- > - > R output > > $Stats > mean sd MC_error val2.5pc median val97.5pc start > sample > D 0.18660 0.033020 3.561e-04 0.13200 0.18300 0.26050 5001 > 10000 > int0 1.19900 0.075700 1.929e-03 1.04900 1.20000 1.34500 5001 > 10000 > int1 1.38100 0.077490 2.523e-03 1.22700 1.38000 1.53300 5001 > 10000 > slope0 -0.02381 0.003932 6.226e-05 -0.03146 -0.02384 -0.01612 5001 > 10000 > slope1 -0.02527 0.004123 7.525e-05 -0.03339 -0.02531 -0.01709 5001 > 10000 > tau 5.52300 0.960900 1.027e-02 3.83900 5.46500 7.57300 5001 > 10000 > > $DIC > Dbar Dhat DIC pD > y 10780 10700 10860 81.33 > total 10780 10700 10860 81.33 > > #------------------------------------------------------------------------- > - > > Does anyone have any idea what is causing the discrepancy in results? > > Best regards > > Everton B. da Rocha > > ------------------------------------------------------------------- > This list is for discussion of modelling issues and the BUGS software. > For help with crashes and error messages, first mail bugs@mrc- > bsu.cam.ac.uk 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 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