reg.ordcomp(var, covar, fact, conf.level = 0.95, coeff.arr = 4, theo = rep(0, nlevels(fact)), p.method = "fdr", p.arr = 4)
p.adjust()
.lm
, residuals
, combinations
, p.adjust
# Tests par defaut
covar<-rep(1:30,3)
var<-c(1:30+rnorm(30,0,3),seq(10,35,25/29)+rnorm(30,0,3),seq(-27,0,27/29)+rnorm(30,0,3))
facteur<-as.factor(rep(LETTERS[1:3],each=30))
reg.ordcomp(var,covar,facteur)
# Intervalles de confiance a 99 pc, ordonnees a l'origine theoriques non nulles
# et correction des p-values par la methode de Bonferroni
reg.ordcomp(var,covar,facteur,conf.level=0.99,theo=c(1,8,-20),p.method="bonferroni")
Run the code above in your browser using DataLab