wglm: Logistic Regression Using IPCW

In tagteam/riskRegression: Risk Regression Models and Prediction Scores for Survival Analysis with Competing Risks

Logistic Regression Using IPCW

Description

Logistic regression over multiple timepoints where right-censoring is handled using inverse probability of censoring weighting (IPCW).

Usage

wglm(
  formula.event,
  times,
  data,
  formula.censor = ~1,
  cause = NA,
  fitter = NULL,
  ties = NULL,
  product.limit = NULL,
  iid = FALSE,
  se = TRUE,
  store = NULL
)

Arguments

formula.event	[formula] a formula with a Surv object on the left hand side and the covariates for the logistic regression on the right hand side.
times	[numeric vector] time points at which to model the probability of experiencing an event.
data	[data.frame] dataset containing the time at which the event occured, the type of event, and regressors used to fit the censoring and logistic models.
formula.censor	[formula] an optional formula indicating on its right hand side the covariates for the censoring model.
cause	[character or numeric] the cause of interest. Defaults to the first cause.
fitter	[character] routine to fit the Cox regression models.
ties	[character] method used to handle ties when using a Cox model ("breslow" or "efron"). Ignored if fitter equals to "prodlim".
product.limit	[logical] if TRUE the survival is computed using the product limit estimator.
iid	[logical] should the influence function of the logistic regression parameters be computed, accounting for the uncertainty of the weights. This can be computationally and memory intensive.
se	[logical] should the variance-covariance matrix of the logistic regression parameters be stored, accounting for the uncertainty of the weights. This can be computationally and memory intensive.
store	[vector] when evaluating the iid, should prediction be only computed for unique covariate sets and mapped back to the original dataset (data="minimal"). Otherwise use data="full".

Details

First, a Cox model is fitted (argument formula.censor) and the censoring probabilities are computed relative to each timepoint (argument times) to obtain the censoring weights. Then, for each timepoint, a logistic regression is fitted with the appropriate censoring weights and where the outcome is the indicator of having experience the event of interest (argument cause) at or before the timepoint.

Value

an object of class "wglm".

See Also

coef.wglm to output the estimated parameters from the logistic regression.
confint.wglm to output the estimated parameters from the logistic regression with their confidence interval.
model.tables.wglm to output a data.frame containing the estimated parameters from the logistic regression with its confidence intervals and p-values.
predictRisk.wglm to evaluate event probabilities (e.g. survival probabilities) conditional on covariates.
summary.wglm for displaying in the console a summary of the results.
weights.wglm to extract the IPCW weights.

Examples

library(survival)

#### simulate data ####
set.seed(10)
n <- 250
tau <- 1:5
d <- sampleData(n, outcome = "competing.risks")
dFull <- d[event!=0] ## (artificially) remove censoring
dSurv <- d[event!=2] ## (artificially) remove competing risk

#### no censoring ####
e.wglm <- wglm(Surv(time,event) ~ X1, 
               times = tau, data = dFull, product.limit = TRUE)
e.wglm ## same as a logistic regression at each timepoint

coef(e.wglm)
confint(e.wglm)
model.tables(e.wglm)

summary(ate(e.wglm, data = dFull, times = tau, treatment = "X1", verbose = FALSE))

#### right-censoring ####
## no covariante in the censoring model (independent censoring)
eC.wglm <- wglm(Surv(time,event) ~ X1,
               times = tau, data = dSurv, product.limit = TRUE)
summary(eC.wglm)

weights(eC.wglm)

## with covariates in the censoring model
eC2.wglm <- wglm(Surv(time,event) ~ X1 + X8, formula.censor = ~ X1*X8,
                 times = tau, data = dSurv)
eC2.wglm

#### Competing risks ####
## here Kaplan-Meier as censoring model
eCR.wglm <- wglm(Surv(time,event) ~ X1, formula.censor = ~X1,
                 times = tau, data = d)
eCR.wglm
summary(eCR.wglm)
eCR.wglm <- wglm(Surv(time,event) ~ X1, formula.censor = ~X1,
                 times = tau, data = d)

tagteam/riskRegression documentation built on April 13, 2025, 9:28 a.m.