66 Handbook of Discrete-Valued Time Series
TABLE 3.1
Results for testing various xed effect multiple GLARMA models for the Road Deaths Series in 17
U.S. States
Model −2log L S G
2
d.f. p-val
FE-I: Unrestricted 5345.69 92 — — —
FE-II: φ
s in 6 groups
FE-III: BAC, ALR, FS same
FE-IV: BAC, ALR, FS, lnOVD same
5347.87
5391.76
5436.57
81
43
27
G
2
II v I
= 2.18
G
2
III v II
= 43.89
G
2
IV v III
= 44.81
G
2
IV v II
= 88.61
11
38
16
54
0.998
0.236
0.00015
0.0021
We next check whether the regression coefcients in W
jt
vary signicantly between indi-
vidual states. We begin with the overall unrestricted t to all 17 states. We refer to this
as Model FE-I, which has −2 log L = 5345.69 with S = 92 parameters. Examination of
the individual estimates φ
ˆ
12
suggested that they could be simplied as follows: Group 1
(State 11, φ
ˆ
12
=−0.081 ± 0.045), Group 2 (States 1, 6, 7, 9, 10, 12:17, φ
ˆ
12
= 0.005 ± 0.013),
Group 3 (States 2, 4, φ
ˆ
12
= 0.066 ± 0.015), Group 4 (State 5, φ
ˆ
12
= 0.212 ± 0.087), Group 5
(State 8, φ
ˆ
12
= 0.401 ± 0.096), Group 6 (State 5, φ
ˆ
12
= 0.545 ± 0.219) in which, at most, 6
φ
12
coefcients are signicant.
The model with the φ
12
restricted to these groups is referred to as Model FE-II in
Table 3.1. Using the likelihood ratio test we obtain G
2
II v I
= 2.18 on 11 d.f.; hence, restriction
of the φ
12
would not be rejected. From this model, we then examined whether or not some
or all of the regression coefcients (other than the intercept which does vary substantially
between states) take common values across all 17 states. Model FE-III restricts the coef-
cients for BAC, ALR, Friday–Saturday to be the same and (see Table 3.1) G
2
=43.89
III v II
on 38 d.f. and associated p-value of 0.24, which is not sufciently strong evidence to sug-
gest that the impact of these variables differs between individual states in a statistically
signicant way. Next, in Model FE-IV, log OMVD was allowed to differ between states.
Compared with Model FE-III or Model FE-II, this risk control variable is strongly statisti-
cally signicant between states with G
2
IV v III
= 44.81 on 16 d.f. and associated p-value of
0.00015 and G
2
IV v II
= 88.61 on 54 d.f. and associated p-value of 0.0021.
Hence, Model FE-III provides a useful summary of the commonality or otherwise of
regression variable impacts on single vehicle night time road deaths across the 17 states. The
tted parameters and associated standard errors are reported in Table 3.2. The six groups
for φ
12
could be reduced to four by removing the nonsignicant cases of Groups 1 and 2.
We did not pursue this here, preferring to move onto the use of a random effects analy-
sis. The impact of lowering the legal BAC level is estimated to be β
ˆ
2
=−0.072 ± 0.022
conrming the statistical signicance of this association found in Bernat et al. (2004).
The xed effects GLARMA model analysis provides a good starting point for the random
effects GLARMA modeling that we turn to in the next section. In particular, it seems plau-
sible from the results of Table 3.1 that random effects will be needed for the intercept term
and the log OMVD term, but not for BAC, ALR, or Friday–Saturday effects. The parame-
ter values reported for Model FE-III in Table 3.2 can provide useful starting values for the
random effects model tting. For xed effects, we use the point estimates of coefcients for
predictors that are common to all series, while for predictors that vary between series, we
use the mean values of point estimates of the coefcients.