Skip to contents

Fits 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


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_obj

A completed Design object with a proportion response.

model_formula

Optional formula for covariate adjustment.

use_rcpp

Whether to use the internal Rcpp GEE solver (TRUE, default) with automatic fallback to geepack::geeglm on failure, or always use geepack::geeglm directly (FALSE, requires geepack to be installed).

verbose

Whether to print progress messages.

smart_cold_start_default

Whether to use smart cold start values.


InferencePropKKGEE$clone()

The objects of this class are cloneable with this method.

Usage

InferencePropKKGEE$clone(deep = FALSE)

Arguments

deep

Whether to make a deep clone.

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
# }