MQL estimation of Random Effects Logistic Regression Model for Odds of Amenorrhea
Clinical Trial of Contracepting Women

The GLIMMIX Procedure

Model Information
Data Set WORK.AMEN
Response Variable y
Response Distribution Binomial
Link Function Logit
Variance Function Default
Variance Matrix Blocked By id
Estimation Technique Residual MPL
Degrees of Freedom Method Containment
Fixed Effects SE Adjustment Sandwich - Classical

Number of Observations Read 4604
Number of Observations Used 3616

Dimensions
G-side Cov. Parameters 1
Columns in X 5
Columns in Z per Subject 1
Subjects (Blocks in V) 1151
Max Obs per Subject 4

Optimization Information
Optimization Technique Dual Quasi-Newton
Parameters in Optimization 1
Lower Boundaries 1
Upper Boundaries 0
Fixed Effects Profiled
Starting From Data

Iteration History
Iteration Restarts Subiterations Objective
Function
Change Max
Gradient
0 0 4 15924.450822 2.00000000 8.839E-6
1 0 4 15749.815 0.07957237 4.066E-6
2 0 2 15753.003353 0.00146408 4.397E-7
3 0 1 15753.012592 0.00000686 0.00037
4 0 1 15753.012419 0.00000312 0.000168
5 0 1 15753.01243 0.00000134 0.000072
6 0 1 15753.012425 0.00000058 0.000031
7 0 1 15753.012427 0.00000025 0.000013
8 0 1 15753.012426 0.00000107 0.00018
9 0 1 15753.012431 0.00000143 0.000077
10 0 1 15753.012425 0.00000062 0.000033
11 0 1 15753.012427 0.00000027 0.000014
12 0 1 15753.012426 0.00000114 0.000192
13 0 1 15753.012431 0.00000152 0.000082
14 0 1 15753.012425 0.00000066 0.000035
15 0 1 15753.012427 0.00000028 0.000015
16 0 1 15753.012426 0.00000122 0.000204
17 0 1 15753.012431 0.00000162 0.000088
18 0 1 15753.012425 0.00000070 0.000038
19 0 1 15753.012427 0.00000030 0.000016
20 0 1 15753.012426 0.00000130 0.000218
21 0 1 15753.012431 0.00000173 0.000093
22 0 1 15753.012425 0.00000075 0.00004
23 0 1 15753.012427 0.00000032 0.000017
24 0 1 15753.012426 0.00000014 7.468E-6
25 0 0 15753.012427 0.00000000 7.531E-6

Convergence criterion (PCONV=1.11022E-8) satisfied.

Fit Statistics
-2 Res Log Pseudo-Likelihood 15753.01
Generalized Chi-Square 2672.50
Gener. Chi-Square / DF 0.74

Estimated G Matrix
Effect Row Col1
Intercept 1 1.2608

Covariance Parameter Estimates
Cov Parm Subject Estimate Standard Error
Intercept id 1.2608 0.1184

Solutions for Fixed Effects
Effect Estimate Standard Error DF t Value Pr > |t|
Intercept -2.2163 0.1770 1150 -12.52 <.0001
time 0.6851 0.1598 2461 4.29 <.0001
time2 -0.03137 0.03223 2461 -0.97 0.3305
time*trt 0.3080 0.1115 2461 2.76 0.0058
time2*trt -0.06111 0.02899 2461 -2.11 0.0351

Type III Tests of Fixed Effects
Effect Num DF Den DF F Value Pr > F
time 1 2461 18.39 <.0001
time2 1 2461 0.95 0.3305
time*trt 1 2461 7.64 0.0058
time2*trt 1 2461 4.44 0.0351

Contrasts
Label Num DF Den DF Chi-Square F Value Pr > ChiSq Pr > F
2 df Test of Trt X Time Interaction 2 2461 10.51 5.26 0.0052 0.0053