Bioequivalence and Bioavailability Forum

FA0(1) for partial replicate in SAS [RSABE / ABEL]

posted by d_labes  – Berlin, Germany, 2013-08-08 11:42 (4297 d 13:05 ago) – Posting: # 11246
Views: 9,970

Dear Helmut!

With John’s nasty data set I got:
- all data

PE= 101.0213% (90% CI: 80.7733–126.3450%) s²WR= 0.4210 (CVWR= 72.35%)

… matching yours and the former evaluation with FA0(2) to full precision.
- complete cases only

PE=  90.5032% (90% CI: 73.5610–111.3474%) s²WR= 0.4001 (CVWR= 70.14%)

… different from yours but matching nearly to full precision the former FA0(2) evaluation (differences in the CI only after the second decimal).

No hint of any convergence problems, but in both cases in the log
NOTE: Estimated G matrix is not positive definite.
Not quite sure if this is a problem.

What puzzles me a little bit is that we have an estimate of s²_WT in the model (via the repeated/grp=tmt sub=subject; statement) also Wolfgang claimed, correctly I think, that "... s²_b,T and s²_w,T cannot be estimated, only their sum." But have to think about that ...

Complete thread:

• RSAB model overspecified Wolfgang 2013-08-05 14:37 
  • RSAB model overspecified jag009 2013-08-05 15:35
  • Proc Mixed FA0(1) for partial replicate d_labes 2013-08-05 15:59
    • Proc Mixed FA0(1) for partial replicate Wolfgang 2013-08-05 16:24
      • Proc Mixed FA0(1) for partial replicate ElMaestro 2013-08-05 17:19
        • Proc Mixed FA0(1) for partial replicate Wolfgang 2013-08-06 08:35
          • Proc Mixed FA0(1) for partial replicate ElMaestro 2013-08-06 10:34
      • Proc Mixed FA0(1) for partial replicate d_labes 2013-08-06 10:30
        • FA0(1) for partial replicate in Phoenix Helmut 2013-08-07 15:13
          • FA0(1) for partial replicate in SASd_labes 2013-08-08 09:42
          • FA0(1) for partial replicate in Phoenix lhughes 2013-08-30 15:17
            • FA0(1) in PHX: works (& more questions) Helmut 2013-08-30 17:07