Introduction to konfound

Sarah Narvaiz, Qinyun Lin, Joshua M. Rosenberg, Kenneth A. Frank, Spiro Maroulis, Wei Wang, Ran Xu

2024-10-17

Note: This is identical to the following open-access paper:

Narvaiz, S., Lin, Q., Rosenberg, J. M., Frank, K. A., Maroulis, S. J., Wang, W., & Xu, R. (2024). konfound: An R Sensitivity Analysis Package to Quantify the Robustness of Causal Inferences. Journal of Open Source Software, 9(95), 1-6. https://joss.theoj.org/papers/10.21105/joss.05779.pdf

Quantifying the Robustness of Inferences

Sensitivity analysis, a statistical method crucial for validating inferences across disciplines, quantifies the conditions that could alter conclusions (Razavi et al. 2021). One line of work is rooted in linear models and foregrounds the sensitivity of inferences to the strength of omitted variables (Frank 2000; Cinelli and Hazlett 2019). A more recent approach is rooted in the potential outcomes framework for causal inference and foregrounds how hypothetical changes in a sample would alter an inference if such cases were otherwise observed (Frank and Min 2007; Frank et al. 2008, 2013; Xu et al. 2019).

One sensitivity measure is the Impact Threshold of a Confounding Variable, or ITCV, which generates statements about the correlation of an omitted, confounding variable with both a predictor of interest and the outcome (Frank 2000). The ITCV index can be calculated for any linear model. The Robustness of an Inference to Replacement, RIR, assesses how replacing a certain percentage of cases with counterfactuals of zero treatment effect could nullify an inference (Frank et al. 2013). The RIR index is more general than the ITCV index.

The sensitivity analysis techniques we describe in this paper and implement in the konfound R package differ from others in several ways. Unlike Linden, Mathur, and VanderWeele (2020), whose approach focuses on dichotomous outcomes and omitted variable sensitivity, our approach extends to continuous outcomes and evaluates both changes in estimates and standard errors. Oster (2019) focuses only on selection into the treatment based on unobservable variables versus observable variables necessary to nullify an estimate. The ITCV index focuses on the relationship of the unobservable to the predictor of interest and to the outcome. More generally, many others used simulation-based approaches, while our approach uses closed-form expressions to generate a single term representing sensitivity. These techniques, along with others, are reviewed and discussed (along with the ITCV and RIR approaches) in Frank et al. (2023).

We have implemented the calculation of both the ITCV and RIR indices in the konfound R package. This package is intended to provide an easy-to-use and principled set of sensitivity techniques that are suitable for a range of model and dependent variable types and use cases. Its audience is broad: primarily social scientists, but also interested individuals in other disciplines (e.g., the health sciences). This paper provides an overview of two core functions within the konfound package, each of which can calculate the ITCV and RIR indices: konfound() and pkonfound(). These functions allow users to calculate the robustness of inferences using a model estimated (with R) or using information about a model from a published study, respectively.

The konfound package is available from the Comprehensive R Archive Network (CRAN) at https://CRAN.R-project.org/package=konfound; it can be installed via the install.packages(“konfound”) function within R.

Functionality

konfound

This function calculates the ITCV and RIR for models fitted in R. This function currently works for linear models fitted with lm(), glm(), and lmer(). The output printed in the R console is the bias that must be present and the number of cases that would have to be replaced with cases for which there is no effect to nullify the inference.

Example for linear models fit with lm()

For this example, we will use the concord1 dataset built into the konfound package. This dataset comes from a study that examines the causal mechanism behind household water conservation in a U.S. city.

We will estimate the effect of the following variables on household water consumption in 1981:

  • household water consumption in 1980 (water80)
  • household income (income)
  • education level of household survey respondent (educat)
  • retirement status of respondent (retire)
  • number of individuals in household in 1980 (peop80)

Here is the code we use to fit a linear model using these variables:

m <- lm(water81 ~ water80 + income + educat + retire + peop80, data = concord1)

The results of the model fitting (which can be obtained by running summary(m) within R) indicate that all of the predictors apart from retire have a statistically significant effect on water consumption. In the example, we focus on the coefficient for peop80 (β = 225.198, SE = 28.704, t = 7.845, p < .001).

ITCV example for linear models fit with lm()

Now, let’s examine the robustness of the peop80 effect by calculating the ITCV:

library(konfound)

konfound(m, peop80, index = "IT")

