139 State Space Models for Count Time Series
The origin of the method and results when applied to the Polio data set is listed in paren-
theses, and additionally, if the results are from application of the method by Nelson and
Leroux (2006) these are also indicated by an additional annotation ‘NL’. The methods can be
roughly partitioned into three groups. Group 1 consists of two implementations of MCEM
(Monte Carlo EM) and a Bayes procedure. Group 2, which is essentially approximate
likelihood based-methods, consists of PQL (penalized quasilikelihood), AL (approximate
likelihood), AL-BC (bias corrected AL), AIS (approximate importance sampling), AIS-BC
(bias corrected AIS), MCNR (Monte Carlo Newton–Raphson), and EIS (efcient impor-
tance sampling). Note that the rst 3 procedures of this group are nonsimulation based,
while the last 4 involve some level of simulation. Group 3 consists of nonlikelihood-based
procedures: GLM (generalized linear model estimates ignoring the latent process), GEE
(generalized estimating equations), CPL
2
(pairwise composite likelihood), and IBC (itera-
tive bias correction using iterative weighted least squares). We exclude from our review the
few studies that have used alternative response distributions or latent process distributions
for these data so that the methods are compared on the same model.
With the exception of the GLM, GEE, CPL
2
, and Bayesian analyses, all other methods
aim to obtain approximations to the likelihood estimates and their standard errors. Clearly
there are both substantial differences and similarities between the results for various meth-
ods, a point also noted in Nelson and Leroux (2006). We now discuss these differences and
similarities in more detail in an attempt to draw some general conclusions about which
methods may be preferred. Of course, this comparison is only for application to a sin-
gle data set and much more research is required before general conclusions can be drawn.
However, this is the only data set for which all the methods listed have been applied. Unfor-
tunately, simulation evidence comparing the variety of methods is rather limited with the
exception of the results in Nelson and Leroux (2006).
6.3.1 Estimate of Trend Coefficient
β
ˆ
2
The GLM, IBC method as implemented by Nelson and Leroux (2006), and CPL
2
give the
most negative trend estimates. It would appear as if the IBC method is not adjusting the bias
of the GLM estimate sufciently well and this may be a result of iterative weighted least
squares being used as the basis for the bias adjustment simulations. It is likely that these
methods are substantially biased. Amongst the remaining methods, there appear to be two
groups of values for the trend coefcient estimates: Group 1, the values for both implemen-
tations of MCEM and the Bayes t; and, Group 2 based on approximations to the likelihood
with and without importance sampling (PQL, AL, AL-BC, AIS, AIS-BC, MCNR, and EIS).
The concordance in Group 2 is perhaps not surprising since they are all aimed at approx-
imating the likelihood. However, it is surprising that the Group 1 do not agree as closely
with the Group 2 results. Turning to comparison of the estimated standard errors, those for
Group 1 appear to be substantially smaller than those for Group 2, and within this latter
group there is considerable agreement. Also note that the MCEM and Bayes methods are
biasing the point estimates towards larger negative values and biasing the associated esti-
mated standard errors downwards. The net effect of these two biases would be to increase
the ratio of estimate to standard error resulting in a higher chance of concluding that there
is a signicant downward trend in Polio cases over the time period of observation. On the
other hand, for Group 2, these test ratios would all be consistent with a conclusion of no
downward trend.