I would like to reassure myself that what I'm doing is right. I am applying an LME model to my data. And I would like to know if the use of EMMEANS to get multiple comparisons is right.
emmeans(mod1.lme, pairwise ~ Status, adjust="bonferroni")
The variable 'Status' has 4 factors, and I would like to know how to adjust my LME for 'Time since event' when only 2 of those factors are affected by that event.
Would be enough to include na.action = na.omit
at the end of the model or I should do anything else to control for that?