SASian potpourri [Study As­sess­ment]

posted by d_labes  – Berlin, Germany, 2016-05-24 14:02 (2855 d 20:12 ago) – Posting: # 16358
Views: 29,056

Dear All!

❝ ... Would be great if one of our SASians could jump in.


Here we go:

Software               Estimate       90% CI
NCSS               B0 1.9414    1.4849    2.3978   
                   B1 0.7617    0.6659    0.8576   
Phoenix/WinNonlin  B0 1.9413858 1.4967592 2.3860125
                   B1 0.7617406 0.6695783 0.8539029
R lmerTest         B0 1.9413858 1.4967584 2.3860132
                   B1 0.7617406 0.6695768 0.8539044
SAS (ML)           B0 1.9413856 1.5174574 2.3675030
                   B1 0.7614780 0.6738351 0.8491210
SAS (REML)         B0 1.9413856 1.4807360 2.4020351
                   B1 0.7617407 0.6664695 0.8570118
Smith et al.       B0 1.94      1.54      2.35     
                   B1 0.7615    0.679     0.844    
SAS (ML/satterth)  B0 1.9425    1.5391    2.3458   
                   B1 0.7615    0.6789    0.8441   

Smith et al. used ML (usual maximum likelihood) as estimation method! Degrees of freedom method is per default "containment" (7 for the intercept, 5 for the regression slope). Astonishing enough I couldn't reproduce their results w.r.t. to the 90% CIs to a sufficient degree of accuracy with this ddfm method, at least sufficient for me.
Choosing ddfm=SATTERTHWAITE gives the desired results.

SAS code:
Proc mixed data=dp method=ML;
  class subject;
  model lnCmax=lnd /s cl alpha=0.1 ddfm=SAT;
  random intercept /subject=subject;
run;



edit: the missing REML+SATTERTH
SAS (REML/satterth) B0 1.9414    1.4968    2.3860     
                    B1 0.7617    0.6696    0.8539   

Bingo! Same as R lmer() and Phoenix/WinNonlin, at least with 4 decimals, SAS default output format. Too lazy to tease out more numbers :cool:.

Regards,

Detlew

Complete thread:

UA Flag
Activity
 Admin contact
22,940 posts in 4,812 threads, 1,639 registered users;
45 visitors (0 registered, 45 guests [including 9 identified bots]).
Forum time: 09:14 CET (Europe/Vienna)

Those people who think they know everything
are a great annoyance to those of us who do.    Isaac Asimov

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