FDA’s model… [Software]

posted by Helmut Homepage – Vienna, Austria, 2014-10-09 20:34 (3910 d 05:33 ago) – Posting: # 13683
Views: 39,138

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,672 registered users;
290 visitors (0 registered, 290 guests [including 21 identified bots]).
Forum time: 02:07 CEST (Europe/Vienna)

I have had my results for a long time:
but I do not yet know how I am to arrive
at them.    Carl Friedrich Gauß

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