Donald’s model [RSABE / ABEL]

posted by Helmut Homepage – Vienna, Austria, 2020-08-14 11:43 (44 d 01:04 ago) – Posting: # 21875
Views: 2,135

Dear zizou,

» in this post I mentioned the oldest study with partial replicate design (sequences RRT,RTR,TRR) I know about.

THX! Funky model by Donald Schuirmann:

PROC MIXED;
CLASS SUBJ SEQ PER GRP TRT;
MODEL Y =
  GRP SEQ GRP*SEQ PER GRP*PER TRT/DDFM=SATTERTH;
RANDOM SUBJ(GRP*SEQ) SUBJ*TRT(GRP*SEQ);
ESTIMATE 'T VS.R' TRT -1 1/CL ALPHA=0.1;
RUN;

Two things are interesting: No REPEATED statement and the covariance structure is not specified. Hence, SAS applies the default TYPE=VC (variance components). Let’s see what happens in PHX/WNL with the EMA’s data set II (of course, group terms deleted from Donald’s setup).
  1. FDA’s code (2001 BE-stats guidance and 2010 progesterone guidance)
    PE 102.26% (90% CI 97.05–107.76%)
    Warning 11091: Newton's algorithm converged with modified Hessian. Output is suspect.

  2. As above but TYPE=FA0(1) instead of TYPE=FA0(2)
    PE 102.26% (90% CI 97.05–107.76%)
  3. Donald
    PE 102.26% (90% CI 94.37–110.82%)
  4. EMA Method B (implicitly DDFM=CONTAIN)
    PE 102.26% (90% CI 97.32–107.46%)
    Same results with DDFM=SATTERTH cause the data set is complete and balanced.
I think that Donald’s model doesn’t make sense cause the replicative structure of the study is ignored (given, like in the EMA’s Method A and B).

Dif-tor heh smusma 🖖
Helmut Schütz
[image]

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

Complete thread:

Activity
 Admin contact
21,077 posts in 4,395 threads, 1,469 registered users;
online 3 (1 registered, 2 guests [including 1 identified bots]).
Forum time: Sunday 12:48 CEST (Europe/Vienna)

Politicians use statistics like drunkards use lampposts:
not for illumination, but for support.    attributed to Hans Kuhn

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