Mixed Effects Negative Binomial Regression Model (Random Intercept and Slope) |
Analysis Excluding Patient 49 (Potential Outlier) |
Clinical Trial of Epileptic Patients |
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 |