R help - Aug 2009 - 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.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]]

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.

[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>