## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----cic-example--------------------------------------------------------------
library(sccic)

# Load workers' comp data
if (requireNamespace("wooldridge", quietly = TRUE)) {
  data("injury", package = "wooldridge")

  result <- cic(
    y_00 = injury$ldurat[injury$highearn == 0 & injury$afchnge == 0],
    y_01 = injury$ldurat[injury$highearn == 0 & injury$afchnge == 1],
    y_10 = injury$ldurat[injury$highearn == 1 & injury$afchnge == 0],
    y_11 = injury$ldurat[injury$highearn == 1 & injury$afchnge == 1]
  )
  print(result)
}

## ----sccic-example------------------------------------------------------------
if (requireNamespace("Synth", quietly = TRUE)) {
  data("basque", package = "Synth")

  # Reshape to wide format
  gdp <- reshape(basque[, c("regionno", "year", "gdpcap")],
                 idvar = "year", timevar = "regionno", direction = "wide")

  y_treated <- gdp[, "gdpcap.17"]  # Basque Country
  donor_cols <- grep("gdpcap\\.", names(gdp), value = TRUE)
  donor_cols <- donor_cols[!donor_cols %in% c("gdpcap.17", "gdpcap.1")]
  donors <- as.matrix(gdp[, donor_cols])

  valid <- complete.cases(y_treated, donors)
  result2 <- sc_cic(y_treated[valid], donors[valid, ],
                    treatment_period = 16, seed = 42)
  print(result2)
}

## ----boot-example, eval = FALSE-----------------------------------------------
#  result <- cic(y_00, y_01, y_10, y_11, boot = TRUE, boot_iters = 500, seed = 42)