Generates balance statistics for CBPS and CBMSM objects from the CBPS package.
Usage
# S3 method for class 'CBPS'
bal.tab(
x,
stats,
int = FALSE,
poly = 1,
distance = NULL,
addl = NULL,
data = NULL,
continuous,
binary,
s.d.denom,
thresholds = NULL,
weights = NULL,
cluster = NULL,
imp = NULL,
pairwise = TRUE,
s.weights = NULL,
abs = FALSE,
subset = NULL,
quick = TRUE,
...
)Arguments
- x
a
CBPSorCBMSMobject; the output of a call toCBPS::CBPS()orCBPS::CBMSM().- stats
character; which statistic(s) should be reported. Seestatsfor allowable options. For binary and multi-category treatments,"mean.diffs"(i.e., mean differences) is the default. For continuous treatments,"correlations"(i.e., treatment-covariate Pearson correlations) is the default. Multiple options are allowed.- int
logicalornumeric; whether or not to include 2-way interactions of covariates included incovsand inaddl. Ifnumeric, will be passed topolyas well.- poly
numeric; the highest polynomial of each continuous covariate to display. For example, if 2, squares of each continuous covariate will be displayed (in addition to the covariate itself); if 3, squares and cubes of each continuous covariate will be displayed, etc. If 1, the default, only the base covariate will be displayed. Ifintis numeric,polywill take on the value ofint.- distance
an optional formula or data frame containing distance values (e.g., propensity scores) or a character vector containing their names. If a formula or variable names are specified,
bal.tab()will look in the argument todata, if specified. Propensity scores generated byCBPS()andCBMSM()are automatically included and named"prop.score". ForCBMSMobjects, can be a list of allowable arguments, one for each time point, but each dataset in the list supplied todistancemust have one row per individual, unlike the data frame in the original call toCBMSM().- addl
an optional formula or data frame containing additional covariates for which to present balance or a character vector containing their names. If a formula or variable names are specified,
bal.tab()will look in the arguments to the input object,covs, anddata, if specified. ForCBMSMobjects, can be a list of allowable arguments, one for each time point, but each dataset in the list supplied toaddlmust have one row per individual, unlike the data frame in the original call toCBMSM().- data
an optional data frame containing variables named in other arguments. For some input object types, this is required.
- continuous
whether mean differences for continuous variables should be standardized (
"std") or raw ("raw"). Default"std". Abbreviations allowed. This option can be set globally usingset.cobalt.options().- binary
whether mean differences for binary variables (i.e., difference in proportion) should be standardized (
"std") or raw ("raw"). Default"raw". Abbreviations allowed. This option can be set globally usingset.cobalt.options().- s.d.denom
character; how the denominator for standardized mean differences should be calculated, if requested. Seecol_w_smd()for allowable options. Abbreviations allowed. If not specified, will be set to"treated"if the estimand of the call toCBPS()is the ATT and "pooled" if the estimand is the ATE.- thresholds
a named vector of balance thresholds, where the name corresponds to the statistic (i.e., in
stats) that the threshold applies to. For example, to request thresholds on mean differences and variance ratios, one can setthresholds = c(m = .05, v = 2). Requesting a threshold automatically requests the display of that statistic. When specified, extra columns are inserted into the Balance table describing whether the requested balance statistics exceeded the threshold or not. Summary tables tallying the number of variables that exceeded and were within the threshold and displaying the variables with the greatest imbalance on that balance measure are added to the output.- weights
a vector, list, or
data.framecontaining weights for each unit, or a string containing the names of the weights variables indata, or an object with aget.w()method or a list thereof. The weights can be, e.g., inverse probability weights or matching weights resulting from a matching algorithm.- cluster
either a vector containing cluster membership for each unit or a string containing the name of the cluster membership variable in
dataor the input object. Seeclass-bal.tab.clusterfor details.- imp
either a vector containing imputation indices for each unit or a string containing the name of the imputation index variable in
dataor the input object. Seeclass-bal.tab.impfor details. Not necessary ifdatais amidsobject.- pairwise
whether balance should be computed for pairs of treatments or for each treatment against all groups combined. See
bal.tab.multi()for details. This can also be used with a binary treatment to assess balance with respect to the full sample.- s.weights
Optional; either a vector containing sampling weights for each unit or a string containing the name of the sampling weight variable in
data. These function like regular weights except that both the adjusted and unadjusted samples will be weighted according to these weights if weights are used. If sampling weights are used inCBPS(), it is critical that they are specified inbal.tab()as well using thes.weightsargument.- abs
logical; whether displayed balance statistics should be in absolute value or not.- subset
a
logicalornumericvector denoting whether each observation should be included or which observations should be included. Iflogical, it should have length equal to the number of units.NAs will be treated asFALSE. This can be used as an alternative toclusterto examine balance on subsets of the data.- quick
logical; ifTRUE, will not compute any values that will not be displayed. Set toFALSEif computed values not displayed will be used later.- ...
for some input types, other arguments that are required or allowed. Otherwise, further arguments to control display of output. See display options for details.
Value
For point treatments, if clusters are not specified, an object of class "bal.tab" containing balance summaries for the CBPS object. See bal.tab() for details.
If clusters are specified, an object of class "bal.tab.cluster" containing balance summaries within each cluster and a summary of balance across clusters. See class-bal.tab.cluster for details.
If CBPS() is used with multi-category treatments, an object of class "bal.tab.multi" containing balance summaries for each pairwise treatment comparison and a summary of balance across pairwise comparisons. See bal.tab.multi() for details.
If CBMSM() is used for longitudinal treatments, an object of class "bal.tab.msm" containing balance summaries for each time period and a summary of balance across time periods. See class-bal.tab.msm for details.
Details
bal.tab.CBPS() and bal.tab.CBMSM() generate a list of balance summaries for the CBPS or CBMSM object given and functions similarly to CBPS::balance().
See also
bal.tab()for details of calculations.class-bal.tab.clusterfor more information on clustered data.bal.tab.multi()for more information on multi-category treatments.class-bal.tab.msmfor more information on longitudinal treatments.
Examples
data("lalonde", package = "cobalt")
## Using CBPS() for generating covariate balancing
## propensity score weights
cbps.out <- CBPS::CBPS(treat ~ age + educ + married + race +
nodegree + re74 + re75,
data = lalonde)
#> [1] "Finding ATT with T=1 as the treatment. Set ATT=2 to find ATT with T=0 as the treatment"
bal.tab(cbps.out)
#> Balance Measures
#> Type Diff.Adj
#> prop.score Distance -0.0057
#> age Contin. -0.0052
#> educ Contin. -0.0017
#> married Binary -0.0029
#> race_black Binary 0.0019
#> race_hispan Binary -0.0002
#> race_white Binary -0.0017
#> nodegree Binary 0.0042
#> re74 Contin. -0.0078
#> re75 Contin. 0.0061
#>
#> Effective sample sizes
#> Control Treated
#> Unadjusted 429. 185
#> Adjusted 99.97 185
