Page 383 - Special Topic Session (STS) - Volume 2
P. 383

STS506 L. Leticia R. et al.
            “memory” of its hidden state about previously seen sequences (Figure 1a).
            These features are combined with the structural properties of the pertaining
                           
                     
            graph ( , … ,  ) and  feed  into  a  dense  layer  that  generates  the  model’s
                     1
                           
            prediction (Figure 1a). is modified to consider the contact network structure
            or network topology, as constant variables (Figure 1b).

            3.  Result
            3.1 Exact simulations
                To assess the ABC and the proposed modifications, we recreate two set of
            synthetic data on epidemics evolving in two different random networks (each
            with  500  vertices)  with  degree  distribution:  (A)  Poisson  (2.42)  (B)
            Polylogarithmic  (0.1,  2).  For  each  network  we  simulate  the  outbreak
            surveillance information using parameter values  = 0.03 and  = 0.01, and
            initial states ((0), (0), (0)) = ( − 2, 2, 0). The generated data is then set “real
            data”. As Dutta, et al. (2018), we assume the knowledge on the number of
            initial cases in  and  but we do not specify the individuals in these states. In
            contrast with Dutta, et al. (2018), we perform the statistical inference only from
            the surveillance-like reports.
                The ABC-MCMC described in Section 2.2 is implemented with discrepancy
            measure
                                                        (   −    ) 2
                                      (, ) = ∑ √   (  >  0) +  (  =  0)  ,
                                               
                                                                      
            where    >  0  is  a  small  constant  introduced  to  define  (,  )  beyond  the
            observed outbreak span.
                We  use  a  Gaussian  kernel  function  and  the  proposal  distribution  
            corresponds to a mixture of Gaussian densities. This  allows to include two
            distributions with different standard deviations (in our case, 0.005 and 0.1) to
            improve the parametric space exploration. For these experiments, we consider
            that  and  have independent prior exponential distribution with parameters
            3, that is equivalent to be approximately 95% confident that each of the real
            values  are  between  0  and  1.  After  removing  the  burn-in  and  thinning  the
            series, the parameters sampled from the posterior produce the statistics in
            Table 1. We observe that for the two networks, the 95% posterior intervals
            contain the true parameter values.

                              A: Poisson                     B: Polylogarithmic
                  Quantile                                                 

                  2.5 %        0.0250        0.0015          0.0176         0.0008

                  50 %         0.0333        0.0127          0.0249         0.0085
                  97.5 %       0.0622        0.0341          0.0411         0.0269

                               Table 1 Statistics of posterior ABC samples.
                                                               372 | I S I   W S C   2 0 1 9
   378   379   380   381   382   383   384   385   386   387   388