# Test par defaut
var1<-c(1:15+rnorm(15,0,4),1:15+rnorm(15,0,1),1:15+rnorm(15,0,8))
var2<-c(-1:-15+rnorm(15,0,4),1:15+rnorm(15,0,1),1:15+rnorm(15,0,8))
facteur<-as.factor(rep(LETTERS[1:3],each=15))
cor.multcomp(var1,var2,facteur)
# Correction des p-values avec la methode de Bonferroni
cor.multcomp(var1,var2,facteur,p.method="bonferroni")
Run the code above in your browser using DataLab