Bioequivalence and Bioavailability Forum

SASian potpourri [Study As­sess­ment]

posted by d_labes  – Berlin, Germany, 2016-05-24 14:02 (3680 d 04:38 ago) – Posting: # 16358
Views: 39,412

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 .

Complete thread:

• Dose Proportionality and Variance AngusMcLean 2016-05-11 16:55 [Study As­sess­ment] 
  • More information, please Helmut 2016-05-12 14:34
    • More information, please AngusMcLean 2016-05-13 16:40
      • Setup in Phoenix/WinNonlin Helmut 2016-05-14 02:26
        • Setup in Phoenix/WinNonlin AngusMcLean 2016-05-14 18:54
          • Setup in Phoenix/WinNonlin Helmut 2016-05-15 14:47
            • Setup in Phoenix/WinNonlin AngusMcLean 2016-05-15 15:17
              • Phoenix 64 Warning Helmut 2016-05-15 15:56
                • Phoenix 64 Warning AngusMcLean 2016-05-15 20:11
                  • OT: imperial vs. metric units Helmut 2016-05-16 16:26
              • Setup in Phoenix/WinNonlin ElMaestro 2016-05-15 20:54
                • Setup in Phoenix/WinNonlin AngusMcLean 2016-05-15 22:30
        • Setup in Phoenix/WinNonlin AngusMcLean 2016-05-16 21:00
          • NCSS vs. PHX/WNL vs. SAS Helmut 2016-05-17 01:50
            • NCSS vs. PHX/WNL vs. SAS - Validation? mittyri 2016-05-18 08:23
              • Diagnostics ElMaestro 2016-05-18 09:20
                • Diagnostics: R and Phoenix Helmut 2016-05-18 15:14
                  • Diagnostics: R zizou 2016-05-22 19:07
                    • Diagnostics: R Helmut 2016-05-23 01:22
                      • SASian potpourrid_labes 2016-05-24 12:02
                        • Compilation Helmut 2016-05-24 14:27
                          • REML or not d_labes 2016-05-24 16:33
                            • complete or not Helmut 2016-05-24 16:57
                          • Compilation AngusMcLean 2016-05-26 16:46
                            • doubts about NCSS Helmut 2016-05-26 19:13
                              • Doubts about NCSS zizou 2016-05-26 23:38
                          • Compilation AngusMcLean 2016-05-28 00:51
                            • Kenward-Roger? Helmut 2016-05-28 15:59
                          • 90% confidence interval for R_dnm Shuanghe 2019-01-04 17:45
                            • 90% confidence interval for R_dnm d_labes 2019-01-05 14:01
                              • Visualizing lmer and limits mittyri 2019-01-06 17:00
                                • Visualizing lmer and limits Shuanghe 2019-01-07 11:05
                                • Visualizing lmer and limits d_labes 2019-01-07 15:08
                                  • Visualizing lmer and limits mittyri 2019-01-13 23:53
                              • 90% confidence interval for R_dnm Shuanghe 2019-01-07 10:53
                                • 90% confidence interval for R_dnm d_labes 2019-01-07 15:17
                                  • 90% confidence interval for R_dnm Shuanghe 2019-01-07 17:11
                                    • 90% confidence interval for R_dnm d_labes 2019-01-07 18:24
                                    • offtop: greek letters and tables mittyri 2019-01-08 00:19
                                      • OT: greek letters and symbols Helmut 2019-02-02 16:04
              • Smith’s paper Helmut 2016-05-18 14:44
                • Smith’s paper d_labes 2019-01-05 15:00