BookmarkSubscribeRSS Feed
Choose Language Hide Translation Bar
Thierry_S
Community Trekker

How to focus LSMeans Comparisons in the Fit Model report for an interaction term with large number of levels?

I have analyzed a biomarker (LOG2 FC P) data set using a Repeat Measure linear fit model as shown below

Model Setup Blocked.png

Now, I would like to compare the treatments (LEAD_TAG) at each time point (Visit_Description) without generating all possible comparison pairs produced by the LSMEans Differences Tukey HSD method (see graph below):

  1. is it possible to do so?
  2. Does it matter? (i.e. is there a "cost" for comparing all pairs [36 comparisons] instead of the desired ones [6 comparisons])

LSMeans Plot for Interaction term.png

Thierry R. Sornasse
0 Kudos
1 ACCEPTED SOLUTION

Accepted Solutions
julian
Staff

Re: How to focus LSMeans Comparisons in the Fit Model report for an interaction term with large number of levels?

Hi @Thierry_S ,

There is a cost in using the Tukey HSD pairwise method since the correction applied to keep your overall alpha at the stated level takes into account the number of comparisons. If you need to look at only 3, you're essentially paying with statistical power by correcting for 36. 

 

I would use one of the LSMeans Contrast options, right next toTukey HSD in the Red Triangle for your interaction between treatment and timepoint. You could either specify the linear contrasts yourself by selecting LSMeans Contrast three times, or by using "Test Slices," which will automatically test each factor at a slice along the other factor. In your case, that means you will get 6 tests: three tests of time, one for each treatment, and three tests of treatment, one for each time. It's the last three you want, and they'll be labeled accordingly (baseline, week 2, week 8), and they test whether you have evidence there is an effect of treatment at the given slice of time. 

 

Note that these tests will not correct for multiple comparisons -- the p-values reported are unadjusted.  With three planned comparisons, a Bonferroni correction will likely be most useful.

 

@julian  

1 REPLY 1
julian
Staff

Re: How to focus LSMeans Comparisons in the Fit Model report for an interaction term with large number of levels?

Hi @Thierry_S ,

There is a cost in using the Tukey HSD pairwise method since the correction applied to keep your overall alpha at the stated level takes into account the number of comparisons. If you need to look at only 3, you're essentially paying with statistical power by correcting for 36. 

 

I would use one of the LSMeans Contrast options, right next toTukey HSD in the Red Triangle for your interaction between treatment and timepoint. You could either specify the linear contrasts yourself by selecting LSMeans Contrast three times, or by using "Test Slices," which will automatically test each factor at a slice along the other factor. In your case, that means you will get 6 tests: three tests of time, one for each treatment, and three tests of treatment, one for each time. It's the last three you want, and they'll be labeled accordingly (baseline, week 2, week 8), and they test whether you have evidence there is an effect of treatment at the given slice of time. 

 

Note that these tests will not correct for multiple comparisons -- the p-values reported are unadjusted.  With three planned comparisons, a Bonferroni correction will likely be most useful.

 

@julian