Treatment-by-Group in PROC MIXED and PROC GLM? [General Sta­tis­tics]

posted by mittyri – Russia, 2024-10-07 00:00 (65 d 15:51 ago) – Posting: # 24219
Views: 1,362

(edited on 2024-10-07 00:31)

Hi BEQool,

❝ So you would assume that Residual Mean Squares (MS) as the denominator was used because nothing is explicitly stated? If Subject MS was used instead as the denominator, you think that they would probably mention it?

yes.
The regulators are not trying to make our life easier. Yes, the model to be used is rarely presented. There are not so many exceptions (i.e. replicate designs for EMA or FDA).

❝ Is there any source or literature what to use as denominator when assessing Group*Treatment interaction? I cant find anything relevant. Nothing is mentioned about it in the new ICH M13A guideline

You are right. There's nothing. Note that when FDA published the code for 2*2*2 studies (a while ago), they did not mention in the text that a sequence should be tested against Subject MS, but just added a RANDOM statement in PROC GLM. Anyway, this decision had some reason.

In any sources with Group*Treatment interaction the denominator is not mentioned directly, but I think that without additional RANDOM statement (or ANOVA postprocessing) which SHOULD be mentioned if in use, the Residual MS is still valid. So you should prove why you modified your ANOVA to get F against Subject MS. The `famous` Sun et al. * paper states:

For each dataset, a Linear Mixed Model is fit with subgroup, sequence, subgroup-by-sequence, period nested within subgroup, treatment, treatment-by-subgroup as the fixed model and subject nested within subgroup and sequence as a random effect model (the same model can also be fit using the General Linear Model (GLM) with fixed effects only but with slight modification).


❝ Additionally, what if the significance between the denominator used differ?

yes, that's why it should be mentioned if the model/ANOVA was modified to use Subject MS




PS: Are you equipped with SAS to check the ANOVA tables for this model (using some dataset without dropouts?) It would be interesting to see what PROC MIXED thinks about this interaction. From this paper I assume they should be the same as PROC GLM. But who knows! I could be wrong with my claim...

Kind regards,
Mittyri

Complete thread:

UA Flag
Activity
 Admin contact
23,336 posts in 4,902 threads, 1,666 registered users;
37 visitors (1 registered, 36 guests [including 8 identified bots]).
Forum time: 14:51 CET (Europe/Vienna)

I’m all in favor of the democratic principle
that one idiot is as good as one genius, but I draw the line
when someone takes the next step and concludes
that two idiots are better than one genius.    Leo Szilard

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC Ing. Helmut Schütz
HTML5