I would like to compare my data using the TukeyHSD() function. I read online it first uses 'model', which is your data and then it is followed by the 'which' parameter that selects data. When I try to run my code I receive an error.
sexe.mult <- TukeyHSD(pois.log, which = "Sex")
sexe.mult$Sexe
The Tuckey HSD test is a post-hoc test: it can only be used after another test (typically, ANOVA). This is why the first argument should be a fitted model object, typically created from aov().
From this error, it looks like pois.log is a vector, not a model (i.e. not the result of a fitting function).