FDA’s model… [Software]

posted by Helmut Homepage – Vienna, Austria, 2014-10-09 20:34 (4270 d 02:04 ago) – Posting: # 13683
Views: 42,601

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,654 posts in 4,992 threads, 1,570 registered users;
195 visitors (0 registered, 195 guests [including 19 identified bots]).
Forum time: 22:38 CEST (Europe/Vienna)

“Data! Data! Data!” he cried impatiently.
“I can’t make bricks without clay!”    Arthur Conan Doyle

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