Bioequivalence and Bioavailability Forum

ANOVA of Model 1 [General Sta­tis­tics]

posted by Helmut  – Vienna, Austria, 2024-10-04 12:11 (58 d 09:32 ago) – Posting: # 24215
Views: 1,354

Hi BEQool,

you have a point!

anova1 <- anova(model1)
num    <- as.numeric(anova1["group:treatment", c(3, 1)])
denom  <- as.numeric(anova1["group:sequence:subject", c(3, 1)])
cat(pf(num[1] / denom[1], num[2], denom[2], lower.tail = FALSE), "\n")

Gives for the first símulation above:

0.1501297

However, I’m puzzled because then p(G×T) of our meta study are not uniformly distributed any more.



Compare them to these plots.

—
Dif-tor heh smusma 🖖🏼 Довге життя Україна! 
Helmut Schütz


The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes

Complete thread:

• Group-by-Treatment Interaction: the final word? Helmut 2024-04-19 01:26 [General Sta­tis­tics] 
  • Group-by-Treatment Interaction: FDA Helmut 2024-09-20 13:30
  • Group-by-Treatment Interaction: the final word? BEQool 2024-10-03 15:09
    • ANOVA of Model 1 Helmut 2024-10-03 15:33
      • ANOVA of Model 1 BEQool 2024-10-03 16:08
        • ANOVA of Model 1Helmut 2024-10-04 10:11
          • Between vs Within mittyri 2024-10-05 00:08
            • Between vs Within BEQool 2024-10-06 20:06
              • Treatment-by-Group in PROC MIXED and PROC GLM? mittyri 2024-10-06 22:00
                • Treatment-by-Group in PROC MIXED and PROC GLM? BEQool 2024-10-07 10:27
                  • Treatment-by-Group in PROC MIXED and PROC GLM? mittyri 2024-10-07 13:11
                    • Treatment-by-Group in PROC MIXED and PROC GLM? BEQool 2024-10-14 11:01