cif_aiptw_pseudo: Augmented Inverse Probability of Treatment Weighted CIFs...

In RobinDenz1/adjustedCurves: Confounder-Adjusted Survival Curves and Cumulative Incidence Functions

Augmented Inverse Probability of Treatment Weighted CIFs using Pseudo-Values

Description

This page explains the details of estimating augmented inverse probability of treatment weighted CIFs using pseudo-values in a competing risks setting (method="aiptw_pseudo" in the adjustedcif function). All regular arguments of the adjustedcif function can be used. Additionally, the outcome_vars argument and the treatment_model argument have to be specified in the adjustedcif call. Further arguments specific to this method are listed below.

Arguments

outcome_vars	[required] A character vector of column names specifying variables to be used when modeling the outcome mechanism using geese. See details and examples.
treatment_model	[required] Must be a glm or multinom model object with variable as response variable. See details and examples.
type_time	A character string specifying how the time should be modeled. Possible values are "factor" (modeling each point in time as a separate variable, the default), "bs" (modeling time using B-Splines) or "ns" (modeling time using natural splines).
spline_df	The number of degrees of freedom used for the natural-spline or B-spline function. Defaults to 5. Ignored if type_time="factor".

Details

• Type of Adjustment: Requires a treatment assignment model (glm or multinom) and a character vector of variable names used to model the outcome mechanism (internally uses geese). In contrast to the "aiptw" method this function does not allow for dependent censoring.

• Doubly-Robust: Estimates are Doubly-Robust.

• Categorical groups: Any number of levels in variable are allowed. Must be a factor variable.

• Approximate Variance: Calculations to approximate the variance and confidence intervals are available.

• Allowed Time Values: Allows both continuous and integer time.

• Bounded Estimates: Estimates are not guaranteed to be bounded in the 0 to 1 probability range.

• Monotone Function: Estimates are not guaranteed to be monotone.

• Dependencies: This method relies on the geepack and prodlim packages.

Instead of only modeling the outcome mechanism or the treatment assignment mechanism, both kind of models are required to use this method. If either of those models are correctly specified, unbiased estimates will be obtained. In contrast to the "aiptw" method, the "aiptw_pseudo" method uses a generalized estimation equation (geese) approach to model the outcome mechanism. The model is fit in the same way as described in the "direct_pseudo" method. Those Direct Standardization based estimates are then transformed using the previously estimated propensity score. This results in the doubly-robust property of the method. More information on this particular method can be found in the original article by Wang (2018). The original article only deals with survival probabilities without competing risks, but the only difference to the CIF estimation with competing risks is the calculation of the pseudo-values. More information on Pseudo-Values is available in Andersen et al. (2017) and Andersen and Perme (2010).

When estimating the geese model the ev_time variable is used as a factor by default. This results in one coefficient being estimated for each unique point in time, which can be very slow computationally if there are a lot of unique points in time and/or the dataset has many rows. In these cases it is recommended to use type_time="bs" or type_time="ns", which results in the ev_time being modeled using B-Splines or Natural Splines. Simulation studies indicate that there is little difference in the estimates when an appropriately large number of spline_df is used.

Value

Adds the following additional objects to the output of the adjustedcif function:

• pseudo_values: The matrix of estimated pseudo-values.

• geese_model: The geese model used to make the predictions.

Author(s)

Jixian Wang supplied the R source code used in the original article, which was used by Robin Denz to create a generalized version of this method with additional functionality and improved performance.

References

Jixian Wang (2018). "A Simple, Doubly Robust, Efficient Estimator for Survival Functions Using Pseudo Observations". In: Pharmaceutical Statistics 17.38-48

James M. Robins and Andrea Rotnitzky (1992). "Recovery of Information and Adjustment for Dependent Censoring Using Surrogate Markers". In: AIDS Epidemiology: Methodological Issues. Ed. by Nicholas P. Jewell, Klaus Dietz, and Vernon T. Farewell. New York: Springer Science + Business Media, pp. 297-331

Per Kragh Andersen, Elisavet Syriopoulou, and Erik T. Parner (2017). "Causal Inference in Survival Analysis using Pseudo-Observations". In: Statistics in Medicine 36, pp. 2669-2681

Per Kragh Andersen and Maja Pohar Perme (2010). "Pseudo-Observations in Survival Analysis". In: Statistical Methods in Medical Research 19, pp. 71-99

Aris Perperoglou, Willi Sauerbrei, Michal Abrahamowicz, and Matthias Schmid (2019). "A Review of Spline Function Procedures in R". in: BMC Medical Research Methodology 19.46, pp. 1-16

See Also

geese, jackknife, ns, bs

Examples

library(adjustedCurves)

if (requireNamespace("geepack") & requireNamespace("prodlim")) {

set.seed(42)

# simulate some data as example
sim_dat <- sim_confounded_crisk(n=30, max_t=5)
sim_dat$group <- as.factor(sim_dat$group)

# estimate a treatment assignment model
glm_mod <- glm(group ~ x1 + x3 + x5 + x6, data=sim_dat, family="binomial")

# use it + pseudo values + geese model to calculate adjusted CIFs
adjcif <- adjustedcif(data=sim_dat,
                      variable="group",
                      ev_time="time",
                      event="event",
                      cause=1,
                      method="aiptw_pseudo",
                      outcome_vars=c("x1", "x2", "x3", "x4", "x5", "x6"),
                      treatment_model=glm_mod,
                      conf_int=FALSE,
                      force_bounds=TRUE,
                      iso_reg=TRUE)

# plot the curves
plot(adjcif)
}

RobinDenz1/adjustedCurves documentation built on Sept. 27, 2024, 7:04 p.m.