inst/doc/sensemakr.R

In sensemakr: Sensitivity Analysis Tools for Regression Models

## ---- include = FALSE---------------------------------------------------------
knitr::opts_chunk$set(
  fig.height = 4.5,
  fig.width = 4.5,
  fig.align = 'center',
  collapse = TRUE,
  comment = "#>"
)

## ---- message=FALSE-----------------------------------------------------------
# loads package
library(sensemakr)

# loads data
data("darfur")

## -----------------------------------------------------------------------------
# runs regression model
darfur.model <- lm(peacefactor ~ directlyharmed  + village +  female +
                     age + farmer_dar + herder_dar + pastvoted + hhsize_darfur, 
                   data = darfur)

## ---- echo=FALSE,  comment = ""-----------------------------------------------
stargazer::stargazer(darfur.model, keep = "directlyharmed", type = "text")

## ---- eval=FALSE--------------------------------------------------------------
#  darfur.complete.model <- lm(peacefactor ~ directlyharmed  + village +  female +
#                                age + farmer_dar + herder_dar + pastvoted + hhsize_darfur +
#                                center*wealth*political_attitudes,
#                              data = darfur)

## ---- results = 'asis'--------------------------------------------------------
# runs sensemakr for sensitivity analysis
# in the darfur example
darfur.sensitivity <- sensemakr(model = darfur.model, 
                                treatment = "directlyharmed",
                                benchmark_covariates = "female",
                                kd = 1:3,
                                ky = 1:3, 
                                q = 1,
                                alpha = 0.05, 
                                reduce = TRUE)

## -----------------------------------------------------------------------------
darfur.sensitivity <- sensemakr(model = darfur.model, 
                                treatment = "directlyharmed",
                                benchmark_covariates = "female",
                                kd = 1:3)

## -----------------------------------------------------------------------------
darfur.sensitivity

## ---- results='asis'----------------------------------------------------------
ovb_minimal_reporting(darfur.sensitivity, format = "html")

## ---- results='hide'----------------------------------------------------------
summary(darfur.sensitivity)

## -----------------------------------------------------------------------------
plot(darfur.sensitivity)

## -----------------------------------------------------------------------------
plot(darfur.sensitivity, sensitivity.of = "t-value")

## ---- fig.width=6-------------------------------------------------------------
plot(darfur.sensitivity, type = "extreme")