I have data that is classified into 3 times (1-3) in each of 4 weeks (1-4) in a repeated measures design (5 subjects). (Total 3 times x 4 weeks = 12 mean values throughout the experiment)
I want to present lsmeans of the three times compared within each week.
For Bonferroni correction, which is correct; dividing p=0.05/3 = 0.0167 (three comparisons within a week; time1 vs. time2; time1 vs. time3; time2 vs. time3). Or, p=0.05/66 = 0.00076 (all possible comparison, which is not my aim).
The code I wrote was:
PROC MIXED DATA=HORMONES; CLASS BUCK WK TIME;
MODEL Trig = WK TIME WK*TIME;
REPEATED / SUBJECT=BUCK TYPE=AR(1);
LSMEANS WK*TIME/PDIFF; RUN;
Thank you.