## Impact Threshold for a Confounding Variable:
## The minimum impact of an omitted variable to invalidate an inference
## for a null hypothesis of 0 effect is based on a correlation of 0.52 with 
## the outcome and at  0.52  with the predictor of interest (conditioning on
## observed covariates) based on a threshold of 0.089 for statistical 
## significance (alpha = 0.05).
## 
## Correspondingly the impact of an omitted variable (as defined in Frank 
## 2000) must be 0.52 X 0.52 = 0.27 to invalidate an inference for a null 
## hypothesis of 0 effect. See Frank (2000) for a description of the method.
## 
## Citation:
## Frank, K. (2000). Impact of a confounding variable on the
## inference of a regression coefficient. Sociological Methods and Research,
## 29(2), 147-194

## For more detailed output, consider setting `to_return` to table

## To consider other predictors of interest, consider setting `test_all` to
## TRUE.

The output indicates that to invalidate the inference that peop80 has an effect on water81 using statistical significance as a threshold (e.g., p = .05), an omitted variable would have to be correlated at 0.520 with peop80 and 0.520 with water81, conditioning on observed covariates.

RIR example for linear models fit with lm()

We can also examine the robustness by calculating the RIR:

konfound(m, peop80, index = "RIR") 

## Robustness of Inference to Replacement (RIR):
## To invalidate an inference,  74.955 % of the 
## estimate would have to be due to bias. 
## This is based on a threshold of 56.4 for statistical 
## significance (alpha = 0.05).
##
## To invalidate an inference,  372  observations would 
## have to be replaced with cases for which the effect is 0 (RIR = 372).
## 
## See Frank et al. (2013) for a description of the method.
## 
## Citation: Frank, K.A., Maroulis, S., Duong, M., and Kelcey, B. (2013).
## What would it take to change an inference?
## Using Rubin's causal model to interpret the robustness of causal inferences.
## Education, Evaluation and Policy Analysis, 35 437-460.

## For more detailed output, consider setting `to_return` to table

## To consider other predictors of interest, 
## consider setting `test_all` to TRUE.

The output presents two interpretations of the RIR. First, 74.956% of the estimated effect of peop80 on water81 would have to be attributed to bias to invalidate the inference. Equivalently, we would expect to have to replace 372 out of the 486 households (about 76%) with cases for which peop80 had no effect to invalidate the inference. We have created guidelines for evaluating the RIR relative to the bias accounted for by observed covariates or published norms (Frank et al. 2013, 2021).

pkonfound

This function quantifies the sensitivity for analyses that have already been conducted, such as in an already-published study or from analysis carried out using other software. This function calculates how much bias there must be in an estimate to invalidate/sustain an inference, which can be interpreted as the percentage of cases that would have to be replaced (e.g., with cases for which the predictor had no effect on the outcome) to invalidate the inference. It also calculates the impact of an omitted variable necessary to invalidate/sustain an inference for a regression coefficient, where impact is defined as the correlation between the omitted variable and the focal predictor multiplied with the correlation between the omitted variable and the outcome.

ITCV example for a regression analysis

For this example, the following estimated quantities from an estimated regression model would be entered into the pkonfound function: unstandardized beta coefficient for the predictor of interest (est_eff = 2), estimated standard error (std_err = .4), sample size (n_obs = 100), and the number of covariates (n_covariates = 3), as follows:

pkonfound(2, .4, 100, 3, index = "IT")

## Impact Threshold for a Confounding Variable:
## The minimum impact of an omitted variable to invalidate an inference for 
## a null hypothesis of 0 effect is based on a correlation of  0.568 with 
## the outcome and at  0.568  with the predictor of interest (conditioning 
## on observed covariates) based on a threshold of 0.201 for statistical 
## significance (alpha = 0.05).
## 
## Correspondingly the impact of an omitted variable (as defined in Frank 
## 2000) must be 0.568 X 0.568 = 0.323 to invalidate an inference for a null
## hypothesis of 0 effect. See Frank (2000) for a description of the method.
## 
## Citation:
## Frank, K. (2000). Impact of a confounding variable on the inference of a 
## regression coefficient. Sociological Methods and Research, 29 (2), 147-194

## For other forms of output, run ?pkonfound and inspect the to_return argument

## For models fit in R, consider use of konfound().

RIR example for a regression analysis

We can also use the same inputs to calculate output for the RIR index:

pkonfound(2, .4, 100, 3, index = "RIR")

## Robustness of Inference to Replacement (RIR):
## To invalidate an inference,  60.29 % of the estimate would have to be 
## due to bias. This is based on a threshold of 0.794 for statistical 
## significance (alpha = 0.05).
## 
## To invalidate an inference,  60  observations would have to be replaced 
## with cases for which the effect is 0 (RIR = 60).
## 
## See Frank et al. (2013) for a description of the method.
## 
## Citation: Frank, K.A., Maroulis, S., Duong, M., and Kelcey, B. (2013).
## What would it take to change an inference?
## Using Rubin's causal model to interpret the robustness of causal inferences.
## Education, Evaluation and Policy Analysis, 35 437-460.

