Greetings: I am doing multiple comparisons of least squares means in proc glimmix, with adjust=tukey, in SAS 9.4. The documentation under 'LSMEANS Statement' for PROC GLIMMIX states the following: "When you specify ADJUST=TUKEY and your data are unbalanced, PROC GLIMMIX uses the approximation described in Kramer (1956) and identifies the adjustment as "Tukey-Kramer" in the results." In the output, however, the adjustment is labeled as "Holm-Tukey" not "Tukey-Kramer". I can not find any reference for a Holm-Tukey test, and do not understand why it is labeled as such. I need to report in a paper with a reference what test I used for the multiple comparison, and now I am not sure what it should be. Does anyone have any helpful information about this discrepancy? The data are proportion cover of plants in plots, analyzed as repeated measures of a factorial design with fixed and random effects. Here are the first few lines of data (there are six forests, with 32 to 40 plots per forest, so it is unbalanced). data cover; input forest $ plot $ fencing $ miviadd $ year cover; LBAL b1 fence addmivi 2012 0.11 LBAL b2 fence nomivi 2012 0.18 LBAL b3 nofence addmivi 2012 0.13 LBAL b4 nofence nomiv 2012 0.10 LBAL b1 fence addmivi 2015 0.13 LBAL b2 fence nomivi 2015 0.22 LBAL b3 nofence addmivi 2015 0.11 LBAL b4 nofence nomiv 2015 0.09 Here is the code for the LSMEANS step: lsmeans forest*fencing*year / slicediff=(forest*year) CL ilink adjust=tukey plot=meanplot; Here is the top of the output for this section: Thanks!
... View more