|
Mixed Effects Log-linear Regression Model (Random Intercept and Slope) |
|
Clinical Trial of Epileptic Patients |
|
The NLMIXED Procedure |
|
Specifications |
|
|
Data Set |
WORK.SEIZURE |
|
Dependent Variable |
Y |
|
Distribution for Dependent Variable |
Poisson |
|
Random Effects |
u1 u2 |
|
Distribution for Random Effects |
Normal |
|
Subject Variable |
Id |
|
Optimization Technique |
Dual Quasi-Newton |
|
Integration Method |
Adaptive Gaussian Quadrature |
|
Dimensions |
|
|
Observations Used |
295 |
|
Observations Not Used |
0 |
|
Total Observations |
295 |
|
Subjects |
59 |
|
Max Obs Per Subject |
5 |
|
Parameters |
7 |
|
Quadrature Points |
50 |
|
Parameters |
|||||||
|
beta1 |
beta2 |
beta3 |
beta4 |
s2u1 |
s2u2 |
cu12 |
NegLogLike |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0 |
0 |
0 |
2389.89383 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.2 |
0 |
0 |
1035.75052 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.4 |
0 |
0 |
1015.0326 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.6 |
0 |
0 |
1012.08586 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.8 |
0 |
0 |
1012.84609 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
1 |
0 |
0 |
1014.69701 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0 |
0.2 |
0 |
1457.05866 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.2 |
0.2 |
0 |
946.672423 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.4 |
0.2 |
0 |
932.319714 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.6 |
0.2 |
0 |
931.236945 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.8 |
0.2 |
0 |
932.857916 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
1 |
0.2 |
0 |
935.198 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0 |
0.4 |
0 |
1427.41194 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.2 |
0.4 |
0 |
946.169637 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.4 |
0.4 |
0 |
932.524815 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.6 |
0.4 |
0 |
931.635202 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.8 |
0.4 |
0 |
933.336091 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
1 |
0.4 |
0 |
935.716799 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0 |
0.6 |
0 |
1419.56258 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.2 |
0.6 |
0 |
949.034202 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.4 |
0.6 |
0 |
935.702247 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.6 |
0.6 |
0 |
934.898051 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.8 |
0.6 |
0 |
936.633504 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
1 |
0.6 |
0 |
939.031209 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0 |
0.8 |
0 |
1417.22116 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.2 |
0.8 |
0 |
952.314329 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.4 |
0.8 |
0 |
939.162429 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.6 |
0.8 |
0 |
938.407697 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.8 |
0.8 |
0 |
940.16304 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
1 |
0.8 |
0 |
942.5704 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0 |
1 |
0 |
1416.91642 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.2 |
1 |
0 |
955.486634 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.4 |
1 |
0 |
942.452802 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.6 |
1 |
0 |
941.73073 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.8 |
1 |
0 |
943.499201 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
1 |
1 |
0 |
945.912902 |
|
Iteration History |
||||||
|
Iter |
|
Calls |
NegLogLike |
Diff |
MaxGrad |
Slope |
|
1 |
|
3 |
926.152102 |
5.084843 |
12.00163 |
-416.31 |
|
2 |
|
5 |
925.933812 |
0.21829 |
13.61071 |
-11.7582 |
|
3 |
|
6 |
925.819654 |
0.114158 |
11.49685 |
-5.02778 |
|
4 |
|
7 |
925.620274 |
0.199381 |
7.70171 |
-3.2169 |
|
5 |
|
8 |
925.458053 |
0.162221 |
8.582898 |
-7.74709 |
|
6 |
|
9 |
925.242355 |
0.215698 |
5.416426 |
-4.34749 |
|
7 |
|
10 |
925.022545 |
0.21981 |
3.43752 |
-1.42986 |
|
8 |
|
12 |
924.972964 |
0.049582 |
1.224833 |
-0.10048 |
|
9 |
|
14 |
924.968495 |
0.004468 |
0.181362 |
-0.01157 |
|
10 |
|
16 |
924.968047 |
0.000448 |
0.137067 |
-0.00168 |
|
11 |
|
18 |
924.967965 |
0.000082 |
0.009213 |
-0.0002 |
|
12 |
|
20 |
924.967965 |
6.796E-7 |
0.000377 |
-1.41E-6 |
|
NOTE: GCONV convergence criterion satisfied. |
|
Fit Statistics |
|
|
-2 Log Likelihood |
1849.9 |
|
AIC (smaller is better) |
1863.9 |
|
AICC (smaller is better) |
1864.3 |
|
BIC (smaller is better) |
1878.5 |
|
Parameter Estimates |
|||||||||
|
Parameter |
Estimate |
Standard Error |
DF |
t Value |
Pr > |t| |
Alpha |
Lower |
Upper |
Gradient |
|
beta1 |
1.0707 |
0.1406 |
57 |
7.62 |
<.0001 |
0.05 |
0.7892 |
1.3521 |
-0.00016 |
|
beta2 |
-0.00039 |
0.1097 |
57 |
-0.00 |
0.9971 |
0.05 |
-0.2201 |
0.2193 |
-0.00038 |
|
beta3 |
0.05131 |
0.1931 |
57 |
0.27 |
0.7914 |
0.05 |
-0.3353 |
0.4379 |
0.000088 |
|
beta4 |
-0.3065 |
0.1513 |
57 |
-2.03 |
0.0474 |
0.05 |
-0.6094 |
-0.00361 |
-0.00022 |
|
s2u1 |
0.5010 |
0.1010 |
57 |
4.96 |
<.0001 |
0.05 |
0.2988 |
0.7032 |
0.00007 |
|
s2u2 |
0.2334 |
0.06075 |
57 |
3.84 |
0.0003 |
0.05 |
0.1117 |
0.3550 |
-0.00003 |
|
cu12 |
0.05408 |
0.05589 |
57 |
0.97 |
0.3373 |
0.05 |
-0.05784 |
0.1660 |
0.000327 |