Mixed Effects Negative Binomial Regression Model (Random Intercept and Slope)
Analysis Excluding Patient 49 (Potential Outlier)
Clinical Trial of Epileptic Patients

The GLIMMIX Procedure

Model Information
Data Set WORK.SEIZURE
Response Variable y
Response Distribution Negative Binomial
Link Function Log
Variance Function Default
Offset Variable ltime
Variance Matrix Blocked By id
Estimation Technique Maximum Likelihood
Likelihood Approximation Gauss-Hermite Quadrature
Degrees of Freedom Method Containment

Class Level Information
Class Levels Values
id 58 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 50 51 52 53 54 55 56 57 58 59

Number of Observations Read 290
Number of Observations Used 290

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

Optimization Information
Optimization Technique Dual Quasi-Newton
Parameters in Optimization 7
Lower Boundaries 3
Upper Boundaries 0
Fixed Effects Not Profiled
Starting From GLM estimates
Quadrature Points 50

Iteration History
Iteration Restarts Evaluations Objective
Function
Change Max
Gradient
0 0 4 1758.3010557 . 186.6201
1 0 5 1750.6327206 7.66833509 184.1917
2 0 11 1720.2691576 30.36356298 58.98885
3 0 5 1716.272248 3.99690959 200.5654
4 0 2 1712.1482567 4.12399137 95.15279
5 0 3 1710.7220829 1.42617372 26.1452
6 0 2 1709.5657251 1.15635782 8.262079
7 0 3 1709.3416724 0.22405275 8.314936
8 0 3 1709.276991 0.06468132 6.572191
9 0 3 1709.2585659 0.01842515 3.145712
10 0 3 1709.2476044 0.01096154 2.358768
11 0 3 1709.2448464 0.00275800 0.671064
12 0 3 1709.2441879 0.00065844 0.026985
13 0 3 1709.2441842 0.00000373 0.001043

Convergence criterion (GCONV=1E-8) satisfied.

Fit Statistics
-2 Log Likelihood 1709.24
AIC (smaller is better) 1723.24
AICC (smaller is better) 1723.64
BIC (smaller is better) 1737.67
CAIC (smaller is better) 1744.67
HQIC (smaller is better) 1728.86

Fit Statistics for Conditional Distribution
-2 log L(y | r. effects) 1485.46
Pearson Chi-Square 207.57
Pearson Chi-Square / DF 0.72

Estimated G Matrix
Effect Row Col1 Col2
Intercept 1 0.4414  
time 2   0.1278

Covariance Parameter Estimates
Cov Parm Subject Estimate Standard Error
Intercept id 0.4414 0.1040
time id 0.1278 0.06783
Scale   0.1173 0.02543

Solutions for Fixed Effects
Effect Estimate Standard Error DF t Value Pr > |t|
Intercept 1.0986 0.1481 56 7.42 <.0001
time 0.005442 0.1160 56 0.05 0.9628
trt 0.003289 0.2055 174 0.02 0.9873
time*trt -0.3585 0.1629 174 -2.20 0.0291

Type III Tests of Fixed Effects
Effect Num DF Den DF F Value Pr > F
time 1 56 0.00 0.9628
trt 1 174 0.00 0.9873
time*trt 1 174 4.84 0.0291