Bioequivalence and Bioavailability Forum • Sample data & results

Sample data & results [Software]

posted by Eva – 2008-05-06 19:42 (6207 d 08:24 ago) – Posting: # 1824
Views: 23,322

Dear Ohlbe and Helmut,

thanks for your fast replies!
Here is a sample data set of 8 subjects, maybe this can help to illustrate my problem...


Subject Period  Sequence Treatment       AUC

------------------------------------------------

1       1       TR       Test            3981.76

2       1       TR       Test            5573.62

3       1       TR       Test            5821.52

4       1       TR       Test            8558.99

5       1       RT       Reference       3958.93

6       1       RT       Reference       5145.91

7       1       RT       Reference       689.6

8       1       RT       Reference       3812.02

1       2       TR       Reference       4446.05

2       2       TR       Reference       4847.12

3       2       TR       Reference       4274.72

4       2       TR       Reference       4839.49

5       2       RT       Test            3219.92

6       2       RT       Test            4587.32

7       2       RT       Test            4436.7

8       2       RT       Test            6168.8

My SAS results (PROC GLM) for these data are


EFFECT             DF    TYPE III SS    MEAN SQUARE     F-VALUE     P-VALUE

---------------------------------------------------------------------------

Sequence           1     0.6105         0.6105          2.46        0.1678

Subject(Sequence)  6     1.5739         0.2623          1.06        0.4739

Treatment          1     0.5365         0.5365          2.16        0.1918

Period             1     0.0774         0.0774          0.31        0.5966

Error              6     1.4887         0.2481

and this is what WinNonlin calculated:


EFFECT             DF    TYPE III SS    MEAN SQUARE     F-VALUE     P-VALUE

---------------------------------------------------------------------------

Sequence           1     0.6105         0.6105          2.33        0.1780

Subject(Sequence)  6     1.5739         0.2623          1.06        0.4739

Treatment          1     0.5365         0.5365          2.16        0.1918

Period             1     0.0774         0.0774          0.31        0.5966

Error              6     1.4887         0.2481

Dear Ohlbe,
I think I have picked the right F- and p-values from my output, I am using PROC GLM (and I just checked, the output is identical with PROC ANOVA) and picked the Sequence, Subject(Sequence), Treatment and Period results from the bottom part and the Error results from the top part.

Dear Helmut,
could this be a problem with my WNL version? (I noticed that you used v5.2, whereas I have v5.1.1)
As you can see above, the sums of squares are not the problem, it is only the F- and p-values that differ.
I understood that the F-value is the ratio of the effect mean square and the mean square error. This is correct for the SAS results (sequence: 0.6105/0.2481=2.46), but not for WinNonLin...

Any ideas?
Thanks again, for your help and for sharing your precious time with my "little" problem...

Kind regards,
Eva

Edit: Treatment in the table changed to Test to avoid confusion with the heading. [Helmut]

Complete thread:

SAS vs. WinNonlin: different sequence effect results Eva 2008-05-06 14:04 [Software]
- SAS vs. WinNonlin: different sequence effect results Ohlbe 2008-05-06 15:56
- imbalanced design? Helmut 2008-05-06 16:13
  - Sample data & resultsEva 2008-05-06 17:42
    - Error factor Ohlbe 2008-05-06 19:25
      - Error tests between/within Helmut 2008-05-06 21:33
        
        Error tests between/within Helmut 2008-05-07 18:18
        The power to know d_labes 2008-05-08 11:16
        
        The power to know Helmut 2008-05-08 17:18
        
        The power to know d_labes 2008-05-09 09:35
        
        The power to know Nirali 2008-05-09 11:00
        
        The power to know d_labes 2008-05-16 08:54
        
        The power to know kevan 2009-05-25 15:46
        
        Bogus statement for 2,2,2-BE ElMaestro 2009-05-25 22:25
        
        Bogus? What? d_labes 2009-05-27 08:57
        
        Linear model on steroids ElMaestro 2009-05-28 19:12
        
        Fixed nowadays what? d_labes 2009-05-27 09:03
        
        Kinetica 5.0 bug Helmut 2008-12-31 16:42