[R] Question regarding finding credible interval using r2winbugs

Ben Bolker bolker at ufl.edu
Mon Aug 17 03:45:38 CEST 2009



coda::HPDinterval(fit,prob=0.9)


SATRAJIT ROYCHOUDHURY-2 wrote:
> 
> Dear
> 
>   I am trying to find a 90% credible interval. I am using the following
> code.
>    fit<-bugs(
>    model.file=BUGScode,
>    data=data,
>    inits = list(geninits1,geninits2),
>    parameters.to.save=keepers,
>    n.chains=nchains,
>    n.iter=runs,
>    n.burnin=burn,
>    n.thin=nthin,
>    DIC= TRUE,
>    bugs.directory="C:/Program Files/WINBUGS.14",
> 
> 
> \
>    )
> 
>    But this is only giving 95% credible interval. How can I get 90%, 80%
>    credible intervals using r2winbugs?
> 
>    Pls help.
> 
>    Waiting for your response.
> 
>    Satrajit
> 
> 	[[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.
> 
> 

-- 
View this message in context: http://www.nabble.com/Question-regarding-finding-credible-interval-using-r2winbugs-tp24998594p24999931.html
Sent from the R help mailing list archive at Nabble.com.




More information about the R-help mailing list