Bioequivalence and Bioavailability Forum • Fixed effects model: changing the F, p values

Fixed effects model: changing the F, p values [Two-Stage / GS Designs]

posted by zizou – Plzeň, Czech Republic, 2016-10-11 03:06 (2747 d 20:02 ago) – Posting: # 16713
Views: 34,989

Dear mittyri,
for the results consistency (R versus Phoenix as you reported), there should be one more term of marvelous effect x) ... you posted it as number 7: (I resorted it for comparison of F and p values, note: int and Trt are missing in the results from R code)


❝   Dependent Hypothesis Numer_DF Denom_DF        F_stat     P_value
❝ 5   Ln(Var)  Group*Per        2       14 3.2110481e+00 0.071153125
❝ 7   Ln(Var)  Trt*Group        1       14 2.3942317e-01 0.632200219
❝ 2   Ln(Var)      Group        1       14 9.9615453e-01 0.335181899
❝ 3   Ln(Var)        Seq        1       14 9.7094813e-01 0.341167964
❝ 4   Ln(Var)  Group*Seq        1       14 5.7364354e-02 0.814182229
❝ 1   Ln(Var)        int        1       14 3.1047428e+03 0.000000000

❝ 6 Ln(Var) Trt 1 14 3.4404496e+00 0.084798088

If you change the muddle:
muddle <- lm(log(Var)~Group+Seq+Seq*Group+Subj%in%(Seq)+Per%in%Group+Trt+Group*Trt, data=data)
and run the mentioned code... Voilà!

Type III sum of squares:

(mimicing SAS/SPSS but with Seq tested against Subj)



log(Var) ~ Group + Seq + Seq * Group + Subj %in% (Seq * Group) + ...

               Df  Sum Sq   Mean Sq  F value     Pr(>F)    

Group:Per       2 0.04514 0.0225681  3.21105   0.071153 .  

Group:Trt       1 0.00168 0.0016827  0.23942   0.632200    

Group           1 0.22546 0.2254640  0.79946   0.386368    

Seq             1 0.27383 0.2738263  0.97095   0.341168    

Group:Seq       1 0.01618 0.0161779  0.05736   0.814182    

Group:Seq:Subj 14 3.94827 0.2820195 40.12641 7.5802e-09 ***

Residuals      14 0.09840 0.0070283                        

---

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

The comparison is telling us that with groups it will be little bit more complicated.

❝ Now Group factor seems to be more suitable.

❝ What bothers me is that Group effect is not equal in R and PHX, cannot figure out the reason

Good luck!
zizou

Complete thread:

Potvin C in the EU Helmut 2013-04-16 17:40 [Two-Stage / GS Designs]
- Potvin C in the EU ElMaestro 2013-04-17 02:31
  - Potvin C in the EU Helmut 2013-04-17 12:23
    - 2 Groups model FDA d_labes 2013-04-17 13:10
      - 2 Groups model FDA Helmut 2013-04-17 16:37
        
        Group effects obsolete? d_labes 2013-04-18 09:55
        
        Group effects FDA/EMA Helmut 2013-04-19 14:34
        
        Group effects FDA/EMA mittyri 2014-08-21 15:34
        
        Group effects FDA/EMA ElMaestro 2014-08-21 17:31
        
        Group effects EMA mittyri 2014-10-01 11:25
        
        Group effects EMA ElMaestro 2014-10-01 11:31
        Group effects EMA Helmut 2014-10-01 13:24
        
        Group effects EMA VStus 2016-10-06 22:06
        
        Group effects EMA ElMaestro 2016-10-07 12:45
        
        Group effects EMA mittyri 2016-10-07 15:06
        
        Fixed effects model with Group term mittyri 2016-10-09 13:27
        
        Fixed effects model: changing the F, p values mittyri 2016-10-10 18:33
        
        Fixed effects model: changing the F, p values ElMaestro 2016-10-10 20:10
        
        Fixed effects model: changing the F, p values Helmut 2016-10-11 00:17
        Holy War of type III d_labes 2016-10-11 13:17
        
        Significance of Group effect in Russia: why the type III is so 'important' mittyri 2016-10-11 15:07
        
        Fixed effects model: changing the F, p valueszizou 2016-10-11 01:06
        
        FDA group model in R mittyri 2016-10-11 13:43
        
        FDA group model in R ElMaestro 2016-10-12 01:07
        FDA group model in R VStus 2016-10-17 15:32
        
        FDA group model in R ElMaestro 2016-10-17 18:58
        
        FDA group model in R VStus 2016-10-20 13:57
        
        Fixed effects model with Group term VStus 2016-10-20 13:54
        
        Group effects FDA/EMA zizou 2016-12-31 02:54
        
        Group effects FDA/EMA Helmut 2017-06-03 14:57
        
        2 Groups model FDA hiren379 2013-08-22 14:32
        
        Homework Helmut 2013-08-23 11:53