Random Effects Logistic Regression Model for Odds of Amenorrhea |
Clinical Trial of Contracepting Women |
The NLMIXED Procedure |
Specifications |
|
Data Set |
WORK.AMEN |
Dependent Variable |
y |
Distribution for Dependent Variable |
Binary |
Random Effects |
u |
Distribution for Random Effects |
Normal |
Subject Variable |
id |
Optimization Technique |
Dual Quasi-Newton |
Integration Method |
Adaptive Gaussian Quadrature |
Dimensions |
|
Observations Used |
3616 |
Observations Not Used |
988 |
Total Observations |
4604 |
Subjects |
1151 |
Max Obs Per Subject |
4 |
Parameters |
6 |
Quadrature Points |
50 |
Parameters |
||||||
b0 |
b1 |
b2 |
b3 |
b4 |
s1 |
NegLogLike |
-2.2461 |
0.703 |
-0.0323 |
0.338 |
-0.0683 |
0 |
2186.2544 |
-2.2461 |
0.703 |
-0.0323 |
0.338 |
-0.0683 |
1 |
2022.12015 |
-2.2461 |
0.703 |
-0.0323 |
0.338 |
-0.0683 |
2 |
1985.21862 |
-2.2461 |
0.703 |
-0.0323 |
0.338 |
-0.0683 |
3 |
1977.43936 |
-2.2461 |
0.703 |
-0.0323 |
0.338 |
-0.0683 |
4 |
1980.13195 |
-2.2461 |
0.703 |
-0.0323 |
0.338 |
-0.0683 |
5 |
1987.30068 |
-2.2461 |
0.703 |
-0.0323 |
0.338 |
-0.0683 |
6 |
1996.50765 |
-2.2461 |
0.703 |
-0.0323 |
0.338 |
-0.0683 |
7 |
2006.62042 |
-2.2461 |
0.703 |
-0.0323 |
0.338 |
-0.0683 |
8 |
2017.0682 |
Iteration History |
||||||
Iter |
|
Calls |
NegLogLike |
Diff |
MaxGrad |
Slope |
1 |
|
4 |
1975.19945 |
2.239914 |
88.74847 |
-2278.13 |
2 |
|
5 |
1959.91912 |
15.28033 |
173.7908 |
-295.465 |
3 |
|
7 |
1950.15246 |
9.766663 |
23.56372 |
-946.798 |
4 |
|
9 |
1941.54662 |
8.605838 |
55.6478 |
-5.32732 |
5 |
|
11 |
1936.37483 |
5.171791 |
64.52357 |
-2.69208 |
6 |
|
13 |
1934.5452 |
1.829628 |
6.353883 |
-1.95532 |
7 |
|
15 |
1934.50882 |
0.036382 |
1.980495 |
-0.04951 |
8 |
|
16 |
1934.47401 |
0.034811 |
4.299917 |
-0.01296 |
9 |
|
18 |
1934.4655 |
0.008513 |
1.29607 |
-0.01213 |
10 |
|
20 |
1934.46498 |
0.000511 |
0.191723 |
-0.00081 |
11 |
|
22 |
1934.46498 |
1.56E-6 |
0.006509 |
-2.52E-6 |
NOTE: GCONV convergence criterion satisfied. |
Fit Statistics |
|
-2 Log Likelihood |
3868.9 |
AIC (smaller is better) |
3880.9 |
AICC (smaller is better) |
3881.0 |
BIC (smaller is better) |
3911.2 |
Parameter Estimates |
|||||||||
Parameter |
Estimate |
Standard Error |
DF |
t Value |
Pr > |t| |
Alpha |
Lower |
Upper |
Gradient |
b0 |
-3.8057 |
0.3050 |
1150 |
-12.48 |
<.0001 |
0.05 |
-4.4040 |
-3.2073 |
0.001134 |
b1 |
1.1332 |
0.2682 |
1150 |
4.22 |
<.0001 |
0.05 |
0.6070 |
1.6595 |
0.002455 |
b2 |
-0.04193 |
0.05481 |
1150 |
-0.76 |
0.4445 |
0.05 |
-0.1495 |
0.06561 |
0.006509 |
b3 |
0.5644 |
0.1922 |
1150 |
2.94 |
0.0034 |
0.05 |
0.1872 |
0.9416 |
0.000525 |
b4 |
-0.1095 |
0.04961 |
1150 |
-2.21 |
0.0274 |
0.05 |
-0.2069 |
-0.01221 |
0.00269 |
s1 |
5.0646 |
0.5840 |
1150 |
8.67 |
<.0001 |
0.05 |
3.9187 |
6.2104 |
0.000103 |
Contrasts |
||||
Label |
Num DF |
Den DF |
F Value |
Pr > F |
Trt X Time Interaction |
2 |
1150 |
6.20 |
0.0021 |