dat <- data.frame(group = c("a", "b", "a", "b"),
t1 = c(9, 4, 8, 10),
t2 = c(8, 8, 12, 10),
t3 = c(8, 5, 10, 11))
dat$group <- factor(dat$group, levels = c("a", "b"))
h <- hypothesis(1:3)
opamod <- opa(dat[,2:4], h, group = dat$group)
z <- compare_groups(opamod, "a", "b")
summary(z)
Run the code above in your browser using DataLab