Une solution est l'utilisation de la fonction lstrends de lsmeans . Exemple : d.f <- data.frame( X=runif( 60), G=c( "a", "b", "c"), P=c( 0, 0.5, 1)) d.f$Y <- d.f$X * d.f$P + rnorm( 60, sd=0.1) d.f$G <- factor( d.f$G) res <- lm( Y ~ X * G, d.f) lstrends( res, pairwise ~ ...