Mixed interest [General Statistics]
❝ ...
❝ The Mixed Procedure
❝
❝ Iteration History
❝
❝ Iteration Evaluations -2 Res Log Like Criterion
❝ 5 1 272.91032497 0.00000001
OK, I end up after 7 iterations at:
-2* REML log(likelihood) 251.724
Same result if use the default convergence criterion (10-10), yours ((10-8), or approach numeric resolution.
❝ Covariance Parameter Estimates
❝
❝ Cov Parm Subject Group Estimate
❝
❝ FA(1,1) Subject 0.5756
❝ FA(2,1) Subject 0.4485
❝ FA(2,2) Subject 0.3981
❝ Residual Subject Treatment R 0.1862 (CV=0.4524)
❝ Residual Subject Treatment T 0.3269
Closer match than we had in our previous example...
Final variance parameter estimates:
lambda(1,1)_11 0.575597
lambda(1,2)_11 0.448489
lambda(2,2)_11 0.355511
Var(PERIOD*TREATMEN*SUBJECT)_21 0.186160
Var(PERIOD*TREATMEN*SUBJECT)_22 0.359059
❝ And just more interesting: Analysis according to D. Brown (EMA)...
❝ ...
❝ The GLM Procedure
❝
❝ Dependent Variable: ln_Cmax ln_Cmax
❝
❝ Sum of
❝ Source DF Squares Mean Square F Value Pr > F
❝
❝ Model 45 38.86413212 0.86364738 4.73 <.0001
❝ Error 42 7.66318204 0.18245672
❝ Corrected Total 87 46.52731416
Confirmed.
Total Observations : 88
Observations Used : 88
Residual SS : 7.66318204
Residual df : 42
Final variance parameter estimates:
Var(Residual) 0.18245672
❝ Very near to the MIXED results. But maybe this coincidence is by chance.
Who knows?
❝ But never trust estimates with questionable convergence in REML (WINNONLIN
❝ definitely warns at least "Output is suspect", R's lme() will throw an error if
❝ not converged and will not give you any parameter estimate or will not give
❝ any CI for the covariance parameters if VarCov is not positive definite, but
❝ the incredible tells all estimates
❝ as if nothing happens ).
Well, if I ask PHX/WNL for intermediate results (besides many pages of Matrazen aka matrices) I got an additional
Warning 11090: Asymptotic covariance matrix not computed. Information matrix is deemed singular.
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:
- SABE Reference Variability and CI Computation preyes323 2010-07-12 17:12 [General Statistics]
- Dataset Helmut 2010-07-12 17:26
- Dataset preyes323 2010-07-12 18:20
- Data format Helmut 2010-07-12 18:38
- Data format preyes323 2010-07-13 01:16
- SAS and/or Phoenix/WinNonlin-experts around? Helmut 2010-07-13 02:29
- Phoenix/WinNonlin-experts around? d_labes 2010-07-13 12:15
- Phoenix/WinNonlin-experts around? Helmut 2010-07-13 13:42
- Mixed interest d_labes 2010-07-13 15:11
- Mixed interestHelmut 2010-07-13 19:19
- lme() answer and beyond ... d_labes 2010-07-14 10:41
- lme() answer and beyond ... Helmut 2010-07-14 13:47
- AIC, BIC and that all ... d_labes 2010-07-15 12:01
- lme() answer and beyond ... Helmut 2010-07-14 13:47
- lme() answer and beyond ... d_labes 2010-07-14 10:41
- Mixed interestHelmut 2010-07-13 19:19
- Mixed interest d_labes 2010-07-13 15:11
- Phoenix/WinNonlin-experts around? Helmut 2010-07-13 13:42
- Phoenix/WinNonlin-experts around? d_labes 2010-07-13 12:15
- To Err is Human d_labes 2010-07-13 11:51
- To Err is Human, but... Helmut 2010-07-13 13:45
- ... to Arr is Pirate d_labes 2010-07-13 15:47
- Brilliant page!!! ElMaestro 2010-07-14 21:03
- ... to Arr is Pirate d_labes 2010-07-13 15:47
- To Err is Human preyes323 2010-07-14 14:45
- Heads up! Helmut 2010-07-14 15:17
- To Err is Teacher d_labes 2010-07-15 11:05
- To Err is Human preyes323 2010-07-16 02:28
- To Err is Human d_labes 2010-07-16 10:42
- To Err is Human preyes323 2010-07-17 11:14
- Regulatory constants d_labes 2010-07-19 10:00
- To Err is Human preyes323 2010-07-17 11:14
- To Err is Human d_labes 2010-07-16 10:42
- To Err is Human, but... Helmut 2010-07-13 13:45
- Data format ElMaestro 2011-02-06 11:54
- SAS System Viewer Helmut 2011-02-06 13:02
- SAS and/or Phoenix/WinNonlin-experts around? Helmut 2010-07-13 02:29
- Data format preyes323 2010-07-13 01:16
- Data format Helmut 2010-07-12 18:38
- Dataset preyes323 2010-07-12 18:20
- Dataset Helmut 2010-07-12 17:26