# NOT RUN {
set.seed(1234)
post <- c(rep(0, 100), rep(1, 100))
tx <- c(rep(0, 100), rbinom(n = 100, prob = 0.27, size = 1))
y <- rbinom(n = 200, prob = 0.1 + 0.02 * post - 0.05 * tx, size = 1)
# Two cohorts, full data
scqe.2cohort.full = scqe(post = post, treatment = tx, outcome = y,
delta=seq(from = -0.1, to = 0.1, by = 0.05))
summary(scqe.2cohort.full)
set.seed(1234)
post = c(rep(0,100), rep(1,100))
tx = c(rep(0, 100), rbinom(n = 100, prob = 0.27, size = 1))
y = rbinom(n= 200, prob = 0.1 + 0.02 * post - 0.05 * tx, size = 1)
# Two cohorts, full data
scqe.2cohort.full = scqe(post = post, treatment = tx, outcome = y,
delta = seq(from = -0.1, to = 0.1, by = 0.05))
scqe.2cohort.full
print(scqe.2cohort.full)
# }
Run the code above in your browser using DataLab