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