Description Usage Arguments Value Author(s) Examples
Function to produce plot based on different diagnostic statistics
1 |
x |
a |
diag |
diagnostic statistics to be used—either the Gelman-Rubin statistic (Rhat) by default or effective sample size (n.eff) |
A plot showing the relevant diagnostic stats for each node in the model
Tao Ding Gianluca Baio
1 2 3 4 5 6 7 8 9 10 11 12 13 14 | ### Read and format the data (binary)
data = read.csv(url("http://www.statistica.it/gianluca/bmeta/Data-bin.csv"))
### List data for binary outcome
data.list <- list(y0=data$y0,y1=data$y1,n0=data$n0,n1=data$n1)
### generate output using bmeta
x <- bmeta(data=data.list,outcome="bin",model="std.norm",type="fix")
### run the diagnostic plot to examine the Gelman-Rubin statistic
diag.plot(x)
### run the diagnostic plot to examine the effective sample size
diag.plot(x,diag="n.eff")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.