tests/testthat/test-two-arm.R

In acoppock/ri2: Randomization Inference for Randomized Experiments

context("Two Arm Trial")


test_that("Basic Two Arm Trial", {
  N <- 100
  declaration <- randomizr::declare_ra(N = N, m = 50)

  Z <- randomizr::conduct_ra(declaration)
  X <- rnorm(N)
  Y <- .9 * X + .2 * Z + rnorm(N)
  W <- runif(N)
  df <- data.frame(Y, X, Z, W)

  ri_out <-
    conduct_ri(
      formula = Y ~ Z,
      declaration = declaration,
      assignment = "Z",
      sharp_hypothesis = 0,
      data = df, sims = 100
    )

  plot(ri_out)
  summary(ri_out)

  ri_out <-
    conduct_ri(
      formula = Y ~ Z + X,
      declaration = declaration,
      assignment = "Z",
      sharp_hypothesis = 0,
      data = df, sims = 100
    )


  plot(ri_out)
  summary(ri_out)

  expect_true(TRUE)
})