View source: R/ordinal_superiority.R
ordinal_superiority.bracl | R Documentation |
ordinal_superiority()
is a method for the estimation and
inference about model-based ordinal superiority scores introduced
in Agresti and Kateri (2017, Section 5) from fitted objects. The
mean bias of the estimates of the ordinal superiority scores can be
corrected.
## S3 method for class 'bracl' ordinal_superiority( object, formula, data, measure = c("gamma", "Delta"), level = 0.95, bc = FALSE )
object |
a fitted object from an ordinal regression
model. Currently only models from class |
formula |
a RHS formula indicating the group variable to use. |
data |
an optional data frame in which to look for variables
with which to compute ordinal superiority measures. If
omitted, an attempt is made to use the data that produced
|
measure |
either |
level |
the confidence level required when computing confidence intervals for the ordinal superiority measures. |
bc |
logical. If |
Agresti, A., Kateri, M. (2017). Ordinal probability effect measures for group comparisons in multinomial cumulative link models. Biometrics, 73 214-219. doi: 10.1111/biom.12565.
data("stemcell", package = "brglm2") # Adjacent category logit (proportional odds) stem <- within(stemcell, {nreligion = as.numeric(religion)}) fit_bracl_p <- bracl(research ~ nreligion + gender, weights = frequency, data = stem, type = "ML", parallel = TRUE) # Estimates and 95% confidence intervals for the probabilities that the response # category for gender "female" is higher than the response category for gender "male", # while adjusting for religion. ordinal_superiority(fit_bracl_p, ~ gender) ## Not run: # And their (very-similar in value here) bias corrected versions # with 99% CIs ordinal_superiority(fit_bracl_p, ~ gender, bc = TRUE, level = 0.99) # Note that the object is refitted with type = "AS_mean" ## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.