I have a rookie question about pairwise comparisons with emmeans following linear mixed effects modeling. I'm planning an experiment and I'm trying to thing of the correct statistical methods to analyze the data.
Let's say I have an experiment with two groups of participants (beginners vs. advanced) who are tested in three conditions (condition A, condition B, condition C). The order of conditions is counterbalanced. I'm interested in the analysis of their performance, so I have the following model formula:
model <- lmer(performance ~ group * condition + condition_order + (1|participant), data)
The model informs me about the main effect of condition (averaged across groups), the main effect of group (averaged across conditions), order effect, and the interaction of group x condition.
What I'd also like to see is whether there are any within-group differences, i.e., whether the performance of beginner and advanced participants differs across the three conditions when we look at both groups separately.
Would a pairwise comparison with emmeans be statistically sound in such a case? If so, is it OK to proceed with it even if there's no interaction detected in the model?