Test post-hoc anova a plusieurs facteurs

Bonjour,
j'ai fait une anova a 3 facteurs avec la fonction Anova du package car.
La variable réponse est : valeur carbone et les facteurs : Taxon et Saison, Année
J'obtiens une différence significative entre les facteurs Taxon x Saison. J'aimerai faire un test post-hoc derrière pour connaitre les différents groupes, mais je ne sais pas comment faire comme c'est avec 2 facteurs.

J'ai essayer avec glht (mod2, linfct=mcp(   ="Tukey")) du package multcomp, mais je peux mettre qu'un facteur.

Pourriez-vous m'aider ?

Merci d'avance
Kikibigoudi

Que veux-tu dire par: "les facteurs Taxon x Saison" ?
Tu veux dire l'interaction des 2 facteurs ?

Utilise le package emmeans, l'argument specs permet de spécifier la formule qui devrait être du type pairwise ~ Taxon + Saison