In mbannick/RoboCar:

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

rm(list=ls())
library(RoboCar)

Adjustment with Continuous Outcomes

The dataset that we will use for this is included in the package.

data <- RoboCar::data_sim
data$A <- as.factor(data$A)

ANOVA Model

fit.anova <- robincar_linear(df = data, 
                             response_col="y",
                             treat_col="A",
                             covariate_cols=c("x1", "x3"),
                             car_scheme="simple",
                             adj_method="ANOVA",
                             vcovHC="HC0")

ANCOVA model with biased coin randomization

fit.ancova <- robincar_linear(df = data, 
                              response_col="y",
                              treat_col="A",
                              strata_cols=c("z1", "z2"),
                              covariate_cols=c("x1", "x3"),
                              car_scheme="biased-coin",
                              adj_method="ANCOVA",
                              vcovHC="HC0")
fit.ancova$result

ANHECOVA model with Pocock-Simon randomization and linear contrast

fit.anhecova <- robincar_linear(df = data, 
                                response_col="y",
                                treat_col="A",
                                strata_cols=c("z1", "z2"),
                                covariate_cols=c("x1", "x3"),
                                car_scheme="biased-coin",
                                adj_method="ANCOVA",
                                vcovHC="HC0",
                                contrast_h="diff")

fit.anhecova$contrast$result