overspecified model? [General Sta­tis­tics]

posted by Helmut Homepage – Vienna, Austria, 2011-06-22 15:36 (5472 d 00:43 ago) – Posting: # 7159
Views: 7,300

Dear RB!

❝ I have a 3 period, 2 treatment crossover (TRR, RTT) design and used the following code (as per FDA guidance) to analyse the data:


FDA 2001 Apendix E


❝ The above analysis resulted in the following note in the SAS log: "Estimated G matrix not positive definite" suggesting that "one or more variance components on the RANDOM statement is/are estimated to be zero and could/should be removed from the model."


Yes, old story. SAS complains about an over-specified model (T is not repeated) and gives a variance for T. ;-) See this thread.

❝ My colleague is suggesting I switch the covariance structure to VC however I'm concerned that this negates the REPEATED statement and is questionable in a 3x2 design. What is your opinion on this point?


Agree.

❝ I used the CSH structure but, as expected, encountered the same notes in the log.


Yep.

❝ Looking through the covariance structures, Compound Symmetry (CS) looks more appropriate as it accounts for correlation & assumes it is constant regardless of the lag between pairs of repeated measurements (sounds OK given the study design).


Interesting idea. At least you don’t get a warning any more. :-D
Tried it with EMA’s example data set II (TRR|RTR|RRT) and got (Phoenix 6.2 = poor man’s SAS):

                         90% CI         CVWR   CVWT
EMA's Model A         97.32  107.46
EMA’s Model B         97.32  107.46
EMA’s crippled model       NA          11.20   NA
FDA FA0(2)            97.05  107.76    11.55   8.65 Warning (SAS: CVWT=3.87)
FDA CSH               not estimable     8.26   7.15 Warning
FDA CS                97.05  107.76    11.55  18.67 no Warning!


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

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

Complete thread:

UA Flag
Activity
 Admin contact
23,654 posts in 4,992 threads, 1,570 registered users;
102 visitors (0 registered, 102 guests [including 14 identified bots]).
Forum time: 16:20 CEST (Europe/Vienna)

The idea is to try and give all the information to help others
to judge the value of your contribution;
not just the information that leads to judgment
in one particular direction or another.    Richard Feynman

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