Re: [R] Question regarding finding credible interval using r2winbugs

[cc'd to list] You may need coda::HPDinterval(as.mcmc(fit),prob=0.9) SATRAJIT ROYCHOUDHURY wrote: > Dear Ben, > > coda::HPDinterval(fit,prob=0.9) > > This command is keep giving me error message as; > Error in HPDinterval(fit) : no applicable method for "HPDinterval" > > Can you pls tell

Re: [R] Question regarding finding credible interval using r2winbugs

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

[R] Question regarding finding credible interval using r2winbugs

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.