
GEE Inference for KK Designs with Proportion Response
Source:R/inference_proportion_KK_combined.R
InferencePropKKGEE.RdFits a Generalized Estimating Equations model with a binomial (quasi-likelihood,
fractional-response) family and logit link, \(\mathrm{logit}\,E[Y_i \mid x_i]
= x_i^\top\beta\), for proportion (continuous values in (0, 1)) responses
under a KK matching-on-the-fly design — the same fractional-logit mean-model
idea as InferencePropFractionalLogit,
extended to jointly account for matched-pair and reservoir clustering via GEE.
Each GEE cluster is either a matched pair (2 members) or a reservoir singleton
(1 member), with an exchangeable working correlation structure — see
$compute_estimate()'s method-level documentation for the full fitting
contract (internal Rcpp solver vs. geepack fallback, hardening/retry
behavior). Inference is quasi-likelihood/estimating-equation based
(likelihood_tier = "quasi"): standard errors are GEE sandwich (robust)
standard errors, not model-likelihood-based.
References
Liang, K.-Y., and Zeger, S. L. (1986). "Longitudinal Data Analysis Using Generalized Linear Models." Biometrika, 73(1), 13-22, doi:10.1093/biomet/73.1.13 , for the GEE estimating-equation framework and sandwich variance estimator used here.
Super class
Inference -> InferencePropKKGEE
Methods
+ 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$is_nonestimable()Inference$set_optimization_alg()Inference$set_seed()Inference$supports()
InferencePropKKGEE$new()
Initialize KK proportion-response GEE inference, validate
the matched/reservoir design, and prepare the exchangeable-working-correlation
fractional-logit GEE fitting machinery used by
InferencePropKKGEE.
Usage
InferencePropKKGEE$new(
des_obj,
model_formula = NULL,
use_rcpp = TRUE,
verbose = FALSE,
smart_cold_start_default = NULL
)Arguments
des_objA completed
Designobject with a proportion response.model_formulaOptional formula for covariate adjustment.
use_rcppWhether to use the internal Rcpp GEE solver (
TRUE, default) with automatic fallback togeepack::geeglmon failure, or always usegeepack::geeglmdirectly (FALSE, requires geepack to be installed).verboseWhether to print progress messages.
smart_cold_start_defaultWhether to use smart cold start values.
Examples
# \donttest{
seq_des = DesignSeqOneByOneKK14$new(n = 10, response_type = 'proportion')
for (i in 1:10) {
seq_des$add_one_subject_to_experiment_and_assign(data.frame(x1 = rnorm(1), x2 = rnorm(1)))
}
seq_des$add_all_subject_responses(runif(10))
inf = InferencePropKKGEE$new(seq_des)
inf$compute_estimate()
#> [1] 0.1663989
# }