# NOT RUN {
library(metan)
cc1 = can_corr(data_ge2,
FG = c(PH, EH, EP),
SG = c(EL, ED, CL, CD, CW, KW, NR))
plot(cc1, 2)
cc2 = can_corr(data_ge2,
FG = c(PH, EH, EP),
SG = c(EL, ED, CL, CD, CW, KW, NR),
means_by = GEN)
plot(cc2, 2, labels = TRUE)
# }
Run the code above in your browser using DataLab