Remove -1 as ref

Author: kylebutts

DESCRIPTION

@@ -1,6 +1,6 @@
 Package: did2s
 Title: Two-Stage Difference-in-Differences Following Gardner (2021)
-Version: 1.2.0
+Version: 1.2.1
 Authors@R: 
     c(person(given = "Kyle",
              family = "Butts",

R/data.R

@@ -64,8 +64,9 @@
 #'
 #' @format A data frame with 550 rows and 5 variables:
 #' \describe{
-#'   \item{sid}{state id, unit of observation}
+#'   \item{state}{state id, unit of observation}
 #'   \item{year}{time in panel data}
+#'   \item{popwt}{population-based weights}
 #'   \item{l_homicide}{log of the number of homicides per capita}
 #'   \item{effyear}{year that castle doctrine is passed}
 #'   \item{post}{0/1 variable for when castle doctrine is active}

R/did2s.R

@@ -62,7 +62,7 @@
 #' es <- did2s(df_hom,
 #'     yname = "dep_var", treatment = "treat", cluster_var = "state",
 #'     first_stage = ~ 0 | unit + year,
-#'     second_stage = ~ i(rel_year, ref=c(-1, Inf)))
+#'     second_stage = ~ i(rel_year, ref=Inf))
 #'
 #' fixest::esttable(es)
 #' ```
@@ -77,15 +77,16 @@
 #' Here's an example using data from Cheng and Hoekstra (2013)
 #' ```{r, comment = "#>", collapse = TRUE}
 #' # Castle Data
-#' castle <- haven::read_dta("https://github.com/scunning1975/mixtape/raw/master/castle.dta")
+#' data(castle, package = "did2s")
 #'
 #' did2s(
 #' 	data = castle,
 #' 	yname = "l_homicide",
-#' 	first_stage = ~ 0 | sid + year,
+#' 	first_stage = ~ 0 | state + year,
 #' 	second_stage = ~ i(post, ref=0),
 #' 	treatment = "post",
-#' 	cluster_var = "state", weights = "popwt"
+#' 	cluster_var = "state",
+#'  weights = "popwt"
 #' )
 #' ```
 #'
@@ -199,6 +200,7 @@ did2s <- function(
     IF <- IF_fs - IF_ss
 
     cl <- data[[cluster_var]]
+
     cov <- Reduce(
       "+",
       lapply(

data-raw/gen_data.R

@@ -40,7 +40,15 @@ castle <- haven::read_dta("data-raw/castle.dta")
 data.table::setDT(castle)
 castle[, time_til := data.table::fifelse(is.na(time_til), -Inf, time_til)]
 
-castle = castle[, .(year, sid, l_homicide, post, effyear, time_til)]
+castle = castle[, .(
+  year,
+  state = sid,
+  popwt,
+  l_homicide,
+  post,
+  effyear,
+  time_til
+)]
 
 df_hom = as.data.frame(df_hom)
 df_het = as.data.frame(df_het)