# NOT RUN {
data(pocmaj)
library(dplyr)
pocmajqt <- as.qtag(pocmaj, id.vars = c("core", "depth"))
pocmajqt %>% group() %>% summarise(mean(Ca))
# equivalent to
pocmaj %>% group_by(core, depth) %>% summarise(mean(Ca))
# }
Run the code above in your browser using DataLab