## For other forms of output, run ?pkonfound and inspect the to_return argument

## For models fit in R, consider use of konfound().

Doing and Learning More

We have created a website including a Shiny interactive web application at http://konfound-it.com that can be applied to linear models, 2x2 tables resulting corresponding to treatment and control by success and failure conditions, and logistic regression models. We are also developing extensions of the sensitivity analysis techniques described in this paper, including preserving the standard error (Frank et al. 2023) and calculating the coefficient of proportionality (Frank et al. 2022) for ITCV analyses. Functionality for designs including mediation, hazard functions, differences in difference, and regression discontinuity are also presently under development. Additional documentation on the R package and its future extensions are available at the http://konfound-it.com website.

Acknowledgements

The research reported here was supported by the Institute of Education Sciences, U.S. Department of Education, through Grant R305D220022 to Michigan State University. The opinions expressed are those of the authors and do not represent views of the Institute or the U.S. Department of Education.

References

Cinelli, C., and C. Hazlett. 2019. Making Sense of Sensitivity: Extending Omitted Variable Bias.” Journal of the Royal Statistical Society Series B: Statistical Methodology 82 (1): 39–67. https://doi.org/10.1111/rssb.12348.
Frank, K. A. 2000. “Impact of a Confounding Variable on a Regression Coefficient.” Sociological Methods & Research 29 (2): 147–94. https://doi.org/10.1177/0049124100029002001.
Frank, K. A., Qinyun Lin, Spiro Maroulis, Shimeng Dai, Nicole Jess, Hung-chang Lin, Sarah Maestrales, Ellen Searle, Jordan Tait, et al. 2022. “Improving Oster’s \(\delta\)*: Exact Calculation for the Coefficient of Proportionality Without Subjective Specification of a Baseline Model.” Ellen and Tait, Jordan, Improving Oster’s \(\delta\)*: Exact Calculation for the Coefficient of Proportionality Without Subjective Specification of a Baseline Model (December 16, 2022). https://doi.org/10.2139/ssrn.4305243.
Frank, K. A., Q. Lin, S. Maroulis, A. S. Mueller, R. Xu, J. M. Rosenberg, C. S. Hayter, et al. 2021. “Hypothetical Case Replacement Can Be Used to Quantify the Robustness of Trial Results.” Journal of Clinical Epidemiology 134: 150–59. https://doi.org/10.1016/j.jclinepi.2021.01.025.
Frank, K. A., Q. Lin, R. Xu, S. Maroulis, and A. Mueller. 2023. “Quantifying the Robustness of Causal Inferences: Sensitivity Analysis for Pragmatic Social Science.” Social Science Research 110 (February). https://doi.org/10.1016/j.ssresearch.2022.102815.
Frank, K. A., S. J. Maroulis, M. Q. Duong, and B. M. Kelcey. 2013. “What Would It Take to Change an Inference? Using Rubin’s Causal Model to Interpret the Robustness of Causal Inferences.” Educational Evaluation and Policy Analysis 35 (4): 437–60. https://doi.org/10.3102/0162373713493129.
Frank, K. A., and K.-S. Min. 2007. “Indices of Robustness for Sample Representation.” Sociological Methodology 37 (1): 349–92. https://doi.org/10.1111/j.1467-9531.2007.00186.x.
Frank, K. A., G. Sykes, D. Anagnostopoulos, M. Cannata, L. Chard, A. Krause, and R. McCrory. 2008. “Does NBPTS Certification Affect the Number of Colleagues a Teacher Helps with Instructional Matters?” Educational Evaluation and Policy Analysis 30 (1): 3–30. https://doi.org/10.3102/0162373707313781.
Linden, Ariel, Maya B Mathur, and Tyler J VanderWeele. 2020. “Conducting Sensitivity Analysis for Unmeasured Confounding in Observational Studies Using e-Values: The Evalue Package.” The Stata Journal 20 (1): 162–75. https://doi.org/10.1177/1536867X20909696.
Oster, Emily. 2019. “Unobservable Selection and Coefficient Stability: Theory and Evidence.” Journal of Business & Economic Statistics 37 (2): 187–204.
Razavi, S., A. Jakeman, A. Saltelli, C. Prieur, B. Iooss, E. Borgonovo, E. Plischke, et al. 2021. “The Future of Sensitivity Analysis: An Essential Discipline for Systems Modeling and Policy Support.” Environmental Modelling & Software 137 (March). https://doi.org/10.1016/j.envsoft.2020.104954.
Xu, R., K. A. Frank, S. J. Maroulis, and J. M. Rosenberg. 2019. “Konfound: Command to Quantify Robustness of Causal Inferences.” The Stata Journal 19 (3): 523–50. https://doi.org/10.1177/1536867X19874223.