# 1. Create a sample data frame
data <- data.frame(
group = c(1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3),
response = c(16, 60, 82, 126, 137, 44, 67, 87, 100, 142, 17, 28, 105, 149, 160),
covariate1 = c(26, 10, 42, 49, 55, 21, 28, 5, 12, 58, 1, 19, 41, 48, 35)
)
# 2. Run the Burnett and Barr method
results <- Burnett_Barr(
formula = response ~ covariate1 + group,
data = data
)
# 3. View the results
print(results)
print(results$anova)
Run the code above in your browser using DataLab