Nothing
## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>",
fig.width = 6,
fig.height = 6
)
## ----echo = TRUE--------------------------------------------------------------
library("brglm2")
data("stemcell", package = "brglm2")
stem <- within(stemcell, religion <- as.numeric(religion))
## ----echo = TRUE--------------------------------------------------------------
stem_formula <- research ~ religion + gender
stemcells_ml <- bracl(stem_formula, weights = frequency, data = stem,
parallel = TRUE, type = "ML")
summary(stemcells_ml)
## ----echo = TRUE--------------------------------------------------------------
class(stemcells_ml)
## ----echo = TRUE--------------------------------------------------------------
stemcells_ml_full <- bracl(stem_formula, weights = frequency, data = stemcell,
parallel = FALSE, type = "ML")
summary(stemcells_ml_full)
## ----echo = TRUE--------------------------------------------------------------
(lrt <- deviance(stemcells_ml) - deviance(stemcells_ml_full))
## ----echo = TRUE--------------------------------------------------------------
(df1 <- df.residual(stemcells_ml) - df.residual(stemcells_ml_full))
## ----echo = TRUE--------------------------------------------------------------
pchisq(lrt, df1, lower.tail = FALSE)
## ----echo = TRUE--------------------------------------------------------------
summary(update(stemcells_ml, type = "AS_mean"))
## ----echo = TRUE--------------------------------------------------------------
summary(update(stemcells_ml, type = "AS_median"))
## ----echo = TRUE--------------------------------------------------------------
predict(stemcells_ml, type = "probs")
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.