Learn R Programming

rqlm (version 4.3-1)

SumStat: Creating summary table for IPTW analysis using stabilized weights

Description

Summary statistics are computed before and after weighting for IPTW analyses. For all covariates included in the weight model, the mean, standard deviation, and standardized mean difference (SMD) are calculated. Weighting can be performed using either stabilized or unstabilized weights, with optional truncation. For binary covariates, summary statistics are computed using unbiased estimators of the binomial mean and variance.

Usage

SumStat(formula, data, trunc=c(0.01,0.99), digits=3)

Value

Summary statistics for unadjusted and IPTW-adjusted populations are provided. For each covariate, the mean and standard deviation, as well as the standardized mean difference (SMD), were calculated. Truncated stabilized weights were used for weighting. Setting the trunc argument to c(0, 1) allows weighting with untruncated stabilized weights.

  • mean0: Mean for group 0.

  • mean1: Mean for group 1.

  • sd0: SD for group 0.

  • sd1: SD for group 1.

  • SMD: SMD for the original population.

  • wmean0: IPTW-weighted mean for group 0.

  • wmean1: IPTW-weighted mean for group 1.

  • wsd0: IPTW-weighted SD for group 0.

  • wsd1: IPTW-weighted SD for group 1.

  • wSMD: IPTW-weighted SMD for the adjusted pseudo-population.

  • type: Type of the variable.

Arguments

formula

An object of class "formula" (or one that can be coerced to that class): a symbolic description of the logistic regression model to be fitted.

data

A data frame, list or environment (or object coercible by as.data.frame to a data frame) containing the variables in the model.

trunc

Quantiles to be truncated the weights (default: 0.01, 0.99).

digits

Number of decimal places in the output (default: 4).

References

Austin, P. C. (2009). Balance diagnostics for comparing the distribution of baseline covariates between treatment groups in propensity-score matched samples. Statistics in Medicine 28, 3083-3107.

Examples

Run this code
data(exdata02)

SumStat(x2 ~ x1 + x3 + x4, data=exdata02)

Run the code above in your browser using DataLab