Compilation [Study Assessment]
Dear all,
in R /
My preference is REML/Satterthwaite because one could reproduce results in three different software packages. I would avoid NCSS; no idea how the calculation is done.
in R /
lmerTest()
one can get maximum likelihood estimation by setting the argument REML=FALSE
in the model statement (default is REML=TRUE
). Satterthwaite’s DF are 10.858 for the intercept and 6.884 for the slope. In lme()
use method="ML"
(default is method="REML"
). Below a compilation (in analogy to BE rounded to four decimal figures and grouped by “similarity”):Software Method Estimate 90% CI
═══════════════════════════════════════════════════
SAS (Smith) ML/Satter. B0 1.94 1.54 2.35
B1 0.7615 0.679 0.844
SAS ML/Satter. B0 1.9425 1.5391 2.3458
B1 0.7615 0.6789 0.8441
R lmerTest ML/Satter. B0 1.9425 1.5391 2.3458
B1 0.7615 0.6789 0.8441
───────────────────────────────────────────────────
R lme ML B0 1.9425 1.5175 2.3675
B1 0.7615 0.6738 0.8491
───────────────────────────────────────────────────
SAS ML/contain. B0 1.9414 1.5175 2.3675
B1 0.7615 0.6738 0.8491
───────────────────────────────────────────────────
SAS REML/Satter. B0 1.9414 1.4968 2.3860
B1 0.7617 0.6696 0.8539
PHX/WNL REML/Satter. B0 1.9414 1.4968 2.3860
B1 0.7617 0.6696 0.8539
R lmerTest REML/Satter. B0 1.9414 1.4968 2.3860
B1 0.7617 0.6696 0.8539
───────────────────────────────────────────────────
NCSS method? B0 1.9414 1.4849 2.3978
B1 0.7617 0.6659 0.8576
───────────────────────────────────────────────────
SAS ML/contain. B0 1.9414 1.4807 2.4020
B1 0.7617 0.6665 0.8570
R lme REML B0 1.9414 1.4807 2.4020
B1 0.7617 0.6665 0.8570
───────────────────────────────────────────────────
PHX/WNL REML/resid. B0 1.9414 1.4515 2.4313
B1 0.7617 0.6665 0.8570
My preference is REML/Satterthwaite because one could reproduce results in three different software packages. I would avoid NCSS; no idea how the calculation is done.
—
Dif-tor heh smusma 🖖🏼 Довге життя Україна!![[image]](https://static.bebac.at/pics/Blue_and_yellow_ribbon_UA.png)
Helmut Schütz
![[image]](https://static.bebac.at/img/CC by.png)
The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
Dif-tor heh smusma 🖖🏼 Довге життя Україна!
![[image]](https://static.bebac.at/pics/Blue_and_yellow_ribbon_UA.png)
Helmut Schütz
![[image]](https://static.bebac.at/img/CC by.png)
The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
Complete thread:
- Dose Proportionality and Variance AngusMcLean 2016-05-11 16:55 [Study Assessment]
- 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
- Phoenix 64 Warning AngusMcLean 2016-05-15 20:11
- Setup in Phoenix/WinNonlin ElMaestro 2016-05-15 20:54
- Setup in Phoenix/WinNonlin AngusMcLean 2016-05-15 22:30
- Phoenix 64 Warning Helmut 2016-05-15 15:56
- Setup in Phoenix/WinNonlin AngusMcLean 2016-05-15 15:17
- Setup in Phoenix/WinNonlin Helmut 2016-05-15 14:47
- 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 potpourri d_labes 2016-05-24 12:02
- CompilationHelmut 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
- doubts about NCSS Helmut 2016-05-26 19:13
- 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
- 90% confidence interval for R_dnm Shuanghe 2019-01-07 17:11
- 90% confidence interval for R_dnm d_labes 2019-01-07 15:17
- Visualizing lmer and limits mittyri 2019-01-06 17:00
- 90% confidence interval for R_dnm d_labes 2019-01-05 14:01
- REML or not d_labes 2016-05-24 16:33
- CompilationHelmut 2016-05-24 14:27
- SASian potpourri d_labes 2016-05-24 12:02
- Diagnostics: R Helmut 2016-05-23 01:22
- Diagnostics: R zizou 2016-05-22 19:07
- Diagnostics: R and Phoenix Helmut 2016-05-18 15:14
- Smith’s paper Helmut 2016-05-18 14:44
- Smith’s paper d_labes 2019-01-05 15:00
- Diagnostics ElMaestro 2016-05-18 09:20
- NCSS vs. PHX/WNL vs. SAS - Validation? mittyri 2016-05-18 08:23
- NCSS vs. PHX/WNL vs. SAS Helmut 2016-05-17 01:50
- Setup in Phoenix/WinNonlin AngusMcLean 2016-05-14 18:54
- Setup in Phoenix/WinNonlin Helmut 2016-05-14 02:26
- More information, please AngusMcLean 2016-05-13 16:40
- More information, please Helmut 2016-05-12 14:34