influences.CumBH {CaseCohortCoxSurvival}R Documentation

influences.CumBH

Description

Computes the influences on the log-relative hazard, baseline hazards at each unique event time, and on the cumulative baseline hazard in a given time interval [Tau1, Tau2]. Can take calibration of the design weights into account.

Usage

influences.CumBH(mod, Tau1 = NULL, Tau2 = NULL, A=NULL, calibrated = NULL)

Arguments

mod

a cox model object, result of function coxph.

Tau1

left bound of the time interval considered for the cumulative baseline hazard and pure risk. Default is the first event time.

Tau2

right bound of the time interval considered for the cumulative baseline hazard and pure risk. Default is the last event time.

A

n \times q matrix with the values of the auxiliary variables used for the calibration of the weights in the whole cohort. Needs to be provided if calibrated = TRUE.

calibrated

are calibrated weights used for the estimation of the parameters? If calibrated = TRUE, the argument below needs to be provided. Default is FALSE.

Details

influences.CumBH works for estimation from a case-cohort with design weights or calibrated weights (case-cohort consisting of the subcohort and cases not in the subcohort, i.e., case-cohort obtained from two phases of sampling).

If covariate information is missing for certain individuals in the phase-two data (i.e., case-cohort obtained from three phases of sampling), use influences.CumBH.missingdata.

influences.CumBH uses the influence formulas provided in Etievant and Gail (2024).

If calibrated = FALSE, the infuences are only provided for the individuals in the case-cohort. If calibrated = TRUE, the influences are provided for all the individuals in the cohort.

Value

infl.beta: matrix with the overall influences on the log-relative hazard estimates.

infl.lambda0.t: matrix with the overall influences on the baseline hazards estimates at each unique event time.

infl.Lambda0.Tau1Tau2.hat: vector with the overall influences on the cumulative baseline hazard estimate in [Tau1, Tau2].

infl2.beta: matrix with the phase-two influences on the log-relative hazard estimates. Returned if calibrated = TRUE.

infl2.lambda0.t: matrix with the phase-two influences on the baseline hazards estimates at each unique event time. Returned if calibrated = TRUE.

infl2.Lambda0.Tau1Tau2.hat: vector with the phase-two influences on the cumulative baseline hazard estimate in [Tau1, Tau2]. Returned if calibrated = TRUE.

beta.hat: vector of length p with log-relative hazard estimates.

lambda0.t.hat: vector with baseline hazards estimates at each unique event time.

Lambda0.Tau1Tau2.hat: cumulative baseline hazard estimate in [Tau1, Tau2].

References

Etievant, L., Gail, M. H. (2024). Cox model inference for relative hazard and pure risk from stratified weight-calibrated case-cohort data. Lifetime Data Analysis, 30, 572-599.

See Also

estimation, estimation.CumBH, estimation.PR, influences, influences.RH, influences.PR, influences.missingdata, influences.RH.missingdata, influences.CumBH.missingdata,
influences.PR.missingdata, robustvariance and variance.

Examples


  data(dataexample.stratified, package="CaseCohortCoxSurvival")
  cohort <- dataexample.stratified$cohort
  casecohort <- cohort[which(cohort$status == 1 |
                       cohort$subcohort == 1),] # the stratified case-cohort
  casecohort$weights <- casecohort$strata.n / casecohort$strata.m
  casecohort$weights[which(casecohort$status == 1)] <- 1

  Tau1 <- 0
  Tau2 <- 8
  x <- c(-1, 1, -0.6) # given covariate profile for the pure risk

  # Estimation using the stratified case cohort with design weights
  mod <- coxph(Surv(event.time, status) ~ X1 + X2 + X3, data = casecohort,
               weight = weights, id = id, robust = TRUE)
  est <- influences(mod, Tau1 = Tau1, Tau2 = Tau2, x = x)

  # print the influences on the cumulative baseline hazard estimate
  # est$infl.Lambda0.Tau1Tau2


[Package CaseCohortCoxSurvival version 0.0.36 Index]