Shuying Yang wrote:
Dear Members,
I have a question about using R2WinBUGS to obtain the WinBUGS results.
By default, when R2WinBUGS returns summary stats, I got mean, sd, 2.5%, 25%,
median, 75% and 97.5%. Could anyone tell me how to modify the code to obtain
5% and 95% summary results?
Hi,
please read the posting guide and use a sensible subject line (I almost
missed your posting).
The used quantiles in the printed output are hardcoded, but you can
calculate arbitrary quantiles from the data.
Ad an example, please consider the schools example from ?bugs, running
the example results in an object called "schools.sim". Now you can
calculate arbitrary quantiles as e.g. in:
apply(schools.sim$sims.matrix, 2, quantile, c(0.05, 0.95))
Uwe Ligges
Many thanks
Alice
_________________________________________________________________
[[elided Hotmail spam]]
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.