## Example: Create a function for creating groupwise t-tests
mydata <- data.frame(y=rnorm(32), x=rnorm(32),
g1=factor(rep(c(1,2),each=16)), g2=factor(rep(c(1,2), each=8)),
g3=factor(rep(c(1,2),each=4)))
t.testBy <- function(formula, data, ...){
createFunBy(formula, data, FUN=t.test, class="t.testBy", ...)
}
t.testBy(y~g1|g2+g3, data=mydata)
Run the code above in your browser using DataLab