# Example data for 2x2 factorial design (Fertilization x Warming)
data <- data.frame(
study_id = 1:2,
control_mean = c(10, 12), control_sd = c(2.0, 2.5), control_n = c(20, 18),
fertilization_mean = c(15, 16), fertilization_sd = c(2.2, 2.8), fertilization_n = c(20, 19),
warming_mean = c(11, 13), warming_sd = c(2.1, 2.6), warming_n = c(21, 17),
fert_warm_mean = c(17, 19), fert_warm_sd = c(2.4, 3.0), fert_warm_n = c(19, 20)
)
# Compute main effect of fertilization
result <- lnRR_main(
data = data,
Ctrl_mean = "control_mean", Ctrl_sd = "control_sd", Ctrl_n = "control_n",
A_mean = "fertilization_mean", A_sd = "fertilization_sd", A_n = "fertilization_n",
B_mean = "warming_mean", B_sd = "warming_sd", B_n = "warming_n",
AB_mean = "fert_warm_mean", AB_sd = "fert_warm_sd", AB_n = "fert_warm_n"
)
Run the code above in your browser using DataLab