ald <- data.frame(trt = c("B","C","B","C"),
variable = c(NA, NA, "y", "y"),
statistic = c("N", "N", "sum", "sum"),
value = c(100, 100, 50, 60))
marginal_variance(ald, ref_trt = "C", comp_trt = "B",
scale = "log_odds", family = "binomial")
Run the code above in your browser using DataLab