Reconstructing the history of individual transmission events between cases is key to understanding what factors facilitate the spread of an infectious disease. Since conducting extended ...contact-tracing investigations can be logistically challenging and costly, statistical inference methods have been developed to reconstruct transmission trees from onset dates and genetic sequences. However, these methods are not as effective if the mutation rate of the virus is very slow, or if sequencing data is sparse.
We developed the package
o2geosocial to combine variables from routinely collected surveillance data with a simple transmission process model. The model reconstructs transmission trees when full genetic sequences are not available, or uninformative. Our model incorporates the reported age-group, onset date, location and genotype of infected cases to infer probabilistic transmission trees. The package also includes functions to summarise and visualise the inferred cluster size distribution.
The results generated by
o2geosocial can highlight regions where importations repeatedly caused large outbreaks, which may indicate a higher regional susceptibility to infections. It can also be used to generate the individual number of secondary transmissions, and show the features associated with individuals involved in high transmission events.
The package is available for download from the Comprehensive R Archive Network (CRAN) and GitHub.
In this paper, the Monte Carlo Markov chain method for solving the modified anomalous fractional sub-diffusion equation is studied. Most of the previous methods are low in temporal and spatial ...accuracy order. Based on the idea of Monte Carlo Markov chain method and compact finite difference schemes, a probability model for solving the modified anomalous fractional sub-diffusion equation is established. Numerical examples are given to show the feasibility of the proposed scheme. Compared with the compact finite difference method, the present method is truly meshless and is easy to be implemented with high temporal and spatial accuracy order. And it is also applied to solve partial differential equation in irregular domains.
•The Monte Carlo Markov chains method is developed to solve the modified anomalous fractional sub-diffusion equation.•The present method is truly meshless and is easy to be implemented with high temporal and spatial accuracy order.•The convergence orders are analyzed and discussed by calculating the numerical experiments.
Abstract
Several strategies have been recently proposed in order to improve Monte Carlo sampling efficiency using machine learning tools. Here, we challenge these methods by considering a class of ...problems that are known to be exponentially hard to sample using conventional local Monte Carlo at low enough temperatures. In particular, we study the antiferromagnetic Potts model on a random graph, which reduces to the coloring of random graphs at zero temperature. We test several machine-learning-assisted Monte Carlo approaches, and we find that they all fail. Our work thus provides good benchmarks for future proposals for smart sampling algorithms.
Random access has been considered for machine-type communication (MTC). In particular, a random access scheme with a pool of preambles is widely studied, which is similar to multichannel ALOHA, to ...support a number of MTC devices. In this paper, we study optimal approaches for user activity detection in random access with preambles over fading channels. Since the computational complexity grows exponentially with the number of preambles, a low-complexity detector is derived using Markov chain Monte Carlo (MCMC) approaches that can approximately solve an optimal maximum a posteriori detection problem. The resulting MCMC detector can enjoy a trade-off between performance and complexity, while its complexity to obtain a sample is linearly proportional to the number of preambles. A performance analysis for optimal detection is also studied to see the optimal performance. Simulation results confirm that the MCMC detector performs better than compressive sensing-based approaches and can provide a near optimal performance under certain conditions with a reasonable computational complexity.
The parametric estimation of stochastic differential equations (SDEs) has been the subject of intense studies already for several decades. The Heston model, for instance, is based on two coupled SDEs ...and is often used in financial mathematics for the dynamics of asset prices and their volatility. Calibrating it to real data would be very useful in many practical scenarios. It is very challenging, however, since the volatility is not directly observable. In this paper, a complete estimation procedure of the Heston model without and with jumps in the asset prices is presented. Bayesian regression combined with the particle filtering method is used as the estimation framework. Within the framework, we propose a novel approach to handle jumps in order to neutralise their negative impact on the estimates of the key parameters of the model. An improvement in the sampling in the particle filtering method is discussed as well. Our analysis is supported by numerical simulations of the Heston model to investigate the performance of the estimators. In addition, a practical follow-along recipe is given to allow finding adequate estimates from any given data.
Novel pathway optimization methods are presented using the 'Global Calculator' model and webtool
1
to goal-seek within a set of optimization constraints. The Global Calculator (GC) is a model used to ...forecast climate-related develop pathways for the world's energy, food and land systems to 2050. The optimization methods enable the GC's user to specify optimization constraints and return a combination of input parameters that satisfy them. The optimization methods evaluated aim to address the challenge of efficiently navigating the GC's ample parameter space (8e
70
parameter combinations) using Monte Carlo Markov Chains and genetic algorithms. The optimization methods are used to calculate an optimal input combination of the 'lever' settings in the GC that satisfy twelve input constraints while minimizing cumulative CO
2
emissions and maximizing GDP output. This optimal development pathway yields a prediction to 2100 of 2,835 GtCO
2
cumulative emissions and a 3.7% increase in GDP with respect to the "business as usual" pathway defined by the International Energy Agency, the IEA's 6DS scenario, a likely extension of current trends. At a similar or lower ambition level as the benchmark scenarios considered to date (distributed effort, consumer reluctance, low action on forests and consumer activism), the optimal pathway shows a significant decrease in CO
2
emissions and increased GDP. The chosen optimization method presented here enables the generation of optimal, user defined/constrained, bespoke pathways to sustainability, relying on the Global Calculator's whole system approach and assumptions.
Morphological models are commonly used to describe microstructures observed in heterogeneous materials. Usually, these models depend upon a set of parameters that must be chosen carefully to match ...experimental observations conducted on the microstructure. A common approach to perform the parameters determination is to try to minimize an objective function, usually taken to be the discrepancy between measurements computed on the simulations and on the experimental observations, respectively. In this article, we present a Bayesian approach for determining the parameters of morphological models, based upon the definition of a posterior distribution for the parameters. A Monte Carlo Markov Chains (MCMC) algorithm is then used to generate samples from the posterior distribution and to identify a set of optimal parameters. We show on several examples that the Bayesian approach allows us to properly identify the optimal parameters of distinct morphological models and to identify potential correlations between the parameters of the models.
Reconstructing the history of individual transmission events between cases is key to understanding what factors facilitate the spread of an infectious disease. Since conducting extended ...contact-tracing investigations can be logistically challenging and costly, statistical inference methods have been developed to reconstruct transmission trees from onset dates and genetic sequences. However, these methods are not as effective if the mutation rate of the virus is very slow, or if sequencing data is sparse.
We developed the package
o2geosocial to combine variables from routinely collected surveillance data with a simple transmission process model. The model reconstructs transmission trees when full genetic sequences are unavailable, or uninformative. Our model incorporates the reported age-group, onset date, location and genotype of infected cases to infer probabilistic transmission trees. The package also includes functions to summarise and visualise the inferred cluster size distribution.
The results generated by
o2geosocial can highlight regions where importations repeatedly caused large outbreaks, which may indicate a higher regional susceptibility to infections. It can also be used to generate the individual number of secondary transmissions, and show the features associated with individuals involved in high transmission events.
The package is available for download from the Comprehensive R Archive Network (CRAN) and GitHub.