update DIDmultiplegt code

Author: graemeblair

DESCRIPTION

@@ -2,4 +2,4 @@ Package: rdssbook
 Title: Research Design: Declare, Diagnose, Redesign
 Version: 0.0.1
 Suggests: pak, ggcorrplot, modelr, ICC, gt, ggspatial, GGally, psych, rdrobust, DIDmultiplegt, car, checkpoint, dagitty, downlit, dplyr, estimatr, fabricatr, forcats, ggdag, ggforce, ggplot2, ggraph, ggrepel, ggridges, ggtext, glmnet, glue, grid, gridExtra, haven, jsonlite, kableExtra, knitr, latex2exp, lme4, magrittr, MASS, metafor, patchwork, prediction, purrr, randomizr, Rcpp, readr, remotes, reshape2, rmarkdown, rprojroot, rr, sf, shiny, stringr, texreg, tibble, tidyr, tidyverse, xml2, broom.mixed, parsermd, rstanarm, CausalQueries, DeclareDesign, DesignLibrary, ragg, margins, MatchIt, cjoint, spdep, grf, bbmle, rdss, vayr, dddag, randnet, interference
-Remotes: acoppock/vayr, DeclareDesign/dddag, cran/randnet, szonszein/interference, cran/DIDmultiplegt, cran/cjoint
+Remotes: acoppock/vayr, DeclareDesign/dddag, cran/randnet, szonszein/interference

errata.qmd

@@ -21,3 +21,4 @@ These errata are errors in the printed book. These changes are implemented in th
 We are committed to maintaining the code in this book as the R ecosystem evolves. The following changes reflect differences between the printed book and the online version. 
 
 - On pg. 103, Declaration 9.3 now reads `.summary = ~tidy(., exponentiate = TRUE)` instead of `.summary = ~tidy_stan(., exponentiate = TRUE)`, reflecting the use of the standard `broom.mixed` package for tidying stan objects. In addition, the library `rdss` is replaced with the library `broom.mixed` in the declaration.
+- On pg. 212-213, Declaration 16.3 now uses updated syntax from the `DIDmultiplegt` package, specifically `mode = 'old'` is added. We also add `Sys.setenv(RGL_USE_NULL = TRUE)` to address package dependency issues.

scripts_declarations/declaration_16.3.R

@@ -1,14 +1,15 @@
 library(rdss) # for helper functions
+Sys.setenv(RGL_USE_NULL = TRUE) # required to use DIDmultiplegt
 library(DIDmultiplegt)
 
 N_units <- 20
 N_time_periods <- 20
 
-declaration_16.3 <- 
+declaration_16.3 <-
   declare_model(
   units = add_level(
-    N = N_units, 
-    U_unit = rnorm(N), 
+    N = N_units,
+    U_unit = rnorm(N),
     D_unit = if_else(U_unit > median(U_unit), 1, 0),
     D_time = sample(1:N_time_periods, N, replace = TRUE)
   ),
@@ -18,21 +19,21 @@ declaration_16.3 <-
     nest = FALSE
   ),
   unit_period = cross_levels(
-    by = join_using(units, periods), 
+    by = join_using(units, periods),
     U = rnorm(N),
-    potential_outcomes(Y ~ U + U_unit + U_time + 
-                         D * (0.2 - 1 * (D_time - as.numeric(periods))), 
+    potential_outcomes(Y ~ U + U_unit + U_time +
+                         D * (0.2 - 1 * (D_time - as.numeric(periods))),
                        conditions = list(D = c(0, 1))),
     D = if_else(D_unit == 1 & as.numeric(periods) >= D_time, 1, 0),
     D_lag = lag_by_group(D, groups = units, n = 1, order_by = periods)
   )
 ) +
   declare_inquiry(
-    ATT = mean(Y_D_1 - Y_D_0), 
+    ATT = mean(Y_D_1 - Y_D_0),
     subset = D == 1
-  ) + 
+  ) +
   declare_inquiry(
-    ATT_switchers = mean(Y_D_1 - Y_D_0), 
+    ATT_switchers = mean(Y_D_1 - Y_D_0),
     subset = D == 1 & D_lag == 0 & !is.na(D_lag)
   ) +
   declare_measurement(Y = reveal_outcomes(Y ~ D)) +
@@ -43,11 +44,12 @@ declaration_16.3 <-
     label = "twoway-fe"
   ) +
   declare_estimator(
-    Y = "Y", 
-    G = "units", 
-    T = "periods", 
+    Y = "Y",
+    G = "units",
+    T = "periods",
     D = "D",
+    mode = "old",
     handler = label_estimator(did_multiplegt_tidy),
     inquiry = c("ATT", "ATT_switchers"),
     label = "chaisemartin"
-  ) 
+  )