My design is in split-split plot arrangement. There are three variables each with three levels and the study has three replications. Its a mixed effects model.
I used agricoale package with R to carry out the analysis. However, I am not able to make a pairwise comparison with Tukey adjustment for different combinations of these factors. When I call the TukeyHSD() function for the output of ssp.plot(), it throws an error that say 'TukeyHSD applied to an object of class list.'
Can I use the lme4 package's lmer() function (used for split plot design) and just add one more factor to make it a split split plot analysis ?
Hi Bhaskar! Have you try using emmeans? I have more more experience with split-plot design using SAS but thinking on your situation, maybe you could try something like this:
- Badges
- Science topic
- Similar topics
- Mathematics
- Statistics
- Sampling Design
- Sample Size