Bogus PROC MIXED/FDA [Software]

posted by d_labes  – Berlin, Germany, 2010-03-25 10:15 (5508 d 16:48 ago) – Posting: # 4970
Views: 11,182

Dear dixit,

I fully agree with EM.

IMHO Proc MIXED stops here with arbitrary values for the inter-inividual variance parameters and also with arbitrary intra-individual variability for Test. An indication is:


                                          Estimated G Matrix


❝ Row    Effect     treatment   subject        Col1        Col2

❝   1    treatment    R             1           0.4529      0.4233   CVinter=75.7%
❝   2    treatment    T             1           0.4233      0.5232   CVinter=82.9%
❝                                                                   
❝         Covariance Parameter Estimates


❝ Cov Parm     Subject    Group          Estimate

❝ ...

❝ Residual     subject    treatment R      0.6289  CVintra=93.6%

❝ Residual     subject    treatment T     0.09081  CVintra=30.8%



Do you really believe that your Test preparation has such a very lower intra-individual variability? If it is the same drug I can't imagine that at all :no:.

Once again said, the FDA model is over-specified for the partial replicate design. And therefore the model fit with Proc MIXED is not reliable.
We had this phenomenon already here were SAS Proc MIXED and WINNONLIN gave totally different values for the covariance parameters.

What to the rescue?

Within a model formulation in Proc MIXED I don't know, as I already said. I have the strong believe (but believe is not The power to know :-D) that it is impossible to get a solution within that because we have a confounding between the subject-by-formulation interaction (an inter-individual term) and the intra-individual variability of the Reference in the design used here.

I would suggest you to go with the so called "Methods of moments".
See

[1]R.J. McNally
Tests for Individual and Population Bioequivalence Using 3-Period Crossover Designs

which can be found here for the necessary statistics within a partial replicate design.
They can be used also in the framework of ABE and scaled ABE. Your intention I guess?

Or neglect any subject-by-formulation interaction and go with the classical model (same as for a 2x2 cross-over) to obtain the ABE 90% confidence intervals and estimate independently the intra-individual variability of the Reference using the intra-subject contrasts (R-R')/sqrt(2).

Good luck.
BTW: Seems no m on your keyboard :-).

Regards,

Detlew

Complete thread:

UA Flag
Activity
 Admin contact
23,424 posts in 4,927 threads, 1,670 registered users;
18 visitors (0 registered, 18 guests [including 1 identified bots]).
Forum time: 04:04 CEST (Europe/Vienna)

Most scientists today are devoid of ideas, full of fear, intent on
producing some paltry result so that they can add to the flood
of inane papers that now constitutes “scientific progress”
in many areas.    Paul Feyerabend

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