416 Handbook of Discrete-Valued Time Series
k
f
1
(k) =
y
1,
j
t
1
−1
k
y
2,
−
t−
j
1
1
α
11
j
1
(1 − α
11
)
y
1,t−1
−j
1
α
12
k−j
1
(1 − α
12
)
y
2,t−1
−k+j
1
,
j
1
=0
s
y
2,t−1
y
1,t−1
j
2
s−j
2
f
2
(s) =
j
2
s − j
2
α
22
(1 − α
22
)
y
2,t−1
−j
2
α
21
(1 − α
21
)
y
1,t−1
−s+j
2
,
j
2
=0
and a bivariate distribution of the form f
3
(r
1
, r
2
) = P(R
1t
= r
1
, R
2t
= r
2
). The functions
f
1
(·) and f
2
(·) are the pmfs of a convolution of two binomial variates. Thus, the conditional
density takes the form
g
1
g
2
f (y
t
|y
t−1
, θ) = f
1
(k)f
2
(s)f
3
(y
1t
− k, y
2t
− s),
k=0
s=0
where g
1
= min(y
1t
, y
1,t−1
) and g
2
= min(y
2t
, y
2,t−1
). Maximum likelihood estimates of
the vector of unknown parameters θ can be obtained by maximization of the conditional
likelihood function
T
L(θ|y) =
f (y
t
|y
t−1
, θ) (19.8)
t=1
for some initial value y
0
. The asymptotic normality of the conditional maximum likelihood
estimate θ
ˆ
has been shown in Franke and Rao (1995) after imposing a set of regular-
ity conditions and applying the results of Billingsley (1961) for the estimation of Markov
processes.
Numerical maximization of (19.8) is straightforward with standard statistical packages.
The binomial convolution implies nite summation and hence it is feasible. Note also that
since the pgf of a binomial distribution is a polynomial, one can derive the pmf of the convo-
lution easily via polynomial multiplication using packages in R. Depending on the choice
for the innovation distribution, the conditional maximum likelihood (CML) approach can
be applied. In Pedeli and Karlis (2013c), a bivariate Poisson and a bivariate negative bino-
mial distribution were used. For the parametric models prediction was discussed. An
interesting result is that for the bivariate Poisson innovations the univariate series have
a Hermite marginal distribution. In Karlis and Pedeli (2013), a copula-based bivariate
innovation distribution was used allowing negative cross-correlation.
When moving to the multivariate case things become more demanding. First of all, a
multivariate discrete distribution is needed for the innovations. As discussed in Section
19.2, such models can be complicated. In Pedeli and Karlis (2013a), a multivariate Pois-
son distribution is assumed with a diagonal matrix A. Even in this case, the pmf of
the multivariate Poisson distribution is demanding since multiple summation is needed.
The conditional likelihood can be derived as in the bivariate case but now this is a con-
volution of several binomials and a multivariate discrete distribution. Alternatively, a
composite likelihood approach can be used. Composite likelihood methods are based
on the idea of constructing lower-dimensional score functions that still contain enough
information about the structure considered but they are computationally more tractable
(Varin, 2008). See also Davis and Yau (2011) for asymptotic properties of composite
likelihood methods applied to linear time series models.