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 me what is the way around of this error? Thanks in advance, Best regards, Satrajit On Sun, Aug 16, 2009 at 9:55 PM, SATRAJIT ROYCHOUDHURY <satrirqa at gmail.com<mailto:satrirqa at gmail.com>> wrote: Thanks a lot! On Sun, Aug 16, 2009 at 9:45 PM, Ben Bolker <bolker at ufl.edu<mailto:bolker at ufl.edu>> wrote: 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<mailto: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<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.
______________________________________________ R-help at r-project.org<mailto: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<http://www.r-project.org/posting-guide.html> and provide commented, minimal, self-contained, reproducible code.
Ben Bolker Associate professor, Biology Dep't, Univ. of Florida bolker at ufl.edu / www.zoology.ufl.edu/bolker GPG key: www.zoology.ufl.edu/bolker/benbolker-publickey.asc -------------- next part -------------- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 261 bytes Desc: OpenPGP digital signature URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20090817/45938e0b/attachment-0002.bin>