
LWA-style Marginal Cox Combined-Likelihood Inference for KK Designs
Source:R/inference_survival_KK_lwa_cox.R
InferenceSurvivalKKLWACoxPHOneLik.RdFits a single combined Cox partial likelihood
\(\lambda(t \mid w, x) = \lambda_0(t)\exp(\beta_T w + \beta_X^\top x)\)
jointly over matched-pair and reservoir subjects for KK matching-on-the-fly
designs with survival responses (a marginal, not stratified, Cox model:
matched pairs do not get pair-specific baseline hazards). This is the
one-likelihood combined-fit analog of
InferenceSurvivalKKLWACoxPHIVWC,
which instead fits and pools two separate estimators. likelihood_tier
= "partial": exposes likelihood-ratio and parametric-likelihood-bootstrap
inference in addition to Wald/asymptotic and Bayesian-bootstrap paths.
References
Cox, D. R. (1972). "Regression Models and Life-Tables." Journal of the Royal Statistical Society, Series B, 34(2), 187-220.
Super class
Inference -> InferenceSurvivalKKLWACoxPHOneLik
Methods
Public methods
InferenceSurvivalKKLWACoxPHOneLik$set_custom_randomization_statistic_function()InferenceSurvivalKKLWACoxPHOneLik$set_custom_randomization_statistic_cpp()InferenceSurvivalKKLWACoxPHOneLik$approximate_randomization_distribution_beta_hat_T()InferenceSurvivalKKLWACoxPHOneLik$supports_rand_pval_for_incidence()InferenceSurvivalKKLWACoxPHOneLik$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 = DesignSeqOneByOneKK14$new(n = 10, response_type = 'survival')
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 = InferenceSurvivalKKLWACoxPHOneLik$new(seq_des)
inf$compute_estimate()
#> [1] -0.3242638
# }