Multiple comparison

[elmaghraby_m]

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.

[PGStats]

It might be simpler to output all uncorrected comparison results with ODS OUTPUT, eliminate unwanted comparisons and obtain the corrected p-values with PROC MULTTEST.