FDA’s model… [Software]

posted by Helmut Homepage – Vienna, Austria, 2014-10-09 20:34 (3862 d 22:10 ago) – Posting: # 13683
Views: 38,637

Dear Yung-jin,

❝ […] it is not caused by "reaching its max. iterations..." as lme().


Exactly.

❝ Is there any random effect or are all effects fixed in model for this situation? Could you try with all effects fixed to see if it can be improved?


Nope. I happens (only for some nasty datasets) in a partial replicate design (full replicates are OK) if following FDA’s mixed effects model for ABE:

PROC MIXED
  data=pk;
  CLASSES SEQ SUBJ PER TRT;
  MODEL LAUCT = SEQ PER TRT/ DDFM=SATTERTH;
  RANDOM TRT/TYPE=FA0(2) SUB=SUBJ G;
  REPEATED/GRP=TRT SUB=SUJ;
  ESTIMATE 'T vs. R' TRT 1 -1/CL ALPHA=0.1;
  ods output Estimates=unsc1;
  title1 'unscaled BE 90% CI - guidance version';
  title2 'AUCt';
run;

data unsc1;
  set unsc1;
  unscabe_lower=exp(lower);
  unscabe_upper=exp(upper);
run;


The problem lies in the random section where the structure of the covariance structure is specified with TYPE=FA0(2). CSH or FAO(1) helps. I don’t see how to go with ‘all fixed effects’ here…

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,424 posts in 4,927 threads, 1,669 registered users;
181 visitors (0 registered, 181 guests [including 4 identified bots]).
Forum time: 18:44 CEST (Europe/Vienna)

[Those] who have an excessive faith in their theories or in their
ideas are not only poorly disposed to make discoveries, but they
also make very poor observations.    Claude Bernard

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