
Adjacent Category Logit Regression Inference for Ordinal Responses
Source:R/inference_ordinal_adj_cat_logit.R
InferenceOrdinalAdjCatLogitRegr.RdFits an adjacent-category logit regression for ordinal responses (via
fast_adjacent_category_logit_cpp — see that page for the full
model, an alternative ordinal parameterization to the cumulative-logit
proportional-odds model) using the treatment indicator and, optionally, all
recorded covariates as predictors. This is a full-likelihood class
(likelihood_tier = "full") supporting score, gradient, and
likelihood-ratio tests, plus parametric likelihood-ratio bootstrap
calibration, in addition to Wald and resampling-based inference.
Bayesian-bootstrap inference is temporarily unavailable because the current
non-uniform weighted hook fits a cumulative-logit surrogate rather than the
adjacent-category likelihood. It will remain disabled until the native
weighted adjacent-category backend described in the package implementation
plan lands.
Super class
Inference -> InferenceOrdinalAdjCatLogitRegr
Methods
Public methods
InferenceOrdinalAdjCatLogitRegr$set_custom_randomization_statistic_function()InferenceOrdinalAdjCatLogitRegr$set_custom_randomization_statistic_cpp()InferenceOrdinalAdjCatLogitRegr$approximate_randomization_distribution_beta_hat_T()InferenceOrdinalAdjCatLogitRegr$supports_rand_pval_for_incidence()InferenceOrdinalAdjCatLogitRegr$compute_rand_two_sided_pval()
+ inherited public methods from Inference
Inference$capabilities()Inference$compute_asymp_confidence_interval()Inference$compute_asymp_two_sided_pval()Inference$compute_estimate()Inference$compute_exact_confidence_interval()Inference$compute_exact_two_sided_pval_for_treatment_effect()Inference$duplicate()Inference$get_analysis_data()Inference$get_covariates()Inference$get_design_object()Inference$get_model_formula()Inference$get_nonestimable_reason()Inference$get_nonestimable_stage()Inference$get_optimization_alg()Inference$get_response()Inference$get_response_type()Inference$get_treatment()Inference$initialize()Inference$is_nonestimable()Inference$set_optimization_alg()Inference$set_seed()Inference$supports()
Examples
# \donttest{
seq_des = DesignSeqOneByOneBernoulli$new(n = 10, response_type = 'ordinal')
for (i in 1:10) {
seq_des$add_one_subject_to_experiment_and_assign(data.frame(x1 = rnorm(1)))
}
seq_des$add_all_subject_responses(sample(1:4, 10, replace = TRUE))
inf = InferenceOrdinalAdjCatLogitRegr$new(seq_des)
inf$compute_estimate()
#> [1] -0.9086307
# }