|
Mixed Effects Log-linear Regression Model (Random Intercept and Slope) |
|
Analysis Excluding Patient 49 (Potential Outlier) |
|
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 |
290 |
|
Observations Not Used |
0 |
|
Total Observations |
290 |
|
Subjects |
58 |
|
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 |
1771.18019 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.2 |
0 |
0 |
975.322584 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.4 |
0 |
0 |
959.914013 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.6 |
0 |
0 |
958.663874 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.8 |
0 |
0 |
960.232556 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
1 |
0 |
0 |
962.544129 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0 |
0.2 |
0 |
1279.83589 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.2 |
0.2 |
0 |
912.380805 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.4 |
0.2 |
0 |
901.394355 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.6 |
0.2 |
0 |
901.372907 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.8 |
0.2 |
0 |
903.487076 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
1 |
0.2 |
0 |
906.099422 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0 |
0.4 |
0 |
1257.0264 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.2 |
0.4 |
0 |
912.184035 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.4 |
0.4 |
0 |
901.857245 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.6 |
0.4 |
0 |
902.013235 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.8 |
0.4 |
0 |
904.19953 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
1 |
0.4 |
0 |
906.847848 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0 |
0.6 |
0 |
1251.39584 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.2 |
0.6 |
0 |
915.070522 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.4 |
0.6 |
0 |
905.039947 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.6 |
0.6 |
0 |
905.275927 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.8 |
0.6 |
0 |
907.494096 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
1 |
0.6 |
0 |
910.157805 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0 |
0.8 |
0 |
1250.12641 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.2 |
0.8 |
0 |
918.320524 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.4 |
0.8 |
0 |
908.461566 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.6 |
0.8 |
0 |
908.74426 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.8 |
0.8 |
0 |
910.980957 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
1 |
0.8 |
0 |
913.653507 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0 |
1 |
0 |
1250.44121 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.2 |
1 |
0 |
921.449939 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.4 |
1 |
0 |
911.703949 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.6 |
1 |
0 |
912.017642 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
0.8 |
1 |
0 |
914.266643 |
|
1.1213 |
0.1244 |
0.0718 |
-0.115 |
1 |
1 |
0 |
916.945041 |
|
Iteration History |
||||||
|
Iter |
|
Calls |
NegLogLike |
Diff |
MaxGrad |
Slope |
|
1 |
|
3 |
895.107816 |
6.265092 |
12.12678 |
-423.342 |
|
2 |
|
5 |
894.353775 |
0.754041 |
15.54505 |
-19.4536 |
|
3 |
|
7 |
894.146094 |
0.207681 |
12.2134 |
-9.97304 |
|
4 |
|
9 |
894.02912 |
0.116974 |
11.30246 |
-3.06491 |
|
5 |
|
10 |
893.879361 |
0.149759 |
9.461355 |
-2.89954 |
|
6 |
|
11 |
893.755246 |
0.124115 |
4.877882 |
-2.33965 |
|
7 |
|
12 |
893.564602 |
0.190644 |
0.762254 |
-0.79132 |
|
8 |
|
14 |
893.550247 |
0.014355 |
0.334841 |
-0.03658 |
|
9 |
|
16 |
893.548565 |
0.001683 |
0.154194 |
-0.00603 |
|
10 |
|
18 |
893.548382 |
0.000183 |
0.057579 |
-0.00049 |
|
11 |
|
20 |
893.548374 |
7.597E-6 |
0.011938 |
-0.00002 |
|
12 |
|
22 |
893.548373 |
7.513E-7 |
0.000343 |
-1.55E-6 |
|
NOTE: GCONV convergence criterion satisfied. |
|
Fit Statistics |
|
|
-2 Log Likelihood |
1787.1 |
|
AIC (smaller is better) |
1801.1 |
|
AICC (smaller is better) |
1801.5 |
|
BIC (smaller is better) |
1815.5 |
|
Parameter Estimates |
|||||||||
|
Parameter |
Estimate |
Standard Error |
DF |
t Value |
Pr > |t| |
Alpha |
Lower |
Upper |
Gradient |
|
beta1 |
1.0692 |
0.1344 |
56 |
7.96 |
<.0001 |
0.05 |
0.8001 |
1.3384 |
0.000237 |
|
beta2 |
0.007805 |
0.1070 |
56 |
0.07 |
0.9421 |
0.05 |
-0.2066 |
0.2222 |
-0.00007 |
|
beta3 |
-0.00794 |
0.1860 |
56 |
-0.04 |
0.9661 |
0.05 |
-0.3806 |
0.3647 |
0.000263 |
|
beta4 |
-0.3461 |
0.1489 |
56 |
-2.33 |
0.0237 |
0.05 |
-0.6443 |
-0.04791 |
0.000217 |
|
s2u1 |
0.4529 |
0.09354 |
56 |
4.84 |
<.0001 |
0.05 |
0.2655 |
0.6403 |
-0.00034 |
|
s2u2 |
0.2163 |
0.05870 |
56 |
3.68 |
0.0005 |
0.05 |
0.09866 |
0.3339 |
-0.00022 |
|
cu12 |
0.01511 |
0.05290 |
56 |
0.29 |
0.7763 |
0.05 |
-0.09087 |
0.1211 |
-0.00018 |