Description Usage Arguments Value Examples
Extract Gelman-Rubin diagnostics as data frame
1 | GelmanDiag(obj, iter)
|
obj |
model object from |
iter |
integer. The number of iterations. |
data frame
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 | ## Not run:
Load packages
library(zoib)
library(zibHelpers)
# Run example of clustered zero-inflated beta regression from
# Liu and Kong (under review, 12-14)
# Load data
data("AlcoholUse", package = "zoib")
AlcoholUse$Grade <- as.factor(AlcoholUse$Grade)
nIter = 50 # Number of iterations, including burn-in
# Estimate
Out <- zoib(Percentage ~ Grade + Days + Gender|1|Grade + Days + Gender|1,
data = AlcoholUse, random = 1, EUID = AlcoholUse$County,
zero.inflation = TRUE, one.inflation = FALSE, joint = FALSE,
n.iter = nIter)
# Gelman-Rubin diagnostics
GelmanDiag(Out, nIter)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.