176 Handbook of Discrete-Valued Time Series
With the unobservable or latent threshold variable vector u
t
= (u
1
, ..., u
T
), equations (8.22)
and (8.23) can be rewritten as
−1/2
u
t
= z
t
β + x
t
+ λ
t
t
, (8.26)
ν
ν
λ
t
∼ G
, , (8.27)
2 2
where the innovations
t
are assumed to be mutually independent and normally dis-
tributed with mean zero and unit variance. A similar data augmentation scheme can be
used if the link function is logistic instead of probit (Polson et al., 2013). Equations (8.24)
and (8.26), conditioned on δ, the vector β, and the mixing variable λ
t
, jointly represent a lin-
ear state space model, so the algorithm of de Jong and Shephard (1995) is used to simulate
the states.
We t an SSM to the binary observations using probit and Student-t links. For each case,
we conducted MCMC simulation for 50,000 iterations. For both cases, the rst 10,000 draws
were discarded as a burn-in period. In order to reduce the autocorrelation between suc-
cessive values of the simulated chain, every 20th value of the chain was stored. Posterior
means and 95% credible intervals were calculated using the resulting 2000 values.
For SMC, we ran the adaptive random walk Metropolis for 300,000 iterations and dis-
carded the rst 100,000 while the particle lter was run with 100 particles only. Note that
there are some differences in the estimates, but not much on the persistence parameter φ.
The R-INLA commands required for the analysis are described in the Appendix.
8.5.1.2 Results
The results obtained with the different approximation schemes are summarized later in
this chapter. They show reasonable agreement between the different schemes for any given
model. Results for the Student-t link are not presented for the INLA analysis because this
option is not yet available in R-INLA.
The main results for the static parameters are summarized in Tables 8.1 and 8.2. In both
models, the posterior means of δ are around 0.93–0.95, showing higher persistence of the
autoregressive parameter for states variables and thus in the binary time series. The heav-
iness of the tails is measured by the shape parameter ν in the BSSM-T model. In Table 8.2,
the posterior mean of ν is around 10–15. This result seems to indicate that the measure-
ment errors of the u
t
threshold variables are better explained by heavy-tailed distributions,
as a consequence the t-links could be more convenient than the probit link. We found
empirically that the inuence of the number of body movements (z
1
) is marginal, since
the corresponding 95% credible intervals for β
1
contain the zero value. On the other hand,
the inuence of the number of body movements not due to sucking (x
2
) is detected to be
statistically signicant. The negative value of the posterior mean for β
2
shows that a higher
number of body movements not due to sucking will reduce the probability of the infant
being in REM sleep.
Figures 8.1 and 8.2 show the posterior smoothed means of the probabilities π
t
for both
links considered. They show substantial agreement between the different schemes for both
models. A more thorough comparison is also possible using t criteria such as BIC or DIC.
We found some differences between the ts from the different models, but in general the
results are in accordance with Czado and Song (2008). Such difference between the methods