Page 117 - Contributed Paper Session (CPS) - Volume 3
P. 117

CPS1956 Areti B. et al.
            a  spatially structured term following an  intrinsic conditional autoregressive
            prior   ∼ ( ,   )  where    is  the  adjacency  matrix  specifying  the
                                  2
                                 
                   
            spatial neighbourhood. The temporal component   and the age component
                                                              
              both  follow  a  Gaussian  random  walk  prior  of  order  1  (RW1)  which  is
             
            implemented as an ICAR prior. The interacton terms   and   follow normal
                                                                ,
                                                                        ,
            iid priors, while the interaction term   follows an age specific random walk
                                                 ,
            prior of order 1 (RW1). All model hyperparameters follow a weakly informative
            half Normal prior. The analysis was carried out using integrated nested Laplace
            approximations (INLA), which is an alternative to the traditional Markov chain
            Monte carlo (McMC) methods for Bayesian inference specifically designed for
            latent Gaussian models (Rue and Held 2005; Blangiardo et al. 2013). In the
            context of spatio-temporal modelling, several Bayesian detection rules exist,
            and a trade off between true and false positives is usually desired. A common
            practice is to use Bayesian exceedance probabilities (Richardson et al. 2004;
            Lawson 2013). These base the detection rule on a cut off value on the posterior
            probability that the space-time interaction is above a reference threshold. The
            choice of the reference threshold value and the cut off value depends on the
            data,  i.e.  number  of  areas,  time  points  and  magnitude  of  expected  cases
            (Ugarte  et  al.  2009).  In  our  work,  we  were  interested  in  detecting  unusual
            behaviour both in terms of increasing or decreasing trend compared to the
            national one and we chose the reference threshold to be 1, and cut off values
            0.05 and 0.95. Unusual areas were detected in terms of space-time interactions
            over  the  last  5  years  of  the  time  period,  i.e.  2013-2017,
            ( Prob(exp( >1))<0.05 or  Prob(exp( >1))>0.95 ).
                                                  ,
                        ,
            National mortality rates were calculated as the population-weighted average
            of age-specific mortality rates. Finally, mortality rates were converted to life
            expectancy using life table methods (Preston, Heuveline and Guillot 2001)

            3.   Result
                Under  the  detection  rule  that  we  specified,  we  found  21  areas  to  be
            unusual  over  the  last  5  years  of  the  study  period  with  an  exceedance
            probability  below  0.05,  suggesting  a  decreased  slope  compared  to  the
            national  one,  and  32  areas  with  an  exceedance  probability  above  0.95,
            suggesting an increased slope respectively (Figure 1).















                                                               106 | I S I   W S C   2 0 1 9
   112   113   114   115   116   117   118   119   120   121   122