Chapter 5 : Mouse feeding example |
Analysis of Variance with PROC MIXED |
Model Information | |
---|---|
Data Set | WORK.MOUSEDIET |
Dependent Variable | LIFETIME |
Covariance Structure | Diagonal |
Estimation Method | REML |
Residual Variance Method | Profile |
Fixed Effects SE Method | Model-Based |
Degrees of Freedom Method | Residual |
Class Level Information | ||
---|---|---|
Class | Levels | Values |
DIET | 6 | N/N85 N/R40 N/R50 NP R/R50 lopro |
Dimensions | |
---|---|
Covariance Parameters | 1 |
Columns in X | 7 |
Columns in Z | 0 |
Subjects | 1 |
Max Obs Per Subject | 349 |
Number of Observations | |
---|---|
Number of Observations Read | 349 |
Number of Observations Used | 349 |
Number of Observations Not Used | 0 |
Covariance Parameter Estimates | |||||||
---|---|---|---|---|---|---|---|
Cov Parm | Estimate | Standard Error | Z Value | Pr Z | Alpha | Lower | Upper |
Residual | 44.5989 | 3.4056 | 13.10 | <.0001 | 0.05 | 38.6102 | 52.1055 |
Fit Statistics | |
---|---|
-2 Res Log Likelihood | 2300.3 |
AIC (smaller is better) | 2302.3 |
AICC (smaller is better) | 2302.4 |
BIC (smaller is better) | 2306.2 |
Type 3 Tests of Fixed Effects | ||||
---|---|---|---|---|
Effect | Num DF | Den DF | F Value | Pr > F |
DIET | 5 | 343 | 57.10 | <.0001 |
Chapter 5 : Mouse feeding example |
Testing homogeniety of Variance with PROC MIXED |
Model Information | |
---|---|
Data Set | WORK.MOUSEDIET |
Dependent Variable | LIFETIME |
Covariance Structure | Variance Components |
Group Effect | DIET |
Estimation Method | REML |
Residual Variance Method | None |
Fixed Effects SE Method | Prasad-Rao-Jeske-Kackar-Harville |
Degrees of Freedom Method | Kenward-Roger |
Class Level Information | ||
---|---|---|
Class | Levels | Values |
DIET | 6 | N/N85 N/R40 N/R50 NP R/R50 lopro |
Dimensions | |
---|---|
Covariance Parameters | 6 |
Columns in X | 7 |
Columns in Z | 0 |
Subjects | 349 |
Max Obs Per Subject | 1 |
Number of Observations | |
---|---|
Number of Observations Read | 349 |
Number of Observations Used | 349 |
Number of Observations Not Used | 0 |
Iteration History | |||
---|---|---|---|
Iteration | Evaluations | -2 Res Log Like | Criterion |
0 | 1 | 2300.34852165 | |
1 | 1 | 2289.27652629 | 0.00000000 |
Convergence criteria met. |
Covariance Parameter Estimates | ||||||||
---|---|---|---|---|---|---|---|---|
Cov Parm | Group | Estimate | Standard Error | Z Value | Pr Z | Alpha | Lower | Upper |
Residual | DIET N/N85 | 26.2687 | 4.9643 | 5.29 | <.0001 | 0.05 | 18.7234 | 39.5319 |
Residual | DIET N/R40 | 44.9356 | 8.2733 | 5.43 | <.0001 | 0.05 | 32.2855 | 66.8452 |
Residual | DIET N/R50 | 60.3448 | 10.2001 | 5.92 | <.0001 | 0.05 | 44.4538 | 86.6356 |
Residual | DIET NP | 37.6223 | 7.6796 | 4.90 | <.0001 | 0.05 | 26.1635 | 58.7189 |
Residual | DIET R/R50 | 44.6645 | 8.5172 | 5.24 | <.0001 | 0.05 | 31.7464 | 67.4911 |
Residual | DIET lopro | 48.8838 | 9.3218 | 5.24 | <.0001 | 0.05 | 34.7453 | 73.8667 |
Fit Statistics | |
---|---|
-2 Res Log Likelihood | 2289.3 |
AIC (smaller is better) | 2301.3 |
AICC (smaller is better) | 2301.5 |
BIC (smaller is better) | 2324.4 |
Null Model Likelihood Ratio Test | ||
---|---|---|
DF | Chi-Square | Pr > ChiSq |
5 | 11.07 | 0.0500 |
Type 3 Tests of Fixed Effects | ||||
---|---|---|---|---|
Effect | Num DF | Den DF | F Value | Pr > F |
DIET | 5 | 155 | 64.71 | <.0001 |
Chapter 5 : Mouse feeding example |
Testing homogeniety of Variance with PROC MIXED |
Model Information | |
---|---|
Data Set | WORK.MOUSEDIET |
Dependent Variable | LIFETIME |
Covariance Structure | Variance Components |
Group Effect | DIET |
Estimation Method | REML |
Residual Variance Method | None |
Fixed Effects SE Method | Prasad-Rao-Jeske-Kackar-Harville |
Degrees of Freedom Method | Kenward-Roger |
Class Level Information | ||
---|---|---|
Class | Levels | Values |
DIET | 6 | N/N85 N/R40 N/R50 NP R/R50 lopro |
Dimensions | |
---|---|
Covariance Parameters | 6 |
Columns in X | 7 |
Columns in Z | 0 |
Subjects | 349 |
Max Obs Per Subject | 1 |
Number of Observations | |
---|---|
Number of Observations Read | 349 |
Number of Observations Used | 349 |
Number of Observations Not Used | 0 |
Iteration History | |||
---|---|---|---|
Iteration | Evaluations | -2 Res Log Like | Criterion |
0 | 1 | 2300.34852165 | |
1 | 1 | 2289.27652629 | 0.00000000 |
Convergence criteria met. |
Covariance Parameter Estimates | ||||||||
---|---|---|---|---|---|---|---|---|
Cov Parm | Group | Estimate | Standard Error | Z Value | Pr Z | Alpha | Lower | Upper |
Residual | DIET N/N85 | 26.2687 | 4.9643 | 5.29 | <.0001 | 0.05 | 18.7234 | 39.5319 |
Residual | DIET N/R40 | 44.9356 | 8.2733 | 5.43 | <.0001 | 0.05 | 32.2855 | 66.8452 |
Residual | DIET N/R50 | 60.3448 | 10.2001 | 5.92 | <.0001 | 0.05 | 44.4538 | 86.6356 |
Residual | DIET NP | 37.6223 | 7.6796 | 4.90 | <.0001 | 0.05 | 26.1635 | 58.7189 |
Residual | DIET R/R50 | 44.6645 | 8.5172 | 5.24 | <.0001 | 0.05 | 31.7464 | 67.4911 |
Residual | DIET lopro | 48.8838 | 9.3218 | 5.24 | <.0001 | 0.05 | 34.7453 | 73.8667 |
Fit Statistics | |
---|---|
-2 Res Log Likelihood | 2289.3 |
AIC (smaller is better) | 2301.3 |
AICC (smaller is better) | 2301.5 |
BIC (smaller is better) | 2324.4 |
Null Model Likelihood Ratio Test | ||
---|---|---|
DF | Chi-Square | Pr > ChiSq |
5 | 11.07 | 0.0500 |
Type 3 Tests of Fixed Effects | ||||
---|---|---|---|---|
Effect | Num DF | Den DF | F Value | Pr > F |
DIET | 5 | 155 | 64.71 | <.0001 |
Contrasts | ||||
---|---|---|---|---|
Label | Num DF | Den DF | F Value | Pr > F |
A: N/N85 N/R50 | 1 | 122 | 70.40 | <.0001 |
B: N/R50 R/R50 | 1 | 124 | 0.21 | 0.6474 |
C: N/R40 N/R50 | 1 | 129 | 4.97 | 0.0275 |
D: N/R50 lopro | 1 | 123 | 3.96 | 0.0489 |
E: N/N85 NP | 1 | 93.9 | 22.77 | <.0001 |
Chapter 5 : Mouse feeding example |
Analysis wtih PROC GLM |
Class Level Information | ||
---|---|---|
Class | Levels | Values |
DIET | 6 | N/N85 N/R40 N/R50 NP R/R50 lopro |
Number of Observations Read | 349 |
---|---|
Number of Observations Used | 349 |
Chapter 5 : Mouse feeding example |
Analysis wtih PROC GLM |
Source | DF | Sum of Squares | Mean Square | F Value | Pr > F |
---|---|---|---|---|---|
Model | 5 | 12733.94181 | 2546.78836 | 57.10 | <.0001 |
Error | 343 | 15297.41532 | 44.59888 | ||
Corrected Total | 348 | 28031.35713 |
R-Square | Coeff Var | Root MSE | LIFETIME Mean |
---|---|---|---|
0.454275 | 17.21323 | 6.678239 | 38.79713 |
Source | DF | Type I SS | Mean Square | F Value | Pr > F |
---|---|---|---|---|---|
DIET | 5 | 12733.94181 | 2546.78836 | 57.10 | <.0001 |
Source | DF | Type III SS | Mean Square | F Value | Pr > F |
---|---|---|---|---|---|
DIET | 5 | 12733.94181 | 2546.78836 | 57.10 | <.0001 |
Contrast | DF | Contrast SS | Mean Square | F Value | Pr > F |
---|---|---|---|---|---|
A: N/N85 N/R50 | 1 | 2917.456121 | 2917.456121 | 65.42 | <.0001 |
B: N/R50 R/R50 | 1 | 10.843803 | 10.843803 | 0.24 | 0.6223 |
C: N/R40 N/R50 | 1 | 258.510054 | 258.510054 | 5.80 | 0.0166 |
D: N/R50 lopro | 1 | 213.507110 | 213.507110 | 4.79 | 0.0293 |
E: N/N85 NP | 1 | 737.128081 | 737.128081 | 16.53 | <.0001 |
Chapter 5 : Mouse feeding example |
Analysis wtih PROC GLM |
Model Information | |
---|---|
Data Set | WORK.MOUSEDIET |
Dependent Variable | LIFETIME |
Covariance Structure | Variance Components |
Group Effect | DIET |
Estimation Method | REML |
Residual Variance Method | None |
Fixed Effects SE Method | Prasad-Rao-Jeske-Kackar-Harville |
Degrees of Freedom Method | Kenward-Roger |
Class Level Information | ||
---|---|---|
Class | Levels | Values |
DIET | 6 | N/N85 N/R40 N/R50 NP R/R50 lopro |
Dimensions | |
---|---|
Covariance Parameters | 6 |
Columns in X | 7 |
Columns in Z | 0 |
Subjects | 349 |
Max Obs Per Subject | 1 |
Number of Observations | |
---|---|
Number of Observations Read | 349 |
Number of Observations Used | 349 |
Number of Observations Not Used | 0 |
Iteration History | |||
---|---|---|---|
Iteration | Evaluations | -2 Res Log Like | Criterion |
0 | 1 | 2300.34852165 | |
1 | 1 | 2289.27652629 | 0.00000000 |
Convergence criteria met. |
Covariance Parameter Estimates | ||||||||
---|---|---|---|---|---|---|---|---|
Cov Parm | Group | Estimate | Standard Error | Z Value | Pr Z | Alpha | Lower | Upper |
Residual | DIET N/N85 | 26.2687 | 4.9643 | 5.29 | <.0001 | 0.05 | 18.7234 | 39.5319 |
Residual | DIET N/R40 | 44.9356 | 8.2733 | 5.43 | <.0001 | 0.05 | 32.2855 | 66.8452 |
Residual | DIET N/R50 | 60.3448 | 10.2001 | 5.92 | <.0001 | 0.05 | 44.4538 | 86.6356 |
Residual | DIET NP | 37.6223 | 7.6796 | 4.90 | <.0001 | 0.05 | 26.1635 | 58.7189 |
Residual | DIET R/R50 | 44.6645 | 8.5172 | 5.24 | <.0001 | 0.05 | 31.7464 | 67.4911 |
Residual | DIET lopro | 48.8838 | 9.3218 | 5.24 | <.0001 | 0.05 | 34.7453 | 73.8667 |
Fit Statistics | |
---|---|
-2 Res Log Likelihood | 2289.3 |
AIC (smaller is better) | 2301.3 |
AICC (smaller is better) | 2301.5 |
BIC (smaller is better) | 2324.4 |
Null Model Likelihood Ratio Test | ||
---|---|---|
DF | Chi-Square | Pr > ChiSq |
5 | 11.07 | 0.0500 |
Type 3 Tests of Fixed Effects | ||||
---|---|---|---|---|
Effect | Num DF | Den DF | F Value | Pr > F |
DIET | 5 | 155 | 64.71 | <.0001 |
Contrasts | ||||
---|---|---|---|---|
Label | Num DF | Den DF | F Value | Pr > F |
1 | 1 | 147 | 141.01 | <.0001 |
2 | 1 | 147 | 141.01 | <.0001 |
3 | 1 | 74.4 | 2.54 | 0.1154 |
4 | 1 | 74.4 | 2.54 | 0.1154 |
5 | 1 | 221 | 303.02 | <.0001 |
6 | 1 | 221 | 303.02 | <.0001 |
7 | . | . | . | . |
8 | 1 | 236 | 312.50 | <.0001 |
Chapter 5 : Mouse feeding example |
Analysis of Variance with PROC MIXED |
Model Information | |
---|---|
Data Set | WORK.MOUSEDIET |
Dependent Variable | LIFETIME |
Covariance Structure | Variance Components |
Group Effect | DIET |
Estimation Method | REML |
Residual Variance Method | None |
Fixed Effects SE Method | Prasad-Rao-Jeske-Kackar-Harville |
Degrees of Freedom Method | Kenward-Roger |
Class Level Information | ||
---|---|---|
Class | Levels | Values |
DIET | 6 | N/N85 N/R40 N/R50 NP R/R50 lopro |
Dimensions | |
---|---|
Covariance Parameters | 6 |
Columns in X | 7 |
Columns in Z | 0 |
Subjects | 349 |
Max Obs Per Subject | 1 |
Number of Observations | |
---|---|
Number of Observations Read | 349 |
Number of Observations Used | 349 |
Number of Observations Not Used | 0 |
Iteration History | |||
---|---|---|---|
Iteration | Evaluations | -2 Res Log Like | Criterion |
0 | 1 | 2300.34852165 | |
1 | 1 | 2289.27652629 | 0.00000000 |
Convergence criteria met. |
Covariance Parameter Estimates | ||||||||
---|---|---|---|---|---|---|---|---|
Cov Parm | Group | Estimate | Standard Error | Z Value | Pr Z | Alpha | Lower | Upper |
Residual | DIET N/N85 | 26.2687 | 4.9643 | 5.29 | <.0001 | 0.05 | 18.7234 | 39.5319 |
Residual | DIET N/R40 | 44.9356 | 8.2733 | 5.43 | <.0001 | 0.05 | 32.2855 | 66.8452 |
Residual | DIET N/R50 | 60.3448 | 10.2001 | 5.92 | <.0001 | 0.05 | 44.4538 | 86.6356 |
Residual | DIET NP | 37.6223 | 7.6796 | 4.90 | <.0001 | 0.05 | 26.1635 | 58.7189 |
Residual | DIET R/R50 | 44.6645 | 8.5172 | 5.24 | <.0001 | 0.05 | 31.7464 | 67.4911 |
Residual | DIET lopro | 48.8838 | 9.3218 | 5.24 | <.0001 | 0.05 | 34.7453 | 73.8667 |
Fit Statistics | |
---|---|
-2 Res Log Likelihood | 2289.3 |
AIC (smaller is better) | 2301.3 |
AICC (smaller is better) | 2301.5 |
BIC (smaller is better) | 2324.4 |
Null Model Likelihood Ratio Test | ||
---|---|---|
DF | Chi-Square | Pr > ChiSq |
5 | 11.07 | 0.0500 |
Type 3 Tests of Fixed Effects | ||||
---|---|---|---|---|
Effect | Num DF | Den DF | F Value | Pr > F |
DIET | 5 | 155 | 64.71 | <.0001 |
Least Squares Means | |||||||||
---|---|---|---|---|---|---|---|---|---|
Effect | DIET | Estimate | Standard Error | DF | t Value | Pr > |t| | Alpha | Lower | Upper |
DIET | N/N85 | 32.6912 | 0.6789 | 56 | 48.16 | <.0001 | 0.05 | 31.3313 | 34.0512 |
DIET | N/R40 | 45.1167 | 0.8654 | 59 | 52.13 | <.0001 | 0.05 | 43.3850 | 46.8483 |
DIET | N/R50 | 42.2972 | 0.9219 | 70 | 45.88 | <.0001 | 0.05 | 40.4585 | 44.1359 |
DIET | NP | 27.4020 | 0.8762 | 48 | 31.27 | <.0001 | 0.05 | 25.6402 | 29.1638 |
DIET | R/R50 | 42.8857 | 0.8931 | 55 | 48.02 | <.0001 | 0.05 | 41.0960 | 44.6755 |
DIET | lopro | 39.6857 | 0.9343 | 55 | 42.48 | <.0001 | 0.05 | 37.8133 | 41.5581 |
Differences of Least Squares Means | ||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Effect | DIET | _DIET | Estimate | Standard Error | DF | t Value | Pr > |t| | Adjustment | Adj P | Alpha | Lower | Upper | Adj Lower | Adj Upper |
DIET | N/N85 | N/R40 | -12.4254 | 1.0999 | 110 | -11.30 | <.0001 | Tukey-Kramer | <.0001 | 0.05 | -14.6052 | -10.2457 | -15.5995 | -9.2514 |
DIET | N/N85 | N/R50 | -9.6060 | 1.1449 | 122 | -8.39 | <.0001 | Tukey-Kramer | <.0001 | 0.05 | -11.8724 | -7.3395 | -12.9099 | -6.3020 |
DIET | N/N85 | NP | 5.2892 | 1.1084 | 93.9 | 4.77 | <.0001 | Tukey-Kramer | <.0001 | 0.05 | 3.0883 | 7.4901 | 2.0905 | 8.4879 |
DIET | N/N85 | R/R50 | -10.1945 | 1.1218 | 103 | -9.09 | <.0001 | Tukey-Kramer | <.0001 | 0.05 | -12.4193 | -7.9697 | -13.4318 | -6.9572 |
DIET | N/N85 | lopro | -6.9945 | 1.1549 | 101 | -6.06 | <.0001 | Tukey-Kramer | <.0001 | 0.05 | -9.2855 | -4.7034 | -10.3273 | -3.6617 |
DIET | N/R40 | N/R50 | 2.8195 | 1.2645 | 129 | 2.23 | 0.0275 | Tukey-Kramer | 0.2302 | 0.05 | 0.3177 | 5.3213 | -0.8295 | 6.4684 |
DIET | N/R40 | NP | 17.7146 | 1.2316 | 106 | 14.38 | <.0001 | Tukey-Kramer | <.0001 | 0.05 | 15.2728 | 20.1564 | 14.1606 | 21.2686 |
DIET | N/R40 | R/R50 | 2.2310 | 1.2436 | 113 | 1.79 | 0.0755 | Tukey-Kramer | 0.4727 | 0.05 | -0.2327 | 4.6946 | -1.3578 | 5.8197 |
DIET | N/R40 | lopro | 5.4310 | 1.2735 | 113 | 4.26 | <.0001 | Tukey-Kramer | 0.0005 | 0.05 | 2.9078 | 7.9541 | 1.7559 | 9.1061 |
DIET | N/R50 | NP | 14.8951 | 1.2719 | 116 | 11.71 | <.0001 | Tukey-Kramer | <.0001 | 0.05 | 12.3759 | 17.4143 | 11.2247 | 18.5656 |
DIET | N/R50 | R/R50 | -0.5885 | 1.2836 | 124 | -0.46 | 0.6474 | Tukey-Kramer | 0.9974 | 0.05 | -3.1290 | 1.9520 | -4.2926 | 3.1155 |
DIET | N/R50 | lopro | 2.6115 | 1.3126 | 123 | 1.99 | 0.0489 | Tukey-Kramer | 0.3530 | 0.05 | 0.01326 | 5.2097 | -1.1763 | 6.3993 |
DIET | NP | R/R50 | -15.4837 | 1.2512 | 103 | -12.38 | <.0001 | Tukey-Kramer | <.0001 | 0.05 | -17.9651 | -13.0022 | -19.0942 | -11.8731 |
DIET | NP | lopro | -12.2837 | 1.2809 | 103 | -9.59 | <.0001 | Tukey-Kramer | <.0001 | 0.05 | -14.8241 | -9.7433 | -15.9801 | -8.5873 |
DIET | R/R50 | lopro | 3.2000 | 1.2925 | 110 | 2.48 | 0.0148 | Tukey-Kramer | 0.1378 | 0.05 | 0.6385 | 5.7615 | -0.5298 | 6.9298 |