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

posted by Helmut Homepage – Vienna, Austria, 2013-08-07 17:13 (4345 d 09:48 ago) – Posting: # 11235
Views: 10,178

Dear Detlew!

❝ Will try it the next time numerically.


Yes please. I tried John’s nasty data set (log-transformed data in full precision) where we ran into convergence problems and got different results for RSABE and ABE (we tried SAS and PHX; complete data set and after exclusion of incomplete data).
When I set up FA0(1) in PHX’ LME I got…
ERROR 11021: Variance matrix is not positive definite in QA.
… which is an indication that the initial variance estimates (by methods of moments) don’t “work”. I tried different starting values to no avail.

Now for the surprise: PHX’ BE-model (which supposedly sits ‘on top’ of LME) converged without a warning. I got

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

… matching Shuanghe’s and your’s in SAS with FDA’s FA0(2)-parameterization.
After exclusion of incomplete data (subjects 8, 14, 17, and 31)

PE  91.2423% (90% CI: 74.0705–112.3950%) s²WR 0.40652 (CVWR 70.82%)

… different to yours’.
Amazing.

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,425 posts in 4,928 threads, 1,681 registered users;
40 visitors (0 registered, 40 guests [including 8 identified bots]).
Forum time: 03:01 CEST (Europe/Vienna)

Fools with Tools Are Still Fools!    Anonymous

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