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
|