R help - Mar 2009 - (no subject)

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?

 

Many thanks

 

Alice

 

 

_________________________________________________________________
[[elided Hotmail spam]]

	[[alternative HTML version deleted]]

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 at 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.