# Example with 3 groups
set.seed(123)
g1 <- rnorm(20, mean = 5, sd = 1.5)
g2 <- rnorm(25, mean = 6.5, sd = 2)
g3 <- rnorm(22, mean = 7.2, sd = 2.5)
means <- c(mean(g1), mean(g2), mean(g3))
vars <- c(var(g1), var(g2), var(g3))
ns <- c(length(g1), length(g2), length(g3))
result <- PPBMT(means, vars, ns, B = 1000, alpha = 0.05)
print(result)
Run the code above in your browser using DataLab