Heterogenous Agents Model with the Worst Out Algorithm

Czech Academy of Sciences Publication Activity Database

Vácha, Lukáš; Vošvrda, Miloslav

-, č. 8 (2006), s. 3-19 ISSN 1801-5999 Institutional research plan: CEZ:AV0Z10750506 Keywords : efficient market hypothesis * fractal market hypothesis * agents' investment horizons * agents' trading strategies * technical trading rules * heterogeneous agent model with stochastic memory * Worst out algorithm Subject RIV: AH - Economics

Stochastic lattice model of synaptic membrane protein domains.

Science.gov (United States)

Li, Yiwei; Kahraman, Osman; Haselwandter, Christoph A

2017-05-01

Neurotransmitter receptor molecules, concentrated in synaptic membrane domains along with scaffolds and other kinds of proteins, are crucial for signal transmission across chemical synapses. In common with other membrane protein domains, synaptic domains are characterized by low protein copy numbers and protein crowding, with rapid stochastic turnover of individual molecules. We study here in detail a stochastic lattice model of the receptor-scaffold reaction-diffusion dynamics at synaptic domains that was found previously to capture, at the mean-field level, the self-assembly, stability, and characteristic size of synaptic domains observed in experiments. We show that our stochastic lattice model yields quantitative agreement with mean-field models of nonlinear diffusion in crowded membranes. Through a combination of analytic and numerical solutions of the master equation governing the reaction dynamics at synaptic domains, together with kinetic Monte Carlo simulations, we find substantial discrepancies between mean-field and stochastic models for the reaction dynamics at synaptic domains. Based on the reaction and diffusion properties of synaptic receptors and scaffolds suggested by previous experiments and mean-field calculations, we show that the stochastic reaction-diffusion dynamics of synaptic receptors and scaffolds provide a simple physical mechanism for collective fluctuations in synaptic domains, the molecular turnover observed at synaptic domains, key features of the observed single-molecule trajectories, and spatial heterogeneity in the effective rates at which receptors and scaffolds are recycled at the cell membrane. Our work sheds light on the physical mechanisms and principles linking the collective properties of membrane protein domains to the stochastic dynamics that rule their molecular components.

Measuring energy efficiency under heterogeneous technologies using a latent class stochastic frontier approach: An application to Chinese energy economy

International Nuclear Information System (INIS)

Lin, Boqiang; Du, Kerui

2014-01-01

The importance of technology heterogeneity in estimating economy-wide energy efficiency has been emphasized by recent literature. Some studies use the metafrontier analysis approach to estimate energy efficiency. However, for such studies, some reliable priori information is needed to divide the sample observations properly, which causes a difficulty in unbiased estimation of energy efficiency. Moreover, separately estimating group-specific frontiers might lose some common information across different groups. In order to overcome these weaknesses, this paper introduces a latent class stochastic frontier approach to measure energy efficiency under heterogeneous technologies. An application of the proposed model to Chinese energy economy is presented. Results show that the overall energy efficiency of China's provinces is not high, with an average score of 0.632 during the period from 1997 to 2010. - Highlights: • We introduce a latent class stochastic frontier approach to measure energy efficiency. • Ignoring technological heterogeneity would cause biased estimates of energy efficiency. • An application of the proposed model to Chinese energy economy is presented. • There is still a long way for China to develop an energy efficient regime

The effect of spatial heterogeneity on the extinction transition in stochastic population dynamics

International Nuclear Information System (INIS)

Kessler, David A; Shnerb, Nadav M

2009-01-01

Stochastic logistic-type growth on a static heterogeneous substrate is studied both above and below the drift-induced delocalization transition. Using agent-based simulations, the delocalization of the highest eigenfunction of the deterministic operator is connected with the large N limit of the stochastic theory. It is seen that the localization length of the deterministic theory controls the divergence of the spatial correlation length with N at the transition. It is argued that, in the presence of a strong wind, the extinction transition belongs to the directed percolation universality class, as any finite colony made of discrete agents is washed away from a heterogeneity with compact support. Some of the difficulties in the analysis of the extinction transition in the presence of a weak wind, where there is a localized active state, are discussed.

Bias and Uncertainty in Regression-Calibrated Models of Groundwater Flow in Heterogeneous Media

DEFF Research Database (Denmark)

Cooley, R.L.; Christensen, Steen

2006-01-01

by a lumped or smoothed m-dimensional approximation γθ*, where γ is an interpolation matrix and θ* is a stochastic vector of parameters. Vector θ* has small enough dimension to allow its estimation with the available data. The consequence of the replacement is that model function f(γθ*) written in terms......Groundwater models need to account for detailed but generally unknown spatial variability (heterogeneity) of the hydrogeologic model inputs. To address this problem we replace the large, m-dimensional stochastic vector β that reflects both small and large scales of heterogeneity in the inputs...... small. Model error is accounted for in the weighted nonlinear regression methodology developed to estimate θ* and assess model uncertainties by incorporating the second-moment matrix of the model errors into the weight matrix. Techniques developed by statisticians to analyze classical nonlinear...

A non-linear dimension reduction methodology for generating data-driven stochastic input models

Science.gov (United States)

Ganapathysubramanian, Baskar; Zabaras, Nicholas

2008-06-01

Stochastic analysis of random heterogeneous media (polycrystalline materials, porous media, functionally graded materials) provides information of significance only if realistic input models of the topology and property variations are used. This paper proposes a framework to construct such input stochastic models for the topology and thermal diffusivity variations in heterogeneous media using a data-driven strategy. Given a set of microstructure realizations (input samples) generated from given statistical information about the medium topology, the framework constructs a reduced-order stochastic representation of the thermal diffusivity. This problem of constructing a low-dimensional stochastic representation of property variations is analogous to the problem of manifold learning and parametric fitting of hyper-surfaces encountered in image processing and psychology. Denote by M the set of microstructures that satisfy the given experimental statistics. A non-linear dimension reduction strategy is utilized to map M to a low-dimensional region, A. We first show that M is a compact manifold embedded in a high-dimensional input space Rn. An isometric mapping F from M to a low-dimensional, compact, connected set A⊂Rd(d≪n) is constructed. Given only a finite set of samples of the data, the methodology uses arguments from graph theory and differential geometry to construct the isometric transformation F:M→A. Asymptotic convergence of the representation of M by A is shown. This mapping F serves as an accurate, low-dimensional, data-driven representation of the property variations. The reduced-order model of the material topology and thermal diffusivity variations is subsequently used as an input in the solution of stochastic partial differential equations that describe the evolution of dependant variables. A sparse grid collocation strategy (Smolyak algorithm) is utilized to solve these stochastic equations efficiently. We showcase the methodology by constructing low

A non-linear dimension reduction methodology for generating data-driven stochastic input models

International Nuclear Information System (INIS)

Ganapathysubramanian, Baskar; Zabaras, Nicholas

2008-01-01

Stochastic analysis of random heterogeneous media (polycrystalline materials, porous media, functionally graded materials) provides information of significance only if realistic input models of the topology and property variations are used. This paper proposes a framework to construct such input stochastic models for the topology and thermal diffusivity variations in heterogeneous media using a data-driven strategy. Given a set of microstructure realizations (input samples) generated from given statistical information about the medium topology, the framework constructs a reduced-order stochastic representation of the thermal diffusivity. This problem of constructing a low-dimensional stochastic representation of property variations is analogous to the problem of manifold learning and parametric fitting of hyper-surfaces encountered in image processing and psychology. Denote by M the set of microstructures that satisfy the given experimental statistics. A non-linear dimension reduction strategy is utilized to map M to a low-dimensional region, A. We first show that M is a compact manifold embedded in a high-dimensional input space R n . An isometric mapping F from M to a low-dimensional, compact, connected set A is contained in R d (d<< n) is constructed. Given only a finite set of samples of the data, the methodology uses arguments from graph theory and differential geometry to construct the isometric transformation F:M→A. Asymptotic convergence of the representation of M by A is shown. This mapping F serves as an accurate, low-dimensional, data-driven representation of the property variations. The reduced-order model of the material topology and thermal diffusivity variations is subsequently used as an input in the solution of stochastic partial differential equations that describe the evolution of dependant variables. A sparse grid collocation strategy (Smolyak algorithm) is utilized to solve these stochastic equations efficiently. We showcase the methodology

Stochastic neuron models

CERN Document Server

Greenwood, Priscilla E

2016-01-01

This book describes a large number of open problems in the theory of stochastic neural systems, with the aim of enticing probabilists to work on them. This includes problems arising from stochastic models of individual neurons as well as those arising from stochastic models of the activities of small and large networks of interconnected neurons. The necessary neuroscience background to these problems is outlined within the text, so readers can grasp the context in which they arise. This book will be useful for graduate students and instructors providing material and references for applying probability to stochastic neuron modeling. Methods and results are presented, but the emphasis is on questions where additional stochastic analysis may contribute neuroscience insight. An extensive bibliography is included. Dr. Priscilla E. Greenwood is a Professor Emerita in the Department of Mathematics at the University of British Columbia. Dr. Lawrence M. Ward is a Professor in the Department of Psychology and the Brain...

Stochastic biomathematical models with applications to neuronal modeling

CERN Document Server

Batzel, Jerry; Ditlevsen, Susanne

2013-01-01

Stochastic biomathematical models are becoming increasingly important as new light is shed on the role of noise in living systems. In certain biological systems, stochastic effects may even enhance a signal, thus providing a biological motivation for the noise observed in living systems. Recent advances in stochastic analysis and increasing computing power facilitate the analysis of more biophysically realistic models, and this book provides researchers in computational neuroscience and stochastic systems with an overview of recent developments. Key concepts are developed in chapters written by experts in their respective fields. Topics include: one-dimensional homogeneous diffusions and their boundary behavior, large deviation theory and its application in stochastic neurobiological models, a review of mathematical methods for stochastic neuronal integrate-and-fire models, stochastic partial differential equation models in neurobiology, and stochastic modeling of spreading cortical depression.

Modeling and forecasting energy consumption for heterogeneous buildings using a physical–statistical approach

International Nuclear Information System (INIS)

Lü, Xiaoshu; Lu, Tao; Kibert, Charles J.; Viljanen, Martti

2015-01-01

Highlights: • This paper presents a new modeling method to forecast energy demands. • The model is based on physical–statistical approach to improving forecast accuracy. • A new method is proposed to address the heterogeneity challenge. • Comparison with measurements shows accurate forecasts of the model. • The first physical–statistical/heterogeneous building energy modeling approach is proposed and validated. - Abstract: Energy consumption forecasting is a critical and necessary input to planning and controlling energy usage in the building sector which accounts for 40% of the world’s energy use and the world’s greatest fraction of greenhouse gas emissions. However, due to the diversity and complexity of buildings as well as the random nature of weather conditions, energy consumption and loads are stochastic and difficult to predict. This paper presents a new methodology for energy demand forecasting that addresses the heterogeneity challenges in energy modeling of buildings. The new method is based on a physical–statistical approach designed to account for building heterogeneity to improve forecast accuracy. The physical model provides a theoretical input to characterize the underlying physical mechanism of energy flows. Then stochastic parameters are introduced into the physical model and the statistical time series model is formulated to reflect model uncertainties and individual heterogeneity in buildings. A new method of model generalization based on a convex hull technique is further derived to parameterize the individual-level model parameters for consistent model coefficients while maintaining satisfactory modeling accuracy for heterogeneous buildings. The proposed method and its validation are presented in detail for four different sports buildings with field measurements. The results show that the proposed methodology and model can provide a considerable improvement in forecasting accuracy

Modeling flow in fractured medium. Uncertainty analysis with stochastic continuum approach

International Nuclear Information System (INIS)

Niemi, A.

1994-01-01

For modeling groundwater flow in formation-scale fractured media, no general method exists for scaling the highly heterogeneous hydraulic conductivity data to model parameters. The deterministic approach is limited in representing the heterogeneity of a medium and the application of fracture network models has both conceptual and practical limitations as far as site-scale studies are concerned. The study investigates the applicability of stochastic continuum modeling at the scale of data support. No scaling of the field data is involved, and the original variability is preserved throughout the modeling. Contributions of various aspects to the total uncertainty in the modeling prediction can also be determined with this approach. Data from five crystalline rock sites in Finland are analyzed. (107 refs., 63 figs., 7 tabs.)

The Limit Behavior of a Stochastic Logistic Model with Individual Time-Dependent Rates

Directory of Open Access Journals (Sweden)

Yilun Shang

2013-01-01

Full Text Available We investigate a variant of the stochastic logistic model that allows individual variation and time-dependent infection and recovery rates. The model is described as a heterogeneous density dependent Markov chain. We show that the process can be approximated by a deterministic process defined by an integral equation as the population size grows.

Heterogeneous Agents Model with the Worst Out Algorithm

Czech Academy of Sciences Publication Activity Database

Vošvrda, Miloslav; Vácha, Lukáš

I, č. 1 (2007), s. 54-66 ISSN 1802-4696 R&D Projects: GA MŠk(CZ) LC06075; GA ČR(CZ) GA402/06/0990 Grant - others:GA UK(CZ) 454/2004/A-EK/FSV Institutional research plan: CEZ:AV0Z10750506 Keywords : Efficient Market s Hypothesis * Fractal Market Hypothesis * agents' investment horizons * agents' trading strategies * technical trading rules * heterogeneous agent model with stochastic memory * Worst out Algorithm Subject RIV: AH - Economics

Stochastic Still Water Response Model

DEFF Research Database (Denmark)

Friis-Hansen, Peter; Ditlevsen, Ove Dalager

2002-01-01

In this study a stochastic field model for the still water loading is formulated where the statistics (mean value, standard deviation, and correlation) of the sectional forces are obtained by integration of the load field over the relevant part of the ship structure. The objective of the model is...... out that an important parameter of the stochastic cargo field model is the mean number of containers delivered by each customer.......In this study a stochastic field model for the still water loading is formulated where the statistics (mean value, standard deviation, and correlation) of the sectional forces are obtained by integration of the load field over the relevant part of the ship structure. The objective of the model...... is to establish the stochastic load field conditional on a given draft and trim of the vessel. The model contributes to a realistic modelling of the stochastic load processes to be used in a reliability evaluation of the ship hull. Emphasis is given to container vessels. The formulation of the model for obtaining...

Stochastic-field cavitation model

International Nuclear Information System (INIS)

Dumond, J.; Magagnato, F.; Class, A.

2013-01-01

Nonlinear phenomena can often be well described using probability density functions (pdf) and pdf transport models. Traditionally, the simulation of pdf transport requires Monte-Carlo codes based on Lagrangian “particles” or prescribed pdf assumptions including binning techniques. Recently, in the field of combustion, a novel formulation called the stochastic-field method solving pdf transport based on Eulerian fields has been proposed which eliminates the necessity to mix Eulerian and Lagrangian techniques or prescribed pdf assumptions. In the present work, for the first time the stochastic-field method is applied to multi-phase flow and, in particular, to cavitating flow. To validate the proposed stochastic-field cavitation model, two applications are considered. First, sheet cavitation is simulated in a Venturi-type nozzle. The second application is an innovative fluidic diode which exhibits coolant flashing. Agreement with experimental results is obtained for both applications with a fixed set of model constants. The stochastic-field cavitation model captures the wide range of pdf shapes present at different locations

Stochastic-field cavitation model

Science.gov (United States)

Dumond, J.; Magagnato, F.; Class, A.

2013-07-01

Nonlinear phenomena can often be well described using probability density functions (pdf) and pdf transport models. Traditionally, the simulation of pdf transport requires Monte-Carlo codes based on Lagrangian "particles" or prescribed pdf assumptions including binning techniques. Recently, in the field of combustion, a novel formulation called the stochastic-field method solving pdf transport based on Eulerian fields has been proposed which eliminates the necessity to mix Eulerian and Lagrangian techniques or prescribed pdf assumptions. In the present work, for the first time the stochastic-field method is applied to multi-phase flow and, in particular, to cavitating flow. To validate the proposed stochastic-field cavitation model, two applications are considered. First, sheet cavitation is simulated in a Venturi-type nozzle. The second application is an innovative fluidic diode which exhibits coolant flashing. Agreement with experimental results is obtained for both applications with a fixed set of model constants. The stochastic-field cavitation model captures the wide range of pdf shapes present at different locations.

An entropy approach to size and variance heterogeneity

NARCIS (Netherlands)

Balasubramanyan, L.; Stefanou, S.E.; Stokes, J.R.

2012-01-01

In this paper, we investigate the effect of bank size differences on cost efficiency heterogeneity using a heteroskedastic stochastic frontier model. This model is implemented by using an information theoretic maximum entropy approach. We explicitly model both bank size and variance heterogeneity

Numerical calculations on heterogeneity of groundwater flow

International Nuclear Information System (INIS)

Follin, S.

1992-01-01

The upscaling of model parameters is a key issue in many research fields concerned with parameter heterogeneity. The upscaling process allows for fewer model blocks and relaxes the numerical problems caused by high contrasts in the hydraulic conductivity. The trade-offs are dependent on the object but the general drawback is an increasing uncertainty about the representativeness. The present study deals with numerical calculations of heterogeneity of groundwater flow and solute transport in hypothetical blocks of fractured hard rock in a '3m scale' and addresses both conceptual and practical problems in numerical simulation. Evidence that the hydraulic conductivity (K) of the rock mass between major fracture zones is highly heterogeneous in a 3m scale is provided by a large number of field investigations. The present uses the documented heterogeneity and investigates flow and transport in a two-dimensional stochastic continuum characterized by a variance in Y = In(K) of σ y 2 = 16, corresponding to about 12 log 10 cycles in K. The study considers anisotropy, channelling, non-Fickian and Fickian transport, and conditional simulation. The major conclusions are: * heterogeneity gives rise to anisotropy in the upscaling process, * the choice of support scale is crucial for the modelling of solute transport. As a consequence of the obtained results, a two-dimensional stochastic discontinuum model is presented, which provides a tool for linking stochastic continuum models to discrete fracture network models. (au) (14 figs., 136 refs.)

Research on nonlinear stochastic dynamical price model

International Nuclear Information System (INIS)

Li Jiaorui; Xu Wei; Xie Wenxian; Ren Zhengzheng

2008-01-01

In consideration of many uncertain factors existing in economic system, nonlinear stochastic dynamical price model which is subjected to Gaussian white noise excitation is proposed based on deterministic model. One-dimensional averaged Ito stochastic differential equation for the model is derived by using the stochastic averaging method, and applied to investigate the stability of the trivial solution and the first-passage failure of the stochastic price model. The stochastic price model and the methods presented in this paper are verified by numerical studies

Alternative Asymmetric Stochastic Volatility Models

NARCIS (Netherlands)

M. Asai (Manabu); M.J. McAleer (Michael)

2010-01-01

textabstractThe stochastic volatility model usually incorporates asymmetric effects by introducing the negative correlation between the innovations in returns and volatility. In this paper, we propose a new asymmetric stochastic volatility model, based on the leverage and size effects. The model is

Challenges in conditioning a stochastic geological model of a heterogeneous glacial aquifer to a comprehensive soft data set

DEFF Research Database (Denmark)

Koch, Julian; He, Xin; Jensen, Karsten Høgh

2014-01-01

In traditional hydrogeological investigations, one geological model is often used based on subjective interpretations and sparse data availability. This deterministic approach usually does not account for any uncertainties. Stochastic simulation methods address this problem and can capture......TEM data set are derived from a histogram probability matching method, where resistivity is paired with the corresponding lithology from the categorized borehole data. This study integrates two distinct data sources into the stochastic modeling process that represent two extremes of the conditioning...

Estimation of heterogeneity in malaria transmission by stochastic modelling of apparent deviations from mass action kinetics

Directory of Open Access Journals (Sweden)

Smith Thomas A

2008-01-01

Full Text Available Abstract Background Quantifying heterogeneity in malaria transmission is a prerequisite for accurate predictive mathematical models, but the variance in field measurements of exposure overestimates true micro-heterogeneity because it is inflated to an uncertain extent by sampling variation. Descriptions of field data also suggest that the rate of Plasmodium falciparum infection is not proportional to the intensity of challenge by infectious vectors. This appears to violate the principle of mass action that is implied by malaria biology. Micro-heterogeneity may be the reason for this anomaly. It is proposed that the level of micro-heterogeneity can be estimated from statistical models that estimate the amount of variation in transmission most compatible with a mass-action model for the relationship of infection to exposure. Methods The relationship between the entomological inoculation rate (EIR for falciparum malaria and infection risk was reanalysed using published data for cohorts of children in Saradidi (western Kenya. Infection risk was treated as binomially distributed, and measurement-error (Poisson and negative binomial models were considered for the EIR. Models were fitted using Bayesian Markov chain Monte Carlo algorithms and model fit compared for models that assume either mass-action kinetics, facilitation, competition or saturation of the infection process with increasing EIR. Results The proportion of inocula that resulted in infection in Saradidi was inversely related to the measured intensity of challenge. Models of facilitation showed, therefore, a poor fit to the data. When sampling error in the EIR was neglected, either competition or saturation needed to be incorporated in the model in order to give a good fit. Negative binomial models for the error in exposure could achieve a comparable fit while incorporating the more parsimonious and biologically plausible mass action assumption. Models that assume negative binomial micro-heterogeneity

Stochastic modeling for the expression of a gene regulated by competing transcription factors.

Directory of Open Access Journals (Sweden)

Hsih-Te Yang

Full Text Available It is widely accepted that gene expression regulation is a stochastic event. The common approach for its computer simulation requires detailed information on the interactions of individual molecules, which is often not available for the analyses of biological experiments. As an alternative approach, we employed a more intuitive model to simulate the experimental result, the Markov-chain model, in which a gene is regulated by activators and repressors, which bind the same site in a mutually exclusive manner. Our stochastic simulation in the presence of both activators and repressors predicted a Hill-coefficient of the dose-response curve closer to the experimentally observed value than the calculated value based on the simple additive effects of activators alone and repressors alone. The simulation also reproduced the heterogeneity of gene expression levels among individual cells observed by Fluorescence Activated Cell Sorting analysis. Therefore, our approach may help to apply stochastic simulations to broader experimental data.

Stochasticity Modeling in Memristors

KAUST Repository

Naous, Rawan

2015-10-26

Diverse models have been proposed over the past years to explain the exhibiting behavior of memristors, the fourth fundamental circuit element. The models varied in complexity ranging from a description of physical mechanisms to a more generalized mathematical modeling. Nonetheless, stochasticity, a widespread observed phenomenon, has been immensely overlooked from the modeling perspective. This inherent variability within the operation of the memristor is a vital feature for the integration of this nonlinear device into the stochastic electronics realm of study. In this paper, experimentally observed innate stochasticity is modeled in a circuit compatible format. The model proposed is generic and could be incorporated into variants of threshold-based memristor models in which apparent variations in the output hysteresis convey the switching threshold shift. Further application as a noise injection alternative paves the way for novel approaches in the fields of neuromorphic engineering circuits design. On the other hand, extra caution needs to be paid to variability intolerant digital designs based on non-deterministic memristor logic.

Stochasticity Modeling in Memristors

KAUST Repository

Naous, Rawan; Al-Shedivat, Maruan; Salama, Khaled N.

2015-01-01

Diverse models have been proposed over the past years to explain the exhibiting behavior of memristors, the fourth fundamental circuit element. The models varied in complexity ranging from a description of physical mechanisms to a more generalized mathematical modeling. Nonetheless, stochasticity, a widespread observed phenomenon, has been immensely overlooked from the modeling perspective. This inherent variability within the operation of the memristor is a vital feature for the integration of this nonlinear device into the stochastic electronics realm of study. In this paper, experimentally observed innate stochasticity is modeled in a circuit compatible format. The model proposed is generic and could be incorporated into variants of threshold-based memristor models in which apparent variations in the output hysteresis convey the switching threshold shift. Further application as a noise injection alternative paves the way for novel approaches in the fields of neuromorphic engineering circuits design. On the other hand, extra caution needs to be paid to variability intolerant digital designs based on non-deterministic memristor logic.

Stochastic equilibria of an asset pricing model with heterogeneous beliefs and random dividends

NARCIS (Netherlands)

Zhu, M.; Wang, D.; Guo, M.

2011-01-01

We investigate dynamical properties of a heterogeneous agent model with random dividends and further study the relationship between dynamical properties of the random model and those of the corresponding deterministic skeleton, which is obtained by setting the random dividends as their constant mean

A Simple, Realistic Stochastic Model of Gastric Emptying.

Directory of Open Access Journals (Sweden)

Jiraphat Yokrattanasak

Full Text Available Several models of Gastric Emptying (GE have been employed in the past to represent the rate of delivery of stomach contents to the duodenum and jejunum. These models have all used a deterministic form (algebraic equations or ordinary differential equations, considering GE as a continuous, smooth process in time. However, GE is known to occur as a sequence of spurts, irregular both in size and in timing. Hence, we formulate a simple stochastic process model, able to represent the irregular decrements of gastric contents after a meal. The model is calibrated on existing literature data and provides consistent predictions of the observed variability in the emptying trajectories. This approach may be useful in metabolic modeling, since it describes well and explains the apparently heterogeneous GE experimental results in situations where common gastric mechanics across subjects would be expected.

Stochastic models: theory and simulation.

Energy Technology Data Exchange (ETDEWEB)

Field, Richard V., Jr.

2008-03-01

Many problems in applied science and engineering involve physical phenomena that behave randomly in time and/or space. Examples are diverse and include turbulent flow over an aircraft wing, Earth climatology, material microstructure, and the financial markets. Mathematical models for these random phenomena are referred to as stochastic processes and/or random fields, and Monte Carlo simulation is the only general-purpose tool for solving problems of this type. The use of Monte Carlo simulation requires methods and algorithms to generate samples of the appropriate stochastic model; these samples then become inputs and/or boundary conditions to established deterministic simulation codes. While numerous algorithms and tools currently exist to generate samples of simple random variables and vectors, no cohesive simulation tool yet exists for generating samples of stochastic processes and/or random fields. There are two objectives of this report. First, we provide some theoretical background on stochastic processes and random fields that can be used to model phenomena that are random in space and/or time. Second, we provide simple algorithms that can be used to generate independent samples of general stochastic models. The theory and simulation of random variables and vectors is also reviewed for completeness.

Stochastic Optimization of Wind Turbine Power Factor Using Stochastic Model of Wind Power

DEFF Research Database (Denmark)

Chen, Peiyuan; Siano, Pierluigi; Bak-Jensen, Birgitte

2010-01-01

This paper proposes a stochastic optimization algorithm that aims to minimize the expectation of the system power losses by controlling wind turbine (WT) power factors. This objective of the optimization is subject to the probability constraints of bus voltage and line current requirements....... The optimization algorithm utilizes the stochastic models of wind power generation (WPG) and load demand to take into account their stochastic variation. The stochastic model of WPG is developed on the basis of a limited autoregressive integrated moving average (LARIMA) model by introducing a crosscorrelation...... structure to the LARIMA model. The proposed stochastic optimization is carried out on a 69-bus distribution system. Simulation results confirm that, under various combinations of WPG and load demand, the system power losses are considerably reduced with the optimal setting of WT power factor as compared...

Stochastic diffusion models for substitutable technological innovations

NARCIS (Netherlands)

Wang, L.; Hu, B.; Yu, X.

2004-01-01

Based on the analysis of firms' stochastic adoption behaviour, this paper first points out the necessity to build more practical stochastic models. And then, stochastic evolutionary models are built for substitutable innovation diffusion system. Finally, through the computer simulation of the

Skeeter Buster: a stochastic, spatially explicit modeling tool for studying Aedes aegypti population replacement and population suppression strategies.

Directory of Open Access Journals (Sweden)

Krisztian Magori

2009-09-01

Full Text Available Dengue is the most important mosquito-borne viral disease affecting humans. The only prevention measure currently available is the control of its vectors, primarily Aedes aegypti. Recent advances in genetic engineering have opened the possibility for a new range of control strategies based on genetically modified mosquitoes. Assessing the potential efficacy of genetic (and conventional strategies requires the availability of modeling tools that accurately describe the dynamics and genetics of Ae. aegypti populations.We describe in this paper a new modeling tool of Ae. aegypti population dynamics and genetics named Skeeter Buster. This model operates at the scale of individual water-filled containers for immature stages and individual properties (houses for adults. The biology of cohorts of mosquitoes is modeled based on the algorithms used in the non-spatial Container Inhabiting Mosquitoes Simulation Model (CIMSiM. Additional features incorporated into Skeeter Buster include stochasticity, spatial structure and detailed population genetics. We observe that the stochastic modeling of individual containers in Skeeter Buster is associated with a strongly reduced temporal variation in stage-specific population densities. We show that heterogeneity in container composition of individual properties has a major impact on spatial heterogeneity in population density between properties. We detail how adult dispersal reduces this spatial heterogeneity. Finally, we present the predicted genetic structure of the population by calculating F(ST values and isolation by distance patterns, and examine the effects of adult dispersal and container movement between properties.We demonstrate that the incorporated stochasticity and level of spatial detail have major impacts on the simulated population dynamics, which could potentially impact predictions in terms of control measures. The capacity to describe population genetics confers the ability to model the outcome

Sequential neural models with stochastic layers

DEFF Research Database (Denmark)

Fraccaro, Marco; Sønderby, Søren Kaae; Paquet, Ulrich

2016-01-01

How can we efficiently propagate uncertainty in a latent state representation with recurrent neural networks? This paper introduces stochastic recurrent neural networks which glue a deterministic recurrent neural network and a state space model together to form a stochastic and sequential neural...... generative model. The clear separation of deterministic and stochastic layers allows a structured variational inference network to track the factorization of the model's posterior distribution. By retaining both the nonlinear recursive structure of a recurrent neural network and averaging over...

Numerical Simulation of the Heston Model under Stochastic Correlation

Directory of Open Access Journals (Sweden)

Long Teng

2017-12-01

Full Text Available Stochastic correlation models have become increasingly important in financial markets. In order to be able to price vanilla options in stochastic volatility and correlation models, in this work, we study the extension of the Heston model by imposing stochastic correlations driven by a stochastic differential equation. We discuss the efficient algorithms for the extended Heston model by incorporating stochastic correlations. Our numerical experiments show that the proposed algorithms can efficiently provide highly accurate results for the extended Heston by including stochastic correlations. By investigating the effect of stochastic correlations on the implied volatility, we find that the performance of the Heston model can be proved by including stochastic correlations.

Transport properties of stochastic Lorentz models

NARCIS (Netherlands)

Beijeren, H. van

Diffusion processes are considered for one-dimensional stochastic Lorentz models, consisting of randomly distributed fixed scatterers and one moving light particle. In waiting time Lorentz models the light particle makes instantaneous jumps between scatterers after a stochastically distributed

Modeling and analysis of stochastic systems

CERN Document Server

Kulkarni, Vidyadhar G

2011-01-01

Based on the author's more than 25 years of teaching experience, Modeling and Analysis of Stochastic Systems, Second Edition covers the most important classes of stochastic processes used in the modeling of diverse systems, from supply chains and inventory systems to genetics and biological systems. For each class of stochastic process, the text includes its definition, characterization, applications, transient and limiting behavior, first passage times, and cost/reward models. Along with reorganizing the material, this edition revises and adds new exercises and examples. New to the second edi

Quantifying geological uncertainty for flow and transport modeling in multi-modal heterogeneous formations

Science.gov (United States)

Feyen, Luc; Caers, Jef

2006-06-01

In this work, we address the problem of characterizing the heterogeneity and uncertainty of hydraulic properties for complex geological settings. Hereby, we distinguish between two scales of heterogeneity, namely the hydrofacies structure and the intrafacies variability of the hydraulic properties. We employ multiple-point geostatistics to characterize the hydrofacies architecture. The multiple-point statistics are borrowed from a training image that is designed to reflect the prior geological conceptualization. The intrafacies variability of the hydraulic properties is represented using conventional two-point correlation methods, more precisely, spatial covariance models under a multi-Gaussian spatial law. We address the different levels and sources of uncertainty in characterizing the subsurface heterogeneity, and explore their effect on groundwater flow and transport predictions. Typically, uncertainty is assessed by way of many images, termed realizations, of a fixed statistical model. However, in many cases, sampling from a fixed stochastic model does not adequately represent the space of uncertainty. It neglects the uncertainty related to the selection of the stochastic model and the estimation of its input parameters. We acknowledge the uncertainty inherent in the definition of the prior conceptual model of aquifer architecture and in the estimation of global statistics, anisotropy, and correlation scales. Spatial bootstrap is used to assess the uncertainty of the unknown statistical parameters. As an illustrative example, we employ a synthetic field that represents a fluvial setting consisting of an interconnected network of channel sands embedded within finer-grained floodplain material. For this highly non-stationary setting we quantify the groundwater flow and transport model prediction uncertainty for various levels of hydrogeological uncertainty. Results indicate the importance of accurately describing the facies geometry, especially for transport

Stochastic variation in telomere shortening rate causes heterogeneity of human fibroblast replicative life span.

Science.gov (United States)

Martin-Ruiz, Carmen; Saretzki, Gabriele; Petrie, Joanne; Ladhoff, Juliane; Jeyapalan, Jessie; Wei, Wenyi; Sedivy, John; von Zglinicki, Thomas

2004-04-23

The replicative life span of human fibroblasts is heterogeneous, with a fraction of cells senescing at every population doubling. To find out whether this heterogeneity is due to premature senescence, i.e. driven by a nontelomeric mechanism, fibroblasts with a senescent phenotype were isolated from growing cultures and clones by flow cytometry. These senescent cells had shorter telomeres than their cycling counterparts at all population doubling levels and both in mass cultures and in individual subclones, indicating heterogeneity in the rate of telomere shortening. Ectopic expression of telomerase stabilized telomere length in the majority of cells and rescued them from early senescence, suggesting a causal role of telomere shortening. Under standard cell culture conditions, there was a minor fraction of cells that showed a senescent phenotype and short telomeres despite active telomerase. This fraction increased under chronic mild oxidative stress, which is known to accelerate telomere shortening. It is possible that even high telomerase activity cannot fully compensate for telomere shortening in all cells. The data show that heterogeneity of the human fibroblast replicative life span can be caused by significant stochastic cell-to-cell variation in telomere shortening.

Modelling Cow Behaviour Using Stochastic Automata

DEFF Research Database (Denmark)

Jónsson, Ragnar Ingi

This report covers an initial study on the modelling of cow behaviour using stochastic automata with the aim of detecting lameness. Lameness in cows is a serious problem that needs to be dealt with because it results in less profitable production units and in reduced quality of life...... for the affected livestock. By featuring training data consisting of measurements of cow activity, three different models are obtained, namely an autonomous stochastic automaton, a stochastic automaton with coinciding state and output and an autonomous stochastic automaton with coinciding state and output, all...... of which describe the cows' activity in the two regarded behavioural scenarios, non-lame and lame. Using the experimental measurement data the different behavioural relations for the two regarded behavioural scenarios are assessed. The three models comprise activity within last hour, activity within last...

Electricity price modeling with stochastic time change

International Nuclear Information System (INIS)

Borovkova, Svetlana; Schmeck, Maren Diane

2017-01-01

In this paper, we develop a novel approach to electricity price modeling, based on the powerful technique of stochastic time change. This technique allows us to incorporate the characteristic features of electricity prices (such as seasonal volatility, time varying mean reversion and seasonally occurring price spikes) into the model in an elegant and economically justifiable way. The stochastic time change introduces stochastic as well as deterministic (e.g., seasonal) features in the price process' volatility and in the jump component. We specify the base process as a mean reverting jump diffusion and the time change as an absolutely continuous stochastic process with seasonal component. The activity rate of the stochastic time change can be related to the factors that influence supply and demand. Here we use the temperature as a proxy for the demand and hence, as the driving factor of the stochastic time change, and show that this choice leads to realistic price paths. We derive properties of the resulting price process and develop the model calibration procedure. We calibrate the model to the historical EEX power prices and apply it to generating realistic price paths by Monte Carlo simulations. We show that the simulated price process matches the distributional characteristics of the observed electricity prices in periods of both high and low demand. - Highlights: • We develop a novel approach to electricity price modeling, based on the powerful technique of stochastic time change. • We incorporate the characteristic features of electricity prices, such as seasonal volatility and spikes into the model. • We use the temperature as a proxy for the demand and hence, as the driving factor of the stochastic time change • We derive properties of the resulting price process and develop the model calibration procedure. • We calibrate the model to the historical EEX power prices and apply it to generating realistic price paths.

Stochastic volatility and stochastic leverage

DEFF Research Database (Denmark)

Veraart, Almut; Veraart, Luitgard A. M.

This paper proposes the new concept of stochastic leverage in stochastic volatility models. Stochastic leverage refers to a stochastic process which replaces the classical constant correlation parameter between the asset return and the stochastic volatility process. We provide a systematic...... treatment of stochastic leverage and propose to model the stochastic leverage effect explicitly, e.g. by means of a linear transformation of a Jacobi process. Such models are both analytically tractable and allow for a direct economic interpretation. In particular, we propose two new stochastic volatility...... models which allow for a stochastic leverage effect: the generalised Heston model and the generalised Barndorff-Nielsen & Shephard model. We investigate the impact of a stochastic leverage effect in the risk neutral world by focusing on implied volatilities generated by option prices derived from our new...

Definition and solution of a stochastic inverse problem for the Manning's n parameter field in hydrodynamic models

Science.gov (United States)

Butler, T.; Graham, L.; Estep, D.; Dawson, C.; Westerink, J. J.

2015-04-01

The uncertainty in spatially heterogeneous Manning's n fields is quantified using a novel formulation and numerical solution of stochastic inverse problems for physics-based models. The uncertainty is quantified in terms of a probability measure and the physics-based model considered here is the state-of-the-art ADCIRC model although the presented methodology applies to other hydrodynamic models. An accessible overview of the formulation and solution of the stochastic inverse problem in a mathematically rigorous framework based on measure theory is presented. Technical details that arise in practice by applying the framework to determine the Manning's n parameter field in a shallow water equation model used for coastal hydrodynamics are presented and an efficient computational algorithm and open source software package are developed. A new notion of "condition" for the stochastic inverse problem is defined and analyzed as it relates to the computation of probabilities. This notion of condition is investigated to determine effective output quantities of interest of maximum water elevations to use for the inverse problem for the Manning's n parameter and the effect on model predictions is analyzed.

Stochastic models of cell motility

DEFF Research Database (Denmark)

Gradinaru, Cristian

2012-01-01

Cell motility and migration are central to the development and maintenance of multicellular organisms, and errors during this process can lead to major diseases. Consequently, the mechanisms and phenomenology of cell motility are currently under intense study. In recent years, a new...... interdisciplinary field focusing on the study of biological processes at the nanoscale level, with a range of technological applications in medicine and biological research, has emerged. The work presented in this thesis is at the interface of cell biology, image processing, and stochastic modeling. The stochastic...... models introduced here are based on persistent random motion, which I apply to real-life studies of cell motility on flat and nanostructured surfaces. These models aim to predict the time-dependent position of cell centroids in a stochastic manner, and conversely determine directly from experimental...

Stochastic Modelling of Hydrologic Systems

DEFF Research Database (Denmark)

Jonsdottir, Harpa

2007-01-01

In this PhD project several stochastic modelling methods are studied and applied on various subjects in hydrology. The research was prepared at Informatics and Mathematical Modelling at the Technical University of Denmark. The thesis is divided into two parts. The first part contains...... an introduction and an overview of the papers published. Then an introduction to basic concepts in hydrology along with a description of hydrological data is given. Finally an introduction to stochastic modelling is given. The second part contains the research papers. In the research papers the stochastic methods...... are described, as at the time of publication these methods represent new contribution to hydrology. The second part also contains additional description of software used and a brief introduction to stiff systems. The system in one of the papers is stiff....

Stochasticity and determinism in models of hematopoiesis.

Science.gov (United States)

Kimmel, Marek

2014-01-01

This chapter represents a novel view of modeling in hematopoiesis, synthesizing both deterministic and stochastic approaches. Whereas the stochastic models work in situations where chance dominates, for example when the number of cells is small, or under random mutations, the deterministic models are more important for large-scale, normal hematopoiesis. New types of models are on the horizon. These models attempt to account for distributed environments such as hematopoietic niches and their impact on dynamics. Mixed effects of such structures and chance events are largely unknown and constitute both a challenge and promise for modeling. Our discussion is presented under the separate headings of deterministic and stochastic modeling; however, the connections between both are frequently mentioned. Four case studies are included to elucidate important examples. We also include a primer of deterministic and stochastic dynamics for the reader's use.

Mitochondrial dysfunction accounts for the stochastic heterogeneity in telomere-dependent senescence.

Directory of Open Access Journals (Sweden)

João F Passos

2007-05-01

Full Text Available Aging is an inherently stochastic process, and its hallmark is heterogeneity between organisms, cell types, and clonal populations, even in identical environments. The replicative lifespan of primary human cells is telomere dependent; however, its heterogeneity is not understood. We show that mitochondrial superoxide production increases with replicative age in human fibroblasts despite an adaptive UCP-2-dependent mitochondrial uncoupling. This mitochondrial dysfunction is accompanied by compromised [Ca(2+]i homeostasis and other indicators of a retrograde response in senescent cells. Replicative senescence of human fibroblasts is delayed by mild mitochondrial uncoupling. Uncoupling reduces mitochondrial superoxide generation, slows down telomere shortening, and delays formation of telomeric gamma-H2A.X foci. This indicates mitochondrial production of reactive oxygen species (ROS as one of the causes of replicative senescence. By sorting early senescent (SES cells from young proliferating fibroblast cultures, we show that SES cells have higher ROS levels, dysfunctional mitochondria, shorter telomeres, and telomeric gamma-H2A.X foci. We propose that mitochondrial ROS is a major determinant of telomere-dependent senescence at the single-cell level that is responsible for cell-to-cell variation in replicative lifespan.

Stochastic Models of Polymer Systems

Science.gov (United States)

2016-01-01

Distribution Unlimited Final Report: Stochastic Models of Polymer Systems The views, opinions and/or findings contained in this report are those of the...ADDRESS. Princeton University PO Box 0036 87 Prospect Avenue - 2nd floor Princeton, NJ 08544 -2020 14-Mar-2014 ABSTRACT Number of Papers published in...peer-reviewed journals: Number of Papers published in non peer-reviewed journals: Final Report: Stochastic Models of Polymer Systems Report Title

CAM Stochastic Volatility Model for Option Pricing

Directory of Open Access Journals (Sweden)

Wanwan Huang

2016-01-01

Full Text Available The coupled additive and multiplicative (CAM noises model is a stochastic volatility model for derivative pricing. Unlike the other stochastic volatility models in the literature, the CAM model uses two Brownian motions, one multiplicative and one additive, to model the volatility process. We provide empirical evidence that suggests a nontrivial relationship between the kurtosis and skewness of asset prices and that the CAM model is able to capture this relationship, whereas the traditional stochastic volatility models cannot. We introduce a control variate method and Monte Carlo estimators for some of the sensitivities (Greeks of the model. We also derive an approximation for the characteristic function of the model.

Modeling Stochastic Route Choice Behaviors with Equivalent Impedance

Directory of Open Access Journals (Sweden)

Jun Li

2015-01-01

Full Text Available A Logit-based route choice model is proposed to address the overlapping and scaling problems in the traditional multinomial Logit model. The nonoverlapping links are defined as a subnetwork, and its equivalent impedance is explicitly calculated in order to simply network analyzing. The overlapping links are repeatedly merged into subnetworks with Logit-based equivalent travel costs. The choice set at each intersection comprises only the virtual equivalent route without overlapping. In order to capture heterogeneity in perception errors of different sizes of networks, different scale parameters are assigned to subnetworks and they are linked to the topological relationships to avoid estimation burden. The proposed model provides an alternative method to model the stochastic route choice behaviors without the overlapping and scaling problems, and it still maintains the simple and closed-form expression from the MNL model. A link-based loading algorithm based on Dial’s algorithm is proposed to obviate route enumeration and it is suitable to be applied on large-scale networks. Finally a comparison between the proposed model and other route choice models is given by numerical examples.

Dynamic heterogeneity in life histories

DEFF Research Database (Denmark)

Tuljapurkar, Shripad; Steiner, Uli; Orzack, Steven Hecht

2009-01-01

or no fixed heterogeneity influences this trait. We propose that dynamic heterogeneity provides a 'neutral' model for assessing the possible role of unobserved 'quality' differences between individuals. We discuss fitness for dynamic life histories, and the implications of dynamic heterogeneity...... generate dynamic heterogeneity: life-history differences produced by stochastic stratum dynamics. We characterize dynamic heterogeneity in a range of species across taxa by properties of the Markov chain: the entropy, which describes the extent of heterogeneity, and the subdominant eigenvalue, which...... distributions of lifetime reproductive success. Dynamic heterogeneity contrasts with fixed heterogeneity: unobserved differences that generate variation between life histories. We show by an example that observed distributions of lifetime reproductive success are often consistent with the claim that little...

Modeling stochasticity and robustness in gene regulatory networks.

Science.gov (United States)

Garg, Abhishek; Mohanram, Kartik; Di Cara, Alessandro; De Micheli, Giovanni; Xenarios, Ioannis

2009-06-15

Understanding gene regulation in biological processes and modeling the robustness of underlying regulatory networks is an important problem that is currently being addressed by computational systems biologists. Lately, there has been a renewed interest in Boolean modeling techniques for gene regulatory networks (GRNs). However, due to their deterministic nature, it is often difficult to identify whether these modeling approaches are robust to the addition of stochastic noise that is widespread in gene regulatory processes. Stochasticity in Boolean models of GRNs has been addressed relatively sparingly in the past, mainly by flipping the expression of genes between different expression levels with a predefined probability. This stochasticity in nodes (SIN) model leads to over representation of noise in GRNs and hence non-correspondence with biological observations. In this article, we introduce the stochasticity in functions (SIF) model for simulating stochasticity in Boolean models of GRNs. By providing biological motivation behind the use of the SIF model and applying it to the T-helper and T-cell activation networks, we show that the SIF model provides more biologically robust results than the existing SIN model of stochasticity in GRNs. Algorithms are made available under our Boolean modeling toolbox, GenYsis. The software binaries can be downloaded from http://si2.epfl.ch/ approximately garg/genysis.html.

Modelling Stochastic Route Choice Behaviours with a Closed-Form Mixed Logit Model

Directory of Open Access Journals (Sweden)

Xinjun Lai

2015-01-01

Full Text Available A closed-form mixed Logit approach is proposed to model the stochastic route choice behaviours. It combines both the advantages of Probit and Logit to provide a flexible form in alternatives correlation and a tractable form in expression; besides, the heterogeneity in alternative variance can also be addressed. Paths are compared by pairs where the superiority of the binary Probit can be fully used. The Probit-based aggregation is also used for a nested Logit structure. Case studies on both numerical and empirical examples demonstrate that the new method is valid and practical. This paper thus provides an operational solution to incorporate the normal distribution in route choice with an analytical expression.

Particle-based modeling of heterogeneous chemical kinetics including mass transfer.

Science.gov (United States)

Sengar, A; Kuipers, J A M; van Santen, Rutger A; Padding, J T

2017-08-01

Connecting the macroscopic world of continuous fields to the microscopic world of discrete molecular events is important for understanding several phenomena occurring at physical boundaries of systems. An important example is heterogeneous catalysis, where reactions take place at active surfaces, but the effective reaction rates are determined by transport limitations in the bulk fluid and reaction limitations on the catalyst surface. In this work we study the macro-micro connection in a model heterogeneous catalytic reactor by means of stochastic rotation dynamics. The model is able to resolve the convective and diffusive interplay between participating species, while including adsorption, desorption, and reaction processes on the catalytic surface. Here we apply the simulation methodology to a simple straight microchannel with a catalytic strip. Dimensionless Damkohler numbers are used to comment on the spatial concentration profiles of reactants and products near the catalyst strip and in the bulk. We end the discussion with an outlook on more complicated geometries and increasingly complex reactions.

Particle-based modeling of heterogeneous chemical kinetics including mass transfer

Science.gov (United States)

Sengar, A.; Kuipers, J. A. M.; van Santen, Rutger A.; Padding, J. T.

2017-08-01

Connecting the macroscopic world of continuous fields to the microscopic world of discrete molecular events is important for understanding several phenomena occurring at physical boundaries of systems. An important example is heterogeneous catalysis, where reactions take place at active surfaces, but the effective reaction rates are determined by transport limitations in the bulk fluid and reaction limitations on the catalyst surface. In this work we study the macro-micro connection in a model heterogeneous catalytic reactor by means of stochastic rotation dynamics. The model is able to resolve the convective and diffusive interplay between participating species, while including adsorption, desorption, and reaction processes on the catalytic surface. Here we apply the simulation methodology to a simple straight microchannel with a catalytic strip. Dimensionless Damkohler numbers are used to comment on the spatial concentration profiles of reactants and products near the catalyst strip and in the bulk. We end the discussion with an outlook on more complicated geometries and increasingly complex reactions.

A stochastic model for density-dependent microwave Snow- and Graupel scattering coefficients of the NOAA JCSDA community radiative transfer model

Science.gov (United States)

Stegmann, Patrick G.; Tang, Guanglin; Yang, Ping; Johnson, Benjamin T.

2018-05-01

A structural model is developed for the single-scattering properties of snow and graupel particles with a strongly heterogeneous morphology and an arbitrary variable mass density. This effort is aimed to provide a mechanism to consider particle mass density variation in the microwave scattering coefficients implemented in the Community Radiative Transfer Model (CRTM). The stochastic model applies a bicontinuous random medium algorithm to a simple base shape and uses the Finite-Difference-Time-Domain (FDTD) method to compute the single-scattering properties of the resulting complex morphology.

Stochastic models in reliability and maintenance

CERN Document Server

2002-01-01

Our daily lives can be maintained by the high-technology systems. Computer systems are typical examples of such systems. We can enjoy our modern lives by using many computer systems. Much more importantly, we have to maintain such systems without failure, but cannot predict when such systems will fail and how to fix such systems without delay. A stochastic process is a set of outcomes of a random experiment indexed by time, and is one of the key tools needed to analyze the future behavior quantitatively. Reliability and maintainability technologies are of great interest and importance to the maintenance of such systems. Many mathematical models have been and will be proposed to describe reliability and maintainability systems by using the stochastic processes. The theme of this book is "Stochastic Models in Reliability and Main­ tainability. " This book consists of 12 chapters on the theme above from the different viewpoints of stochastic modeling. Chapter 1 is devoted to "Renewal Processes," under which cla...

The multivariate supOU stochastic volatility model

DEFF Research Database (Denmark)

Barndorff-Nielsen, Ole; Stelzer, Robert

Using positive semidefinite supOU (superposition of Ornstein-Uhlenbeck type) processes to describe the volatility, we introduce a multivariate stochastic volatility model for financial data which is capable of modelling long range dependence effects. The finiteness of moments and the second order...... structure of the volatility, the log returns, as well as their "squares" are discussed in detail. Moreover, we give several examples in which long memory effects occur and study how the model as well as the simple Ornstein-Uhlenbeck type stochastic volatility model behave under linear transformations....... In particular, the models are shown to be preserved under invertible linear transformations. Finally, we discuss how (sup)OU stochastic volatility models can be combined with a factor modelling approach....

A Fractionally Integrated Wishart Stochastic Volatility Model

NARCIS (Netherlands)

M. Asai (Manabu); M.J. McAleer (Michael)

2013-01-01

textabstractThere has recently been growing interest in modeling and estimating alternative continuous time multivariate stochastic volatility models. We propose a continuous time fractionally integrated Wishart stochastic volatility (FIWSV) process. We derive the conditional Laplace transform of

A stochastic SIS epidemic model with vaccination

Science.gov (United States)

Cao, Boqiang; Shan, Meijing; Zhang, Qimin; Wang, Weiming

2017-11-01

In this paper, we investigate the basic features of an SIS type infectious disease model with varying population size and vaccinations in presence of environment noise. By applying the Markov semigroup theory, we propose a stochastic reproduction number R0s which can be seen as a threshold parameter to utilize in identifying the stochastic extinction and persistence: If R0s disease-free absorbing set for the stochastic epidemic model, which implies that disease dies out with probability one; while if R0s > 1, under some mild extra conditions, the SDE model has an endemic stationary distribution which results in the stochastic persistence of the infectious disease. The most interesting finding is that large environmental noise can suppress the outbreak of the disease.

Population stochastic modelling (PSM)--an R package for mixed-effects models based on stochastic differential equations.

Science.gov (United States)

Klim, Søren; Mortensen, Stig Bousgaard; Kristensen, Niels Rode; Overgaard, Rune Viig; Madsen, Henrik

2009-06-01

The extension from ordinary to stochastic differential equations (SDEs) in pharmacokinetic and pharmacodynamic (PK/PD) modelling is an emerging field and has been motivated in a number of articles [N.R. Kristensen, H. Madsen, S.H. Ingwersen, Using stochastic differential equations for PK/PD model development, J. Pharmacokinet. Pharmacodyn. 32 (February(1)) (2005) 109-141; C.W. Tornøe, R.V. Overgaard, H. Agersø, H.A. Nielsen, H. Madsen, E.N. Jonsson, Stochastic differential equations in NONMEM: implementation, application, and comparison with ordinary differential equations, Pharm. Res. 22 (August(8)) (2005) 1247-1258; R.V. Overgaard, N. Jonsson, C.W. Tornøe, H. Madsen, Non-linear mixed-effects models with stochastic differential equations: implementation of an estimation algorithm, J. Pharmacokinet. Pharmacodyn. 32 (February(1)) (2005) 85-107; U. Picchini, S. Ditlevsen, A. De Gaetano, Maximum likelihood estimation of a time-inhomogeneous stochastic differential model of glucose dynamics, Math. Med. Biol. 25 (June(2)) (2008) 141-155]. PK/PD models are traditionally based ordinary differential equations (ODEs) with an observation link that incorporates noise. This state-space formulation only allows for observation noise and not for system noise. Extending to SDEs allows for a Wiener noise component in the system equations. This additional noise component enables handling of autocorrelated residuals originating from natural variation or systematic model error. Autocorrelated residuals are often partly ignored in PK/PD modelling although violating the hypothesis for many standard statistical tests. This article presents a package for the statistical program R that is able to handle SDEs in a mixed-effects setting. The estimation method implemented is the FOCE(1) approximation to the population likelihood which is generated from the individual likelihoods that are approximated using the Extended Kalman Filter's one-step predictions.

Conceptualizing a tool to optimize therapy based on dynamic heterogeneity

International Nuclear Information System (INIS)

Liao, David; Estévez-Salmerón, Luis; Tlsty, Thea D

2012-01-01

Complex biological systems often display a randomness paralleled in processes studied in fundamental physics. This simple stochasticity emerges owing to the complexity of the system and underlies a fundamental aspect of biology called phenotypic stochasticity. Ongoing stochastic fluctuations in phenotype at the single-unit level can contribute to two emergent population phenotypes. Phenotypic stochasticity not only generates heterogeneity within a cell population, but also allows reversible transitions back and forth between multiple states. This phenotypic interconversion tends to restore a population to a previous composition after that population has been depleted of specific members. We call this tendency homeostatic heterogeneity. These concepts of dynamic heterogeneity can be applied to populations composed of molecules, cells, individuals, etc. Here we discuss the concept that phenotypic stochasticity both underlies the generation of heterogeneity within a cell population and can be used to control population composition, contributing, in particular, to both the ongoing emergence of drug resistance and an opportunity for depleting drug-resistant cells. Using notions of both ‘large’ and ‘small’ numbers of biomolecular components, we rationalize our use of Markov processes to model the generation and eradication of drug-resistant cells. Using these insights, we have developed a graphical tool, called a metronomogram, that we propose will allow us to optimize dosing frequencies and total course durations for clinical benefit. (paper)

Consistent Stochastic Modelling of Meteocean Design Parameters

DEFF Research Database (Denmark)

Sørensen, John Dalsgaard; Sterndorff, M. J.

2000-01-01

Consistent stochastic models of metocean design parameters and their directional dependencies are essential for reliability assessment of offshore structures. In this paper a stochastic model for the annual maximum values of the significant wave height, and the associated wind velocity, current...

Stochastic linear programming models, theory, and computation

CERN Document Server

Kall, Peter

2011-01-01

This new edition of Stochastic Linear Programming: Models, Theory and Computation has been brought completely up to date, either dealing with or at least referring to new material on models and methods, including DEA with stochastic outputs modeled via constraints on special risk functions (generalizing chance constraints, ICC’s and CVaR constraints), material on Sharpe-ratio, and Asset Liability Management models involving CVaR in a multi-stage setup. To facilitate use as a text, exercises are included throughout the book, and web access is provided to a student version of the authors’ SLP-IOR software. Additionally, the authors have updated the Guide to Available Software, and they have included newer algorithms and modeling systems for SLP. The book is thus suitable as a text for advanced courses in stochastic optimization, and as a reference to the field. From Reviews of the First Edition: "The book presents a comprehensive study of stochastic linear optimization problems and their applications. … T...

Stochastic Wake Modelling Based on POD Analysis

Directory of Open Access Journals (Sweden)

David Bastine

2018-03-01

Full Text Available In this work, large eddy simulation data is analysed to investigate a new stochastic modeling approach for the wake of a wind turbine. The data is generated by the large eddy simulation (LES model PALM combined with an actuator disk with rotation representing the turbine. After applying a proper orthogonal decomposition (POD, three different stochastic models for the weighting coefficients of the POD modes are deduced resulting in three different wake models. Their performance is investigated mainly on the basis of aeroelastic simulations of a wind turbine in the wake. Three different load cases and their statistical characteristics are compared for the original LES, truncated PODs and the stochastic wake models including different numbers of POD modes. It is shown that approximately six POD modes are enough to capture the load dynamics on large temporal scales. Modeling the weighting coefficients as independent stochastic processes leads to similar load characteristics as in the case of the truncated POD. To complete this simplified wake description, we show evidence that the small-scale dynamics can be captured by adding to our model a homogeneous turbulent field. In this way, we present a procedure to derive stochastic wake models from costly computational fluid dynamics (CFD calculations or elaborated experimental investigations. These numerically efficient models provide the added value of possible long-term studies. Depending on the aspects of interest, different minimalized models may be obtained.

An improved multi-value cellular automata model for heterogeneous bicycle traffic flow

Energy Technology Data Exchange (ETDEWEB)

Jin, Sheng [College of Civil Engineering and Architecture, Zhejiang University, Hangzhou, 310058 China (China); Qu, Xiaobo [Griffith School of Engineering, Griffith University, Gold Coast, 4222 Australia (Australia); Xu, Cheng [Department of Transportation Management Engineering, Zhejiang Police College, Hangzhou, 310053 China (China); College of Transportation, Jilin University, Changchun, 130022 China (China); Ma, Dongfang, E-mail: mdf2004@zju.edu.cn [Ocean College, Zhejiang University, Hangzhou, 310058 China (China); Wang, Dianhai [College of Civil Engineering and Architecture, Zhejiang University, Hangzhou, 310058 China (China)

2015-10-16

This letter develops an improved multi-value cellular automata model for heterogeneous bicycle traffic flow taking the higher maximum speed of electric bicycles into consideration. The update rules of both regular and electric bicycles are improved, with maximum speeds of two and three cells per second respectively. Numerical simulation results for deterministic and stochastic cases are obtained. The fundamental diagrams and multiple states effects under different model parameters are analyzed and discussed. Field observations were made to calibrate the slowdown probabilities. The results imply that the improved extended Burgers cellular automata (IEBCA) model is more consistent with the field observations than previous models and greatly enhances the realism of the bicycle traffic model. - Highlights: • We proposed an improved multi-value CA model with higher maximum speed. • Update rules are introduced for heterogeneous bicycle traffic with maximum speed 2 and 3 cells/s. • Simulation results of the proposed model are consistent with field bicycle data. • Slowdown probabilities of both regular and electric bicycles are calibrated.

An improved multi-value cellular automata model for heterogeneous bicycle traffic flow

International Nuclear Information System (INIS)

Jin, Sheng; Qu, Xiaobo; Xu, Cheng; Ma, Dongfang; Wang, Dianhai

2015-01-01

This letter develops an improved multi-value cellular automata model for heterogeneous bicycle traffic flow taking the higher maximum speed of electric bicycles into consideration. The update rules of both regular and electric bicycles are improved, with maximum speeds of two and three cells per second respectively. Numerical simulation results for deterministic and stochastic cases are obtained. The fundamental diagrams and multiple states effects under different model parameters are analyzed and discussed. Field observations were made to calibrate the slowdown probabilities. The results imply that the improved extended Burgers cellular automata (IEBCA) model is more consistent with the field observations than previous models and greatly enhances the realism of the bicycle traffic model. - Highlights: • We proposed an improved multi-value CA model with higher maximum speed. • Update rules are introduced for heterogeneous bicycle traffic with maximum speed 2 and 3 cells/s. • Simulation results of the proposed model are consistent with field bicycle data. • Slowdown probabilities of both regular and electric bicycles are calibrated

Stochastic stability and bifurcation in a macroeconomic model

International Nuclear Information System (INIS)

Li Wei; Xu Wei; Zhao Junfeng; Jin Yanfei

2007-01-01

On the basis of the work of Goodwin and Puu, a new business cycle model subject to a stochastically parametric excitation is derived in this paper. At first, we reduce the model to a one-dimensional diffusion process by applying the stochastic averaging method of quasi-nonintegrable Hamiltonian system. Secondly, we utilize the methods of Lyapunov exponent and boundary classification associated with diffusion process respectively to analyze the stochastic stability of the trivial solution of system. The numerical results obtained illustrate that the trivial solution of system must be globally stable if it is locally stable in the state space. Thirdly, we explore the stochastic Hopf bifurcation of the business cycle model according to the qualitative changes in stationary probability density of system response. It is concluded that the stochastic Hopf bifurcation occurs at two critical parametric values. Finally, some explanations are given in a simply way on the potential applications of stochastic stability and bifurcation analysis

Population stochastic modelling (PSM)-An R package for mixed-effects models based on stochastic differential equations

DEFF Research Database (Denmark)

Klim, Søren; Mortensen, Stig Bousgaard; Kristensen, Niels Rode

2009-01-01

are often partly ignored in PK/PD modelling although violating the hypothesis for many standard statistical tests. This article presents a package for the statistical program R that is able to handle SDEs in a mixed-effects setting. The estimation method implemented is the FOCE1 approximation......The extension from ordinary to stochastic differential equations (SDEs) in pharmacokinetic and pharmacodynamic (PK/PD) modelling is an emerging field and has been motivated in a number of articles [N.R. Kristensen, H. Madsen, S.H. Ingwersen, Using stochastic differential equations for PK/PD model...... development, J. Pharmacokinet. Pharmacodyn. 32 (February(l)) (2005) 109-141; C.W. Tornoe, R.V Overgaard, H. Agerso, H.A. Nielsen, H. Madsen, E.N. Jonsson, Stochastic differential equations in NONMEM: implementation, application, and comparison with ordinary differential equations, Pharm. Res. 22 (August(8...

Stochastic volatility models and Kelvin waves

Science.gov (United States)

Lipton, Alex; Sepp, Artur

2008-08-01

We use stochastic volatility models to describe the evolution of an asset price, its instantaneous volatility and its realized volatility. In particular, we concentrate on the Stein and Stein model (SSM) (1991) for the stochastic asset volatility and the Heston model (HM) (1993) for the stochastic asset variance. By construction, the volatility is not sign definite in SSM and is non-negative in HM. It is well known that both models produce closed-form expressions for the prices of vanilla option via the Lewis-Lipton formula. However, the numerical pricing of exotic options by means of the finite difference and Monte Carlo methods is much more complex for HM than for SSM. Until now, this complexity was considered to be an acceptable price to pay for ensuring that the asset volatility is non-negative. We argue that having negative stochastic volatility is a psychological rather than financial or mathematical problem, and advocate using SSM rather than HM in most applications. We extend SSM by adding volatility jumps and obtain a closed-form expression for the density of the asset price and its realized volatility. We also show that the current method of choice for solving pricing problems with stochastic volatility (via the affine ansatz for the Fourier-transformed density function) can be traced back to the Kelvin method designed in the 19th century for studying wave motion problems arising in fluid dynamics.

Stochastic volatility models and Kelvin waves

International Nuclear Information System (INIS)

Lipton, Alex; Sepp, Artur

2008-01-01

We use stochastic volatility models to describe the evolution of an asset price, its instantaneous volatility and its realized volatility. In particular, we concentrate on the Stein and Stein model (SSM) (1991) for the stochastic asset volatility and the Heston model (HM) (1993) for the stochastic asset variance. By construction, the volatility is not sign definite in SSM and is non-negative in HM. It is well known that both models produce closed-form expressions for the prices of vanilla option via the Lewis-Lipton formula. However, the numerical pricing of exotic options by means of the finite difference and Monte Carlo methods is much more complex for HM than for SSM. Until now, this complexity was considered to be an acceptable price to pay for ensuring that the asset volatility is non-negative. We argue that having negative stochastic volatility is a psychological rather than financial or mathematical problem, and advocate using SSM rather than HM in most applications. We extend SSM by adding volatility jumps and obtain a closed-form expression for the density of the asset price and its realized volatility. We also show that the current method of choice for solving pricing problems with stochastic volatility (via the affine ansatz for the Fourier-transformed density function) can be traced back to the Kelvin method designed in the 19th century for studying wave motion problems arising in fluid dynamics

Stochastic volatility models and Kelvin waves

Energy Technology Data Exchange (ETDEWEB)

Lipton, Alex [Merrill Lynch, Mlfc Main, 2 King Edward Street, London EC1A 1HQ (United Kingdom); Sepp, Artur [Merrill Lynch, 4 World Financial Center, New York, NY 10080 (United States)], E-mail: Alex_Lipton@ml.com, E-mail: Artur_Sepp@ml.com

2008-08-29

We use stochastic volatility models to describe the evolution of an asset price, its instantaneous volatility and its realized volatility. In particular, we concentrate on the Stein and Stein model (SSM) (1991) for the stochastic asset volatility and the Heston model (HM) (1993) for the stochastic asset variance. By construction, the volatility is not sign definite in SSM and is non-negative in HM. It is well known that both models produce closed-form expressions for the prices of vanilla option via the Lewis-Lipton formula. However, the numerical pricing of exotic options by means of the finite difference and Monte Carlo methods is much more complex for HM than for SSM. Until now, this complexity was considered to be an acceptable price to pay for ensuring that the asset volatility is non-negative. We argue that having negative stochastic volatility is a psychological rather than financial or mathematical problem, and advocate using SSM rather than HM in most applications. We extend SSM by adding volatility jumps and obtain a closed-form expression for the density of the asset price and its realized volatility. We also show that the current method of choice for solving pricing problems with stochastic volatility (via the affine ansatz for the Fourier-transformed density function) can be traced back to the Kelvin method designed in the 19th century for studying wave motion problems arising in fluid dynamics.

Stochastic Modelling Of The Repairable System

Directory of Open Access Journals (Sweden)

Andrzejczak Karol

2015-11-01

Full Text Available All reliability models consisting of random time factors form stochastic processes. In this paper we recall the definitions of the most common point processes which are used for modelling of repairable systems. Particularly this paper presents stochastic processes as examples of reliability systems for the support of the maintenance related decisions. We consider the simplest one-unit system with a negligible repair or replacement time, i.e., the unit is operating and is repaired or replaced at failure, where the time required for repair and replacement is negligible. When the repair or replacement is completed, the unit becomes as good as new and resumes operation. The stochastic modelling of recoverable systems constitutes an excellent method of supporting maintenance related decision-making processes and enables their more rational use.

Hybrid approaches for multiple-species stochastic reaction–diffusion models

International Nuclear Information System (INIS)

Spill, Fabian; Guerrero, Pilar; Alarcon, Tomas; Maini, Philip K.; Byrne, Helen

2015-01-01

Reaction–diffusion models are used to describe systems in fields as diverse as physics, chemistry, ecology and biology. The fundamental quantities in such models are individual entities such as atoms and molecules, bacteria, cells or animals, which move and/or react in a stochastic manner. If the number of entities is large, accounting for each individual is inefficient, and often partial differential equation (PDE) models are used in which the stochastic behaviour of individuals is replaced by a description of the averaged, or mean behaviour of the system. In some situations the number of individuals is large in certain regions and small in others. In such cases, a stochastic model may be inefficient in one region, and a PDE model inaccurate in another. To overcome this problem, we develop a scheme which couples a stochastic reaction–diffusion system in one part of the domain with its mean field analogue, i.e. a discretised PDE model, in the other part of the domain. The interface in between the two domains occupies exactly one lattice site and is chosen such that the mean field description is still accurate there. In this way errors due to the flux between the domains are small. Our scheme can account for multiple dynamic interfaces separating multiple stochastic and deterministic domains, and the coupling between the domains conserves the total number of particles. The method preserves stochastic features such as extinction not observable in the mean field description, and is significantly faster to simulate on a computer than the pure stochastic model. - Highlights: • A novel hybrid stochastic/deterministic reaction–diffusion simulation method is given. • Can massively speed up stochastic simulations while preserving stochastic effects. • Can handle multiple reacting species. • Can handle moving boundaries

Hybrid approaches for multiple-species stochastic reaction–diffusion models

Energy Technology Data Exchange (ETDEWEB)

Spill, Fabian, E-mail: fspill@bu.edu [Department of Biomedical Engineering, Boston University, 44 Cummington Street, Boston, MA 02215 (United States); Department of Mechanical Engineering, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139 (United States); Guerrero, Pilar [Department of Mathematics, University College London, Gower Street, London WC1E 6BT (United Kingdom); Alarcon, Tomas [Centre de Recerca Matematica, Campus de Bellaterra, Edifici C, 08193 Bellaterra (Barcelona) (Spain); Departament de Matemàtiques, Universitat Atonòma de Barcelona, 08193 Bellaterra (Barcelona) (Spain); Maini, Philip K. [Wolfson Centre for Mathematical Biology, Mathematical Institute, University of Oxford, Oxford OX2 6GG (United Kingdom); Byrne, Helen [Wolfson Centre for Mathematical Biology, Mathematical Institute, University of Oxford, Oxford OX2 6GG (United Kingdom); Computational Biology Group, Department of Computer Science, University of Oxford, Oxford OX1 3QD (United Kingdom)

2015-10-15

Reaction–diffusion models are used to describe systems in fields as diverse as physics, chemistry, ecology and biology. The fundamental quantities in such models are individual entities such as atoms and molecules, bacteria, cells or animals, which move and/or react in a stochastic manner. If the number of entities is large, accounting for each individual is inefficient, and often partial differential equation (PDE) models are used in which the stochastic behaviour of individuals is replaced by a description of the averaged, or mean behaviour of the system. In some situations the number of individuals is large in certain regions and small in others. In such cases, a stochastic model may be inefficient in one region, and a PDE model inaccurate in another. To overcome this problem, we develop a scheme which couples a stochastic reaction–diffusion system in one part of the domain with its mean field analogue, i.e. a discretised PDE model, in the other part of the domain. The interface in between the two domains occupies exactly one lattice site and is chosen such that the mean field description is still accurate there. In this way errors due to the flux between the domains are small. Our scheme can account for multiple dynamic interfaces separating multiple stochastic and deterministic domains, and the coupling between the domains conserves the total number of particles. The method preserves stochastic features such as extinction not observable in the mean field description, and is significantly faster to simulate on a computer than the pure stochastic model. - Highlights: • A novel hybrid stochastic/deterministic reaction–diffusion simulation method is given. • Can massively speed up stochastic simulations while preserving stochastic effects. • Can handle multiple reacting species. • Can handle moving boundaries.

Dynamic heterogeneity and life history variability in the kittiwake

DEFF Research Database (Denmark)

Steiner, Uli; Tuljapurkar, Shripad; Orzack, Steven Hecht

2010-01-01

1. Understanding the evolution of life histories requires an assessment of the process that generates variation in life histories. Within-population heterogeneity of life histories can be dynamically generated by stochastic variation of reproduction and survival or be generated by individual...... differences that are fixed at birth. 2. We show for the kittiwake that dynamic heterogeneity is a sufficient explanation of observed variation of life histories. 3. The total heterogeneity in life histories has a small contribution from reproductive stage dynamics and a large contribution from survival...... differences. We quantify the diversity in life histories by metrics computed from the generating stochastic process. 4. We show how dynamic heterogeneity can be used as a null model and also how it can lead to positive associations between reproduction and survival across the life span. 5. We believe our...

Model selection for integrated pest management with stochasticity.

Science.gov (United States)

Akman, Olcay; Comar, Timothy D; Hrozencik, Daniel

2018-04-07

In Song and Xiang (2006), an integrated pest management model with periodically varying climatic conditions was introduced. In order to address a wider range of environmental effects, the authors here have embarked upon a series of studies resulting in a more flexible modeling approach. In Akman et al. (2013), the impact of randomly changing environmental conditions is examined by incorporating stochasticity into the birth pulse of the prey species. In Akman et al. (2014), the authors introduce a class of models via a mixture of two birth-pulse terms and determined conditions for the global and local asymptotic stability of the pest eradication solution. With this work, the authors unify the stochastic and mixture model components to create further flexibility in modeling the impacts of random environmental changes on an integrated pest management system. In particular, we first determine the conditions under which solutions of our deterministic mixture model are permanent. We then analyze the stochastic model to find the optimal value of the mixing parameter that minimizes the variance in the efficacy of the pesticide. Additionally, we perform a sensitivity analysis to show that the corresponding pesticide efficacy determined by this optimization technique is indeed robust. Through numerical simulations we show that permanence can be preserved in our stochastic model. Our study of the stochastic version of the model indicates that our results on the deterministic model provide informative conclusions about the behavior of the stochastic model. Copyright © 2017 Elsevier Ltd. All rights reserved.

Infinite-degree-corrected stochastic block model

DEFF Research Database (Denmark)

Herlau, Tue; Schmidt, Mikkel Nørgaard; Mørup, Morten

2014-01-01

In stochastic block models, which are among the most prominent statistical models for cluster analysis of complex networks, clusters are defined as groups of nodes with statistically similar link probabilities within and between groups. A recent extension by Karrer and Newman [Karrer and Newman...... corrected stochastic block model as a nonparametric Bayesian model, incorporating a parameter to control the amount of degree correction that can then be inferred from data. Additionally, our formulation yields principled ways of inferring the number of groups as well as predicting missing links...

Modeling in biopharmaceutics, pharmacokinetics and pharmacodynamics homogeneous and heterogeneous approaches

CERN Document Server

Macheras, Panos

2016-01-01

The state of the art in Biopharmaceutics, Pharmacokinetics, and Pharmacodynamics Modeling is presented in this new second edition book. It shows how advanced physical and mathematical methods can expand classical models in order to cover heterogeneous drug-biological processes and therapeutic effects in the body. The book is divided into four parts; the first deals with the fundamental principles of fractals, diffusion and nonlinear dynamics; the second with drug dissolution, release, and absorption; the third with epirical, compartmental, and stochastic pharmacokinetic models, with two new chapters, one on fractional pharmacokinetics and one on bioequivalence; and the fourth mainly with classical and nonclassical aspects of pharmacodynamics. The classical models that have relevance and application to these sciences are also considered throughout. This second edition has new information on reaction limited models of dissolution, non binary biopharmaceutic classification system, time varying models, and interf...

A Stochastic Model for Malaria Transmission Dynamics

Directory of Open Access Journals (Sweden)

Rachel Waema Mbogo

2018-01-01

Full Text Available Malaria is one of the three most dangerous infectious diseases worldwide (along with HIV/AIDS and tuberculosis. In this paper we compare the disease dynamics of the deterministic and stochastic models in order to determine the effect of randomness in malaria transmission dynamics. Relationships between the basic reproduction number for malaria transmission dynamics between humans and mosquitoes and the extinction thresholds of corresponding continuous-time Markov chain models are derived under certain assumptions. The stochastic model is formulated using the continuous-time discrete state Galton-Watson branching process (CTDSGWbp. The reproduction number of deterministic models is an essential quantity to predict whether an epidemic will spread or die out. Thresholds for disease extinction from stochastic models contribute crucial knowledge on disease control and elimination and mitigation of infectious diseases. Analytical and numerical results show some significant differences in model predictions between the stochastic and deterministic models. In particular, we find that malaria outbreak is more likely if the disease is introduced by infected mosquitoes as opposed to infected humans. These insights demonstrate the importance of a policy or intervention focusing on controlling the infected mosquito population if the control of malaria is to be realized.

Approximate reduction of linear population models governed by stochastic differential equations: application to multiregional models.

Science.gov (United States)

Sanz, Luis; Alonso, Juan Antonio

2017-12-01

In this work we develop approximate aggregation techniques in the context of slow-fast linear population models governed by stochastic differential equations and apply the results to the treatment of populations with spatial heterogeneity. Approximate aggregation techniques allow one to transform a complex system involving many coupled variables and in which there are processes with different time scales, by a simpler reduced model with a fewer number of 'global' variables, in such a way that the dynamics of the former can be approximated by that of the latter. In our model we contemplate a linear fast deterministic process together with a linear slow process in which the parameters are affected by additive noise, and give conditions for the solutions corresponding to positive initial conditions to remain positive for all times. By letting the fast process reach equilibrium we build a reduced system with a lesser number of variables, and provide results relating the asymptotic behaviour of the first- and second-order moments of the population vector for the original and the reduced system. The general technique is illustrated by analysing a multiregional stochastic system in which dispersal is deterministic and the rate growth of the populations in each patch is affected by additive noise.

Stochastic modeling and analysis of telecoms networks

CERN Document Server

Decreusefond, Laurent

2012-01-01

This book addresses the stochastic modeling of telecommunication networks, introducing the main mathematical tools for that purpose, such as Markov processes, real and spatial point processes and stochastic recursions, and presenting a wide list of results on stability, performances and comparison of systems.The authors propose a comprehensive mathematical construction of the foundations of stochastic network theory: Markov chains, continuous time Markov chains are extensively studied using an original martingale-based approach. A complete presentation of stochastic recursions from an

Stochastic line motion and stochastic flux conservation for nonideal hydromagnetic models

International Nuclear Information System (INIS)

Eyink, Gregory L.

2009-01-01

We prove that smooth solutions of nonideal (viscous and resistive) incompressible magnetohydrodynamic (MHD) equations satisfy a stochastic law of flux conservation. This property implies that the magnetic flux through a surface is equal to the average of the magnetic fluxes through an ensemble of surfaces advected backward in time by the plasma velocity perturbed with a random white noise. Our result is an analog of the well-known Alfven theorem of ideal MHD and is valid for any value of the magnetic Prandtl number. A second stochastic conservation law is shown to hold at unit Prandtl number, a random version of the generalized Kelvin theorem derived by Bekenstein and Oron for ideal MHD. These stochastic conservation laws are not only shown to be consequences of the nonideal MHD equations but are proved in fact to be equivalent to those equations. We derive similar results for two more refined hydromagnetic models, Hall MHD and the two-fluid plasma model, still assuming incompressible velocities and isotropic transport coefficients. Finally, we use these results to discuss briefly the infinite-Reynolds-number limit of hydromagnetic turbulence and to support the conjecture that flux conservation remains stochastic in that limit.

Stochastic cellular automata model of cell migration, proliferation and differentiation: validation with in vitro cultures of muscle satellite cells.

Science.gov (United States)

Garijo, N; Manzano, R; Osta, R; Perez, M A

2012-12-07

Cell migration and proliferation has been modelled in the literature as a process similar to diffusion. However, using diffusion models to simulate the proliferation and migration of cells tends to create a homogeneous distribution in the cell density that does not correlate to empirical observations. In fact, the mechanism of cell dispersal is not diffusion. Cells disperse by crawling or proliferation, or are transported in a moving fluid. The use of cellular automata, particle models or cell-based models can overcome this limitation. This paper presents a stochastic cellular automata model to simulate the proliferation, migration and differentiation of cells. These processes are considered as completely stochastic as well as discrete. The model developed was applied to predict the behaviour of in vitro cell cultures performed with adult muscle satellite cells. Moreover, non homogeneous distribution of cells has been observed inside the culture well and, using the above mentioned stochastic cellular automata model, we have been able to predict this heterogeneous cell distribution and compute accurate quantitative results. Differentiation was also incorporated into the computational simulation. The results predicted the myotube formation that typically occurs with adult muscle satellite cells. In conclusion, we have shown how a stochastic cellular automata model can be implemented and is capable of reproducing the in vitro behaviour of adult muscle satellite cells. Copyright © 2012 Elsevier Ltd. All rights reserved.

Modeling in biopharmaceutics, pharmacokinetics, and pharmacodynamics homogeneous and heterogeneous approaches

CERN Document Server

Macheras, Panos

2006-01-01

The state of the art in Biopharmaceutics, Pharmacokinetics, and Pharmacodynamics Modeling is presented in this book. It shows how advanced physical and mathematical methods can expand classical models in order to cover heterogeneous drug-biological processes and therapeutic effects in the body. The book is divided into four parts; the first deals with the fundamental principles of fractals, diffusion and nonlinear dynamics; the second with drug dissolution, release, and absorption; the third with empirical, compartmental, and stochastic pharmacokinetic models, and the fourth mainly with nonclassical aspects of pharmacodynamics. The classical models that have relevance and application to these sciences are also considered throughout. Many examples are used to illustrate the intrinsic complexity of drug administration related phenomena in the human, justifying the use of advanced modeling methods. This timely and useful book will appeal to graduate students and researchers in pharmacology, pharmaceutical scienc...

A stochastic model of enzyme kinetics

Science.gov (United States)

Stefanini, Marianne; Newman, Timothy; McKane, Alan

2003-10-01

Enzyme kinetics is generally modeled by deterministic rate equations, and in the simplest case leads to the well-known Michaelis-Menten equation. It is plausible that stochastic effects will play an important role at low enzyme concentrations. We have addressed this by constructing a simple stochastic model which can be exactly solved in the steady-state. Throughout a wide range of parameter values Michaelis-Menten dynamics is replaced by a new and simple theoretical result.

Stochastic Control - External Models

DEFF Research Database (Denmark)

Poulsen, Niels Kjølstad

2005-01-01

This note is devoted to control of stochastic systems described in discrete time. We are concerned with external descriptions or transfer function model, where we have a dynamic model for the input output relation only (i.e.. no direct internal information). The methods are based on LTI systems...

Accounting for heterogeneity in the measurement of hospital performance

OpenAIRE

Philippe K. Widmer; Peter Zweifel; Mehdi Farsi

2011-01-01

With prospective payment of hospitals becoming more common, measuring their performance is gaining in importance. However, the standard cost frontier model yields biased efficiency scores because it ignores technological heterogeneity between hospitals. In this paper, efficiency scores are derived from a random intercept and an extended random parameter frontier model, designed to overcome the problem of unobserved heterogeneity in stochastic frontier analysis. Using a sample of 100 Swiss ...

Pricing Liquidity Risk with Heterogeneous Investment Horizons

NARCIS (Netherlands)

Beber, Alessandro; Driessen, Joost; Neuberger, A.; Tuijp, P

We develop an asset pricing model with stochastic transaction costs and investors with heterogeneous horizons. Depending on their horizon, investors hold different sets of assets in equilibrium. This generates segmentation and spillover effects for expected returns, where the liquidity (risk)

Stochastic simulation of regional groundwater flow in Beishan area

International Nuclear Information System (INIS)

Dong Yanhui; Li Guomin

2010-01-01

Because of the hydrogeological complexity, traditional thinking of aquifer characteristics is not appropriate for groundwater system in Beishan area. Uncertainty analysis of groundwater models is needed to examine the hydrologic effects of spatial heterogeneity. In this study, fast Fourier transform spectral method (FFTS) was used to generate the random horizontal permeability parameters. Depth decay and vertical anisotropy of hydraulic conductivity were included to build random permeability models. Based on high-performance computers, hundreds of groundwater flow models were simulated. Through stochastic simulations, the effect of heterogeneity to groundwater flow pattern was analyzed. (authors)

Stochastic Watershed Models for Risk Based Decision Making

Science.gov (United States)

Vogel, R. M.

2017-12-01

Over half a century ago, the Harvard Water Program introduced the field of operational or synthetic hydrology providing stochastic streamflow models (SSMs), which could generate ensembles of synthetic streamflow traces useful for hydrologic risk management. The application of SSMs, based on streamflow observations alone, revolutionized water resources planning activities, yet has fallen out of favor due, in part, to their inability to account for the now nearly ubiquitous anthropogenic influences on streamflow. This commentary advances the modern equivalent of SSMs, termed `stochastic watershed models' (SWMs) useful as input to nearly all modern risk based water resource decision making approaches. SWMs are deterministic watershed models implemented using stochastic meteorological series, model parameters and model errors, to generate ensembles of streamflow traces that represent the variability in possible future streamflows. SWMs combine deterministic watershed models, which are ideally suited to accounting for anthropogenic influences, with recent developments in uncertainty analysis and principles of stochastic simulation

Dynamics of a Stochastic Intraguild Predation Model

Directory of Open Access Journals (Sweden)

Zejing Xing

2016-04-01

Full Text Available Intraguild predation (IGP is a widespread ecological phenomenon which occurs when one predator species attacks another predator species with which it competes for a shared prey species. The objective of this paper is to study the dynamical properties of a stochastic intraguild predation model. We analyze stochastic persistence and extinction of the stochastic IGP model containing five cases and establish the sufficient criteria for global asymptotic stability of the positive solutions. This study shows that it is possible for the coexistence of three species under the influence of environmental noise, and that the noise may have a positive effect for IGP species. A stationary distribution of the stochastic IGP model is established and it has the ergodic property, suggesting that the time average of population size with the development of time is equal to the stationary distribution in space. Finally, we show that our results may be extended to two well-known biological systems: food chains and exploitative competition.

Stochastic Modeling of Traffic Air Pollution

DEFF Research Database (Denmark)

Thoft-Christensen, Palle

2014-01-01

In this paper, modeling of traffic air pollution is discussed with special reference to infrastructures. A number of subjects related to health effects of air pollution and the different types of pollutants are briefly presented. A simple model for estimating the social cost of traffic related air...... and using simple Monte Carlo techniques to obtain a stochastic estimate of the costs of traffic air pollution for infrastructures....... pollution is derived. Several authors have published papers on this very complicated subject, but no stochastic modelling procedure have obtained general acceptance. The subject is discussed basis of a deterministic model. However, it is straightforward to modify this model to include uncertain parameters...

Stochastic modelling of turbulence

DEFF Research Database (Denmark)

Sørensen, Emil Hedevang Lohse

previously been shown to be closely connected to the energy dissipation. The incorporation of the small scale dynamics into the spatial model opens the door to a fully fledged stochastic model of turbulence. Concerning the interaction of wind and wind turbine, a new method is proposed to extract wind turbine...

Synthetic Sediments and Stochastic Groundwater Hydrology

Science.gov (United States)

Wilson, J. L.

2002-12-01

For over twenty years the groundwater community has pursued the somewhat elusive goal of describing the effects of aquifer heterogeneity on subsurface flow and chemical transport. While small perturbation stochastic moment methods have significantly advanced theoretical understanding, why is it that stochastic applications use instead simulations of flow and transport through multiple realizations of synthetic geology? Allan Gutjahr was a principle proponent of the Fast Fourier Transform method for the synthetic generation of aquifer properties and recently explored new, more geologically sound, synthetic methods based on multi-scale Markov random fields. Focusing on sedimentary aquifers, how has the state-of-the-art of synthetic generation changed and what new developments can be expected, for example, to deal with issues like conceptual model uncertainty, the differences between measurement and modeling scales, and subgrid scale variability? What will it take to get stochastic methods, whether based on moments, multiple realizations, or some other approach, into widespread application?

Pricing Liquidity Risk with Heterogeneous Investment Horizons

NARCIS (Netherlands)

Beber, A.; Driessen, J.; Tuijp, P.F.A.

2012-01-01

We develop a new asset pricing model with stochastic transaction costs and investors with heterogenous horizons. Short-term investors hold only liquid assets in equilibrium. This generates segmentation effects in the pricing of liquid versus illiquid assets. Specifically, the liquidity (risk) premia

Approximate models for broken clouds in stochastic radiative transfer theory

International Nuclear Information System (INIS)

Doicu, Adrian; Efremenko, Dmitry S.; Loyola, Diego; Trautmann, Thomas

2014-01-01

This paper presents approximate models in stochastic radiative transfer theory. The independent column approximation and its modified version with a solar source computed in a full three-dimensional atmosphere are formulated in a stochastic framework and for arbitrary cloud statistics. The nth-order stochastic models describing the independent column approximations are equivalent to the nth-order stochastic models for the original radiance fields in which the gradient vectors are neglected. Fast approximate models are further derived on the basis of zeroth-order stochastic models and the independent column approximation. The so-called “internal mixing” models assume a combination of the optical properties of the cloud and the clear sky, while the “external mixing” models assume a combination of the radiances corresponding to completely overcast and clear skies. A consistent treatment of internal and external mixing models is provided, and a new parameterization of the closure coefficient in the effective thickness approximation is given. An efficient computation of the closure coefficient for internal mixing models, using a previously derived vector stochastic model as a reference, is also presented. Equipped with appropriate look-up tables for the closure coefficient, these models can easily be integrated into operational trace gas retrieval systems that exploit absorption features in the near-IR solar spectrum. - Highlights: • Independent column approximation in a stochastic setting. • Fast internal and external mixing models for total and diffuse radiances. • Efficient optimization of internal mixing models to match reference models

Parameter estimation in stochastic rainfall-runoff models

DEFF Research Database (Denmark)

Jonsdottir, Harpa; Madsen, Henrik; Palsson, Olafur Petur

2006-01-01

A parameter estimation method for stochastic rainfall-runoff models is presented. The model considered in the paper is a conceptual stochastic model, formulated in continuous-discrete state space form. The model is small and a fully automatic optimization is, therefore, possible for estimating all...... the parameter values are optimal for simulation or prediction. The data originates from Iceland and the model is designed for Icelandic conditions, including a snow routine for mountainous areas. The model demands only two input data series, precipitation and temperature and one output data series...

Test models for improving filtering with model errors through stochastic parameter estimation

International Nuclear Information System (INIS)

Gershgorin, B.; Harlim, J.; Majda, A.J.

2010-01-01

The filtering skill for turbulent signals from nature is often limited by model errors created by utilizing an imperfect model for filtering. Updating the parameters in the imperfect model through stochastic parameter estimation is one way to increase filtering skill and model performance. Here a suite of stringent test models for filtering with stochastic parameter estimation is developed based on the Stochastic Parameterization Extended Kalman Filter (SPEKF). These new SPEKF-algorithms systematically correct both multiplicative and additive biases and involve exact formulas for propagating the mean and covariance including the parameters in the test model. A comprehensive study is presented of robust parameter regimes for increasing filtering skill through stochastic parameter estimation for turbulent signals as the observation time and observation noise are varied and even when the forcing is incorrectly specified. The results here provide useful guidelines for filtering turbulent signals in more complex systems with significant model errors.

Modelling and application of stochastic processes

CERN Document Server

1986-01-01

The subject of modelling and application of stochastic processes is too vast to be exhausted in a single volume. In this book, attention is focused on a small subset of this vast subject. The primary emphasis is on realization and approximation of stochastic systems. Recently there has been considerable interest in the stochastic realization problem, and hence, an attempt has been made here to collect in one place some of the more recent approaches and algorithms for solving the stochastic realiza­ tion problem. Various different approaches for realizing linear minimum-phase systems, linear nonminimum-phase systems, and bilinear systems are presented. These approaches range from time-domain methods to spectral-domain methods. An overview of the chapter contents briefly describes these approaches. Also, in most of these chapters special attention is given to the problem of developing numerically ef­ ficient algorithms for obtaining reduced-order (approximate) stochastic realizations. On the application side,...

12th Workshop on Stochastic Models, Statistics and Their Applications

CERN Document Server

Rafajłowicz, Ewaryst; Szajowski, Krzysztof

2015-01-01

This volume presents the latest advances and trends in stochastic models and related statistical procedures. Selected peer-reviewed contributions focus on statistical inference, quality control, change-point analysis and detection, empirical processes, time series analysis, survival analysis and reliability, statistics for stochastic processes, big data in technology and the sciences, statistical genetics, experiment design, and stochastic models in engineering. Stochastic models and related statistical procedures play an important part in furthering our understanding of the challenging problems currently arising in areas of application such as the natural sciences, information technology, engineering, image analysis, genetics, energy and finance, to name but a few. This collection arises from the 12th Workshop on Stochastic Models, Statistics and Their Applications, Wroclaw, Poland.

Improving sub-pixel imperviousness change prediction by ensembling heterogeneous non-linear regression models

Directory of Open Access Journals (Sweden)

Drzewiecki Wojciech

2016-12-01

Full Text Available In this work nine non-linear regression models were compared for sub-pixel impervious surface area mapping from Landsat images. The comparison was done in three study areas both for accuracy of imperviousness coverage evaluation in individual points in time and accuracy of imperviousness change assessment. The performance of individual machine learning algorithms (Cubist, Random Forest, stochastic gradient boosting of regression trees, k-nearest neighbors regression, random k-nearest neighbors regression, Multivariate Adaptive Regression Splines, averaged neural networks, and support vector machines with polynomial and radial kernels was also compared with the performance of heterogeneous model ensembles constructed from the best models trained using particular techniques.

Constraining Stochastic Parametrisation Schemes Using High-Resolution Model Simulations

Science.gov (United States)

Christensen, H. M.; Dawson, A.; Palmer, T.

2017-12-01

Stochastic parametrisations are used in weather and climate models as a physically motivated way to represent model error due to unresolved processes. Designing new stochastic schemes has been the target of much innovative research over the last decade. While a focus has been on developing physically motivated approaches, many successful stochastic parametrisation schemes are very simple, such as the European Centre for Medium-Range Weather Forecasts (ECMWF) multiplicative scheme `Stochastically Perturbed Parametrisation Tendencies' (SPPT). The SPPT scheme improves the skill of probabilistic weather and seasonal forecasts, and so is widely used. However, little work has focused on assessing the physical basis of the SPPT scheme. We address this matter by using high-resolution model simulations to explicitly measure the `error' in the parametrised tendency that SPPT seeks to represent. The high resolution simulations are first coarse-grained to the desired forecast model resolution before they are used to produce initial conditions and forcing data needed to drive the ECMWF Single Column Model (SCM). By comparing SCM forecast tendencies with the evolution of the high resolution model, we can measure the `error' in the forecast tendencies. In this way, we provide justification for the multiplicative nature of SPPT, and for the temporal and spatial scales of the stochastic perturbations. However, we also identify issues with the SPPT scheme. It is therefore hoped these measurements will improve both holistic and process based approaches to stochastic parametrisation. Figure caption: Instantaneous snapshot of the optimal SPPT stochastic perturbation, derived by comparing high-resolution simulations with a low resolution forecast model.

Hybrid approaches for multiple-species stochastic reaction-diffusion models

Science.gov (United States)

Spill, Fabian; Guerrero, Pilar; Alarcon, Tomas; Maini, Philip K.; Byrne, Helen

2015-10-01

Reaction-diffusion models are used to describe systems in fields as diverse as physics, chemistry, ecology and biology. The fundamental quantities in such models are individual entities such as atoms and molecules, bacteria, cells or animals, which move and/or react in a stochastic manner. If the number of entities is large, accounting for each individual is inefficient, and often partial differential equation (PDE) models are used in which the stochastic behaviour of individuals is replaced by a description of the averaged, or mean behaviour of the system. In some situations the number of individuals is large in certain regions and small in others. In such cases, a stochastic model may be inefficient in one region, and a PDE model inaccurate in another. To overcome this problem, we develop a scheme which couples a stochastic reaction-diffusion system in one part of the domain with its mean field analogue, i.e. a discretised PDE model, in the other part of the domain. The interface in between the two domains occupies exactly one lattice site and is chosen such that the mean field description is still accurate there. In this way errors due to the flux between the domains are small. Our scheme can account for multiple dynamic interfaces separating multiple stochastic and deterministic domains, and the coupling between the domains conserves the total number of particles. The method preserves stochastic features such as extinction not observable in the mean field description, and is significantly faster to simulate on a computer than the pure stochastic model.

Hybrid approaches for multiple-species stochastic reaction-diffusion models.

KAUST Repository

Spill, Fabian; Guerrero, Pilar; Alarcon, Tomas; Maini, Philip K; Byrne, Helen

2015-01-01

Reaction-diffusion models are used to describe systems in fields as diverse as physics, chemistry, ecology and biology. The fundamental quantities in such models are individual entities such as atoms and molecules, bacteria, cells or animals, which move and/or react in a stochastic manner. If the number of entities is large, accounting for each individual is inefficient, and often partial differential equation (PDE) models are used in which the stochastic behaviour of individuals is replaced by a description of the averaged, or mean behaviour of the system. In some situations the number of individuals is large in certain regions and small in others. In such cases, a stochastic model may be inefficient in one region, and a PDE model inaccurate in another. To overcome this problem, we develop a scheme which couples a stochastic reaction-diffusion system in one part of the domain with its mean field analogue, i.e. a discretised PDE model, in the other part of the domain. The interface in between the two domains occupies exactly one lattice site and is chosen such that the mean field description is still accurate there. In this way errors due to the flux between the domains are small. Our scheme can account for multiple dynamic interfaces separating multiple stochastic and deterministic domains, and the coupling between the domains conserves the total number of particles. The method preserves stochastic features such as extinction not observable in the mean field description, and is significantly faster to simulate on a computer than the pure stochastic model.

Hybrid approaches for multiple-species stochastic reaction-diffusion models.

KAUST Repository

Spill, Fabian

2015-10-01

Reaction-diffusion models are used to describe systems in fields as diverse as physics, chemistry, ecology and biology. The fundamental quantities in such models are individual entities such as atoms and molecules, bacteria, cells or animals, which move and/or react in a stochastic manner. If the number of entities is large, accounting for each individual is inefficient, and often partial differential equation (PDE) models are used in which the stochastic behaviour of individuals is replaced by a description of the averaged, or mean behaviour of the system. In some situations the number of individuals is large in certain regions and small in others. In such cases, a stochastic model may be inefficient in one region, and a PDE model inaccurate in another. To overcome this problem, we develop a scheme which couples a stochastic reaction-diffusion system in one part of the domain with its mean field analogue, i.e. a discretised PDE model, in the other part of the domain. The interface in between the two domains occupies exactly one lattice site and is chosen such that the mean field description is still accurate there. In this way errors due to the flux between the domains are small. Our scheme can account for multiple dynamic interfaces separating multiple stochastic and deterministic domains, and the coupling between the domains conserves the total number of particles. The method preserves stochastic features such as extinction not observable in the mean field description, and is significantly faster to simulate on a computer than the pure stochastic model.

Overview of medium heterogeneity and transport processes

International Nuclear Information System (INIS)

Tsang, Y.; Tsang, C.F.

1993-11-01

Medium heterogeneity can have significant impact on the behavior of solute transport. Tracer breakthrough curves from transport in a heterogeneous medium are distinctly different from that in a homogeneous porous medium. Usually the shape of the breakthrough curves are highly non-symmetrical with a fast rise at early times and very long tail at late times, and often, they consist of multiple peaks. Moreover, unlike transport in a homogeneous medium where the same transport parameters describe the entire medium, transport through heterogeneous media gives rise to breakthrough curves which have strong spatial dependence. These inherent characteristics of transport in heterogeneous medium present special challenge to the performance assessment of a potential high level nuclear waste repository with respect to the possible release of radio nuclides to the accessible environment. Since an inherently desirable site characteristic for a waste repository is that flow and transport should be slow, then transport measurements in site characterization efforts will necessarily be spatially small and temporally short compare to the scales which are of relevance to performance assessment predictions. In this paper we discuss the role of medium heterogeneity in site characterization and performance assessment. Our discussion will be based on a specific example of a 3D heterogeneous stochastic model of a site generally similar to, the Aespoe Island, the site of the Hard Rock Laboratory in Southern Sweden. For our study, alternative 3D stochastic fields of hydraulic conductivities conditioned on ''point'' measurements shall be generated. Results of stochastic flow and transport simulations would be used to address the issues of (1) the relationship of tracer breakthrough with the structure of heterogeneity, and (2) the inference from small scale testing results to large scale and long term predictions

A computationally efficient description of heterogeneous freezing: A simplified version of the Soccer ball model

Science.gov (United States)

Niedermeier, Dennis; Ervens, Barbara; Clauss, Tina; Voigtländer, Jens; Wex, Heike; Hartmann, Susan; Stratmann, Frank

2014-01-01

In a recent study, the Soccer ball model (SBM) was introduced for modeling and/or parameterizing heterogeneous ice nucleation processes. The model applies classical nucleation theory. It allows for a consistent description of both apparently singular and stochastic ice nucleation behavior, by distributing contact angles over the nucleation sites of a particle population assuming a Gaussian probability density function. The original SBM utilizes the Monte Carlo technique, which hampers its usage in atmospheric models, as fairly time-consuming calculations must be performed to obtain statistically significant results. Thus, we have developed a simplified and computationally more efficient version of the SBM. We successfully used the new SBM to parameterize experimental nucleation data of, e.g., bacterial ice nucleation. Both SBMs give identical results; however, the new model is computationally less expensive as confirmed by cloud parcel simulations. Therefore, it is a suitable tool for describing heterogeneous ice nucleation processes in atmospheric models.

A stochastic modeling of recurrent measles epidemic | Kassem ...

African Journals Online (AJOL)

A simple stochastic mathematical model is developed and investigated for the dynamics of measles epidemic. The model, which is a multi-dimensional diffusion process, includes susceptible individuals, latent (exposed), infected and removed individuals. Stochastic effects are assumed to arise in the process of infection of ...

A cavitation model based on Eulerian stochastic fields

Science.gov (United States)

Magagnato, F.; Dumond, J.

2013-12-01

Non-linear phenomena can often be described using probability density functions (pdf) and pdf transport models. Traditionally the simulation of pdf transport requires Monte-Carlo codes based on Lagrangian "particles" or prescribed pdf assumptions including binning techniques. Recently, in the field of combustion, a novel formulation called the stochastic-field method solving pdf transport based on Eulerian fields has been proposed which eliminates the necessity to mix Eulerian and Lagrangian techniques or prescribed pdf assumptions. In the present work, for the first time the stochastic-field method is applied to multi-phase flow and in particular to cavitating flow. To validate the proposed stochastic-field cavitation model, two applications are considered. Firstly, sheet cavitation is simulated in a Venturi-type nozzle. The second application is an innovative fluidic diode which exhibits coolant flashing. Agreement with experimental results is obtained for both applications with a fixed set of model constants. The stochastic-field cavitation model captures the wide range of pdf shapes present at different locations.

Stochastic quantization for the axial model

International Nuclear Information System (INIS)

Farina, C.; Montani, H.; Albuquerque, L.C.

1991-01-01

We use bosonization ideas to solve the axial model in the stochastic quantization framework. We obtain the fermion propagator of the theory decoupling directly the Langevin equation, instead of the Fokker-Planck equation. In the Appendix we calculate explicitly the anomalous divergence of the axial-vector current by using a regularization that does not break the Markovian character of the stochastic process

Stochastic model simulation using Kronecker product analysis and Zassenhaus formula approximation.

Science.gov (United States)

Caglar, Mehmet Umut; Pal, Ranadip

2013-01-01

Probabilistic Models are regularly applied in Genetic Regulatory Network modeling to capture the stochastic behavior observed in the generation of biological entities such as mRNA or proteins. Several approaches including Stochastic Master Equations and Probabilistic Boolean Networks have been proposed to model the stochastic behavior in genetic regulatory networks. It is generally accepted that Stochastic Master Equation is a fundamental model that can describe the system being investigated in fine detail, but the application of this model is computationally enormously expensive. On the other hand, Probabilistic Boolean Network captures only the coarse-scale stochastic properties of the system without modeling the detailed interactions. We propose a new approximation of the stochastic master equation model that is able to capture the finer details of the modeled system including bistabilities and oscillatory behavior, and yet has a significantly lower computational complexity. In this new method, we represent the system using tensors and derive an identity to exploit the sparse connectivity of regulatory targets for complexity reduction. The algorithm involves an approximation based on Zassenhaus formula to represent the exponential of a sum of matrices as product of matrices. We derive upper bounds on the expected error of the proposed model distribution as compared to the stochastic master equation model distribution. Simulation results of the application of the model to four different biological benchmark systems illustrate performance comparable to detailed stochastic master equation models but with considerably lower computational complexity. The results also demonstrate the reduced complexity of the new approach as compared to commonly used Stochastic Simulation Algorithm for equivalent accuracy.

in Heterogeneous Media

Directory of Open Access Journals (Sweden)

Saeed Balouchi

2013-01-01

Full Text Available Fractured reservoirs contain about 85 and 90 percent of oil and gas resources respectively in Iran. A comprehensive study and investigation of fractures as the main factor affecting fluid flow or perhaps barrier seems necessary for reservoir development studies. High degrees of heterogeneity and sparseness of data have incapacitated conventional deterministic methods in fracture network modeling. Recently, simulated annealing (SA has been applied to generate stochastic realizations of spatially correlated fracture networks by assuming that the elastic energy of fractures follows Boltzmann distribution. Although SA honors local variability, the objective function of geometrical fracture modeling is defined for homogeneous conditions. In this study, after the introduction of SA and the derivation of the energy function, a novel technique is presented to adjust the model with highly heterogeneous data for a fractured field from the southwest of Iran. To this end, the regular object-based model is combined with a grid-based technique to cover the heterogeneity of reservoir properties. The original SA algorithm is also modified by being constrained in different directions and weighting the energy function to make it appropriate for heterogeneous conditions. The simulation results of the presented approach are in good agreement with the observed field data.

Analysis of stochastic effects in Kaldor-type business cycle discrete model

Science.gov (United States)

Bashkirtseva, Irina; Ryashko, Lev; Sysolyatina, Anna

2016-07-01

We study nonlinear stochastic phenomena in the discrete Kaldor model of business cycles. A numerical parametric analysis of stochastically forced attractors (equilibria, closed invariant curves, discrete cycles) of this model is performed using the stochastic sensitivity functions technique. A spatial arrangement of random states in stochastic attractors is modeled by confidence domains. The phenomenon of noise-induced transitions ;chaos-order; is discussed.

Stochastic Spectral Descent for Discrete Graphical Models

International Nuclear Information System (INIS)

Carlson, David; Hsieh, Ya-Ping; Collins, Edo; Carin, Lawrence; Cevher, Volkan

2015-01-01

Interest in deep probabilistic graphical models has in-creased in recent years, due to their state-of-the-art performance on many machine learning applications. Such models are typically trained with the stochastic gradient method, which can take a significant number of iterations to converge. Since the computational cost of gradient estimation is prohibitive even for modestly sized models, training becomes slow and practically usable models are kept small. In this paper we propose a new, largely tuning-free algorithm to address this problem. Our approach derives novel majorization bounds based on the Schatten- norm. Intriguingly, the minimizers of these bounds can be interpreted as gradient methods in a non-Euclidean space. We thus propose using a stochastic gradient method in non-Euclidean space. We both provide simple conditions under which our algorithm is guaranteed to converge, and demonstrate empirically that our algorithm leads to dramatically faster training and improved predictive ability compared to stochastic gradient descent for both directed and undirected graphical models.

Noise-Driven Phenotypic Heterogeneity with Finite Correlation Time in Clonal Populations.

Directory of Open Access Journals (Sweden)

UnJin Lee

Full Text Available There has been increasing awareness in the wider biological community of the role of clonal phenotypic heterogeneity in playing key roles in phenomena such as cellular bet-hedging and decision making, as in the case of the phage-λ lysis/lysogeny and B. Subtilis competence/vegetative pathways. Here, we report on the effect of stochasticity in growth rate, cellular memory/intermittency, and its relation to phenotypic heterogeneity. We first present a linear stochastic differential model with finite auto-correlation time, where a randomly fluctuating growth rate with a negative average is shown to result in exponential growth for sufficiently large fluctuations in growth rate. We then present a non-linear stochastic self-regulation model where the loss of coherent self-regulation and an increase in noise can induce a shift from bounded to unbounded growth. An important consequence of these models is that while the average change in phenotype may not differ for various parameter sets, the variance of the resulting distributions may considerably change. This demonstrates the necessity of understanding the influence of variance and heterogeneity within seemingly identical clonal populations, while providing a mechanism for varying functional consequences of such heterogeneity. Our results highlight the importance of a paradigm shift from a deterministic to a probabilistic view of clonality in understanding selection as an optimization problem on noise-driven processes, resulting in a wide range of biological implications, from robustness to environmental stress to the development of drug resistance.

From complex to simple: interdisciplinary stochastic models

International Nuclear Information System (INIS)

Mazilu, D A; Zamora, G; Mazilu, I

2012-01-01

We present two simple, one-dimensional, stochastic models that lead to a qualitative understanding of very complex systems from biology, nanoscience and social sciences. The first model explains the complicated dynamics of microtubules, stochastic cellular highways. Using the theory of random walks in one dimension, we find analytical expressions for certain physical quantities, such as the time dependence of the length of the microtubules, and diffusion coefficients. The second one is a stochastic adsorption model with applications in surface deposition, epidemics and voter systems. We introduce the ‘empty interval method’ and show sample calculations for the time-dependent particle density. These models can serve as an introduction to the field of non-equilibrium statistical physics, and can also be used as a pedagogical tool to exemplify standard statistical physics concepts, such as random walks or the kinetic approach of the master equation. (paper)

Seismic modeling of multidimensional heterogeneity scales of Mallik gas hydrate reservoirs, Northwest Territories of Canada

Science.gov (United States)

Huang, Jun-Wei; Bellefleur, Gilles; Milkereit, Bernd

2009-07-01

In hydrate-bearing sediments, the velocity and attenuation of compressional and shear waves depend primarily on the spatial distribution of hydrates in the pore space of the subsurface lithologies. Recent characterizations of gas hydrate accumulations based on seismic velocity and attenuation generally assume homogeneous sedimentary layers and neglect effects from large- and small-scale heterogeneities of hydrate-bearing sediments. We present an algorithm, based on stochastic medium theory, to construct heterogeneous multivariable models that mimic heterogeneities of hydrate-bearing sediments at the level of detail provided by borehole logging data. Using this algorithm, we model some key petrophysical properties of gas hydrates within heterogeneous sediments near the Mallik well site, Northwest Territories, Canada. The modeled density, and P and S wave velocities used in combination with a modified Biot-Gassmann theory provide a first-order estimate of the in situ volume of gas hydrate near the Mallik 5L-38 borehole. Our results suggest a range of 528 to 768 × 106 m3/km2 of natural gas trapped within hydrates, nearly an order of magnitude lower than earlier estimates which did not include effects of small-scale heterogeneities. Further, the petrophysical models are combined with a 3-D finite difference modeling algorithm to study seismic attenuation due to scattering and leaky mode propagation. Simulations of a near-offset vertical seismic profile and cross-borehole numerical surveys demonstrate that attenuation of seismic energy may not be directly related to the intrinsic attenuation of hydrate-bearing sediments but, instead, may be largely attributed to scattering from small-scale heterogeneities and highly attenuate leaky mode propagation of seismic waves through larger-scale heterogeneities in sediments.

Stochastic fractional differential equations: Modeling, method and analysis

International Nuclear Information System (INIS)

Pedjeu, Jean-C.; Ladde, Gangaram S.

2012-01-01

By introducing a concept of dynamic process operating under multi-time scales in sciences and engineering, a mathematical model described by a system of multi-time scale stochastic differential equations is formulated. The classical Picard–Lindelöf successive approximations scheme is applied to the model validation problem, namely, existence and uniqueness of solution process. Naturally, this leads to the problem of finding closed form solutions of both linear and nonlinear multi-time scale stochastic differential equations of Itô–Doob type. Finally, to illustrate the scope of ideas and presented results, multi-time scale stochastic models for ecological and epidemiological processes in population dynamic are outlined.

A medium-term coalition-forming model of heterogeneous DERs for a commercial virtual power plant

International Nuclear Information System (INIS)

Shabanzadeh, Morteza; Sheikh-El-Eslami, Mohammad-Kazem; Haghifam, Mahmoud-Reza

2016-01-01

Highlights: • A medium-term coalition-forming scheme is proposed for commercial VPPs. • Decision making on the optimal selection of VPP coalition members, bilateral and forward contracting, and pool involvement. • VPP acts as an arbitrageur by exercising arbitrage between diverse energy trading floors. • Stochastic programming approach applied to characterize the uncertainty and to derive informed decisions. - Abstract: Within a medium-term market horizon, this research work provides a methodology that allows a commercial virtual power plant (CVPP) to form an optimal coalition of heterogeneous distributed energy resources (DERs) based on weekly bilateral contracting, futures-market involvement, and pool participation. The established model aims at composing an optimal portfolio of available DERs and jointly takes into account the risk associated with the energy trading strategy of the CVPP. Perceiving the fact that pool prices have highly uncertain nature, a framework based on stochastic programming approach is utilized to model this decision-making problem. The proposed framework consists of two stages. The first stage deals with decisions regarding DERs optimal selection for the VPP coalition, the amount of agreed quantity in the bilateral negotiation, and the type and quantity selection of futures-market contracts as well. In the second stage, decisions are made based on the most plausible realizations of the stochastic prices in the day-ahead market. For a given pre-specified risk level on profit volatility, the main objective is to maximize the expected profit for the VPP manager over the planning horizon. The efficiency and applicability of the developed model is illustrated and analyzed by its implementation in a system with few heterogeneous DERs and through different scenarios, and finally thereby meaningful conclusions are duly drawn.

Tsunamis: stochastic models of occurrence and generation mechanisms

Science.gov (United States)

Geist, Eric L.; Oglesby, David D.

2014-01-01

The devastating consequences of the 2004 Indian Ocean and 2011 Japan tsunamis have led to increased research into many different aspects of the tsunami phenomenon. In this entry, we review research related to the observed complexity and uncertainty associated with tsunami generation, propagation, and occurrence described and analyzed using a variety of stochastic methods. In each case, seismogenic tsunamis are primarily considered. Stochastic models are developed from the physical theories that govern tsunami evolution combined with empirical models fitted to seismic and tsunami observations, as well as tsunami catalogs. These stochastic methods are key to providing probabilistic forecasts and hazard assessments for tsunamis. The stochastic methods described here are similar to those described for earthquakes (Vere-Jones 2013) and volcanoes (Bebbington 2013) in this encyclopedia.

ARIMA-Based Time Series Model of Stochastic Wind Power Generation

DEFF Research Database (Denmark)

Chen, Peiyuan; Pedersen, Troels; Bak-Jensen, Birgitte

2010-01-01

This paper proposes a stochastic wind power model based on an autoregressive integrated moving average (ARIMA) process. The model takes into account the nonstationarity and physical limits of stochastic wind power generation. The model is constructed based on wind power measurement of one year from...... the Nysted offshore wind farm in Denmark. The proposed limited-ARIMA (LARIMA) model introduces a limiter and characterizes the stochastic wind power generation by mean level, temporal correlation and driving noise. The model is validated against the measurement in terms of temporal correlation...... and probability distribution. The LARIMA model outperforms a first-order transition matrix based discrete Markov model in terms of temporal correlation, probability distribution and model parameter number. The proposed LARIMA model is further extended to include the monthly variation of the stochastic wind power...

Modeling stochasticity in biochemical reaction networks

International Nuclear Information System (INIS)

Constantino, P H; Vlysidis, M; Smadbeck, P; Kaznessis, Y N

2016-01-01

Small biomolecular systems are inherently stochastic. Indeed, fluctuations of molecular species are substantial in living organisms and may result in significant variation in cellular phenotypes. The chemical master equation (CME) is the most detailed mathematical model that can describe stochastic behaviors. However, because of its complexity the CME has been solved for only few, very small reaction networks. As a result, the contribution of CME-based approaches to biology has been very limited. In this review we discuss the approach of solving CME by a set of differential equations of probability moments, called moment equations. We present different approaches to produce and to solve these equations, emphasizing the use of factorial moments and the zero information entropy closure scheme. We also provide information on the stability analysis of stochastic systems. Finally, we speculate on the utility of CME-based modeling formalisms, especially in the context of synthetic biology efforts. (topical review)

STOCHASTIC CHARACTERISTICS AND MODELING OF RELATIVE ...

African Journals Online (AJOL)

Test

Results are highly accurate and promising for all models based on Lewis' criteria. ... hydrological cycle. Future increases in ... STOCHASTIC CHARACTERISTICS AND MODELING OF RELATIVE HUMIDITY OF OGUN BASIN, NIGERIA. 71 ...

A probabilistic graphical model based stochastic input model construction

International Nuclear Information System (INIS)

Wan, Jiang; Zabaras, Nicholas

2014-01-01

Model reduction techniques have been widely used in modeling of high-dimensional stochastic input in uncertainty quantification tasks. However, the probabilistic modeling of random variables projected into reduced-order spaces presents a number of computational challenges. Due to the curse of dimensionality, the underlying dependence relationships between these random variables are difficult to capture. In this work, a probabilistic graphical model based approach is employed to learn the dependence by running a number of conditional independence tests using observation data. Thus a probabilistic model of the joint PDF is obtained and the PDF is factorized into a set of conditional distributions based on the dependence structure of the variables. The estimation of the joint PDF from data is then transformed to estimating conditional distributions under reduced dimensions. To improve the computational efficiency, a polynomial chaos expansion is further applied to represent the random field in terms of a set of standard random variables. This technique is combined with both linear and nonlinear model reduction methods. Numerical examples are presented to demonstrate the accuracy and efficiency of the probabilistic graphical model based stochastic input models. - Highlights: • Data-driven stochastic input models without the assumption of independence of the reduced random variables. • The problem is transformed to a Bayesian network structure learning problem. • Examples are given in flows in random media

Deterministic and stochastic CTMC models from Zika disease transmission

Science.gov (United States)

Zevika, Mona; Soewono, Edy

2018-03-01

Zika infection is one of the most important mosquito-borne diseases in the world. Zika virus (ZIKV) is transmitted by many Aedes-type mosquitoes including Aedes aegypti. Pregnant women with the Zika virus are at risk of having a fetus or infant with a congenital defect and suffering from microcephaly. Here, we formulate a Zika disease transmission model using two approaches, a deterministic model and a continuous-time Markov chain stochastic model. The basic reproduction ratio is constructed from a deterministic model. Meanwhile, the CTMC stochastic model yields an estimate of the probability of extinction and outbreaks of Zika disease. Dynamical simulations and analysis of the disease transmission are shown for the deterministic and stochastic models.

Weather Derivatives and Stochastic Modelling of Temperature

Directory of Open Access Journals (Sweden)

Fred Espen Benth

2011-01-01

Full Text Available We propose a continuous-time autoregressive model for the temperature dynamics with volatility being the product of a seasonal function and a stochastic process. We use the Barndorff-Nielsen and Shephard model for the stochastic volatility. The proposed temperature dynamics is flexible enough to model temperature data accurately, and at the same time being analytically tractable. Futures prices for commonly traded contracts at the Chicago Mercantile Exchange on indices like cooling- and heating-degree days and cumulative average temperatures are computed, as well as option prices on them.

Stochastic Modeling Of Wind Turbine Drivetrain Components

DEFF Research Database (Denmark)

Rafsanjani, Hesam Mirzaei; Sørensen, John Dalsgaard

2014-01-01

reliable components are needed for wind turbine. In this paper focus is on reliability of critical components in drivetrain such as bearings and shafts. High failure rates of these components imply a need for more reliable components. To estimate the reliability of these components, stochastic models...... are needed for initial defects and damage accumulation. In this paper, stochastic models are formulated considering some of the failure modes observed in these components. The models are based on theoretical considerations, manufacturing uncertainties, size effects of different scales. It is illustrated how...

Applied stochastic modelling

CERN Document Server

Morgan, Byron JT; Tanner, Martin Abba; Carlin, Bradley P

2008-01-01

Introduction and Examples Introduction Examples of data sets Basic Model Fitting Introduction Maximum-likelihood estimation for a geometric model Maximum-likelihood for the beta-geometric model Modelling polyspermy Which model? What is a model for? Mechanistic models Function Optimisation Introduction MATLAB: graphs and finite differences Deterministic search methods Stochastic search methods Accuracy and a hybrid approach Basic Likelihood ToolsIntroduction Estimating standard errors and correlations Looking at surfaces: profile log-likelihoods Confidence regions from profiles Hypothesis testing in model selectionScore and Wald tests Classical goodness of fit Model selection biasGeneral Principles Introduction Parameterisation Parameter redundancy Boundary estimates Regression and influence The EM algorithm Alternative methods of model fitting Non-regular problemsSimulation Techniques Introduction Simulating random variables Integral estimation Verification Monte Carlo inference Estimating sampling distributi...

On changes of measure in stochastic volatility models

Directory of Open Access Journals (Sweden)

Bernard Wong

2006-01-01

models. This had led many researchers to “assume the condition away,” even though the condition is not innocuous, and nonsensical results can occur if it is in fact not satisfied. We provide an applicable theorem to check the conditions for a general class of Markovian stochastic volatility models. As an example we will also provide a detailed analysis of the Stein and Stein and Heston stochastic volatility models.

Gompertzian stochastic model with delay effect to cervical cancer growth

International Nuclear Information System (INIS)

Mazlan, Mazma Syahidatul Ayuni binti; Rosli, Norhayati binti; Bahar, Arifah

2015-01-01

In this paper, a Gompertzian stochastic model with time delay is introduced to describe the cervical cancer growth. The parameters values of the mathematical model are estimated via Levenberg-Marquardt optimization method of non-linear least squares. We apply Milstein scheme for solving the stochastic model numerically. The efficiency of mathematical model is measured by comparing the simulated result and the clinical data of cervical cancer growth. Low values of Mean-Square Error (MSE) of Gompertzian stochastic model with delay effect indicate good fits

Gompertzian stochastic model with delay effect to cervical cancer growth

Energy Technology Data Exchange (ETDEWEB)

Mazlan, Mazma Syahidatul Ayuni binti; Rosli, Norhayati binti [Faculty of Industrial Sciences and Technology, Universiti Malaysia Pahang, Lebuhraya Tun Razak, 26300 Gambang, Pahang (Malaysia); Bahar, Arifah [Department of Mathematical Sciences, Faculty of Science, Universiti Teknologi Malaysia, 81310 Johor Bahru, Johor and UTM Centre for Industrial and Applied Mathematics (UTM-CIAM), Universiti Teknologi Malaysia, 81310 Johor Bahru, Johor (Malaysia)

2015-02-03

In this paper, a Gompertzian stochastic model with time delay is introduced to describe the cervical cancer growth. The parameters values of the mathematical model are estimated via Levenberg-Marquardt optimization method of non-linear least squares. We apply Milstein scheme for solving the stochastic model numerically. The efficiency of mathematical model is measured by comparing the simulated result and the clinical data of cervical cancer growth. Low values of Mean-Square Error (MSE) of Gompertzian stochastic model with delay effect indicate good fits.

Stochastic Modelling, Analysis, and Simulations of the Solar Cycle Dynamic Process

Science.gov (United States)

Turner, Douglas C.; Ladde, Gangaram S.

2018-03-01

Analytical solutions, discretization schemes and simulation results are presented for the time delay deterministic differential equation model of the solar dynamo presented by Wilmot-Smith et al. In addition, this model is extended under stochastic Gaussian white noise parametric fluctuations. The introduction of stochastic fluctuations incorporates variables affecting the dynamo process in the solar interior, estimation error of parameters, and uncertainty of the α-effect mechanism. Simulation results are presented and analyzed to exhibit the effects of stochastic parametric volatility-dependent perturbations. The results generalize and extend the work of Hazra et al. In fact, some of these results exhibit the oscillatory dynamic behavior generated by the stochastic parametric additative perturbations in the absence of time delay. In addition, the simulation results of the modified stochastic models influence the change in behavior of the very recently developed stochastic model of Hazra et al.

Stochastic models of the Social Security trust funds.

Science.gov (United States)

Burdick, Clark; Manchester, Joyce

Each year in March, the Board of Trustees of the Social Security trust funds reports on the current and projected financial condition of the Social Security programs. Those programs, which pay monthly benefits to retired workers and their families, to the survivors of deceased workers, and to disabled workers and their families, are financed through the Old-Age, Survivors, and Disability Insurance (OASDI) Trust Funds. In their 2003 report, the Trustees present, for the first time, results from a stochastic model of the combined OASDI trust funds. Stochastic modeling is an important new tool for Social Security policy analysis and offers the promise of valuable new insights into the financial status of the OASDI trust funds and the effects of policy changes. The results presented in this article demonstrate that several stochastic models deliver broadly consistent results even though they use very different approaches and assumptions. However, they also show that the variation in trust fund outcomes differs as the approach and assumptions are varied. Which approach and assumptions are best suited for Social Security policy analysis remains an open question. Further research is needed before the promise of stochastic modeling is fully realized. For example, neither parameter uncertainty nor variability in ultimate assumption values is recognized explicitly in the analyses. Despite this caveat, stochastic modeling results are already shedding new light on the range and distribution of trust fund outcomes that might occur in the future.

A Backward-Lagrangian-Stochastic Footprint Model for the Urban Environment

Science.gov (United States)

Wang, Chenghao; Wang, Zhi-Hua; Yang, Jiachuan; Li, Qi

2018-02-01

Built terrains, with their complexity in morphology, high heterogeneity, and anthropogenic impact, impose substantial challenges in Earth-system modelling. In particular, estimation of the source areas and footprints of atmospheric measurements in cities requires realistic representation of the landscape characteristics and flow physics in urban areas, but has hitherto been heavily reliant on large-eddy simulations. In this study, we developed physical parametrization schemes for estimating urban footprints based on the backward-Lagrangian-stochastic algorithm, with the built environment represented by street canyons. The vertical profile of mean streamwise velocity is parametrized for the urban canopy and boundary layer. Flux footprints estimated by the proposed model show reasonable agreement with analytical predictions over flat surfaces without roughness elements, and with experimental observations over sparse plant canopies. Furthermore, comparisons of canyon flow and turbulence profiles and the subsequent footprints were made between the proposed model and large-eddy simulation data. The results suggest that the parametrized canyon wind and turbulence statistics, based on the simple similarity theory used, need to be further improved to yield more realistic urban footprint modelling.

Managerial performance and cost efficiency of Japanese local public hospitals: a latent class stochastic frontier model.

Science.gov (United States)

Besstremyannaya, Galina

2011-09-01

The paper explores the link between managerial performance and cost efficiency of 617 Japanese general local public hospitals in 1999-2007. Treating managerial performance as unobservable heterogeneity, the paper employs a panel data stochastic cost frontier model with latent classes. Financial parameters associated with better managerial performance are found to be positively significant in explaining the probability of belonging to the more efficient latent class. The analysis of latent class membership was consistent with the conjecture that unobservable technological heterogeneity reflected in the existence of the latent classes is related to managerial performance. The findings may support the cause for raising efficiency of Japanese local public hospitals by enhancing the quality of management. Copyright © 2011 John Wiley & Sons, Ltd.

Biochemical Network Stochastic Simulator (BioNetS: software for stochastic modeling of biochemical networks

Directory of Open Access Journals (Sweden)

Elston Timothy C

2004-03-01

Full Text Available Abstract Background Intrinsic fluctuations due to the stochastic nature of biochemical reactions can have large effects on the response of biochemical networks. This is particularly true for pathways that involve transcriptional regulation, where generally there are two copies of each gene and the number of messenger RNA (mRNA molecules can be small. Therefore, there is a need for computational tools for developing and investigating stochastic models of biochemical networks. Results We have developed the software package Biochemical Network Stochastic Simulator (BioNetS for efficientlyand accurately simulating stochastic models of biochemical networks. BioNetS has a graphical user interface that allows models to be entered in a straightforward manner, and allows the user to specify the type of random variable (discrete or continuous for each chemical species in the network. The discrete variables are simulated using an efficient implementation of the Gillespie algorithm. For the continuous random variables, BioNetS constructs and numerically solvesthe appropriate chemical Langevin equations. The software package has been developed to scale efficiently with network size, thereby allowing large systems to be studied. BioNetS runs as a BioSpice agent and can be downloaded from http://www.biospice.org. BioNetS also can be run as a stand alone package. All the required files are accessible from http://x.amath.unc.edu/BioNetS. Conclusions We have developed BioNetS to be a reliable tool for studying the stochastic dynamics of large biochemical networks. Important features of BioNetS are its ability to handle hybrid models that consist of both continuous and discrete random variables and its ability to model cell growth and division. We have verified the accuracy and efficiency of the numerical methods by considering several test systems.

Electricity Market Stochastic Dynamic Model and Its Mean Stability Analysis

Directory of Open Access Journals (Sweden)

Zhanhui Lu

2014-01-01

Full Text Available Based on the deterministic dynamic model of electricity market proposed by Alvarado, a stochastic electricity market model, considering the random nature of demand sides, is presented in this paper on the assumption that generator cost function and consumer utility function are quadratic functions. The stochastic electricity market model is a generalization of the deterministic dynamic model. Using the theory of stochastic differential equations, stochastic process theory, and eigenvalue techniques, the determining conditions of the mean stability for this electricity market model under small Gauss type random excitation are provided and testified theoretically. That is, if the demand elasticity of suppliers is nonnegative and the demand elasticity of consumers is negative, then the stochastic electricity market model is mean stable. It implies that the stability can be judged directly by initial data without any computation. Taking deterministic electricity market data combined with small Gauss type random excitation as numerical samples to interpret random phenomena from a statistical perspective, the results indicate the conclusions above are correct, valid, and practical.

Aspects if stochastic models for short-term hydropower scheduling and bidding

Energy Technology Data Exchange (ETDEWEB)

Belsnes, Michael Martin [Sintef Energy, Trondheim (Norway); Follestad, Turid [Sintef Energy, Trondheim (Norway); Wolfgang, Ove [Sintef Energy, Trondheim (Norway); Fosso, Olav B. [Dep. of electric power engineering NTNU, Trondheim (Norway)

2012-07-01

This report discusses challenges met when turning from deterministic to stochastic decision support models for short-term hydropower scheduling and bidding. The report describes characteristics of the short-term scheduling and bidding problem, different market and bidding strategies, and how a stochastic optimization model can be formulated. A review of approaches for stochastic short-term modelling and stochastic modelling for the input variables inflow and market prices is given. The report discusses methods for approximating the predictive distribution of uncertain variables by scenario trees. Benefits of using a stochastic over a deterministic model are illustrated by a case study, where increased profit is obtained to a varying degree depending on the reservoir filling and price structure. Finally, an approach for assessing the effect of using a size restricted scenario tree to approximate the predictive distribution for stochastic input variables is described. The report is a summary of the findings of Work package 1 of the research project #Left Double Quotation Mark#Optimal short-term scheduling of wind and hydro resources#Right Double Quotation Mark#. The project aims at developing a prototype for an operational stochastic short-term scheduling model. Based on the investigations summarized in the report, it is concluded that using a deterministic equivalent formulation of the stochastic optimization problem is convenient and sufficient for obtaining a working prototype. (author)

Characterizing economic trends by Bayesian stochastic model specifi cation search

OpenAIRE

Grassi, Stefano; Proietti, Tommaso

2010-01-01

We apply a recently proposed Bayesian model selection technique, known as stochastic model specification search, for characterising the nature of the trend in macroeconomic time series. We illustrate that the methodology can be quite successfully applied to discriminate between stochastic and deterministic trends. In particular, we formulate autoregressive models with stochastic trends components and decide on whether a specific feature of the series, i.e. the underlying level and/or the rate...

Index Option Pricing Models with Stochastic Volatility and Stochastic Interest Rates

NARCIS (Netherlands)

Jiang, G.J.; van der Sluis, P.J.

2000-01-01

This paper specifies a multivariate stochastic volatility (SV) model for the S&P500 index and spot interest rate processes. We first estimate the multivariate SV model via the efficient method of moments (EMM) technique based on observations of underlying state variables, and then investigate the

Dynamic optimization deterministic and stochastic models

CERN Document Server

Hinderer, Karl; Stieglitz, Michael

2016-01-01

This book explores discrete-time dynamic optimization and provides a detailed introduction to both deterministic and stochastic models. Covering problems with finite and infinite horizon, as well as Markov renewal programs, Bayesian control models and partially observable processes, the book focuses on the precise modelling of applications in a variety of areas, including operations research, computer science, mathematics, statistics, engineering, economics and finance. Dynamic Optimization is a carefully presented textbook which starts with discrete-time deterministic dynamic optimization problems, providing readers with the tools for sequential decision-making, before proceeding to the more complicated stochastic models. The authors present complete and simple proofs and illustrate the main results with numerous examples and exercises (without solutions). With relevant material covered in four appendices, this book is completely self-contained.

Stochastic models for atmospheric dispersion

DEFF Research Database (Denmark)

Ditlevsen, Ove Dalager

2003-01-01

Simple stochastic differential equation models have been applied by several researchers to describe the dispersion of tracer particles in the planetary atmospheric boundary layer and to form the basis for computer simulations of particle paths. To obtain the drift coefficient, empirical vertical...... positions close to the boundaries. Different rules have been suggested in the literature with justifications based on simulation studies. Herein the relevant stochastic differential equation model is formulated in a particular way. The formulation is based on the marginal transformation of the position...... velocity distributions that depend on height above the ground both with respect to standard deviation and skewness are substituted into the stationary Fokker/Planck equation. The particle position distribution is taken to be uniform *the well/mixed condition( and also a given dispersion coefficient...

Hopf bifurcation of the stochastic model on business cycle

International Nuclear Information System (INIS)

Xu, J; Wang, H; Ge, G

2008-01-01

A stochastic model on business cycle was presented in thas paper. Simplifying the model through the quasi Hamiltonian theory, the Ito diffusion process was obtained. According to Oseledec multiplicative ergodic theory and singular boundary theory, the conditions of local and global stability were acquired. Solving the stationary FPK equation and analyzing the stationary probability density, the stochastic Hopf bifurcation was explained. The result indicated that the change of parameter awas the key factor to the appearance of the stochastic Hopf bifurcation

Integrated modeling and up-scaling of landfill processes and heterogeneity using stochastic approach

NARCIS (Netherlands)

Bun, A.; Heimovaara, T.J.; Baviskar, S.M.; van Turnhout, A.G.; Konstantaki, L.A.

2012-01-01

Municipal solid waste landfills are a very complex and heterogeneous systems. The waste in a landfill body is a heterogeneous mixture of a wide range of materials containing high levels of organic matter, high amounts of salts and a wide range of different organic and inorganic substances, such as

Stochastic models for predicting pitting corrosion damage of HLRW containers

International Nuclear Information System (INIS)

Henshall, G.A.

1991-10-01

Stochastic models for predicting aqueous pitting corrosion damage of high-level radioactive-waste containers are described. These models could be used to predict the time required for the first pit to penetrate a container and the increase in the number of breaches at later times, both of which would be useful in the repository system performance analysis. Monte Carlo implementations of the stochastic models are described, and predictions of induction time, survival probability and pit depth distributions are presented. These results suggest that the pit nucleation probability decreases with exposure time and that pit growth may be a stochastic process. The advantages and disadvantages of the stochastic approach, methods for modeling the effects of environment, and plans for future work are discussed

Brain-inspired Stochastic Models and Implementations

KAUST Repository

Al-Shedivat, Maruan

2015-05-12

One of the approaches to building artificial intelligence (AI) is to decipher the princi- ples of the brain function and to employ similar mechanisms for solving cognitive tasks, such as visual perception or natural language understanding, using machines. The recent breakthrough, named deep learning, demonstrated that large multi-layer networks of arti- ficial neural-like computing units attain remarkable performance on some of these tasks. Nevertheless, such artificial networks remain to be very loosely inspired by the brain, which rich structures and mechanisms may further suggest new algorithms or even new paradigms of computation. In this thesis, we explore brain-inspired probabilistic mechanisms, such as neural and synaptic stochasticity, in the context of generative models. The two questions we ask here are: (i) what kind of models can describe a neural learning system built of stochastic components? and (ii) how can we implement such systems e ̆ciently? To give specific answers, we consider two well known models and the corresponding neural architectures: the Naive Bayes model implemented with a winner-take-all spiking neural network and the Boltzmann machine implemented in a spiking or non-spiking fashion. We propose and analyze an e ̆cient neuromorphic implementation of the stochastic neu- ral firing mechanism and study the e ̄ects of synaptic unreliability on learning generative energy-based models implemented with neural networks.

Stochastic models to simulate paratuberculosis in dairy herds

DEFF Research Database (Denmark)

Nielsen, Søren Saxmose; Weber, M.F.; Kudahl, Anne Margrethe Braad

2011-01-01

Stochastic simulation models are widely accepted as a means of assessing the impact of changes in daily management and the control of different diseases, such as paratuberculosis, in dairy herds. This paper summarises and discusses the assumptions of four stochastic simulation models and their use...... the models are somewhat different in their underlying principles and do put slightly different values on the different strategies, their overall findings are similar. Therefore, simulation models may be useful in planning paratuberculosis strategies in dairy herds, although as with all models caution...

Stochastic population oscillations in spatial predator-prey models

International Nuclear Information System (INIS)

Taeuber, Uwe C

2011-01-01

It is well-established that including spatial structure and stochastic noise in models for predator-prey interactions invalidates the classical deterministic Lotka-Volterra picture of neutral population cycles. In contrast, stochastic models yield long-lived, but ultimately decaying erratic population oscillations, which can be understood through a resonant amplification mechanism for density fluctuations. In Monte Carlo simulations of spatial stochastic predator-prey systems, one observes striking complex spatio-temporal structures. These spreading activity fronts induce persistent correlations between predators and prey. In the presence of local particle density restrictions (finite prey carrying capacity), there exists an extinction threshold for the predator population. The accompanying continuous non-equilibrium phase transition is governed by the directed-percolation universality class. We employ field-theoretic methods based on the Doi-Peliti representation of the master equation for stochastic particle interaction models to (i) map the ensuing action in the vicinity of the absorbing state phase transition to Reggeon field theory, and (ii) to quantitatively address fluctuation-induced renormalizations of the population oscillation frequency, damping, and diffusion coefficients in the species coexistence phase.

Linking agent-based models and stochastic models of financial markets.

Science.gov (United States)

Feng, Ling; Li, Baowen; Podobnik, Boris; Preis, Tobias; Stanley, H Eugene

2012-05-29

It is well-known that financial asset returns exhibit fat-tailed distributions and long-term memory. These empirical features are the main objectives of modeling efforts using (i) stochastic processes to quantitatively reproduce these features and (ii) agent-based simulations to understand the underlying microscopic interactions. After reviewing selected empirical and theoretical evidence documenting the behavior of traders, we construct an agent-based model to quantitatively demonstrate that "fat" tails in return distributions arise when traders share similar technical trading strategies and decisions. Extending our behavioral model to a stochastic model, we derive and explain a set of quantitative scaling relations of long-term memory from the empirical behavior of individual market participants. Our analysis provides a behavioral interpretation of the long-term memory of absolute and squared price returns: They are directly linked to the way investors evaluate their investments by applying technical strategies at different investment horizons, and this quantitative relationship is in agreement with empirical findings. Our approach provides a possible behavioral explanation for stochastic models for financial systems in general and provides a method to parameterize such models from market data rather than from statistical fitting.

Stochastic differential equation model to Prendiville processes

International Nuclear Information System (INIS)

Granita; Bahar, Arifah

2015-01-01

The Prendiville process is another variation of the logistic model which assumes linearly decreasing population growth rate. It is a continuous time Markov chain (CTMC) taking integer values in the finite interval. The continuous time Markov chain can be approximated by stochastic differential equation (SDE). This paper discusses the stochastic differential equation of Prendiville process. The work started with the forward Kolmogorov equation in continuous time Markov chain of Prendiville process. Then it was formulated in the form of a central-difference approximation. The approximation was then used in Fokker-Planck equation in relation to the stochastic differential equation of the Prendiville process. The explicit solution of the Prendiville process was obtained from the stochastic differential equation. Therefore, the mean and variance function of the Prendiville process could be easily found from the explicit solution

Stochastic differential equation model to Prendiville processes

Energy Technology Data Exchange (ETDEWEB)

Granita, E-mail: granitafc@gmail.com [Dept. of Mathematical Science, Universiti Teknologi Malaysia, 81310, Johor Malaysia (Malaysia); Bahar, Arifah [Dept. of Mathematical Science, Universiti Teknologi Malaysia, 81310, Johor Malaysia (Malaysia); UTM Center for Industrial & Applied Mathematics (UTM-CIAM) (Malaysia)

2015-10-22

The Prendiville process is another variation of the logistic model which assumes linearly decreasing population growth rate. It is a continuous time Markov chain (CTMC) taking integer values in the finite interval. The continuous time Markov chain can be approximated by stochastic differential equation (SDE). This paper discusses the stochastic differential equation of Prendiville process. The work started with the forward Kolmogorov equation in continuous time Markov chain of Prendiville process. Then it was formulated in the form of a central-difference approximation. The approximation was then used in Fokker-Planck equation in relation to the stochastic differential equation of the Prendiville process. The explicit solution of the Prendiville process was obtained from the stochastic differential equation. Therefore, the mean and variance function of the Prendiville process could be easily found from the explicit solution.

Stochastic growth logistic model with aftereffect for batch fermentation process

Energy Technology Data Exchange (ETDEWEB)

Rosli, Norhayati; Ayoubi, Tawfiqullah [Faculty of Industrial Sciences and Technology, Universiti Malaysia Pahang, Lebuhraya Tun Razak, 26300 Gambang, Pahang (Malaysia); Bahar, Arifah; Rahman, Haliza Abdul [Department of Mathematical Sciences, Faculty of Science, Universiti Teknologi Malaysia, 81310 Johor Bahru, Johor (Malaysia); Salleh, Madihah Md [Department of Biotechnology Industry, Faculty of Biosciences and Bioengineering, Universiti Teknologi Malaysia, 81310 Johor Bahru, Johor (Malaysia)

2014-06-19

In this paper, the stochastic growth logistic model with aftereffect for the cell growth of C. acetobutylicum P262 and Luedeking-Piret equations for solvent production in batch fermentation system is introduced. The parameters values of the mathematical models are estimated via Levenberg-Marquardt optimization method of non-linear least squares. We apply Milstein scheme for solving the stochastic models numerically. The effciency of mathematical models is measured by comparing the simulated result and the experimental data of the microbial growth and solvent production in batch system. Low values of Root Mean-Square Error (RMSE) of stochastic models with aftereffect indicate good fits.

Stochastic growth logistic model with aftereffect for batch fermentation process

Science.gov (United States)

Rosli, Norhayati; Ayoubi, Tawfiqullah; Bahar, Arifah; Rahman, Haliza Abdul; Salleh, Madihah Md

2014-06-01

In this paper, the stochastic growth logistic model with aftereffect for the cell growth of C. acetobutylicum P262 and Luedeking-Piret equations for solvent production in batch fermentation system is introduced. The parameters values of the mathematical models are estimated via Levenberg-Marquardt optimization method of non-linear least squares. We apply Milstein scheme for solving the stochastic models numerically. The effciency of mathematical models is measured by comparing the simulated result and the experimental data of the microbial growth and solvent production in batch system. Low values of Root Mean-Square Error (RMSE) of stochastic models with aftereffect indicate good fits.

Stochastic growth logistic model with aftereffect for batch fermentation process

International Nuclear Information System (INIS)

Rosli, Norhayati; Ayoubi, Tawfiqullah; Bahar, Arifah; Rahman, Haliza Abdul; Salleh, Madihah Md

2014-01-01

In this paper, the stochastic growth logistic model with aftereffect for the cell growth of C. acetobutylicum P262 and Luedeking-Piret equations for solvent production in batch fermentation system is introduced. The parameters values of the mathematical models are estimated via Levenberg-Marquardt optimization method of non-linear least squares. We apply Milstein scheme for solving the stochastic models numerically. The effciency of mathematical models is measured by comparing the simulated result and the experimental data of the microbial growth and solvent production in batch system. Low values of Root Mean-Square Error (RMSE) of stochastic models with aftereffect indicate good fits

Transport in Stochastic Media

International Nuclear Information System (INIS)

Haran, O.; Shvarts, D.; Thieberger, R.

1998-01-01

Classical transport of neutral particles in a binary, scattering, stochastic media is discussed. It is assumed that the cross-sections of the constituent materials and their volume fractions are known. The inner structure of the media is stochastic, but there exist a statistical knowledge about the lump sizes, shapes and arrangement. The transmission through the composite media depends on the specific heterogeneous realization of the media. The current research focuses on the averaged transmission through an ensemble of realizations, frm which an effective cross-section for the media can be derived. The problem of one dimensional transport in stochastic media has been studied extensively [1]. In the one dimensional description of the problem, particles are transported along a line populated with alternating material segments of random lengths. The current work discusses transport in two-dimensional stochastic media. The phenomenon that is unique to the multi-dimensional description of the problem is obstacle bypassing. Obstacle bypassing tends to reduce the opacity of the media, thereby reducing its effective cross-section. The importance of this phenomenon depends on the manner in which the obstacles are arranged in the media. Results of transport simulations in multi-dimensional stochastic media are presented. Effective cross-sections derived from the simulations are compared against those obtained for the one-dimensional problem, and against those obtained from effective multi-dimensional models, which are partially based on a Markovian assumption

Towards Model Checking Stochastic Process Algebra

NARCIS (Netherlands)

Hermanns, H.; Grieskamp, W.; Santen, T.; Katoen, Joost P.; Stoddart, B.; Meyer-Kayser, J.; Siegle, M.

2000-01-01

Stochastic process algebras have been proven useful because they allow behaviour-oriented performance and reliability modelling. As opposed to traditional performance modelling techniques, the behaviour- oriented style supports composition and abstraction in a natural way. However, analysis of

Microstructural Quantification, Property Prediction, and Stochastic Reconstruction of Heterogeneous Materials Using Limited X-Ray Tomography Data

Science.gov (United States)

Li, Hechao

An accurate knowledge of the complex microstructure of a heterogeneous material is crucial for quantitative structure-property relations establishment and its performance prediction and optimization. X-ray tomography has provided a non-destructive means for microstructure characterization in both 3D and 4D (i.e., structural evolution over time). Traditional reconstruction algorithms like filtered-back-projection (FBP) method or algebraic reconstruction techniques (ART) require huge number of tomographic projections and segmentation process before conducting microstructural quantification. This can be quite time consuming and computationally intensive. In this thesis, a novel procedure is first presented that allows one to directly extract key structural information in forms of spatial correlation functions from limited x-ray tomography data. The key component of the procedure is the computation of a "probability map", which provides the probability of an arbitrary point in the material system belonging to specific phase. The correlation functions of interest are then readily computed from the probability map. Using effective medium theory, accurate predictions of physical properties (e.g., elastic moduli) can be obtained. Secondly, a stochastic optimization procedure that enables one to accurately reconstruct material microstructure from a small number of x-ray tomographic projections (e.g., 20 - 40) is presented. Moreover, a stochastic procedure for multi-modal data fusion is proposed, where both X-ray projections and correlation functions computed from limited 2D optical images are fused to accurately reconstruct complex heterogeneous materials in 3D. This multi-modal reconstruction algorithm is proved to be able to integrate the complementary data to perform an excellent optimization procedure, which indicates its high efficiency in using limited structural information. Finally, the accuracy of the stochastic reconstruction procedure using limited X

Stochastic forward and inverse groundwater flow and solute transport modeling

NARCIS (Netherlands)

Janssen, G.M.C.M.

2008-01-01

Keywords: calibration, inverse modeling, stochastic modeling, nonlinear biodegradation, stochastic-convective, advective-dispersive, travel time, network design, non-Gaussian distribution, multimodal distribution, representers

This thesis offers three new approaches that contribute

Stochastic models for predicting environmental impact in aquatic ecosystems

International Nuclear Information System (INIS)

Stewart-Oaten, A.

1986-01-01

The purpose of stochastic predictions are discussed in relation to the environmental impacts of nuclear power plants on aquatic ecosystems. One purpose is to aid in making rational decisions about whether a power plant should be built, where, and how it should be designed. The other purpose is to check on the models themselves in the light of what eventually happens. The author discusses the role or statistical decision theory in the decision-making problem. Various types of stochastic models and their problems are presented. In addition some suggestions are made for generating usable stochastic models, and checking and improving on them. 12 references

Introduction to modeling and analysis of stochastic systems

CERN Document Server

Kulkarni, V G

2011-01-01

This is an introductory-level text on stochastic modeling. It is suited for undergraduate students in engineering, operations research, statistics, mathematics, actuarial science, business management, computer science, and public policy. It employs a large number of examples to teach the students to use stochastic models of real-life systems to predict their performance, and use this analysis to design better systems. The book is devoted to the study of important classes of stochastic processes: discrete and continuous time Markov processes, Poisson processes, renewal and regenerative processes, semi-Markov processes, queueing models, and diffusion processes. The book systematically studies the short-term and the long-term behavior, cost/reward models, and first passage times. All the material is illustrated with many examples, and case studies. The book provides a concise review of probability in the appendix. The book emphasizes numerical answers to the problems. A collection of MATLAB programs to accompany...

Stochastic Averaging and Stochastic Extremum Seeking

CERN Document Server

Liu, Shu-Jun

2012-01-01

Stochastic Averaging and Stochastic Extremum Seeking develops methods of mathematical analysis inspired by the interest in reverse engineering  and analysis of bacterial  convergence by chemotaxis and to apply similar stochastic optimization techniques in other environments. The first half of the text presents significant advances in stochastic averaging theory, necessitated by the fact that existing theorems are restricted to systems with linear growth, globally exponentially stable average models, vanishing stochastic perturbations, and prevent analysis over infinite time horizon. The second half of the text introduces stochastic extremum seeking algorithms for model-free optimization of systems in real time using stochastic perturbations for estimation of their gradients. Both gradient- and Newton-based algorithms are presented, offering the user the choice between the simplicity of implementation (gradient) and the ability to achieve a known, arbitrary convergence rate (Newton). The design of algorithms...

A matrix approach to the statistics of longevity in heterogeneous frailty models

Directory of Open Access Journals (Sweden)

Hal Caswell

2014-09-01

Full Text Available Background: The gamma-Gompertz model is a fixed frailty model in which baseline mortality increasesexponentially with age, frailty has a proportional effect on mortality, and frailty at birth follows a gamma distribution. Mortality selects against the more frail, so the marginal mortality rate decelerates, eventually reaching an asymptote. The gamma-Gompertz is one of a wider class of frailty models, characterized by the choice of baseline mortality, effects of frailty, distributions of frailty, and assumptions about the dynamics of frailty. Objective: To develop a matrix model to compute all the statistical properties of longevity from thegamma-Gompertz and related models. Methods: I use the vec-permutation matrix formulation to develop a model in which individuals are jointly classified by age and frailty. The matrix is used to project the age and frailty dynamicsof a cohort and the fundamental matrix is used to obtain the statistics of longevity. Results: The model permits calculation of the mean, variance, coefficient of variation, skewness and all moments of longevity, the marginal mortality and survivorship functions, the dynamics of the frailty distribution, and other quantities. The matrix formulation extends naturally to other frailty models. I apply the analysis to the gamma-Gompertz model (for humans and laboratory animals, the gamma-Makeham model, and the gamma-Siler model, and to a hypothetical dynamic frailty model characterized by diffusion of frailty with reflecting boundaries.The matrix model permits partitioning the variance in longevity into components due to heterogeneity and to individual stochasticity. In several published human data sets, heterogeneity accounts for less than 10Š of the variance in longevity. In laboratory populations of five invertebrate animal species, heterogeneity accounts for 46Š to 83Š ofthe total variance in longevity.

Predicting Footbridge Response using Stochastic Load Models

DEFF Research Database (Denmark)

Pedersen, Lars; Frier, Christian

2013-01-01

Walking parameters such as step frequency, pedestrian mass, dynamic load factor, etc. are basically stochastic, although it is quite common to adapt deterministic models for these parameters. The present paper considers a stochastic approach to modeling the action of pedestrians, but when doing so...... decisions need to be made in terms of statistical distributions of walking parameters and in terms of the parameters describing the statistical distributions. The paper explores how sensitive computations of bridge response are to some of the decisions to be made in this respect. This is useful...

Threshold Dynamics of a Stochastic Chemostat Model with Two Nutrients and One Microorganism

Directory of Open Access Journals (Sweden)

Jian Zhang

2017-01-01

Full Text Available A new stochastic chemostat model with two substitutable nutrients and one microorganism is proposed and investigated. Firstly, for the corresponding deterministic model, the threshold for extinction and permanence of the microorganism is obtained by analyzing the stability of the equilibria. Then, for the stochastic model, the threshold of the stochastic chemostat for extinction and permanence of the microorganism is explored. Difference of the threshold of the deterministic model and the stochastic model shows that a large stochastic disturbance can affect the persistence of the microorganism and is harmful to the cultivation of the microorganism. To illustrate this phenomenon, we give some computer simulations with different intensity of stochastic noise disturbance.

Transfers in heterogeneous environments; Transferts en milieux heterogenes

Energy Technology Data Exchange (ETDEWEB)

Flesselles, J M [Saint-Gobain Recherche, 93 - Aubervilliers (France); Gouesbet, G; Mees, L; Roze, C; Girasole, Th; Grehan, G [Laboratoire d' Electromagnetisme et Systemes Particulaires (LESP), UMR CNRS 6614, CORIA. Universite de Rouen et INSA de Rouen, 76 - Saint-Etienne du Rouvray (France); Goyheneche, J M; Vignoles, G; Coindreau, O [Laboratoire des Composites Thermostructuraux (LCTS), UMR 5801, 33 - Pessac (France); Moyne, Ch [LEMTA (UMR 7563) CNRS-INPL-UHP, 54 - Vandoeuvre les Nancy (France); Coussy, O [Institut Navier - ENPC, 77 - Marne-la-Vallee (France); Lassabatere, Th [Electricite de France Les Renardieres, Dept. Materiaux Mecanique des Composants, 77 - Moret sur Loing (France); Tadrist, L [IUSTI - UMR 6595, 13 - Marseille (France)

2004-07-01

This document gathers the articles and transparencies of the invited talks given at the 2004 French congress of thermal engineering about transfers in heterogeneous environment. Content: transfer phenomena in industrial glass furnaces; simple and multiple scattering diagnosis by femto-second pulsed laser: application to particulate diagnoses; thermal modeling of thermo-structural composites; hybrid mixtures theory, average volumic measurement, periodical or stochastic homogenization: advance in scale change processes; thermo-hydro-chemical-mechanical coupling in porous medium: application to young concrete structures and to clay barriers of disposal facilities; transfers and flows in fluidization: recent advances and future challenges. (J.S.)

The Asymptotic Behaviour of a Stochastic 3D LANS-α Model

International Nuclear Information System (INIS)

Caraballo, Tomas; Marquez-Duran, Antonio M.; Real, Jose

2006-01-01

The long-time behaviour of a stochastic 3D LANS-α model on a bounded domain is analysed. First, we reformulate the model as an abstract problem. Next, we establish sufficient conditions ensuring the existence of stationary (steady state) solutions of this abstract nonlinear stochastic evolution equation, and study the stability properties of the model. Finally, we analyse the effects produced by stochastic perturbations in the deterministic version of the system (persistence of exponential stability as well as possible stabilisation effects produced by the noise). The general results are applied to our stochastic LANS-α system throughout the paper

Fitting PAC spectra with stochastic models: PolyPacFit

Energy Technology Data Exchange (ETDEWEB)

Zacate, M. O., E-mail: zacatem1@nku.edu [Northern Kentucky University, Department of Physics and Geology (United States); Evenson, W. E. [Utah Valley University, College of Science and Health (United States); Newhouse, R.; Collins, G. S. [Washington State University, Department of Physics and Astronomy (United States)

2010-04-15

PolyPacFit is an advanced fitting program for time-differential perturbed angular correlation (PAC) spectroscopy. It incorporates stochastic models and provides robust options for customization of fits. Notable features of the program include platform independence and support for (1) fits to stochastic models of hyperfine interactions, (2) user-defined constraints among model parameters, (3) fits to multiple spectra simultaneously, and (4) any spin nuclear probe.

Stochastic modelling of two-phase flows including phase change

International Nuclear Information System (INIS)

Hurisse, O.; Minier, J.P.

2011-01-01

Stochastic modelling has already been developed and applied for single-phase flows and incompressible two-phase flows. In this article, we propose an extension of this modelling approach to two-phase flows including phase change (e.g. for steam-water flows). Two aspects are emphasised: a stochastic model accounting for phase transition and a modelling constraint which arises from volume conservation. To illustrate the whole approach, some remarks are eventually proposed for two-fluid models. (authors)

Improving sub-pixel imperviousness change prediction by ensembling heterogeneous non-linear regression models

Science.gov (United States)

Drzewiecki, Wojciech

2016-12-01

In this work nine non-linear regression models were compared for sub-pixel impervious surface area mapping from Landsat images. The comparison was done in three study areas both for accuracy of imperviousness coverage evaluation in individual points in time and accuracy of imperviousness change assessment. The performance of individual machine learning algorithms (Cubist, Random Forest, stochastic gradient boosting of regression trees, k-nearest neighbors regression, random k-nearest neighbors regression, Multivariate Adaptive Regression Splines, averaged neural networks, and support vector machines with polynomial and radial kernels) was also compared with the performance of heterogeneous model ensembles constructed from the best models trained using particular techniques. The results proved that in case of sub-pixel evaluation the most accurate prediction of change may not necessarily be based on the most accurate individual assessments. When single methods are considered, based on obtained results Cubist algorithm may be advised for Landsat based mapping of imperviousness for single dates. However, Random Forest may be endorsed when the most reliable evaluation of imperviousness change is the primary goal. It gave lower accuracies for individual assessments, but better prediction of change due to more correlated errors of individual predictions. Heterogeneous model ensembles performed for individual time points assessments at least as well as the best individual models. In case of imperviousness change assessment the ensembles always outperformed single model approaches. It means that it is possible to improve the accuracy of sub-pixel imperviousness change assessment using ensembles of heterogeneous non-linear regression models.

A stochastic model for quantum measurement

International Nuclear Information System (INIS)

Budiyono, Agung

2013-01-01

We develop a statistical model of microscopic stochastic deviation from classical mechanics based on a stochastic process with a transition probability that is assumed to be given by an exponential distribution of infinitesimal stationary action. We apply the statistical model to stochastically modify a classical mechanical model for the measurement of physical quantities reproducing the prediction of quantum mechanics. The system+apparatus always has a definite configuration at all times, as in classical mechanics, fluctuating randomly following a continuous trajectory. On the other hand, the wavefunction and quantum mechanical Hermitian operator corresponding to the physical quantity arise formally as artificial mathematical constructs. During a single measurement, the wavefunction of the whole system+apparatus evolves according to a Schrödinger equation and the configuration of the apparatus acts as the pointer of the measurement so that there is no wavefunction collapse. We will also show that while the outcome of each single measurement event does not reveal the actual value of the physical quantity prior to measurement, its average in an ensemble of identical measurements is equal to the average of the actual value of the physical quantity prior to measurement over the distribution of the configuration of the system. (paper)

The critical domain size of stochastic population models.

Science.gov (United States)

Reimer, Jody R; Bonsall, Michael B; Maini, Philip K

2017-02-01

Identifying the critical domain size necessary for a population to persist is an important question in ecology. Both demographic and environmental stochasticity impact a population's ability to persist. Here we explore ways of including this variability. We study populations with distinct dispersal and sedentary stages, which have traditionally been modelled using a deterministic integrodifference equation (IDE) framework. Individual-based models (IBMs) are the most intuitive stochastic analogues to IDEs but yield few analytic insights. We explore two alternate approaches; one is a scaling up to the population level using the Central Limit Theorem, and the other a variation on both Galton-Watson branching processes and branching processes in random environments. These branching process models closely approximate the IBM and yield insight into the factors determining the critical domain size for a given population subject to stochasticity.

Alternative Approaches to Technical Efficiency Estimation in the Stochastic Frontier Model

OpenAIRE

Acquah, H. de-Graft; Onumah, E. E.

2014-01-01

Estimating the stochastic frontier model and calculating technical efficiency of decision making units are of great importance in applied production economic works. This paper estimates technical efficiency from the stochastic frontier model using Jondrow, and Battese and Coelli approaches. In order to compare alternative methods, simulated data with sample sizes of 60 and 200 are generated from stochastic frontier model commonly applied to agricultural firms. Simulated data is employed to co...

The effect of stochasticity on the lac operon: an evolutionary perspective.

Directory of Open Access Journals (Sweden)

Milan van Hoek

2007-06-01

Full Text Available The role of stochasticity on gene expression is widely discussed. Both potential advantages and disadvantages have been revealed. In some systems, noise in gene expression has been quantified, in among others the lac operon of Escherichia coli. Whether stochastic gene expression in this system is detrimental or beneficial for the cells is, however, still unclear. We are interested in the effects of stochasticity from an evolutionary point of view. We study this question in the lac operon, taking a computational approach: using a detailed, quantitative, spatial model, we evolve through a mutation-selection process the shape of the promoter function and therewith the effective amount of stochasticity. We find that noise values for lactose, the natural inducer, are much lower than for artificial, nonmetabolizable inducers, because these artificial inducers experience a stronger positive feedback. In the evolved promoter functions, noise due to stochasticity in gene expression, when induced by lactose, only plays a very minor role in short-term physiological adaptation, because other sources of population heterogeneity dominate. Finally, promoter functions evolved in the stochastic model evolve to higher repressed transcription rates than those evolved in a deterministic version of the model. This causes these promoter functions to experience less stochasticity in gene expression. We show that a high repression rate and hence high stochasticity increases the delay in lactose uptake in a variable environment. We conclude that the lac operon evolved such that the impact of stochastic gene expression is minor in its natural environment, but happens to respond with much stronger stochasticity when confronted with artificial inducers. In this particular system, we have shown that stochasticity is detrimental. Moreover, we demonstrate that in silico evolution in a quantitative model, by mutating the parameters of interest, is a promising way to unravel

Stochastic Modelling of Shiroro River Stream flow Process

OpenAIRE

Musa, J. J

2013-01-01

Economists, social scientists and engineers provide insights into the drivers of anthropogenic climate change and the options for adaptation and mitigation, and yet other scientists, including geographers and biologists, study the impacts of climate change. This project concentrates mainly on the discharge from the Shiroro River. A stochastic approach is presented for modeling a time series by an Autoregressive Moving Average model (ARMA). The development and use of a stochastic stream flow m...

TIME-DEPENDENT STOCHASTIC ACCELERATION MODEL FOR FERMI BUBBLES

Energy Technology Data Exchange (ETDEWEB)

Sasaki, Kento; Asano, Katsuaki; Terasawa, Toshio, E-mail: kentos@icrr.u-tokyo.ac.jp, E-mail: asanok@icrr.u-tokyo.ac.jp, E-mail: terasawa@icrr.u-tokyo.ac.jp [Institute for Cosmic Ray Research, The University of Tokyo, 5-1-5 Kashiwanoha, Kashiwa, Chiba 277-8582 (Japan)

2015-12-01

We study stochastic acceleration models for the Fermi bubbles. Turbulence is excited just behind the shock front via Kelvin–Helmholtz, Rayleigh–Taylor, or Richtmyer–Meshkov instabilities, and plasma particles are continuously accelerated by the interaction with the turbulence. The turbulence gradually decays as it goes away from the shock fronts. Adopting a phenomenological model for the stochastic acceleration, we explicitly solve the temporal evolution of the particle energy distribution in the turbulence. Our results show that the spatial distribution of high-energy particles is different from those for a steady solution. We also show that the contribution of electrons that escaped from the acceleration regions significantly softens the photon spectrum. The photon spectrum and surface brightness profile are reproduced by our models. If the escape efficiency is very high, the radio flux from the escaped low-energy electrons can be comparable to that of the WMAP haze. We also demonstrate hadronic models with the stochastic acceleration, but they are unlikely in the viewpoint of the energy budget.

Computer Aided Continuous Time Stochastic Process Modelling

DEFF Research Database (Denmark)

Kristensen, N.R.; Madsen, Henrik; Jørgensen, Sten Bay

2001-01-01

A grey-box approach to process modelling that combines deterministic and stochastic modelling is advocated for identification of models for model-based control of batch and semi-batch processes. A computer-aided tool designed for supporting decision-making within the corresponding modelling cycle...

Review of "Stochastic Modelling for Systems Biology" by Darren Wilkinson

Directory of Open Access Journals (Sweden)

Bullinger Eric

2006-12-01

Full Text Available Abstract "Stochastic Modelling for Systems Biology" by Darren Wilkinson introduces the peculiarities of stochastic modelling in biology. This book is particularly suited to as a textbook or for self-study, and for readers with a theoretical background.

Numerical Approach to Spatial Deterministic-Stochastic Models Arising in Cell Biology.

Science.gov (United States)

Schaff, James C; Gao, Fei; Li, Ye; Novak, Igor L; Slepchenko, Boris M

2016-12-01

Hybrid deterministic-stochastic methods provide an efficient alternative to a fully stochastic treatment of models which include components with disparate levels of stochasticity. However, general-purpose hybrid solvers for spatially resolved simulations of reaction-diffusion systems are not widely available. Here we describe fundamentals of a general-purpose spatial hybrid method. The method generates realizations of a spatially inhomogeneous hybrid system by appropriately integrating capabilities of a deterministic partial differential equation solver with a popular particle-based stochastic simulator, Smoldyn. Rigorous validation of the algorithm is detailed, using a simple model of calcium 'sparks' as a testbed. The solver is then applied to a deterministic-stochastic model of spontaneous emergence of cell polarity. The approach is general enough to be implemented within biologist-friendly software frameworks such as Virtual Cell.

Stochastic Modelling of River Geometry

DEFF Research Database (Denmark)

Sørensen, John Dalsgaard; Schaarup-Jensen, K.

1996-01-01

Numerical hydrodynamic river models are used in a large number of applications to estimate critical events for rivers. These estimates are subject to a number of uncertainties. In this paper, the problem to evaluate these estimates using probabilistic methods is considered. Stochastic models for ...... for river geometries are formulated and a coupling between hydraulic computational methods and numerical reliability methods is presented....

A stochastic agent-based model of pathogen propagation in dynamic multi-relational social networks

Science.gov (United States)

Khan, Bilal; Dombrowski, Kirk; Saad, Mohamed

2015-01-01

We describe a general framework for modeling and stochastic simulation of epidemics in realistic dynamic social networks, which incorporates heterogeneity in the types of individuals, types of interconnecting risk-bearing relationships, and types of pathogens transmitted across them. Dynamism is supported through arrival and departure processes, continuous restructuring of risk relationships, and changes to pathogen infectiousness, as mandated by natural history; dynamism is regulated through constraints on the local agency of individual nodes and their risk behaviors, while simulation trajectories are validated using system-wide metrics. To illustrate its utility, we present a case study that applies the proposed framework towards a simulation of HIV in artificial networks of intravenous drug users (IDUs) modeled using data collected in the Social Factors for HIV Risk survey. PMID:25859056

Accounting for the Decreasing Reaction Potential of Heterogeneous Aquifers in a Stochastic Framework of Aquifer-Scale Reactive Transport

Science.gov (United States)

Loschko, Matthias; Wöhling, Thomas; Rudolph, David L.; Cirpka, Olaf A.

2018-01-01

Many groundwater contaminants react with components of the aquifer matrix, causing a depletion of the aquifer's reactivity with time. We discuss conceptual simplifications of reactive transport that allow the implementation of a decreasing reaction potential in reactive-transport simulations in chemically and hydraulically heterogeneous aquifers without relying on a fully explicit description. We replace spatial coordinates by travel-times and use the concept of relative reactivity, which represents the reaction-partner supply from the matrix relative to a reference. Microorganisms facilitating the reactions are not explicitly modeled. Solute mixing is neglected. Streamlines, obtained by particle tracking, are discretized in travel-time increments with variable content of reaction partners in the matrix. As exemplary reactive system, we consider aerobic respiration and denitrification with simplified reaction equations: Dissolved oxygen undergoes conditional zero-order decay, nitrate follows first-order decay, which is inhibited in the presence of dissolved oxygen. Both reactions deplete the bioavailable organic carbon of the matrix, which in turn determines the relative reactivity. These simplifications reduce the computational effort, facilitating stochastic simulations of reactive transport on the aquifer scale. In a one-dimensional test case with a more detailed description of the reactions, we derive a potential relationship between the bioavailable organic-carbon content and the relative reactivity. In a three-dimensional steady-state test case, we use the simplified model to calculate the decreasing denitrification potential of an artificial aquifer over 200 years in an ensemble of 200 members. We demonstrate that the uncertainty in predicting the nitrate breakthrough in a heterogeneous aquifer decreases with increasing scale of observation.

Stochastic Cell Fate Progression in Embryonic Stem Cells

Science.gov (United States)

Zou, Ling-Nan; Doyle, Adele; Jang, Sumin; Ramanathan, Sharad

2013-03-01

Studies on the directed differentiation of embryonic stem (ES) cells suggest that some early developmental decisions may be stochastic in nature. To identify the sources of this stochasticity, we analyzed the heterogeneous expression of key transcription factors in single ES cells as they adopt distinct germ layer fates. We find that under sufficiently stringent signaling conditions, the choice of lineage is unambiguous. ES cells flow into differentiated fates via diverging paths, defined by sequences of transitional states that exhibit characteristic co-expression of multiple transcription factors. These transitional states have distinct responses to morphogenic stimuli; by sequential exposure to multiple signaling conditions, ES cells are steered towards specific fates. However, the rate at which cells travel down a developmental path is stochastic: cells exposed to the same signaling condition for the same amount of time can populate different states along the same path. The heterogeneity of cell states seen in our experiments therefore does not reflect the stochastic selection of germ layer fates, but the stochastic rate of progression along a chosen developmental path. Supported in part by the Jane Coffin Childs Fund

Stochastic Volatility and DSGE Models

DEFF Research Database (Denmark)

Andreasen, Martin Møller

This paper argues that a specification of stochastic volatility commonly used to analyze the Great Moderation in DSGE models may not be appropriate, because the level of a process with this specification does not have conditional or unconditional moments. This is unfortunate because agents may...

The finite state projection approach to analyze dynamics of heterogeneous populations

Science.gov (United States)

Johnson, Rob; Munsky, Brian

2017-06-01

Population modeling aims to capture and predict the dynamics of cell populations in constant or fluctuating environments. At the elementary level, population growth proceeds through sequential divisions of individual cells. Due to stochastic effects, populations of cells are inherently heterogeneous in phenotype, and some phenotypic variables have an effect on division or survival rates, as can be seen in partial drug resistance. Therefore, when modeling population dynamics where the control of growth and division is phenotype dependent, the corresponding model must take account of the underlying cellular heterogeneity. The finite state projection (FSP) approach has often been used to analyze the statistics of independent cells. Here, we extend the FSP analysis to explore the coupling of cell dynamics and biomolecule dynamics within a population. This extension allows a general framework with which to model the state occupations of a heterogeneous, isogenic population of dividing and expiring cells. The method is demonstrated with a simple model of cell-cycle progression, which we use to explore possible dynamics of drug resistance phenotypes in dividing cells. We use this method to show how stochastic single-cell behaviors affect population level efficacy of drug treatments, and we illustrate how slight modifications to treatment regimens may have dramatic effects on drug efficacy.

Stability analysis of multi-group deterministic and stochastic epidemic models with vaccination rate

International Nuclear Information System (INIS)

Wang Zhi-Gang; Gao Rui-Mei; Fan Xiao-Ming; Han Qi-Xing

2014-01-01

We discuss in this paper a deterministic multi-group MSIR epidemic model with a vaccination rate, the basic reproduction number ℛ 0 , a key parameter in epidemiology, is a threshold which determines the persistence or extinction of the disease. By using Lyapunov function techniques, we show if ℛ 0 is greater than 1 and the deterministic model obeys some conditions, then the disease will prevail, the infective persists and the endemic state is asymptotically stable in a feasible region. If ℛ 0 is less than or equal to 1, then the infective disappear so the disease dies out. In addition, stochastic noises around the endemic equilibrium will be added to the deterministic MSIR model in order that the deterministic model is extended to a system of stochastic ordinary differential equations. In the stochastic version, we carry out a detailed analysis on the asymptotic behavior of the stochastic model. In addition, regarding the value of ℛ 0 , when the stochastic system obeys some conditions and ℛ 0 is greater than 1, we deduce the stochastic system is stochastically asymptotically stable. Finally, the deterministic and stochastic model dynamics are illustrated through computer simulations. (general)

Stochastic differential equations used to model conjugation

DEFF Research Database (Denmark)

Philipsen, Kirsten Riber; Christiansen, Lasse Engbo

Stochastic differential equations (SDEs) are used to model horizontal transfer of antibiotic resis- tance by conjugation. The model describes the concentration of donor, recipient, transconjugants and substrate. The strength of the SDE model over the traditional ODE models is that the noise can...

Analysis of radon-induced lung cancer risk by a stochastic state-vector model of radiation carcinogenesis

International Nuclear Information System (INIS)

Crawford-Brown, Douglas J.; Hofmann, Werner

2002-01-01

A biologically based state-vector model (SVM) of radiation carcinogenesis has been extended to incorporate stochasticity of cellular transitions and specific in vivo irradiation conditions in the lungs. Dose-rate-dependent cellular transitions related to the formation of double-stranded DNA breaks, repair of breaks, interactions (translocations) between breaks, fixation of breaks, cellular inactivation, stimulated mitosis and promotion through loss of intercellular communication are simulated by Monte Carlo methods. The stochastic SVM has been applied to the analysis of lung cancer incidence in uranium miners exposed to alpha-emitting radon progeny. When incorporating in vivo features of cell differentiation, stimulated cell division and heterogeneity of cellular doses into the model, excellent agreement between epidemiological data and modelling results could be obtained. At low doses, the model predicts a non-linear dose-response relationship; e.g., computed lung cancer risk at 20 WLM is about half of current lung cancer estimates based on the linear hypothesis. The model also predicts a slight dose rate effect; e.g., at a cumulative exposure of 20 WLM, calculated lung cancer incidence for an exposure rate 0.27 WLM/year (assuming an exposure time of 73 years) is smaller by a factor of 1.2 than that for an exposure rate of 10 WLM/year. (author)

Stochastic Differential Equation-Based Flexible Software Reliability Growth Model

Directory of Open Access Journals (Sweden)

P. K. Kapur

2009-01-01

Full Text Available Several software reliability growth models (SRGMs have been developed by software developers in tracking and measuring the growth of reliability. As the size of software system is large and the number of faults detected during the testing phase becomes large, so the change of the number of faults that are detected and removed through each debugging becomes sufficiently small compared with the initial fault content at the beginning of the testing phase. In such a situation, we can model the software fault detection process as a stochastic process with continuous state space. In this paper, we propose a new software reliability growth model based on Itô type of stochastic differential equation. We consider an SDE-based generalized Erlang model with logistic error detection function. The model is estimated and validated on real-life data sets cited in literature to show its flexibility. The proposed model integrated with the concept of stochastic differential equation performs comparatively better than the existing NHPP-based models.

A primer on stochastic epidemic models: Formulation, numerical simulation, and analysis

Directory of Open Access Journals (Sweden)

Linda J.S. Allen

2017-05-01

Full Text Available Some mathematical methods for formulation and numerical simulation of stochastic epidemic models are presented. Specifically, models are formulated for continuous-time Markov chains and stochastic differential equations. Some well-known examples are used for illustration such as an SIR epidemic model and a host-vector malaria model. Analytical methods for approximating the probability of a disease outbreak are also discussed. Keywords: Branching process, Continuous-time Markov chain, Minor outbreak, Stochastic differential equation, 2000 MSC: 60H10, 60J28, 92D30

Some Remarks on Stochastic Versions of the Ramsey Growth Model

Czech Academy of Sciences Publication Activity Database

Sladký, Karel

2012-01-01

Roč. 19, č. 29 (2012), s. 139-152 ISSN 1212-074X R&D Projects: GA ČR GAP402/10/1610; GA ČR GAP402/10/0956; GA ČR GAP402/11/0150 Institutional support: RVO:67985556 Keywords : Economic dynamics * Ramsey growth model with disturbance * stochastic dynamic programming * multistage stochastic programs Subject RIV: BB - Applied Statistics, Operational Research http://library.utia.cas.cz/separaty/2013/E/sladky-some remarks on stochastic versions of the ramsey growth model.pdf

Stochastic modeling of consumer preferences for health care institutions.

Science.gov (United States)

Malhotra, N K

1983-01-01

This paper proposes a stochastic procedure for modeling consumer preferences via LOGIT analysis. First, a simple, non-technical exposition of the use of a stochastic approach in health care marketing is presented. Second, a study illustrating the application of the LOGIT model in assessing consumer preferences for hospitals is given. The paper concludes with several implications of the proposed approach.

Simulating biological processes: stochastic physics from whole cells to colonies

Science.gov (United States)

Earnest, Tyler M.; Cole, John A.; Luthey-Schulten, Zaida

2018-05-01

The last few decades have revealed the living cell to be a crowded spatially heterogeneous space teeming with biomolecules whose concentrations and activities are governed by intrinsically random forces. It is from this randomness, however, that a vast array of precisely timed and intricately coordinated biological functions emerge that give rise to the complex forms and behaviors we see in the biosphere around us. This seemingly paradoxical nature of life has drawn the interest of an increasing number of physicists, and recent years have seen stochastic modeling grow into a major subdiscipline within biological physics. Here we review some of the major advances that have shaped our understanding of stochasticity in biology. We begin with some historical context, outlining a string of important experimental results that motivated the development of stochastic modeling. We then embark upon a fairly rigorous treatment of the simulation methods that are currently available for the treatment of stochastic biological models, with an eye toward comparing and contrasting their realms of applicability, and the care that must be taken when parameterizing them. Following that, we describe how stochasticity impacts several key biological functions, including transcription, translation, ribosome biogenesis, chromosome replication, and metabolism, before considering how the functions may be coupled into a comprehensive model of a ‘minimal cell’. Finally, we close with our expectation for the future of the field, focusing on how mesoscopic stochastic methods may be augmented with atomic-scale molecular modeling approaches in order to understand life across a range of length and time scales.

Stochastic ontogenetic growth model

Science.gov (United States)

West, B. J.; West, D.

2012-02-01

An ontogenetic growth model (OGM) for a thermodynamically closed system is generalized to satisfy both the first and second law of thermodynamics. The hypothesized stochastic ontogenetic growth model (SOGM) is shown to entail the interspecies allometry relation by explicitly averaging the basal metabolic rate and the total body mass over the steady-state probability density for the total body mass (TBM). This is the first derivation of the interspecies metabolic allometric relation from a dynamical model and the asymptotic steady-state distribution of the TBM is fit to data and shown to be inverse power law.

Multivariate moment closure techniques for stochastic kinetic models

International Nuclear Information System (INIS)

Lakatos, Eszter; Ale, Angelique; Kirk, Paul D. W.; Stumpf, Michael P. H.

2015-01-01

Stochastic effects dominate many chemical and biochemical processes. Their analysis, however, can be computationally prohibitively expensive and a range of approximation schemes have been proposed to lighten the computational burden. These, notably the increasingly popular linear noise approximation and the more general moment expansion methods, perform well for many dynamical regimes, especially linear systems. At higher levels of nonlinearity, it comes to an interplay between the nonlinearities and the stochastic dynamics, which is much harder to capture correctly by such approximations to the true stochastic processes. Moment-closure approaches promise to address this problem by capturing higher-order terms of the temporally evolving probability distribution. Here, we develop a set of multivariate moment-closures that allows us to describe the stochastic dynamics of nonlinear systems. Multivariate closure captures the way that correlations between different molecular species, induced by the reaction dynamics, interact with stochastic effects. We use multivariate Gaussian, gamma, and lognormal closure and illustrate their use in the context of two models that have proved challenging to the previous attempts at approximating stochastic dynamics: oscillations in p53 and Hes1. In addition, we consider a larger system, Erk-mediated mitogen-activated protein kinases signalling, where conventional stochastic simulation approaches incur unacceptably high computational costs

Multivariate moment closure techniques for stochastic kinetic models

Energy Technology Data Exchange (ETDEWEB)

Lakatos, Eszter, E-mail: e.lakatos13@imperial.ac.uk; Ale, Angelique; Kirk, Paul D. W.; Stumpf, Michael P. H., E-mail: m.stumpf@imperial.ac.uk [Department of Life Sciences, Centre for Integrative Systems Biology and Bioinformatics, Imperial College London, London SW7 2AZ (United Kingdom)

2015-09-07

Stochastic effects dominate many chemical and biochemical processes. Their analysis, however, can be computationally prohibitively expensive and a range of approximation schemes have been proposed to lighten the computational burden. These, notably the increasingly popular linear noise approximation and the more general moment expansion methods, perform well for many dynamical regimes, especially linear systems. At higher levels of nonlinearity, it comes to an interplay between the nonlinearities and the stochastic dynamics, which is much harder to capture correctly by such approximations to the true stochastic processes. Moment-closure approaches promise to address this problem by capturing higher-order terms of the temporally evolving probability distribution. Here, we develop a set of multivariate moment-closures that allows us to describe the stochastic dynamics of nonlinear systems. Multivariate closure captures the way that correlations between different molecular species, induced by the reaction dynamics, interact with stochastic effects. We use multivariate Gaussian, gamma, and lognormal closure and illustrate their use in the context of two models that have proved challenging to the previous attempts at approximating stochastic dynamics: oscillations in p53 and Hes1. In addition, we consider a larger system, Erk-mediated mitogen-activated protein kinases signalling, where conventional stochastic simulation approaches incur unacceptably high computational costs.

Compositional Modelling of Stochastic Hybrid Systems

NARCIS (Netherlands)

Strubbe, S.N.

2005-01-01

In this thesis we present a modelling framework for compositional modelling of stochastic hybrid systems. Hybrid systems consist of a combination of continuous and discrete dynamics. The state space of a hybrid system is hybrid in the sense that it consists of a continuous component and a discrete

Stochastic models for turbulent reacting flows

Energy Technology Data Exchange (ETDEWEB)

Kerstein, A. [Sandia National Laboratories, Livermore, CA (United States)

1993-12-01

The goal of this program is to develop and apply stochastic models of various processes occurring within turbulent reacting flows in order to identify the fundamental mechanisms governing these flows, to support experimental studies of these flows, and to further the development of comprehensive turbulent reacting flow models.

Multi-scenario modelling of uncertainty in stochastic chemical systems

International Nuclear Information System (INIS)

Evans, R. David; Ricardez-Sandoval, Luis A.

2014-01-01

Uncertainty analysis has not been well studied at the molecular scale, despite extensive knowledge of uncertainty in macroscale systems. The ability to predict the effect of uncertainty allows for robust control of small scale systems such as nanoreactors, surface reactions, and gene toggle switches. However, it is difficult to model uncertainty in such chemical systems as they are stochastic in nature, and require a large computational cost. To address this issue, a new model of uncertainty propagation in stochastic chemical systems, based on the Chemical Master Equation, is proposed in the present study. The uncertain solution is approximated by a composite state comprised of the averaged effect of samples from the uncertain parameter distributions. This model is then used to study the effect of uncertainty on an isomerization system and a two gene regulation network called a repressilator. The results of this model show that uncertainty in stochastic systems is dependent on both the uncertain distribution, and the system under investigation. -- Highlights: •A method to model uncertainty on stochastic systems was developed. •The method is based on the Chemical Master Equation. •Uncertainty in an isomerization reaction and a gene regulation network was modelled. •Effects were significant and dependent on the uncertain input and reaction system. •The model was computationally more efficient than Kinetic Monte Carlo

Identifiability in stochastic models

CERN Document Server

1992-01-01

The problem of identifiability is basic to all statistical methods and data analysis, occurring in such diverse areas as Reliability Theory, Survival Analysis, and Econometrics, where stochastic modeling is widely used. Mathematics dealing with identifiability per se is closely related to the so-called branch of ""characterization problems"" in Probability Theory. This book brings together relevant material on identifiability as it occurs in these diverse fields.

Study on individual stochastic model of GNSS observations for precise kinematic applications

Science.gov (United States)

Próchniewicz, Dominik; Szpunar, Ryszard

2015-04-01

The proper definition of mathematical positioning model, which is defined by functional and stochastic models, is a prerequisite to obtain the optimal estimation of unknown parameters. Especially important in this definition is realistic modelling of stochastic properties of observations, which are more receiver-dependent and time-varying than deterministic relationships. This is particularly true with respect to precise kinematic applications which are characterized by weakening model strength. In this case, incorrect or simplified definition of stochastic model causes that the performance of ambiguity resolution and accuracy of position estimation can be limited. In this study we investigate the methods of describing the measurement noise of GNSS observations and its impact to derive precise kinematic positioning model. In particular stochastic modelling of individual components of the variance-covariance matrix of observation noise performed using observations from a very short baseline and laboratory GNSS signal generator, is analyzed. Experimental test results indicate that the utilizing the individual stochastic model of observations including elevation dependency and cross-correlation instead of assumption that raw measurements are independent with the same variance improves the performance of ambiguity resolution as well as rover positioning accuracy. This shows that the proposed stochastic assessment method could be a important part in complex calibration procedure of GNSS equipment.

Stochastic Management of Non-Point Source Contamination: Joint Impact of Aquifer Heterogeneity and Well Characteristics

Science.gov (United States)

Henri, C. V.; Harter, T.

2017-12-01

Agricultural activities are recognized as the preeminent origin of non-point source (NPS) contamination of water bodies through the leakage of nitrate, salt and agrochemicals. A large fraction of world agricultural activities and therefore NPS contamination occurs over unconsolidated alluvial deposit basins offering soil composition and topography favorable to productive farming. These basins represent also important groundwater reservoirs. The over-exploitation of aquifers coupled with groundwater pollution by agriculture-related NPS contaminant has led to a rapid deterioration of the quality of these groundwater basins. The management of groundwater contamination from NPS is challenged by the inherent complexity of aquifers systems. Contaminant transport dynamics are highly uncertain due to the heterogeneity of hydraulic parameters controlling groundwater flow. Well characteristics are also key uncertain elements affecting pollutant transport and NPS management but quantifying uncertainty in NPS management under these conditions is not well documented. Our work focuses on better understanding the joint impact of aquifer heterogeneity and pumping well characteristics (extraction rate and depth) on (1) the transport of contaminants from NPS and (2) the spatio-temporal extension of the capture zone. To do so, we generate a series of geostatistically equivalent 3D heterogeneous aquifers and simulate the flow and non-reactive solute transport from NPS to extraction wells within a stochastic framework. The propagation of the uncertainty on the hydraulic conductivity field is systematically analyzed. A sensitivity analysis of the impact of extraction well characteristics (pumping rate and screen depth) is also conducted. Results highlight the significant role that heterogeneity and well characteristics plays on management metrics. We finally show that, in case of NPS contamination, the joint impact of regional longitudinal and transverse vertical hydraulic gradients and

Estimation of Stochastic Volatility Models by Nonparametric Filtering

DEFF Research Database (Denmark)

Kanaya, Shin; Kristensen, Dennis

2016-01-01

/estimated volatility process replacing the latent process. Our estimation strategy is applicable to both parametric and nonparametric stochastic volatility models, and can handle both jumps and market microstructure noise. The resulting estimators of the stochastic volatility model will carry additional biases...... and variances due to the first-step estimation, but under regularity conditions we show that these vanish asymptotically and our estimators inherit the asymptotic properties of the infeasible estimators based on observations of the volatility process. A simulation study examines the finite-sample properties...

A Simulation-Based Dynamic Stochastic Route Choice Model for Evacuation

Directory of Open Access Journals (Sweden)

Xing Zhao

2012-01-01

Full Text Available This paper establishes a dynamic stochastic route choice model for evacuation to simulate the propagation process of traffic flow and estimate the stochastic route choice under evacuation situations. The model contains a lane-group-based cell transmission model (CTM which sets different traffic capacities for links with different turning movements to flow out in an evacuation situation, an actual impedance model which is to obtain the impedance of each route in time units at each time interval and a stochastic route choice model according to the probit-based stochastic user equilibrium. In this model, vehicles loading at each origin at each time interval are assumed to choose an evacuation route under determinate road network, signal design, and OD demand. As a case study, the proposed model is validated on the network nearby Nanjing Olympic Center after the opening ceremony of the 10th National Games of the People's Republic of China. The traffic volumes and clearing time at five exit points of the evacuation zone are calculated by the model to compare with survey data. The results show that this model can appropriately simulate the dynamic route choice and evolution process of the traffic flow on the network in an evacuation situation.

On the small-time behavior of stochastic logistic models

Directory of Open Access Journals (Sweden)

Dung Tien Nguyen

2017-09-01

Full Text Available In this paper we investigate the small-time behaviors of the solution to  a stochastic logistic model. The obtained results allow us to estimate the number of individuals in the population and can be used to study stochastic prey-predator systems.

Combining deterministic and stochastic velocity fields in the analysis of deep crustal seismic data

Science.gov (United States)

Larkin, Steven Paul

Standard crustal seismic modeling obtains deterministic velocity models which ignore the effects of wavelength-scale heterogeneity, known to exist within the Earth's crust. Stochastic velocity models are a means to include wavelength-scale heterogeneity in the modeling. These models are defined by statistical parameters obtained from geologic maps of exposed crystalline rock, and are thus tied to actual geologic structures. Combining both deterministic and stochastic velocity models into a single model allows a realistic full wavefield (2-D) to be computed. By comparing these simulations to recorded seismic data, the effects of wavelength-scale heterogeneity can be investigated. Combined deterministic and stochastic velocity models are created for two datasets, the 1992 RISC seismic experiment in southeastern California and the 1986 PASSCAL seismic experiment in northern Nevada. The RISC experiment was located in the transition zone between the Salton Trough and the southern Basin and Range province. A high-velocity body previously identified beneath the Salton Trough is constrained to pinch out beneath the Chocolate Mountains to the northeast. The lateral extent of this body is evidence for the ephemeral nature of rifting loci as a continent is initially rifted. Stochastic modeling of wavelength-scale structures above this body indicate that little more than 5% mafic intrusion into a more felsic continental crust is responsible for the observed reflectivity. Modeling of the wide-angle RISC data indicates that coda waves following PmP are initially dominated by diffusion of energy out of the near-surface basin as the wavefield reverberates within this low-velocity layer. At later times, this coda consists of scattered body waves and P to S conversions. Surface waves do not play a significant role in this coda. Modeling of the PASSCAL dataset indicates that a high-gradient crust-mantle transition zone or a rough Moho interface is necessary to reduce precritical Pm

Spatial stochastic regression modelling of urban land use

International Nuclear Information System (INIS)

Arshad, S H M; Jaafar, J; Abiden, M Z Z; Latif, Z A; Rasam, A R A

2014-01-01

Urbanization is very closely linked to industrialization, commercialization or overall economic growth and development. This results in innumerable benefits of the quantity and quality of the urban environment and lifestyle but on the other hand contributes to unbounded development, urban sprawl, overcrowding and decreasing standard of living. Regulation and observation of urban development activities is crucial. The understanding of urban systems that promotes urban growth are also essential for the purpose of policy making, formulating development strategies as well as development plan preparation. This study aims to compare two different stochastic regression modeling techniques for spatial structure models of urban growth in the same specific study area. Both techniques will utilize the same datasets and their results will be analyzed. The work starts by producing an urban growth model by using stochastic regression modeling techniques namely the Ordinary Least Square (OLS) and Geographically Weighted Regression (GWR). The two techniques are compared to and it is found that, GWR seems to be a more significant stochastic regression model compared to OLS, it gives a smaller AICc (Akaike's Information Corrected Criterion) value and its output is more spatially explainable

Stochastic Geometric Models with Non-stationary Spatial Correlations in Lagrangian Fluid Flows

Science.gov (United States)

Gay-Balmaz, François; Holm, Darryl D.

2018-01-01

Inspired by spatiotemporal observations from satellites of the trajectories of objects drifting near the surface of the ocean in the National Oceanic and Atmospheric Administration's "Global Drifter Program", this paper develops data-driven stochastic models of geophysical fluid dynamics (GFD) with non-stationary spatial correlations representing the dynamical behaviour of oceanic currents. Three models are considered. Model 1 from Holm (Proc R Soc A 471:20140963, 2015) is reviewed, in which the spatial correlations are time independent. Two new models, called Model 2 and Model 3, introduce two different symmetry breaking mechanisms by which the spatial correlations may be advected by the flow. These models are derived using reduction by symmetry of stochastic variational principles, leading to stochastic Hamiltonian systems, whose momentum maps, conservation laws and Lie-Poisson bracket structures are used in developing the new stochastic Hamiltonian models of GFD.

Stochastic Geometric Models with Non-stationary Spatial Correlations in Lagrangian Fluid Flows

Science.gov (United States)

Gay-Balmaz, François; Holm, Darryl D.

2018-06-01

Inspired by spatiotemporal observations from satellites of the trajectories of objects drifting near the surface of the ocean in the National Oceanic and Atmospheric Administration's "Global Drifter Program", this paper develops data-driven stochastic models of geophysical fluid dynamics (GFD) with non-stationary spatial correlations representing the dynamical behaviour of oceanic currents. Three models are considered. Model 1 from Holm (Proc R Soc A 471:20140963, 2015) is reviewed, in which the spatial correlations are time independent. Two new models, called Model 2 and Model 3, introduce two different symmetry breaking mechanisms by which the spatial correlations may be advected by the flow. These models are derived using reduction by symmetry of stochastic variational principles, leading to stochastic Hamiltonian systems, whose momentum maps, conservation laws and Lie-Poisson bracket structures are used in developing the new stochastic Hamiltonian models of GFD.

Markov Chain Models for the Stochastic Modeling of Pitting Corrosion

Directory of Open Access Journals (Sweden)

A. Valor

2013-01-01

Full Text Available The stochastic nature of pitting corrosion of metallic structures has been widely recognized. It is assumed that this kind of deterioration retains no memory of the past, so only the current state of the damage influences its future development. This characteristic allows pitting corrosion to be categorized as a Markov process. In this paper, two different models of pitting corrosion, developed using Markov chains, are presented. Firstly, a continuous-time, nonhomogeneous linear growth (pure birth Markov process is used to model external pitting corrosion in underground pipelines. A closed-form solution of the system of Kolmogorov's forward equations is used to describe the transition probability function in a discrete pit depth space. The transition probability function is identified by correlating the stochastic pit depth mean with the empirical deterministic mean. In the second model, the distribution of maximum pit depths in a pitting experiment is successfully modeled after the combination of two stochastic processes: pit initiation and pit growth. Pit generation is modeled as a nonhomogeneous Poisson process, in which induction time is simulated as the realization of a Weibull process. Pit growth is simulated using a nonhomogeneous Markov process. An analytical solution of Kolmogorov's system of equations is also found for the transition probabilities from the first Markov state. Extreme value statistics is employed to find the distribution of maximum pit depths.

A stochastic model for the financial market with discontinuous prices

Directory of Open Access Journals (Sweden)

Leda D. Minkova

1996-01-01

Full Text Available This paper models some situations occurring in the financial market. The asset prices evolve according to a stochastic integral equation driven by a Gaussian martingale. A portfolio process is constrained in such a way that the wealth process covers some obligation. A solution to a linear stochastic integral equation is obtained in a class of cadlag stochastic processes.

Stochastic higher spin six vertex model and Macdonald measures

Science.gov (United States)

Borodin, Alexei

2018-02-01

We prove an identity that relates the q-Laplace transform of the height function of a (higher spin inhomogeneous) stochastic six vertex model in a quadrant on one side and a multiplicative functional of a Macdonald measure on the other. The identity is used to prove the GUE Tracy-Widom asymptotics for two instances of the stochastic six vertex model via asymptotic analysis of the corresponding Schur measures.

Integrated reservoir characterization: Improvement in heterogeneities stochastic modelling by integration of additional external constraints

Energy Technology Data Exchange (ETDEWEB)

Doligez, B.; Eschard, R. [Institut Francais du Petrole, Rueil Malmaison (France); Geffroy, F. [Centre de Geostatistique, Fontainebleau (France)] [and others

1997-08-01

The classical approach to construct reservoir models is to start with a fine scale geological model which is informed with petrophysical properties. Then scaling-up techniques allow to obtain a reservoir model which is compatible with the fluid flow simulators. Geostatistical modelling techniques are widely used to build the geological models before scaling-up. These methods provide equiprobable images of the area under investigation, which honor the well data, and which variability is the same than the variability computed from the data. At an appraisal phase, when few data are available, or when the wells are insufficient to describe all the heterogeneities and the behavior of the field, additional constraints are needed to obtain a more realistic geological model. For example, seismic data or stratigraphic models can provide average reservoir information with an excellent areal coverage, but with a poor vertical resolution. New advances in modelisation techniques allow now to integrate this type of additional external information in order to constrain the simulations. In particular, 2D or 3D seismic derived information grids, or sand-shale ratios maps coming from stratigraphic models can be used as external drifts to compute the geological image of the reservoir at the fine scale. Examples are presented to illustrate the use of these new tools, their impact on the final reservoir model, and their sensitivity to some key parameters.

Improved ensemble-mean forecast skills of ENSO events by a zero-mean stochastic model-error model of an intermediate coupled model

Science.gov (United States)

Zheng, F.; Zhu, J.

2015-12-01

To perform an ensemble-based ENSO probabilistic forecast, the crucial issue is to design a reliable ensemble prediction strategy that should include the major uncertainties of a forecast system. In this study, we developed a new general ensemble perturbation technique to improve the ensemble-mean predictive skill of forecasting ENSO using an intermediate coupled model (ICM). The model uncertainties are first estimated and analyzed from EnKF analysis results through assimilating observed SST. Then, based on the pre-analyzed properties of the model errors, a zero-mean stochastic model-error model is developed to mainly represent the model uncertainties induced by some important physical processes missed in the coupled model (i.e., stochastic atmospheric forcing/MJO, extra-tropical cooling and warming, Indian Ocean Dipole mode, etc.). Each member of an ensemble forecast is perturbed by the stochastic model-error model at each step during the 12-month forecast process, and the stochastical perturbations are added into the modeled physical fields to mimic the presence of these high-frequency stochastic noises and model biases and their effect on the predictability of the coupled system. The impacts of stochastic model-error perturbations on ENSO deterministic predictions are examined by performing two sets of 21-yr retrospective forecast experiments. The two forecast schemes are differentiated by whether they considered the model stochastic perturbations, with both initialized by the ensemble-mean analysis states from EnKF. The comparison results suggest that the stochastic model-error perturbations have significant and positive impacts on improving the ensemble-mean prediction skills during the entire 12-month forecast process. Because the nonlinear feature of the coupled model can induce the nonlinear growth of the added stochastic model errors with model integration, especially through the nonlinear heating mechanism with the vertical advection term of the model, the

Stochastic hyperfine interactions modeling library

Science.gov (United States)

Zacate, Matthew O.; Evenson, William E.

2011-04-01

The stochastic hyperfine interactions modeling library (SHIML) provides a set of routines to assist in the development and application of stochastic models of hyperfine interactions. The library provides routines written in the C programming language that (1) read a text description of a model for fluctuating hyperfine fields, (2) set up the Blume matrix, upon which the evolution operator of the system depends, and (3) find the eigenvalues and eigenvectors of the Blume matrix so that theoretical spectra of experimental techniques that measure hyperfine interactions can be calculated. The optimized vector and matrix operations of the BLAS and LAPACK libraries are utilized; however, there was a need to develop supplementary code to find an orthonormal set of (left and right) eigenvectors of complex, non-Hermitian matrices. In addition, example code is provided to illustrate the use of SHIML to generate perturbed angular correlation spectra for the special case of polycrystalline samples when anisotropy terms of higher order than A can be neglected. Program summaryProgram title: SHIML Catalogue identifier: AEIF_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEIF_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU GPL 3 No. of lines in distributed program, including test data, etc.: 8224 No. of bytes in distributed program, including test data, etc.: 312 348 Distribution format: tar.gz Programming language: C Computer: Any Operating system: LINUX, OS X RAM: Varies Classification: 7.4 External routines: TAPP [1], BLAS [2], a C-interface to BLAS [3], and LAPACK [4] Nature of problem: In condensed matter systems, hyperfine methods such as nuclear magnetic resonance (NMR), Mössbauer effect (ME), muon spin rotation (μSR), and perturbed angular correlation spectroscopy (PAC) measure electronic and magnetic structure within Angstroms of nuclear probes through the hyperfine interaction. When

Scalable inference for stochastic block models

KAUST Repository

Peng, Chengbin; Zhang, Zhihua; Wong, Ka-Chun; Zhang, Xiangliang; Keyes, David E.

2017-01-01

Community detection in graphs is widely used in social and biological networks, and the stochastic block model is a powerful probabilistic tool for describing graphs with community structures. However, in the era of "big data," traditional inference

Stochastic Parametrisations and Regime Behaviour of Atmospheric Models

Science.gov (United States)

Arnold, Hannah; Moroz, Irene; Palmer, Tim

2013-04-01

The presence of regimes is a characteristic of non-linear, chaotic systems (Lorenz, 2006). In the atmosphere, regimes emerge as familiar circulation patterns such as the El-Nino Southern Oscillation (ENSO), the North Atlantic Oscillation (NAO) and Scandinavian Blocking events. In recent years there has been much interest in the problem of identifying and studying atmospheric regimes (Solomon et al, 2007). In particular, how do these regimes respond to an external forcing such as anthropogenic greenhouse gas emissions? The importance of regimes in observed trends over the past 50-100 years indicates that in order to predict anthropogenic climate change, our climate models must be able to represent accurately natural circulation regimes, their statistics and variability. It is well established that representing model uncertainty as well as initial condition uncertainty is important for reliable weather forecasts (Palmer, 2001). In particular, stochastic parametrisation schemes have been shown to improve the skill of weather forecast models (e.g. Berner et al., 2009; Frenkel et al., 2012; Palmer et al., 2009). It is possible that including stochastic physics as a representation of model uncertainty could also be beneficial in climate modelling, enabling the simulator to explore larger regions of the climate attractor including other flow regimes. An alternative representation of model uncertainty is a perturbed parameter scheme, whereby physical parameters in subgrid parametrisation schemes are perturbed about their optimal value. Perturbing parameters gives a greater control over the ensemble than multi-model or multiparametrisation ensembles, and has been used as a representation of model uncertainty in climate prediction (Stainforth et al., 2005; Rougier et al., 2009). We investigate the effect of including representations of model uncertainty on the regime behaviour of a simulator. A simple chaotic model of the atmosphere, the Lorenz '96 system, is used to study

Model tracking dual stochastic controller design under irregular internal noises

International Nuclear Information System (INIS)

Lee, Jong Bok; Heo, Hoon; Cho, Yun Hyun; Ji, Tae Young

2006-01-01

Although many methods about the control of irregular external noise have been introduced and implemented, it is still necessary to design a controller that will be more effective and efficient methods to exclude for various noises. Accumulation of errors due to model tracking, internal noises (thermal noise, shot noise and l/f noise) that come from elements such as resistor, diode and transistor etc. in the circuit system and numerical errors due to digital process often destabilize the system and reduce the system performance. New stochastic controller is adopted to remove those noises using conventional controller simultaneously. Design method of a model tracking dual controller is proposed to improve the stability of system while removing external and internal noises. In the study, design process of the model tracking dual stochastic controller is introduced that improves system performance and guarantees robustness under irregular internal noises which can be created internally. The model tracking dual stochastic controller utilizing F-P-K stochastic control technique developed earlier is implemented to reveal its performance via simulation

Simulation of nuclear plant operation into a stochastic energy production model

International Nuclear Information System (INIS)

Pacheco, R.L.

1983-04-01

A simulation model of nuclear plant operation is developed to fit into a stochastic energy production model. In order to improve the stochastic model used, and also reduce its computational time burdened by the aggregation of the model of nuclear plant operation, a study of tail truncation of the unsupplied demand distribution function has been performed. (E.G.) [pt

A complementarity model for solving stochastic natural gas market equilibria

International Nuclear Information System (INIS)

Jifang Zhuang; Gabriel, S.A.

2008-01-01

This paper presents a stochastic equilibrium model for deregulated natural gas markets. Each market participant (pipeline operators, producers, etc.) solves a stochastic optimization problem whose optimality conditions, when combined with market-clearing conditions give rise to a certain mixed complementarity problem (MiCP). The stochastic aspects are depicted by a recourse problem for each player in which the first-stage decisions relate to long-term contracts and the second-stage decisions relate to spot market activities for three seasons. Besides showing that such a market model is an instance of a MiCP, we provide theoretical results concerning long-term and spot market prices and solve the resulting MiCP for a small yet representative market. We also note an interesting observation for the value of the stochastic solution for non-optimization problems. (author)

A complementarity model for solving stochastic natural gas market equilibria

International Nuclear Information System (INIS)

Zhuang Jifang; Gabriel, Steven A.

2008-01-01

This paper presents a stochastic equilibrium model for deregulated natural gas markets. Each market participant (pipeline operators, producers, etc.) solves a stochastic optimization problem whose optimality conditions, when combined with market-clearing conditions give rise to a certain mixed complementarity problem (MiCP). The stochastic aspects are depicted by a recourse problem for each player in which the first-stage decisions relate to long-term contracts and the second-stage decisions relate to spot market activities for three seasons. Besides showing that such a market model is an instance of a MiCP, we provide theoretical results concerning long-term and spot market prices and solve the resulting MiCP for a small yet representative market. We also note an interesting observation for the value of the stochastic solution for non-optimization problems

Pseudodynamic Source Characterization for Strike-Slip Faulting Including Stress Heterogeneity and Super-Shear Ruptures

KAUST Repository

Mena, B.; Dalguer, L. A.; Mai, Paul Martin

2012-01-01

. (2004), we propose new relationships for PD models for moderate‐to‐large strike‐slip earthquakes that include local supershear rupture speed due to stress heterogeneities. We conduct dynamic rupture simulations using stochastic initial stress

Deterministic and stochastic models for middle east respiratory syndrome (MERS)

Science.gov (United States)

Suryani, Dessy Rizki; Zevika, Mona; Nuraini, Nuning

2018-03-01

World Health Organization (WHO) data stated that since September 2012, there were 1,733 cases of Middle East Respiratory Syndrome (MERS) with 628 death cases that occurred in 27 countries. MERS was first identified in Saudi Arabia in 2012 and the largest cases of MERS outside Saudi Arabia occurred in South Korea in 2015. MERS is a disease that attacks the respiratory system caused by infection of MERS-CoV. MERS-CoV transmission occurs directly through direct contact between infected individual with non-infected individual or indirectly through contaminated object by the free virus. Suspected, MERS can spread quickly because of the free virus in environment. Mathematical modeling is used to illustrate the transmission of MERS disease using deterministic model and stochastic model. Deterministic model is used to investigate the temporal dynamic from the system to analyze the steady state condition. Stochastic model approach using Continuous Time Markov Chain (CTMC) is used to predict the future states by using random variables. From the models that were built, the threshold value for deterministic models and stochastic models obtained in the same form and the probability of disease extinction can be computed by stochastic model. Simulations for both models using several of different parameters are shown, and the probability of disease extinction will be compared with several initial conditions.

Considerations when ranking stochastically modeled oil sands resource models for mining applications

Energy Technology Data Exchange (ETDEWEB)

Etris, E.L. [Society of Petroleum Engineers, Canadian Section, Calgary, AB (Canada)]|[Petro-Canada, Calgary, AB (Canada); Idris, Y.; Hunter, A.C. [Petro-Canada, Calgary, AB (Canada)

2008-10-15

Alberta's Athabasca oil sands deposit has been targeted as a major resource for development. Bitumen recovery operations fall into 2 categories, namely mining and in situ operations. Mining recovery is done above ground level and consists of open pit digging, disaggregation of the bitumen-saturated sediment through crushing followed by pipeline transport in a water-based slurry and then separation of oil, water and sediment. In situ recovery consists of drilling wells and stimulating the oil sands in the subsurface with a thermal treatment to reduce the viscosity of the bitumen and allow it to come to the surface. Steam assisted gravity drainage (SAGD) is the most popular thermal treatment currently in use. Resource models that simulate the recovery process are needed for both mining and in situ recovery operations. Both types can benefit from the advantages of a stochastic modeling process for resource model building and uncertainty evaluation. Stochastic modeling provides a realistic geology and allows for multiple realizations, which mining operations can use to evaluate the variability of recoverable bitumen volumes and develop mine plans accordingly. This paper described the processes of stochastic modelling and of determining the appropriate single realization for mine planning as applied to the Fort Hills oil sands mine which is currently in the early planning stage. The modeling exercise was used to estimate the in-place resource and quantify the uncertainty in resource volumes. The stochastic models were checked against those generated from conventional methods to identify any differences and to make the appropriate adaptations. 13 refs., 3 tabs., 16 figs.

A computer model of the biosphere, to estimate stochastic and non-stochastic effects of radionuclides on humans

International Nuclear Information System (INIS)

Laurens, J.M.

1985-01-01

A computer code was written to model food chains in order to estimate the internal and external doses, for stochastic and non-stochastic effects, on humans (adults and infants). Results are given for 67 radionuclides, for unit concentration in water (1 Bq/L) and in atmosphere (1 Bq/m 3 )

A coupled stochastic rainfall-evapotranspiration model for hydrological impact analysis

Science.gov (United States)

Pham, Minh Tu; Vernieuwe, Hilde; De Baets, Bernard; Verhoest, Niko E. C.

2018-02-01

A hydrological impact analysis concerns the study of the consequences of certain scenarios on one or more variables or fluxes in the hydrological cycle. In such an exercise, discharge is often considered, as floods originating from extremely high discharges often cause damage. Investigating the impact of extreme discharges generally requires long time series of precipitation and evapotranspiration to be used to force a rainfall-runoff model. However, such kinds of data may not be available and one should resort to stochastically generated time series, even though the impact of using such data on the overall discharge, and especially on the extreme discharge events, is not well studied. In this paper, stochastically generated rainfall and corresponding evapotranspiration time series, generated by means of vine copulas, are used to force a simple conceptual hydrological model. The results obtained are comparable to the modelled discharge using observed forcing data. Yet, uncertainties in the modelled discharge increase with an increasing number of stochastically generated time series used. Notwithstanding this finding, it can be concluded that using a coupled stochastic rainfall-evapotranspiration model has great potential for hydrological impact analysis.

Stochastic Modeling and Analysis of Power System with Renewable Generation

DEFF Research Database (Denmark)

Chen, Peiyuan

Unlike traditional fossil-fuel based power generation, renewable generation such as wind power relies on uncontrollable prime sources such as wind speed. Wind speed varies stochastically, which to a large extent determines the stochastic behavior of power generation from wind farms...... that such a stochastic model can be used to simulate the effect of load management on the load duration curve. As CHP units are turned on and off by regulating power, CHP generation has discrete output and thus can be modeled by a transition matrix based discrete Markov chain. As the CHP generation has a strong diurnal...

Extinction in neutrally stable stochastic Lotka-Volterra models

Science.gov (United States)

Dobrinevski, Alexander; Frey, Erwin

2012-05-01

Populations of competing biological species exhibit a fascinating interplay between the nonlinear dynamics of evolutionary selection forces and random fluctuations arising from the stochastic nature of the interactions. The processes leading to extinction of species, whose understanding is a key component in the study of evolution and biodiversity, are influenced by both of these factors. Here, we investigate a class of stochastic population dynamics models based on generalized Lotka-Volterra systems. In the case of neutral stability of the underlying deterministic model, the impact of intrinsic noise on the survival of species is dramatic: It destroys coexistence of interacting species on a time scale proportional to the population size. We introduce a new method based on stochastic averaging which allows one to understand this extinction process quantitatively by reduction to a lower-dimensional effective dynamics. This is performed analytically for two highly symmetrical models and can be generalized numerically to more complex situations. The extinction probability distributions and other quantities of interest we obtain show excellent agreement with simulations.

Model predictive control classical, robust and stochastic

CERN Document Server

Kouvaritakis, Basil

2016-01-01

For the first time, a textbook that brings together classical predictive control with treatment of up-to-date robust and stochastic techniques. Model Predictive Control describes the development of tractable algorithms for uncertain, stochastic, constrained systems. The starting point is classical predictive control and the appropriate formulation of performance objectives and constraints to provide guarantees of closed-loop stability and performance. Moving on to robust predictive control, the text explains how similar guarantees may be obtained for cases in which the model describing the system dynamics is subject to additive disturbances and parametric uncertainties. Open- and closed-loop optimization are considered and the state of the art in computationally tractable methods based on uncertainty tubes presented for systems with additive model uncertainty. Finally, the tube framework is also applied to model predictive control problems involving hard or probabilistic constraints for the cases of multiplic...

A low-bias simulation scheme for the SABR stochastic volatility model

NARCIS (Netherlands)

B. Chen (Bin); C.W. Oosterlee (Cornelis); J.A.M. van der Weide

2012-01-01

htmlabstractThe Stochastic Alpha Beta Rho Stochastic Volatility (SABR-SV) model is widely used in the financial industry for the pricing of fixed income instruments. In this paper we develop an lowbias simulation scheme for the SABR-SV model, which deals efficiently with (undesired)

Stochastic mixed-mode oscillations in a three-species predator-prey model

Science.gov (United States)

Sadhu, Susmita; Kuehn, Christian

2018-03-01

The effect of demographic stochasticity, in the form of Gaussian white noise, in a predator-prey model with one fast and two slow variables is studied. We derive the stochastic differential equations (SDEs) from a discrete model. For suitable parameter values, the deterministic drift part of the model admits a folded node singularity and exhibits a singular Hopf bifurcation. We focus on the parameter regime near the Hopf bifurcation, where small amplitude oscillations exist as stable dynamics in the absence of noise. In this regime, the stochastic model admits noise-driven mixed-mode oscillations (MMOs), which capture the intermediate dynamics between two cycles of population outbreaks. We perform numerical simulations to calculate the distribution of the random number of small oscillations between successive spikes for varying noise intensities and distance to the Hopf bifurcation. We also study the effect of noise on a suitable Poincaré map. Finally, we prove that the stochastic model can be transformed into a normal form near the folded node, which can be linked to recent results on the interplay between deterministic and stochastic small amplitude oscillations. The normal form can also be used to study the parameter influence on the noise level near folded singularities.

A Stochastic Model of the Yeast Cell Cycle Reveals Roles for Feedback Regulation in Limiting Cellular Variability.

Science.gov (United States)

Barik, Debashis; Ball, David A; Peccoud, Jean; Tyson, John J

2016-12-01

The cell division cycle of eukaryotes is governed by a complex network of cyclin-dependent protein kinases (CDKs) and auxiliary proteins that govern CDK activities. The control system must function reliably in the context of molecular noise that is inevitable in tiny yeast cells, because mistakes in sequencing cell cycle events are detrimental or fatal to the cell or its progeny. To assess the effects of noise on cell cycle progression requires not only extensive, quantitative, experimental measurements of cellular heterogeneity but also comprehensive, accurate, mathematical models of stochastic fluctuations in the CDK control system. In this paper we provide a stochastic model of the budding yeast cell cycle that accurately accounts for the variable phenotypes of wild-type cells and more than 20 mutant yeast strains simulated in different growth conditions. We specifically tested the role of feedback regulations mediated by G1- and SG2M-phase cyclins to minimize the noise in cell cycle progression. Details of the model are informed and tested by quantitative measurements (by fluorescence in situ hybridization) of the joint distributions of mRNA populations in yeast cells. We use the model to predict the phenotypes of ~30 mutant yeast strains that have not yet been characterized experimentally.

A Stochastic Model of the Yeast Cell Cycle Reveals Roles for Feedback Regulation in Limiting Cellular Variability.

Directory of Open Access Journals (Sweden)

Debashis Barik

2016-12-01

Full Text Available The cell division cycle of eukaryotes is governed by a complex network of cyclin-dependent protein kinases (CDKs and auxiliary proteins that govern CDK activities. The control system must function reliably in the context of molecular noise that is inevitable in tiny yeast cells, because mistakes in sequencing cell cycle events are detrimental or fatal to the cell or its progeny. To assess the effects of noise on cell cycle progression requires not only extensive, quantitative, experimental measurements of cellular heterogeneity but also comprehensive, accurate, mathematical models of stochastic fluctuations in the CDK control system. In this paper we provide a stochastic model of the budding yeast cell cycle that accurately accounts for the variable phenotypes of wild-type cells and more than 20 mutant yeast strains simulated in different growth conditions. We specifically tested the role of feedback regulations mediated by G1- and SG2M-phase cyclins to minimize the noise in cell cycle progression. Details of the model are informed and tested by quantitative measurements (by fluorescence in situ hybridization of the joint distributions of mRNA populations in yeast cells. We use the model to predict the phenotypes of ~30 mutant yeast strains that have not yet been characterized experimentally.

Stochastic resonance in models of neuronal ensembles

International Nuclear Information System (INIS)

Chialvo, D.R.; Longtin, A.; Mueller-Gerkin, J.

1997-01-01

Two recently suggested mechanisms for the neuronal encoding of sensory information involving the effect of stochastic resonance with aperiodic time-varying inputs are considered. It is shown, using theoretical arguments and numerical simulations, that the nonmonotonic behavior with increasing noise of the correlation measures used for the so-called aperiodic stochastic resonance (ASR) scenario does not rely on the cooperative effect typical of stochastic resonance in bistable and excitable systems. Rather, ASR with slowly varying signals is more properly interpreted as linearization by noise. Consequently, the broadening of the open-quotes resonance curveclose quotes in the multineuron stochastic resonance without tuning scenario can also be explained by this linearization. Computation of the input-output correlation as a function of both signal frequency and noise for the model system further reveals conditions where noise-induced firing with aperiodic inputs will benefit from stochastic resonance rather than linearization by noise. Thus, our study clarifies the tuning requirements for the optimal transduction of subthreshold aperiodic signals. It also shows that a single deterministic neuron can perform as well as a network when biased into a suprathreshold regime. Finally, we show that the inclusion of a refractory period in the spike-detection scheme produces a better correlation between instantaneous firing rate and input signal. copyright 1997 The American Physical Society

A Rigorous Temperature-Dependent Stochastic Modelling and Testing for MEMS-Based Inertial Sensor Errors

Directory of Open Access Journals (Sweden)

Spiros Pagiatakis

2009-10-01

Full Text Available In this paper, we examine the effect of changing the temperature points on MEMS-based inertial sensor random error. We collect static data under different temperature points using a MEMS-based inertial sensor mounted inside a thermal chamber. Rigorous stochastic models, namely Autoregressive-based Gauss-Markov (AR-based GM models are developed to describe the random error behaviour. The proposed AR-based GM model is initially applied to short stationary inertial data to develop the stochastic model parameters (correlation times. It is shown that the stochastic model parameters of a MEMS-based inertial unit, namely the ADIS16364, are temperature dependent. In addition, field kinematic test data collected at about 17 °C are used to test the performance of the stochastic models at different temperature points in the filtering stage using Unscented Kalman Filter (UKF. It is shown that the stochastic model developed at 20 °C provides a more accurate inertial navigation solution than the ones obtained from the stochastic models developed at −40 °C, −20 °C, 0 °C, +40 °C, and +60 °C. The temperature dependence of the stochastic model is significant and should be considered at all times to obtain optimal navigation solution for MEMS-based INS/GPS integration.

A Rigorous Temperature-Dependent Stochastic Modelling and Testing for MEMS-Based Inertial Sensor Errors.

Science.gov (United States)

El-Diasty, Mohammed; Pagiatakis, Spiros

2009-01-01

In this paper, we examine the effect of changing the temperature points on MEMS-based inertial sensor random error. We collect static data under different temperature points using a MEMS-based inertial sensor mounted inside a thermal chamber. Rigorous stochastic models, namely Autoregressive-based Gauss-Markov (AR-based GM) models are developed to describe the random error behaviour. The proposed AR-based GM model is initially applied to short stationary inertial data to develop the stochastic model parameters (correlation times). It is shown that the stochastic model parameters of a MEMS-based inertial unit, namely the ADIS16364, are temperature dependent. In addition, field kinematic test data collected at about 17 °C are used to test the performance of the stochastic models at different temperature points in the filtering stage using Unscented Kalman Filter (UKF). It is shown that the stochastic model developed at 20 °C provides a more accurate inertial navigation solution than the ones obtained from the stochastic models developed at -40 °C, -20 °C, 0 °C, +40 °C, and +60 °C. The temperature dependence of the stochastic model is significant and should be considered at all times to obtain optimal navigation solution for MEMS-based INS/GPS integration.

Stochastic volatility and multi-dimensional modeling in the European energy market

Energy Technology Data Exchange (ETDEWEB)

Vos, Linda

2012-07-01

In energy prices there is evidence for stochastic volatility. Stochastic volatility has effect on the price of path-dependent options and therefore has to be modeled properly. We introduced a multi-dimensional non-Gaussian stochastic volatility model with leverage which can be used in energy pricing. It captures special features of energy prices like price spikes, mean-reversion, stochastic volatility and inverse leverage. Moreover it allows modeling dependencies between different commodities.The derived forward price dynamics based on this multi-variate spot price model, provides a very flexible structure. It includes cotango, backwardation and hump shape forward curves.Alternatively energy prices could be modeled by a 2-factor model consisting of a non-Gaussian stable CARMA process and a non-stationary trend models by a Levy process. Also this model is able to capture special features like price spikes, mean reversion and the low frequency dynamics in the market. An robust L1-filter is introduced to filter out the states of the CARMA process. When applying to German electricity EEX exchange data an overall negative risk-premium is found. However close to delivery a positive risk-premium is observed.(Author)

Stochastic Load Models and Footbridge Response

DEFF Research Database (Denmark)

Pedersen, Lars; Frier, Christian

2015-01-01

Pedestrians may cause vibrations in footbridges and these vibrations may potentially be annoying. This calls for predictions of footbridge vibration levels and the paper considers a stochastic approach to modeling the action of pedestrians assuming walking parameters such as step frequency, pedes...

On cross-currency models with stochastic volatility and correlated interest rates

NARCIS (Netherlands)

Grzelak, L.A.; Oosterlee, C.W.

2010-01-01

We construct multi-currency models with stochastic volatility and correlated stochastic interest rates with a full matrix of correlations. We first deal with a foreign exchange (FX) model of Heston-type, in which the domestic and foreign interest rates are generated by the short-rate process of

On the Realistic Stochastic Model of GPS Observables: Implementation and Performance

Science.gov (United States)

Zangeneh-Nejad, F.; Amiri-Simkooei, A. R.; Sharifi, M. A.; Asgari, J.

2015-12-01

High-precision GPS positioning requires a realistic stochastic model of observables. A realistic GPS stochastic model of observables should take into account different variances for different observation types, correlations among different observables, the satellite elevation dependence of observables precision, and the temporal correlation of observables. Least-squares variance component estimation (LS-VCE) is applied to GPS observables using the geometry-based observation model (GBOM). To model the satellite elevation dependent of GPS observables precision, an exponential model depending on the elevation angles of the satellites are also employed. Temporal correlation of the GPS observables is modelled by using a first-order autoregressive noise model. An important step in the high-precision GPS positioning is double difference integer ambiguity resolution (IAR). The fraction or percentage of success among a number of integer ambiguity fixing is called the success rate. A realistic estimation of the GNSS observables covariance matrix plays an important role in the IAR. We consider the ambiguity resolution success rate for two cases, namely a nominal and a realistic stochastic model of the GPS observables using two GPS data sets collected by the Trimble R8 receiver. The results confirm that applying a more realistic stochastic model can significantly improve the IAR success rate on individual frequencies, either on L1 or on L2. An improvement of 20% was achieved to the empirical success rate results. The results also indicate that introducing the realistic stochastic model leads to a larger standard deviation for the baseline components by a factor of about 2.6 on the data sets considered.

Inexact Multistage Stochastic Chance Constrained Programming Model for Water Resources Management under Uncertainties

Directory of Open Access Journals (Sweden)

Hong Zhang

2017-01-01

Full Text Available In order to formulate water allocation schemes under uncertainties in the water resources management systems, an inexact multistage stochastic chance constrained programming (IMSCCP model is proposed. The model integrates stochastic chance constrained programming, multistage stochastic programming, and inexact stochastic programming within a general optimization framework to handle the uncertainties occurring in both constraints and objective. These uncertainties are expressed as probability distributions, interval with multiply distributed stochastic boundaries, dynamic features of the long-term water allocation plans, and so on. Compared with the existing inexact multistage stochastic programming, the IMSCCP can be used to assess more system risks and handle more complicated uncertainties in water resources management systems. The IMSCCP model is applied to a hypothetical case study of water resources management. In order to construct an approximate solution for the model, a hybrid algorithm, which incorporates stochastic simulation, back propagation neural network, and genetic algorithm, is proposed. The results show that the optimal value represents the maximal net system benefit achieved with a given confidence level under chance constraints, and the solutions provide optimal water allocation schemes to multiple users over a multiperiod planning horizon.

The Role of Stochastic Models in Interpreting the Origins of Biological Chirality

Directory of Open Access Journals (Sweden)

Gábor Lente

2010-04-01

Full Text Available This review summarizes recent stochastic modeling efforts in the theoretical research aimed at interpreting the origins of biological chirality. Stochastic kinetic models, especially those based on the continuous time discrete state approach, have great potential in modeling absolute asymmetric reactions, experimental examples of which have been reported in the past decade. An overview of the relevant mathematical background is given and several examples are presented to show how the significant numerical problems characteristic of the use of stochastic models can be overcome by non-trivial, but elementary algebra. In these stochastic models, a particulate view of matter is used rather than the concentration-based view of traditional chemical kinetics using continuous functions to describe the properties system. This has the advantage of giving adequate description of single-molecule events, which were probably important in the origin of biological chirality. The presented models can interpret and predict the random distribution of enantiomeric excess among repetitive experiments, which is the most striking feature of absolute asymmetric reactions. It is argued that the use of the stochastic kinetic approach should be much more widespread in the relevant literature.

Moment Closure for the Stochastic Logistic Model

National Research Council Canada - National Science Library

Singh, Abhyudai; Hespanha, Joao P

2006-01-01

..., which we refer to as the moment closure function. In this paper, a systematic procedure for constructing moment closure functions of arbitrary order is presented for the stochastic logistic model...

Stochastic persistence and stationary distribution in an SIS epidemic model with media coverage

Science.gov (United States)

Guo, Wenjuan; Cai, Yongli; Zhang, Qimin; Wang, Weiming

2018-02-01

This paper aims to study an SIS epidemic model with media coverage from a general deterministic model to a stochastic differential equation with environment fluctuation. Mathematically, we use the Markov semigroup theory to prove that the basic reproduction number R0s can be used to control the dynamics of stochastic system. Epidemiologically, we show that environment fluctuation can inhibit the occurrence of the disease, namely, in the case of disease persistence for the deterministic model, the disease still dies out with probability one for the stochastic model. So to a great extent the stochastic perturbation under media coverage affects the outbreak of the disease.

Effects of demographic structure on key properties of stochastic density-independent population dynamics.

Science.gov (United States)

Vindenes, Yngvild; Sæther, Bernt-Erik; Engen, Steinar

2012-12-01

The development of stochastic demography has largely been based on age structured populations, although other types of demographic structure, especially permanent and dynamic heterogeneity, are likely common in natural populations. The combination of stochasticity and demographic structure is a challenge for analyses of population dynamics and extinction risk, because the population structure will fluctuate around the stable structure and the population size shows transient fluctuations. However, by using a diffusion approximation for the total reproductive value, density-independent dynamics of structured populations can be described with only three population parameters: the expected population growth rate, the environmental variance and the demographic variance. These parameters depend on population structure via the state-specific vital rates and transition rates. Once they are found, the diffusion approximation represents a substantial reduction in model complexity. Here, we review and compare the key population parameters across a wide range of demographic structure, from the case of no structure to the most general case of dynamic heterogeneity, and for both discrete and continuous types. We focus on the demographic variance, but also show how environmental stochasticity can be included. This study brings together results from recent models, each considering a specific type of population structure, and places them in a general framework for structured populations. Comparison across different types of demographic structure reveals that the reproductive value is an essential concept for understanding how population structure affects stochastic dynamics and extinction risk. Copyright © 2011 Elsevier Inc. All rights reserved.

Persistence and extinction for a stochastic logistic model with infinite delay

Directory of Open Access Journals (Sweden)

Chun Lu

2013-11-01

Full Text Available This article, studies a stochastic logistic model with infinite delay. Using a phase space, we establish sufficient conditions for the extinction, nonpersistence in the mean, weak persistence, and stochastic permanence. A threshold between weak persistence and extinction is obtained. Our results state that different types of environmental noises have different effects on the persistence and extinction, and that the delay has no impact on the persistence and extinction for the stochastic model in the autonomous case. Numerical simulations illustrate the theoretical results.

Up scaling two-phase flow in heterogeneous porous media; Mise a l'echelle des ecoulements diphasiques dans les milieux poreux heterogenes

Energy Technology Data Exchange (ETDEWEB)

Artus, V.

2003-11-01

For two-phase flow in heterogeneous media, the emergence of different flow regimes at large-scale is driven by local interactions between the viscous coupling and the heterogeneity. In particular, when the viscosity ratio is favorable, viscous effects induce a transverse flow that stabilizes the front while flooding. However, most of recent stochastic models neglect the influence of the viscous coupling. We developed a stochastic model for the dynamics of the front, taking the viscous coupling into account. For stable cases, this model relates the statistical properties of the front to the statistical properties of the permeability field. For stable flow in stratified media, we show that the front is stationary by parts in the reservoir. These parts can be identified as large-scale hydrodynamic layers and separately coarsened in the large-scale simulation model. For flows with favorable viscosity ratios in isotropic reservoirs, we show that a stationary front occurs, in a statistical sense. For unfavorable viscosity ratios, the flow is driven by the development of viscous fingering. These different regimes lead to different large-scale saturation profiles that can be matched with a macro-dispersion equation, if the effective convective flux is modified to take into account stabilizing or destabilizing viscous effects. (author)

Some recent developments in stochastic volatility modelling

DEFF Research Database (Denmark)

Barndorff-Nielsen, Ole Eiler; Nicolato, Elisa; Shephard, N.

2002-01-01

This paper reviews and puts in context some of our recent work on stochastic volatility (SV) modelling for financial economics. Here our main focus is on: (i) the relationship between subordination and SV, (ii) OU based volatility models, (iii) exact option pricing, (iv) realized power variation...

Scalable inference for stochastic block models

KAUST Repository

Peng, Chengbin

2017-12-08

Community detection in graphs is widely used in social and biological networks, and the stochastic block model is a powerful probabilistic tool for describing graphs with community structures. However, in the era of "big data," traditional inference algorithms for such a model are increasingly limited due to their high time complexity and poor scalability. In this paper, we propose a multi-stage maximum likelihood approach to recover the latent parameters of the stochastic block model, in time linear with respect to the number of edges. We also propose a parallel algorithm based on message passing. Our algorithm can overlap communication and computation, providing speedup without compromising accuracy as the number of processors grows. For example, to process a real-world graph with about 1.3 million nodes and 10 million edges, our algorithm requires about 6 seconds on 64 cores of a contemporary commodity Linux cluster. Experiments demonstrate that the algorithm can produce high quality results on both benchmark and real-world graphs. An example of finding more meaningful communities is illustrated consequently in comparison with a popular modularity maximization algorithm.

Fiber Bundle Model Under Heterogeneous Loading

Science.gov (United States)

Roy, Subhadeep; Goswami, Sanchari

2018-03-01

The present work deals with the behavior of fiber bundle model under heterogeneous loading condition. The model is explored both in the mean-field limit as well as with local stress concentration. In the mean field limit, the failure abruptness decreases with increasing order k of heterogeneous loading. In this limit, a brittle to quasi-brittle transition is observed at a particular strength of disorder which changes with k. On the other hand, the model is hardly affected by such heterogeneity in the limit where local stress concentration plays a crucial role. The continuous limit of the heterogeneous loading is also studied and discussed in this paper. Some of the important results related to fiber bundle model are reviewed and their responses to our new scheme of heterogeneous loading are studied in details. Our findings are universal with respect to the nature of the threshold distribution adopted to assign strength to an individual fiber.

A stochastic analysis for a phytoplankton-zooplankton model

International Nuclear Information System (INIS)

Ge, G; Wang, H-L; Xu, J

2008-01-01

A simple phytoplankton-zooplankton nonlinear dynamical model was proposed to study the coexistence of all the species and a Hopf bifurcation was observed. In order to study the effect of environmental robustness on this system, we have stochastically perturbed the system with respect to white noise around its positive interior equilibrium. We have observed that the system remains stochastically stable around the positive equilibrium for same parametric values in the deterministic situation

A hybrid small-world network/semi-physical model for predicting wildfire spread in heterogeneous landscapes

International Nuclear Information System (INIS)

Billaud, Y; Kaiss, A; Drissi, M; Pizzo, Y; Porterie, B; Zekri, N; Acem, Z; Collin, A; Boulet, P; Santoni, P-A; Bosseur, F

2012-01-01

This paper presents the latest developments and validation results of a hybrid model which combines a broad-scale stochastic small-world network model with a macroscopic deterministic approach, to simulate the effects of large fires burning in heterogeneous landscapes. In the extended version of the model, vegetation is depicted as an amorphous network of combustible cells, and both radiation and convection from the flaming zone are considered in the preheating process of unburned cells. Examples are given to illustrate small-world effects and fire behavior near the percolation threshold. The model is applied to a Mediterranean fire that occurred in Corsica in 2009 showing a good agreement in terms of rate of spread, and area and shape of the burn. A study, based on a fractional factorial plan, is conducted to evaluate the influence of variations of model parameters on fire propagation.

ARMA modeling of stochastic processes in nuclear reactor with significant detection noise

International Nuclear Information System (INIS)

Zavaljevski, N.

1992-01-01

The theoretical basis of ARMA modelling of stochastic processes in nuclear reactor was presented in a previous paper, neglecting observational noise. The identification of real reactor data indicated that in some experiments the detection noise is significant. Thus a more rigorous theoretical modelling of stochastic processes in nuclear reactor is performed. Starting from the fundamental stochastic differential equations of the Langevin type for the interaction of the detector with neutron field, a new theoretical ARMA model is developed. preliminary identification results confirm the theoretical expectations. (author)

INFERENCE AND SENSITIVITY IN STOCHASTIC WIND POWER FORECAST MODELS.

KAUST Repository

Elkantassi, Soumaya

2017-10-03

Reliable forecasting of wind power generation is crucial to optimal control of costs in generation of electricity with respect to the electricity demand. Here, we propose and analyze stochastic wind power forecast models described by parametrized stochastic differential equations, which introduce appropriate fluctuations in numerical forecast outputs. We use an approximate maximum likelihood method to infer the model parameters taking into account the time correlated sets of data. Furthermore, we study the validity and sensitivity of the parameters for each model. We applied our models to Uruguayan wind power production as determined by historical data and corresponding numerical forecasts for the period of March 1 to May 31, 2016.

INFERENCE AND SENSITIVITY IN STOCHASTIC WIND POWER FORECAST MODELS.

KAUST Repository

Elkantassi, Soumaya; Kalligiannaki, Evangelia; Tempone, Raul

2017-01-01

Reliable forecasting of wind power generation is crucial to optimal control of costs in generation of electricity with respect to the electricity demand. Here, we propose and analyze stochastic wind power forecast models described by parametrized stochastic differential equations, which introduce appropriate fluctuations in numerical forecast outputs. We use an approximate maximum likelihood method to infer the model parameters taking into account the time correlated sets of data. Furthermore, we study the validity and sensitivity of the parameters for each model. We applied our models to Uruguayan wind power production as determined by historical data and corresponding numerical forecasts for the period of March 1 to May 31, 2016.

Parameter discovery in stochastic biological models using simulated annealing and statistical model checking.

Science.gov (United States)

Hussain, Faraz; Jha, Sumit K; Jha, Susmit; Langmead, Christopher J

2014-01-01

Stochastic models are increasingly used to study the behaviour of biochemical systems. While the structure of such models is often readily available from first principles, unknown quantitative features of the model are incorporated into the model as parameters. Algorithmic discovery of parameter values from experimentally observed facts remains a challenge for the computational systems biology community. We present a new parameter discovery algorithm that uses simulated annealing, sequential hypothesis testing, and statistical model checking to learn the parameters in a stochastic model. We apply our technique to a model of glucose and insulin metabolism used for in-silico validation of artificial pancreata and demonstrate its effectiveness by developing parallel CUDA-based implementation for parameter synthesis in this model.

Stochastic modelling in design of mechanical properties of nanometals

International Nuclear Information System (INIS)

Tengen, T.B.; Wejrzanowski, T.; Iwankiewicz, R.; Kurzydlowski, K.J.

2010-01-01

Polycrystalline nanometals are being fabricated through different processing routes and conditions. The consequence is that nanometals having the same mean grain size may have different grain size dispersion and, hence, may have different material properties. This has often led to conflicting reports from both theoretical and experimental findings about the evolutions of the mechanical properties of nanomaterials. The present paper employs stochastic model to study the impact of microstructure evolution during grain growth on the mechanical properties of polycrystalline nanometals. The stochastic model for grain growth and the stochastic model for changes in mechanical properties of nanomaterials are proposed. The model for the mechanical properties developed is tested on aluminium samples.Many salient features of the mechanical properties of the aluminium samples are revealed. The results show that the different mechanisms of grain growth impart different nature of response to the material mechanical properties. The conventional, homologous and anomalous temperature dependences of the yield stress have also been revealed to be due to different nature of interactions of the microstructures during evolution.

Bayesian inference for hybrid discrete-continuous stochastic kinetic models

International Nuclear Information System (INIS)

Sherlock, Chris; Golightly, Andrew; Gillespie, Colin S

2014-01-01

We consider the problem of efficiently performing simulation and inference for stochastic kinetic models. Whilst it is possible to work directly with the resulting Markov jump process (MJP), computational cost can be prohibitive for networks of realistic size and complexity. In this paper, we consider an inference scheme based on a novel hybrid simulator that classifies reactions as either ‘fast’ or ‘slow’ with fast reactions evolving as a continuous Markov process whilst the remaining slow reaction occurrences are modelled through a MJP with time-dependent hazards. A linear noise approximation (LNA) of fast reaction dynamics is employed and slow reaction events are captured by exploiting the ability to solve the stochastic differential equation driving the LNA. This simulation procedure is used as a proposal mechanism inside a particle MCMC scheme, thus allowing Bayesian inference for the model parameters. We apply the scheme to a simple application and compare the output with an existing hybrid approach and also a scheme for performing inference for the underlying discrete stochastic model. (paper)

Fitting Social Network Models Using Varying Truncation Stochastic Approximation MCMC Algorithm

KAUST Repository

Jin, Ick Hoon

2013-10-01

The exponential random graph model (ERGM) plays a major role in social network analysis. However, parameter estimation for the ERGM is a hard problem due to the intractability of its normalizing constant and the model degeneracy. The existing algorithms, such as Monte Carlo maximum likelihood estimation (MCMLE) and stochastic approximation, often fail for this problem in the presence of model degeneracy. In this article, we introduce the varying truncation stochastic approximation Markov chain Monte Carlo (SAMCMC) algorithm to tackle this problem. The varying truncation mechanism enables the algorithm to choose an appropriate starting point and an appropriate gain factor sequence, and thus to produce a reasonable parameter estimate for the ERGM even in the presence of model degeneracy. The numerical results indicate that the varying truncation SAMCMC algorithm can significantly outperform the MCMLE and stochastic approximation algorithms: for degenerate ERGMs, MCMLE and stochastic approximation often fail to produce any reasonable parameter estimates, while SAMCMC can do; for nondegenerate ERGMs, SAMCMC can work as well as or better than MCMLE and stochastic approximation. The data and source codes used for this article are available online as supplementary materials. © 2013 American Statistical Association, Institute of Mathematical Statistics, and Interface Foundation of North America.

Discontinuous finite element and characteristics methods for neutrons transport equation solution in heterogeneous grids

International Nuclear Information System (INIS)

Masiello, E.

2006-01-01

The principal goal of this manuscript is devoted to the investigation of a new type of heterogeneous mesh adapted to the shape of the fuel pins (fuel-clad-moderator). The new heterogeneous mesh guarantees the spatial modelling of the pin-cell with a minimum of regions. Two methods are investigated for the spatial discretization of the transport equation: the discontinuous finite element method and the method of characteristics for structured cells. These methods together with the new representation of the pin-cell result in an appreciable reduction of calculation points. They allow an exact modelling of the fuel pin-cell without spatial homogenization. A new synthetic acceleration technique based on an angular multigrid is also presented for the speed up of the inner iterations. These methods are good candidates for transport calculations for a nuclear reactor core. A second objective of this work is the application of method of characteristics for non-structured geometries to the study of double heterogeneity problem. The letters is characterized by fuel material with a stochastic dispersion of heterogeneous grains, and until now was solved with a model based on collision probabilities. We propose a new statistical model based on renewal-Markovian theory, which makes possible to take into account the stochastic nature of the problem and to avoid the approximations of the collision probability model. The numerical solution of this model is guaranteed by the method of characteristics. (author)

New Results on a Stochastic Duel Game with Each Force Consisting of Heterogeneous Units

Science.gov (United States)

2013-02-01

study of duel models dates back to the 1910s, when Lanchester (1916) proposed differential equations that govern the strength of each force through time...which gave rise to what later became known as Lanchester models. A stream of works extended the Lanchester models—which are deterministic in nature...869. Kress, M. and Talmor, I. (1999). A new look at the 3:1 rule of combat through markov stochastic lanchester models. The Journal of the Operational

High-Resolution Replication Profiles Define the Stochastic Nature of Genome Replication Initiation and Termination

Directory of Open Access Journals (Sweden)

Michelle Hawkins

2013-11-01

Full Text Available Eukaryotic genome replication is stochastic, and each cell uses a different cohort of replication origins. We demonstrate that interpreting high-resolution Saccharomyces cerevisiae genome replication data with a mathematical model allows quantification of the stochastic nature of genome replication, including the efficiency of each origin and the distribution of termination events. Single-cell measurements support the inferred values for stochastic origin activation time. A strain, in which three origins were inactivated, confirmed that the distribution of termination events is primarily dictated by the stochastic activation time of origins. Cell-to-cell variability in origin activity ensures that termination events are widely distributed across virtually the whole genome. We propose that the heterogeneity in origin usage contributes to genome stability by limiting potentially deleterious events from accumulating at particular loci.

Double diffusivity model under stochastic forcing

Science.gov (United States)

Chattopadhyay, Amit K.; Aifantis, Elias C.

2017-05-01

The "double diffusivity" model was proposed in the late 1970s, and reworked in the early 1980s, as a continuum counterpart to existing discrete models of diffusion corresponding to high diffusivity paths, such as grain boundaries and dislocation lines. It was later rejuvenated in the 1990s to interpret experimental results on diffusion in polycrystalline and nanocrystalline specimens where grain boundaries and triple grain boundary junctions act as high diffusivity paths. Technically, the model pans out as a system of coupled Fick-type diffusion equations to represent "regular" and "high" diffusivity paths with "source terms" accounting for the mass exchange between the two paths. The model remit was extended by analogy to describe flow in porous media with double porosity, as well as to model heat conduction in media with two nonequilibrium local temperature baths, e.g., ion and electron baths. Uncoupling of the two partial differential equations leads to a higher-ordered diffusion equation, solutions of which could be obtained in terms of classical diffusion equation solutions. Similar equations could also be derived within an "internal length" gradient (ILG) mechanics formulation applied to diffusion problems, i.e., by introducing nonlocal effects, together with inertia and viscosity, in a mechanics based formulation of diffusion theory. While being remarkably successful in studies related to various aspects of transport in inhomogeneous media with deterministic microstructures and nanostructures, its implications in the presence of stochasticity have not yet been considered. This issue becomes particularly important in the case of diffusion in nanopolycrystals whose deterministic ILG-based theoretical calculations predict a relaxation time that is only about one-tenth of the actual experimentally verified time scale. This article provides the "missing link" in this estimation by adding a vital element in the ILG structure, that of stochasticity, that takes into

Stochastic modeling concepts in groundwater and risk assessment: potential application to marine problems

International Nuclear Information System (INIS)

Hamed, Maged M.

2000-01-01

Parameter uncertainty is ubiquitous in marine environmental processes. Failure to account for this uncertainty may lead to erroneous results, and may have significant environmental and economic ramifications. Stochastic modeling of oil spill transport and fate is, therefore, central in the development of an oil spill contingency plan for new oil and gas projects. Over the past twenty years, several stochastic modeling tools have been developed for modeling parameter uncertainty, including the spectral, perturbation, and simulation methods. In this work we explore the application of a new stochastic methodology, the first-order reliability method (FORM), in oil spill modeling. FORM was originally developed in the structural reliability field and has been recently applied to various environmental problems. The method has many appealing features that makes it a powerful tool for modeling complex environmental systems. The theory of FORM is presented, identifying the features that distinguish the method from other stochastic tools. Different formulations to the reliability-based stochastic oil spill modeling are presented in a decision-analytic context. (Author)

Stochastic modeling analysis and simulation

CERN Document Server

Nelson, Barry L

1995-01-01

A coherent introduction to the techniques for modeling dynamic stochastic systems, this volume also offers a guide to the mathematical, numerical, and simulation tools of systems analysis. Suitable for advanced undergraduates and graduate-level industrial engineers and management science majors, it proposes modeling systems in terms of their simulation, regardless of whether simulation is employed for analysis. Beginning with a view of the conditions that permit a mathematical-numerical analysis, the text explores Poisson and renewal processes, Markov chains in discrete and continuous time, se

Digital hardware implementation of a stochastic two-dimensional neuron model.

Science.gov (United States)

Grassia, F; Kohno, T; Levi, T

2016-11-01

This study explores the feasibility of stochastic neuron simulation in digital systems (FPGA), which realizes an implementation of a two-dimensional neuron model. The stochasticity is added by a source of current noise in the silicon neuron using an Ornstein-Uhlenbeck process. This approach uses digital computation to emulate individual neuron behavior using fixed point arithmetic operation. The neuron model's computations are performed in arithmetic pipelines. It was designed in VHDL language and simulated prior to mapping in the FPGA. The experimental results confirmed the validity of the developed stochastic FPGA implementation, which makes the implementation of the silicon neuron more biologically plausible for future hybrid experiments. Copyright © 2017 Elsevier Ltd. All rights reserved.

Stochastic modelling of conjugate heat transfer in near-wall turbulence

International Nuclear Information System (INIS)

Pozorski, Jacek; Minier, Jean-Pierre

2006-01-01

The paper addresses the conjugate heat transfer in turbulent flows with temperature assumed to be a passive scalar. The Lagrangian approach is applied and the heat transfer is modelled with the use of stochastic particles. The intensity of thermal fluctuations in near-wall turbulence is determined from the scalar probability density function (PDF) with externally provided dynamical statistics. A stochastic model for the temperature field in the wall material is proposed and boundary conditions for stochastic particles at the solid-fluid interface are formulated. The heated channel flow with finite-thickness walls is considered as a validation case. Computation results for the mean temperature profiles and the variance of thermal fluctuations are presented and compared with available DNS data

Stochastic modelling of conjugate heat transfer in near-wall turbulence

Energy Technology Data Exchange (ETDEWEB)

Pozorski, Jacek [Institute of Fluid-Flow Machinery, Polish Academy of Sciences, Fiszera 14, 80952 Gdansk (Poland)]. E-mail: jp@imp.gda.pl; Minier, Jean-Pierre [Research and Development Division, Electricite de France, 6 quai Watier, 78400 Chatou (France)

2006-10-15

The paper addresses the conjugate heat transfer in turbulent flows with temperature assumed to be a passive scalar. The Lagrangian approach is applied and the heat transfer is modelled with the use of stochastic particles. The intensity of thermal fluctuations in near-wall turbulence is determined from the scalar probability density function (PDF) with externally provided dynamical statistics. A stochastic model for the temperature field in the wall material is proposed and boundary conditions for stochastic particles at the solid-fluid interface are formulated. The heated channel flow with finite-thickness walls is considered as a validation case. Computation results for the mean temperature profiles and the variance of thermal fluctuations are presented and compared with available DNS data.

Dynamic heterogeneity and DNA methylation in embryonic stem cells.

KAUST Repository

Singer, Zakary S

2014-07-01

Cell populations can be strikingly heterogeneous, composed of multiple cellular states, each exhibiting stochastic noise in its gene expression. A major challenge is to disentangle these two types of variability and to understand the dynamic processes and mechanisms that control them. Embryonic stem cells (ESCs) provide an ideal model system to address this issue because they exhibit heterogeneous and dynamic expression of functionally important regulatory factors. We analyzed gene expression in individual ESCs using single-molecule RNA-FISH and quantitative time-lapse movies. These data discriminated stochastic switching between two coherent (correlated) gene expression states and burst-like transcriptional noise. We further showed that the "2i" signaling pathway inhibitors modulate both types of variation. Finally, we found that DNA methylation plays a key role in maintaining these metastable states. Together, these results show how ESC gene expression states and dynamics arise from a combination of intrinsic noise, coherent cellular states, and epigenetic regulation.

Analysis and reconstruction of stochastic coupled map lattice models

International Nuclear Information System (INIS)

Coca, Daniel; Billings, Stephen A.

2003-01-01

The Letter introduces a general stochastic coupled lattice map model together with an algorithm to estimate the nodal equations involved based only on a small set of observable variables and in the presence of stochastic perturbations. More general forms of the Frobenius-Perron and the transfer operators, which describe the evolution of densities under the action of the CML transformation, are derived

Stochastic bifurcation in a model of love with colored noise

Science.gov (United States)

Yue, Xiaokui; Dai, Honghua; Yuan, Jianping

2015-07-01

In this paper, we wish to examine the stochastic bifurcation induced by multiplicative Gaussian colored noise in a dynamical model of love where the random factor is used to describe the complexity and unpredictability of psychological systems. First, the dynamics in deterministic love-triangle model are considered briefly including equilibrium points and their stability, chaotic behaviors and chaotic attractors. Then, the influences of Gaussian colored noise with different parameters are explored such as the phase plots, top Lyapunov exponents, stationary probability density function (PDF) and stochastic bifurcation. The stochastic P-bifurcation through a qualitative change of the stationary PDF will be observed and bifurcation diagram on parameter plane of correlation time and noise intensity is presented to find the bifurcation behaviors in detail. Finally, the top Lyapunov exponent is computed to determine the D-bifurcation when the noise intensity achieves to a critical value. By comparison, we find there is no connection between two kinds of stochastic bifurcation.

Inefficiency persistence and heterogeneity in Colombian electricity utilities

International Nuclear Information System (INIS)

Galán, Jorge E.; Pollitt, Michael G.

2014-01-01

The electricity reform in Colombia has exhibited gains in terms of reliability but its effects on firm efficiency and service quality have not been clear. Previous studies evaluating the performance of distribution companies after the reform have not found evidence of improvements, although large differences in efficiency have been found among firms. This suggests high inefficiency persistence and heterogeneity in the Colombian distribution sector. In this paper, we propose an extension of dynamic stochastic frontier models that accounts for unobserved heterogeneity in the inefficiency persistence and in the technology. The model incorporates total expenses, service quality and energy losses in an efficiency analysis of Colombian distributors over fifteen years after the reform. We identify the presence of high inefficiency persistence in the sector, and important differences between firms. In particular, rural companies and firms with small customers present low persistence and evidence the largest gains in efficiency during the period. However, increases in efficiency are only manifested during thelast five years when the main improvements in service quality and energy losses are presented. Overall, inefficiency persistence, customer density and consumption density are found to be important criteria to be considered for regulatory purposes. - Highlights: • We evaluate efficiency of Colombian electricity distributors after the reform. • We use a stochastic frontier model with dynamic effects and heterogeneity. • We find high inefficiency persistence but important differences among firms. • High persistent and low efficient firms should draw the attention of the regulator. • Recent regulation in quality has increased not only efficiency but also tariffs

Analysis of novel stochastic switched SILI epidemic models with continuous and impulsive control

Science.gov (United States)

Gao, Shujing; Zhong, Deming; Zhang, Yan

2018-04-01

In this paper, we establish two new stochastic switched epidemic models with continuous and impulsive control. The stochastic perturbations are considered for the natural death rate in each equation of the models. Firstly, a stochastic switched SILI model with continuous control schemes is investigated. By using Lyapunov-Razumikhin method, the sufficient conditions for extinction in mean are established. Our result shows that the disease could be die out theoretically if threshold value R is less than one, regardless of whether the disease-free solutions of the corresponding subsystems are stable or unstable. Then, a stochastic switched SILI model with continuous control schemes and pulse vaccination is studied. The threshold value R is derived. The global attractivity of the model is also obtained. At last, numerical simulations are carried out to support our results.

Heterogeneous traffic flow modelling using second-order macroscopic continuum model

Science.gov (United States)

Mohan, Ranju; Ramadurai, Gitakrishnan

2017-01-01

Modelling heterogeneous traffic flow lacking in lane discipline is one of the emerging research areas in the past few years. The two main challenges in modelling are: capturing the effect of varying size of vehicles, and the lack in lane discipline, both of which together lead to the 'gap filling' behaviour of vehicles. The same section length of the road can be occupied by different types of vehicles at the same time, and the conventional measure of traffic concentration, density (vehicles per lane per unit length), is not a good measure for heterogeneous traffic modelling. First aim of this paper is to have a parsimonious model of heterogeneous traffic that can capture the unique phenomena of gap filling. Second aim is to emphasize the suitability of higher-order models for modelling heterogeneous traffic. Third, the paper aims to suggest area occupancy as concentration measure of heterogeneous traffic lacking in lane discipline. The above mentioned two main challenges of heterogeneous traffic flow are addressed by extending an existing second-order continuum model of traffic flow, using area occupancy for traffic concentration instead of density. The extended model is calibrated and validated with field data from an arterial road in Chennai city, and the results are compared with those from few existing generalized multi-class models.

Stochastic Subspace Modelling of Turbulence

DEFF Research Database (Denmark)

Sichani, Mahdi Teimouri; Pedersen, B. J.; Nielsen, Søren R.K.

2009-01-01

positive definite cross-spectral density matrix a frequency response matrix is constructed which determines the turbulence vector as a linear filtration of Gaussian white noise. Finally, an accurate state space modelling method is proposed which allows selection of an appropriate model order......, and estimation of a state space model for the vector turbulence process incorporating its phase spectrum in one stage, and its results are compared with a conventional ARMA modelling method.......Turbulence of the incoming wind field is of paramount importance to the dynamic response of civil engineering structures. Hence reliable stochastic models of the turbulence should be available from which time series can be generated for dynamic response and structural safety analysis. In the paper...

Reflected stochastic differential equation models for constrained animal movement

Science.gov (United States)

Hanks, Ephraim M.; Johnson, Devin S.; Hooten, Mevin B.

2017-01-01

Movement for many animal species is constrained in space by barriers such as rivers, shorelines, or impassable cliffs. We develop an approach for modeling animal movement constrained in space by considering a class of constrained stochastic processes, reflected stochastic differential equations. Our approach generalizes existing methods for modeling unconstrained animal movement. We present methods for simulation and inference based on augmenting the constrained movement path with a latent unconstrained path and illustrate this augmentation with a simulation example and an analysis of telemetry data from a Steller sea lion (Eumatopias jubatus) in southeast Alaska.

Stochastic interest rates model in compounding | Galadima ...

African Journals Online (AJOL)

Stochastic interest rates model in compounding. ... in finance, real estate, insurance, accounting and other areas of business administration. The assumption that future rates are fixed and known with certainty at the beginning of an investment, ...

A stochastic cloud model for cloud and ozone retrievals from UV measurements

International Nuclear Information System (INIS)

Efremenko, Dmitry S.; Schüssler, Olena; Doicu, Adrian; Loyola, Diego

2016-01-01

The new generation of satellite instruments provides measurements in and around the Oxygen A-band on a global basis and with a relatively high spatial resolution. These data are commonly used for the determination of cloud properties. A stochastic model and radiative transfer model, previously developed by the authors, is used as the forward model component in retrievals of cloud parameters and ozone total and partial columns. The cloud retrieval algorithm combines local and global optimization routines, and yields a retrieval accuracy of about 1% and a fast computational time. Retrieved parameters are the cloud optical thickness and the cloud-top height. It was found that the use of the independent pixel approximation instead of the stochastic cloud model leads to large errors in the retrieved cloud parameters, as well as, in the retrieved ozone height resolved partial columns. The latter can be reduced by using the stochastic cloud model to compute the optimal value of the regularization parameter in the framework of Tikhonov regularization. - Highlights: • A stochastic radiative transfer model for retrieving clouds/ozone is designed. • Errors of independent pixel approximation (IPA) for O3 total column are small. • The error of IPA for ozone profile retrieval may become large. • The use of stochastic model reduces the error of ozone profile retrieval.

Threshold Dynamics in Stochastic SIRS Epidemic Models with Nonlinear Incidence and Vaccination

Directory of Open Access Journals (Sweden)

Lei Wang

2017-01-01

Full Text Available In this paper, the dynamical behaviors for a stochastic SIRS epidemic model with nonlinear incidence and vaccination are investigated. In the models, the disease transmission coefficient and the removal rates are all affected by noise. Some new basic properties of the models are found. Applying these properties, we establish a series of new threshold conditions on the stochastically exponential extinction, stochastic persistence, and permanence in the mean of the disease with probability one for the models. Furthermore, we obtain a sufficient condition on the existence of unique stationary distribution for the model. Finally, a series of numerical examples are introduced to illustrate our main theoretical results and some conjectures are further proposed.

Rate response of neurons subject to fast or frozen noise: from stochastic and homogeneous to deterministic and heterogeneous populations.

Science.gov (United States)

Alijani, Azadeh Khajeh; Richardson, Magnus J E

2011-07-01

The response of a neuronal population to afferent drive can be expected to be sensitive to both the distribution and dynamics of membrane voltages within the population. Voltage fluctuations can be driven by synaptic noise, neuromodulators, or cellular inhomogeneities: processes ranging from millisecond autocorrelation times to effectively static or "frozen" noise. Here we extend previous studies of filtered fluctuations to the experimentally verified exponential integrate-and-fire model. How fast or frozen fluctuations affect the steady-state rate and firing-rate response are both examined using perturbative solutions and limits of a 1 + 2 dimensional Fokker-Planck equation. The central finding is that, under conditions of a more-or-less constant population voltage variance, the firing-rate response is only weakly dependent on the fluctuation filter constant: The voltage distribution is the principal determinant of the population response. This result is unexpected given the nature of the systems underlying the extreme limits of fast and frozen fluctuations; the first limit represents a homogeneous population of neurons firing stochastically, whereas the second limit is equivalent to a heterogeneous population of neurons firing deterministically.

Models of the stochastic activity of neurones

CERN Document Server

Holden, Arun Vivian

1976-01-01

These notes have grown from a series of seminars given at Leeds between 1972 and 1975. They represent an attempt to gather together the different kinds of model which have been proposed to account for the stochastic activity of neurones, and to provide an introduction to this area of mathematical biology. A striking feature of the electrical activity of the nervous system is that it appears stochastic: this is apparent at all levels of recording, ranging from intracellular recordings to the electroencephalogram. The chapters start with fluctuations in membrane potential, proceed through single unit and synaptic activity and end with the behaviour of large aggregates of neurones: L have chgaen this seque~~e\\/~~';uggest that the interesting behaviourr~f :the nervous system - its individuality, variability and dynamic forms - may in part result from the stochastic behaviour of its components. I would like to thank Dr. Julio Rubio for reading and commenting on the drafts, Mrs. Doris Beighton for producing the fin...

Stochastic modelling of avascular tumour growth and therapy

International Nuclear Information System (INIS)

Sahoo, S; Sahoo, A; Shearer, S F C

2011-01-01

In this paper, a generalized stochastic model for the growth of avascular tumours is presented. This model captures the dynamical evolution of avascular tumour cell subpopulations by incorporating Gaussian white noise into the growth rate of the mitotic function. This work generalizes the deterministic model proposed by Sherratt and Chaplain (2001 J. Math. Biol. 43 291) where they formulated a tumour model in an in vivo setting, in terms of continuum densities of proliferating, quiescent and necrotic cells. Detailed simulations of our model show that the inclusion of Gaussian noise in the original model of Sherratt and Chaplain substantially distorts the overall structure of the density profiles in addition to reducing the speed of tumour growth. Within this stochastic carcinogenesis framework the action of therapy is also investigated by replacing Gaussian white noise with a therapy term. We compare a constant therapy protocol with a logarithmic time-dependent protocol. Our results predict that a logarithmic therapy is more effective than the constant therapy protocol.

Stochastic dynamics modeling solute transport in porous media modeling solute transport in porous media

CERN Document Server

Kulasiri, Don

2002-01-01

Most of the natural and biological phenomena such as solute transport in porous media exhibit variability which can not be modeled by using deterministic approaches. There is evidence in natural phenomena to suggest that some of the observations can not be explained by using the models which give deterministic solutions. Stochastic processes have a rich repository of objects which can be used to express the randomness inherent in the system and the evolution of the system over time. The attractiveness of the stochastic differential equations (SDE) and stochastic partial differential equations (SPDE) come from the fact that we can integrate the variability of the system along with the scientific knowledge pertaining to the system. One of the aims of this book is to explaim some useufl concepts in stochastic dynamics so that the scientists and engineers with a background in undergraduate differential calculus could appreciate the applicability and appropriateness of these developments in mathematics. The ideas ...

The effect of soil heterogeneity on ATES performance

Science.gov (United States)

Sommer, W.; Rijnaarts, H.; Grotenhuis, T.; van Gaans, P.

2012-04-01

Due to an increasing demand for sustainable energy, application of Aquifer Thermal Energy Storage (ATES) is growing rapidly. Large-scale application of ATES is limited by the space that is available in the subsurface. Especially in urban areas, suboptimal performance is expected due to thermal interference between individual wells of a single system, or interference with other ATES systems or groundwater abstractions. To avoid thermal interference there are guidelines on well spacing. However, these guidelines, and also design calculations, are based on the assumption of a homogeneous subsurface, while studies report a standard deviation in logpermeability of 1 to 2 for unconsolidated aquifers (Gelhar, 1993). Such heterogeneity may create preferential pathways, reducing ATES performance due to increased advective heat loss or interference between ATES wells. The role of hydraulic heterogeneity of the subsurface related to ATES performance has received little attention in literature. Previous research shows that even small amounts of heterogeneity can result in considerable uncertainty in the distribution of thermal energy in the subsurface and an increased radius of influence (Ferguson, 2007). This is supported by subsurface temperature measurements around ATES wells, which suggest heterogeneity gives rise to preferential pathways and short-circuiting between ATES wells (Bridger and Allen, 2010). Using 3-dimensional stochastic heat transport modeling, we quantified the influence of heterogeneity on the performance of a doublet well energy storage system. The following key parameters are varied to study their influence on thermal recovery and thermal balance: 1) regional flow velocity, 2) distance between wells and 3) characteristics of the heterogeneity. Results show that heterogeneity at the scale of a doublet ATES system introduces an uncertainty up to 18% in expected thermal recovery. The uncertainty increases with decreasing distance between ATES wells. The

Methods and models in mathematical biology deterministic and stochastic approaches

CERN Document Server

Müller, Johannes

2015-01-01

This book developed from classes in mathematical biology taught by the authors over several years at the Technische Universität München. The main themes are modeling principles, mathematical principles for the analysis of these models, and model-based analysis of data. The key topics of modern biomathematics are covered: ecology, epidemiology, biochemistry, regulatory networks, neuronal networks, and population genetics. A variety of mathematical methods are introduced, ranging from ordinary and partial differential equations to stochastic graph theory and  branching processes. A special emphasis is placed on the interplay between stochastic and deterministic models.

Quantifying hidden individual heterogeneity

DEFF Research Database (Denmark)

Steiner, Ulrich; Lenart, Adam; Vaupel, James W.

Aging is assumed to be driven by the accumulation of damage or some other aging factor which shapes demographic patterns, including the classical late age mortality plateaus. However to date, heterogeneity in these damage stages is not observed. Here, we estimate underlying stage distributions...... and stage dynamics, based on observed survival patterns of isoclonal bacteria. Our results reveal demographic dynamics being dominated by low damage stages and transmission of damage from mother to daughters is low. Still, our models are too simplistic and deterministic. Explaining the observed data...... requires more stochastic processes as our current models includes. We are only at the beginning of understanding the diverse mechanism behind aging and the shaping of senescence....

Doubly stochastic Poisson process models for precipitation at fine time-scales

Science.gov (United States)

Ramesh, Nadarajah I.; Onof, Christian; Xie, Dichao

2012-09-01

This paper considers a class of stochastic point process models, based on doubly stochastic Poisson processes, in the modelling of rainfall. We examine the application of this class of models, a neglected alternative to the widely-known Poisson cluster models, in the analysis of fine time-scale rainfall intensity. These models are mainly used to analyse tipping-bucket raingauge data from a single site but an extension to multiple sites is illustrated which reveals the potential of this class of models to study the temporal and spatial variability of precipitation at fine time-scales.

Stochastic models, estimation, and control

CERN Document Server

Maybeck, Peter S

1982-01-01

This volume builds upon the foundations set in Volumes 1 and 2. Chapter 13 introduces the basic concepts of stochastic control and dynamic programming as the fundamental means of synthesizing optimal stochastic control laws.

GillesPy: A Python Package for Stochastic Model Building and Simulation

OpenAIRE

Abel, John H.; Drawert, Brian; Hellander, Andreas; Petzold, Linda R.

2016-01-01

GillesPy is an open-source Python package for model construction and simulation of stochastic biochemical systems. GillesPy consists of a Python framework for model building and an interface to the StochKit2 suite of efficient simulation algorithms based on the Gillespie stochastic simulation algorithms (SSA). To enable intuitive model construction and seamless integration into the scientific Python stack, we present an easy to understand, action-oriented programming interface. Here, we descr...

Stochastic Analysis 2010

CERN Document Server

Crisan, Dan

2011-01-01

"Stochastic Analysis" aims to provide mathematical tools to describe and model high dimensional random systems. Such tools arise in the study of Stochastic Differential Equations and Stochastic Partial Differential Equations, Infinite Dimensional Stochastic Geometry, Random Media and Interacting Particle Systems, Super-processes, Stochastic Filtering, Mathematical Finance, etc. Stochastic Analysis has emerged as a core area of late 20th century Mathematics and is currently undergoing a rapid scientific development. The special volume "Stochastic Analysis 2010" provides a sa

Predicting population extinction or disease outbreaks with stochastic models

Directory of Open Access Journals (Sweden)

Linda J. S. Allen

2017-01-01

Full Text Available Models of exponential growth, logistic growth and epidemics are common applications in undergraduate differential equation courses. The corresponding stochastic models are not part of these courses, although when population sizes are small their behaviour is often more realistic and distinctly different from deterministic models. For example, the randomness associated with births and deaths may lead to population extinction even in an exponentially growing population. Some background in continuous-time Markov chains and applications to populations, epidemics and cancer are presented with a goal to introduce this topic into the undergraduate mathematics curriculum that will encourage further investigation into problems on conservation, infectious diseases and cancer therapy. MATLAB programs for graphing sample paths of stochastic models are provided in the Appendix.

Predicting the Stochastic Properties of the Shallow Subsurface for Improved Geophysical Modeling

Science.gov (United States)

Stroujkova, A.; Vynne, J.; Bonner, J.; Lewkowicz, J.

2005-12-01

Strong ground motion data from numerous explosive field experiments and from moderate to large earthquakes show significant variations in amplitude and waveform shape with respect to both azimuth and range. Attempts to model these variations using deterministic models have often been unsuccessful. It has been hypothesized that a stochastic description of the geological medium is a more realistic approach. To estimate the stochastic properties of the shallow subsurface, we use Measurement While Drilling (MWD) data, which are routinely collected by mines in order to facilitate design of blast patterns. The parameters, such as rotation speed of the drill, torque, and penetration rate, are used to compute the rock's Specific Energy (SE), which is then related to a blastability index. We use values of SE measured at two different mines and calibrated to laboratory measurements of rock properties to determine correlation lengths of the subsurface rocks in 2D, needed to obtain 2D and 3D stochastic models. The stochastic models are then combined with the deterministic models and used to compute synthetic seismic waveforms.

Stochastic Growth Models with No Discounting

Czech Academy of Sciences Publication Activity Database

Sladký, Karel

2007-01-01

Roč. 15, č. 4 (2007), s. 88-98 ISSN 0572-3043 R&D Projects: GA ČR(CZ) GA402/06/0990; GA ČR GA402/05/0115 Institutional research plan: CEZ:AV0Z10750506 Keywords : economic dynamics * stochastic version of the Ramsey growth model * Markov decision processes Subject RIV: AH - Economics

Application of Stochastic Partial Differential Equations to Reservoir Property Modelling

KAUST Repository

Potsepaev, R.

2010-09-06

Existing algorithms of geostatistics for stochastic modelling of reservoir parameters require a mapping (the \\'uvt-transform\\') into the parametric space and reconstruction of a stratigraphic co-ordinate system. The parametric space can be considered to represent a pre-deformed and pre-faulted depositional environment. Existing approximations of this mapping in many cases cause significant distortions to the correlation distances. In this work we propose a coordinate free approach for modelling stochastic textures through the application of stochastic partial differential equations. By avoiding the construction of a uvt-transform and stratigraphic coordinates, one can generate realizations directly in the physical space in the presence of deformations and faults. In particular the solution of the modified Helmholtz equation driven by Gaussian white noise is a zero mean Gaussian stationary random field with exponential correlation function (in 3-D). This equation can be used to generate realizations in parametric space. In order to sample in physical space we introduce a stochastic elliptic PDE with tensor coefficients, where the tensor is related to correlation anisotropy and its variation is physical space.

Survival Analysis of a Nonautonomous Logistic Model with Stochastic Perturbation

Directory of Open Access Journals (Sweden)

Chun Lu

2012-01-01

Full Text Available Taking white noise into account, a stochastic nonautonomous logistic model is proposed and investigated. Sufficient conditions for extinction, nonpersistence in the mean, weak persistence, stochastic permanence, and global asymptotic stability are established. Moreover, the threshold between weak persistence and extinction is obtained. Finally, we introduce some numerical simulink graphics to illustrate our main results.

Exercise effects in a virtual type 1 diabetes patient: Using stochastic differential equations for model extension

DEFF Research Database (Denmark)

Duun-Henriksen, Anne Katrine; Schmidt, S.; Nørgaard, K.

2013-01-01

extension incorporating exercise effects on insulin and glucose dynamics. Our model is constructed as a stochastic state space model consisting of a set of stochastic differential equations (SDEs). In a stochastic state space model, the residual error is split into random measurement error...

Expansion or extinction: deterministic and stochastic two-patch models with Allee effects.

Science.gov (United States)

Kang, Yun; Lanchier, Nicolas

2011-06-01

We investigate the impact of Allee effect and dispersal on the long-term evolution of a population in a patchy environment. Our main focus is on whether a population already established in one patch either successfully invades an adjacent empty patch or undergoes a global extinction. Our study is based on the combination of analytical and numerical results for both a deterministic two-patch model and a stochastic counterpart. The deterministic model has either two, three or four attractors. The existence of a regime with exactly three attractors only appears when patches have distinct Allee thresholds. In the presence of weak dispersal, the analysis of the deterministic model shows that a high-density and a low-density populations can coexist at equilibrium in nearby patches, whereas the analysis of the stochastic model indicates that this equilibrium is metastable, thus leading after a large random time to either a global expansion or a global extinction. Up to some critical dispersal, increasing the intensity of the interactions leads to an increase of both the basin of attraction of the global extinction and the basin of attraction of the global expansion. Above this threshold, for both the deterministic and the stochastic models, the patches tend to synchronize as the intensity of the dispersal increases. This results in either a global expansion or a global extinction. For the deterministic model, there are only two attractors, while the stochastic model no longer exhibits a metastable behavior. In the presence of strong dispersal, the limiting behavior is entirely determined by the value of the Allee thresholds as the global population size in the deterministic and the stochastic models evolves as dictated by their single-patch counterparts. For all values of the dispersal parameter, Allee effects promote global extinction in terms of an expansion of the basin of attraction of the extinction equilibrium for the deterministic model and an increase of the

Dynamic stochastic accumulation model with application to pension savings management

Directory of Open Access Journals (Sweden)

Melicherčik Igor

2010-01-01

Full Text Available We propose a dynamic stochastic accumulation model for determining optimal decision between stock and bond investments during accumulation of pension savings. Stock prices are assumed to be driven by the geometric Brownian motion. Interest rates are modeled by means of the Cox-Ingersoll-Ross model. The optimal decision as a solution to the corresponding dynamic stochastic program is a function of the duration of saving, the level of savings and the short rate. Qualitative and quantitative properties of the optimal solution are analyzed. The model is tested on the funded pillar of the Slovak pension system. The results are calculated for various risk preferences of a saver.

Stochastic modeling of financial electricity contracts

International Nuclear Information System (INIS)

Benth, Fred Espen; Koekebakker, Steen

2008-01-01

We discuss the modeling of electricity contracts traded in many deregulated power markets. These forward/futures type contracts deliver (either physically or financially) electricity over a specified time period, and is frequently referred to as swaps since they in effect represent an exchange of fixed for floating electricity price. We propose to use the Heath-Jarrow-Morton approach to model swap prices since the notion of a spot price is not easily defined in these markets. For general stochastic dynamical models, we connect the spot price, the instantaneous-delivery forward price and the swap price, and analyze two different ways to apply the Heath-Jarrow-Morton approach to swap pricing: Either one specifies a dynamics for the non-existing instantaneous-delivery forwards and derives the implied swap dynamics, or one models directly on the swaps. The former is shown to lead to quite complicated stochastic models for the swap price, even when the forward dynamics is simple. The latter has some theoretical problems due to a no-arbitrage condition that has to be satisfied for swaps with overlapping delivery periods. To overcome this problem, a practical modeling approach is analyzed. The market is supposed only to consist of non-overlapping swaps, and these are modelled directly. A thorough empirical study is performed using data collected from Nord Pool. Our investigations demonstrate that it is possible to state reasonable models for the swap price dynamics which is analytically tractable for risk management and option pricing purposes, however, this is an area of further research. (author)

Real-time forecasting of an epidemic using a discrete time stochastic model: a case study of pandemic influenza (H1N1-2009

Directory of Open Access Journals (Sweden)

Nishiura Hiroshi

2011-02-01

Full Text Available Abstract Background Real-time forecasting of epidemics, especially those based on a likelihood-based approach, is understudied. This study aimed to develop a simple method that can be used for the real-time epidemic forecasting. Methods A discrete time stochastic model, accounting for demographic stochasticity and conditional measurement, was developed and applied as a case study to the weekly incidence of pandemic influenza (H1N1-2009 in Japan. By imposing a branching process approximation and by assuming the linear growth of cases within each reporting interval, the epidemic curve is predicted using only two parameters. The uncertainty bounds of the forecasts are computed using chains of conditional offspring distributions. Results The quality of the forecasts made before the epidemic peak appears largely to depend on obtaining valid parameter estimates. The forecasts of both weekly incidence and final epidemic size greatly improved at and after the epidemic peak with all the observed data points falling within the uncertainty bounds. Conclusions Real-time forecasting using the discrete time stochastic model with its simple computation of the uncertainty bounds was successful. Because of the simplistic model structure, the proposed model has the potential to additionally account for various types of heterogeneity, time-dependent transmission dynamics and epidemiological details. The impact of such complexities on forecasting should be explored when the data become available as part of the disease surveillance.

Backward-stochastic-differential-equation approach to modeling of gene expression.

Science.gov (United States)

Shamarova, Evelina; Chertovskih, Roman; Ramos, Alexandre F; Aguiar, Paulo

2017-03-01

In this article, we introduce a backward method to model stochastic gene expression and protein-level dynamics. The protein amount is regarded as a diffusion process and is described by a backward stochastic differential equation (BSDE). Unlike many other SDE techniques proposed in the literature, the BSDE method is backward in time; that is, instead of initial conditions it requires the specification of end-point ("final") conditions, in addition to the model parametrization. To validate our approach we employ Gillespie's stochastic simulation algorithm (SSA) to generate (forward) benchmark data, according to predefined gene network models. Numerical simulations show that the BSDE method is able to correctly infer the protein-level distributions that preceded a known final condition, obtained originally from the forward SSA. This makes the BSDE method a powerful systems biology tool for time-reversed simulations, allowing, for example, the assessment of the biological conditions (e.g., protein concentrations) that preceded an experimentally measured event of interest (e.g., mitosis, apoptosis, etc.).

Stochastic inequalities and applications to dynamics analysis of a novel SIVS epidemic model with jumps

Directory of Open Access Journals (Sweden)

Xiaona Leng

2017-06-01

Full Text Available Abstract This paper proposes a new nonlinear stochastic SIVS epidemic model with double epidemic hypothesis and Lévy jumps. The main purpose of this paper is to investigate the threshold dynamics of the stochastic SIVS epidemic model. By using the technique of a series of stochastic inequalities, we obtain sufficient conditions for the persistence in mean and extinction of the stochastic system and the threshold which governs the extinction and the spread of the epidemic diseases. Finally, this paper describes the results of numerical simulations investigating the dynamical effects of stochastic disturbance. Our results significantly improve and generalize the corresponding results in recent literatures. The developed theoretical methods and stochastic inequalities technique can be used to investigate the high-dimensional nonlinear stochastic differential systems.

On the interpretation of double-packer tests in heterogeneous porous media: Numerical simulations using the stochastic continuum analogue

International Nuclear Information System (INIS)

Follin, S.

1992-12-01

Flow in fractured crystalline (hard) rocks is of interest in Sweden for assessing the postclosure radiological safety of a deep repository for high-level nuclear waste. For simulation of flow and mass transport in the far field different porous media concepts are often used, whereas discrete fracture/channel network concepts are often used for near-field simulations. Due to lack of data, it is generally necessary to have resort to single-hole double-packer test data for the far-field simulations, i.e., test data on a small scale are regularized in order to fit a comparatively coarser numerical discretization, which is governed by various computational constraints. In the present study the Monte Carlo method is used to investigate the relationship between the transmissivity value interpreted and the corresponding radius of influence in conjunction with single-hole double-packer tests in heterogeneous formations. The numerical flow domain is treated as a two-dimensional heterogeneous porous medium with a spatially varying diffusivity on 3 m scale. The Monte Carlo simulations demonstrate the sensitivity to the correlation range of a spatially varying diffusivity field. In contradiction to what is tacitly assumed in stochastic subsurface hydrology, the results show that the lateral support scale (e.g., the radius of influence) of transmissivity measurements in heterogeneous porous media is a random variable, which is affected by both the hydraulic and statistical characteristics. If these results are general, the traditional methods for scaling-up, assuming a constant lateral scale of support and a multi normal distribution, may lead to an underestimation of the persistence and connectivity of transmissive zones, particularly in highly heterogeneous porous media

Pricing long-dated insurance contracts with stochastic interest rates and stochastic volatility

NARCIS (Netherlands)

van Haastrecht, A.; Lord, R.; Pelsser, A.; Schrager, D.

2009-01-01

We consider the pricing of long-dated insurance contracts under stochastic interest rates and stochastic volatility. In particular, we focus on the valuation of insurance options with long-term equity or foreign exchange exposures. Our modeling framework extends the stochastic volatility model of

Stochastic Modelling of the Diffusion Coefficient for Concrete

DEFF Research Database (Denmark)

Thoft-Christensen, Palle

In the paper, a new stochastic modelling of the diffusion coefficient D is presented. The modelling is based on physical understanding of the diffusion process and on some recent experimental results. The diffusion coefficients D is strongly dependent on the w/c ratio and the temperature....

Stochastic models for transport in a fluidized bed

NARCIS (Netherlands)

Dehling, H.G; Hoffmann, A.C; Stuut, H.W.

1999-01-01

In this paper we study stochastic models for the transport of particles in a fluidized bed reactor and compute the associated residence time distribution (RTD). Our main model is basically a diffusion process in [0;A] with reflecting/absorbing boundary conditions, modified by allowing jumps to the

Individualism in plant populations: using stochastic differential equations to model individual neighbourhood-dependent plant growth.

Science.gov (United States)

Lv, Qiming; Schneider, Manuel K; Pitchford, Jonathan W

2008-08-01

We study individual plant growth and size hierarchy formation in an experimental population of Arabidopsis thaliana, within an integrated analysis that explicitly accounts for size-dependent growth, size- and space-dependent competition, and environmental stochasticity. It is shown that a Gompertz-type stochastic differential equation (SDE) model, involving asymmetric competition kernels and a stochastic term which decreases with the logarithm of plant weight, efficiently describes individual plant growth, competition, and variability in the studied population. The model is evaluated within a Bayesian framework and compared to its deterministic counterpart, and to several simplified stochastic models, using distributional validation. We show that stochasticity is an important determinant of size hierarchy and that SDE models outperform the deterministic model if and only if structural components of competition (asymmetry; size- and space-dependence) are accounted for. Implications of these results are discussed in the context of plant ecology and in more general modelling situations.

Stochastic simulation of karst conduit networks

Science.gov (United States)

Pardo-Igúzquiza, Eulogio; Dowd, Peter A.; Xu, Chaoshui; Durán-Valsero, Juan José

2012-01-01

Karst aquifers have very high spatial heterogeneity. Essentially, they comprise a system of pipes (i.e., the network of conduits) superimposed on rock porosity and on a network of stratigraphic surfaces and fractures. This heterogeneity strongly influences the hydraulic behavior of the karst and it must be reproduced in any realistic numerical model of the karst system that is used as input to flow and transport modeling. However, the directly observed karst conduits are only a small part of the complete karst conduit system and knowledge of the complete conduit geometry and topology remains spatially limited and uncertain. Thus, there is a special interest in the stochastic simulation of networks of conduits that can be combined with fracture and rock porosity models to provide a realistic numerical model of the karst system. Furthermore, the simulated model may be of interest per se and other uses could be envisaged. The purpose of this paper is to present an efficient method for conditional and non-conditional stochastic simulation of karst conduit networks. The method comprises two stages: generation of conduit geometry and generation of topology. The approach adopted is a combination of a resampling method for generating conduit geometries from templates and a modified diffusion-limited aggregation method for generating the network topology. The authors show that the 3D karst conduit networks generated by the proposed method are statistically similar to observed karst conduit networks or to a hypothesized network model. The statistical similarity is in the sense of reproducing the tortuosity index of conduits, the fractal dimension of the network, the direction rose of directions, the Z-histogram and Ripley's K-function of the bifurcation points (which differs from a random allocation of those bifurcation points). The proposed method (1) is very flexible, (2) incorporates any experimental data (conditioning information) and (3) can easily be modified when

Stochastic dynamical models for ecological regime shifts

DEFF Research Database (Denmark)

Møller, Jan Kloppenborg; Carstensen, Jacob; Madsen, Henrik

the physical and biological knowledge of the system, and nonlinearities introduced here can generate regime shifts or enhance the probability of regime shifts in the case of stochastic models, typically characterized by a threshold value for the known driver. A simple model for light competition between...... definition and stability of regimes become less subtle. Ecological regime shifts and their modeling must be viewed in a probabilistic manner, particularly if such model results are to be used in ecosystem management....

Reliability models for a nonrepairable system with heterogeneous components having a phase-type time-to-failure distribution

International Nuclear Information System (INIS)

Kim, Heungseob; Kim, Pansoo

2017-01-01

This research paper presents practical stochastic models for designing and analyzing the time-dependent reliability of nonrepairable systems. The models are formulated for nonrepairable systems with heterogeneous components having phase-type time-to-failure distributions by a structured continuous time Markov chain (CTMC). The versatility of the phase-type distributions enhances the flexibility and practicality of the systems. By virtue of these benefits, studies in reliability engineering can be more advanced than the previous studies. This study attempts to solve a redundancy allocation problem (RAP) by using these new models. The implications of mixing components, redundancy levels, and redundancy strategies are simultaneously considered to maximize the reliability of a system. An imperfect switching case in a standby redundant system is also considered. Furthermore, the experimental results for a well-known RAP benchmark problem are presented to demonstrate the approximating error of the previous reliability function for a standby redundant system and the usefulness of the current research. - Highlights: • Phase-type time-to-failure distribution is used for components. • Reliability model for nonrepairable system is developed using Markov chain. • System is composed of heterogeneous components. • Model provides the real value of standby system reliability not an approximation. • Redundancy allocation problem is used to show usefulness of this model.

Stochastic processes, multiscale modeling, and numerical methods for computational cellular biology

CERN Document Server

2017-01-01

This book focuses on the modeling and mathematical analysis of stochastic dynamical systems along with their simulations. The collected chapters will review fundamental and current topics and approaches to dynamical systems in cellular biology. This text aims to develop improved mathematical and computational methods with which to study biological processes. At the scale of a single cell, stochasticity becomes important due to low copy numbers of biological molecules, such as mRNA and proteins that take part in biochemical reactions driving cellular processes. When trying to describe such biological processes, the traditional deterministic models are often inadequate, precisely because of these low copy numbers. This book presents stochastic models, which are necessary to account for small particle numbers and extrinsic noise sources. The complexity of these models depend upon whether the biochemical reactions are diffusion-limited or reaction-limited. In the former case, one needs to adopt the framework of s...

Stochastic Modeling of Past Volcanic Crises

Science.gov (United States)

Woo, Gordon

2018-01-01

The statistical foundation of disaster risk analysis is past experience. From a scientific perspective, history is just one realization of what might have happened, given the randomness and chaotic dynamics of Nature. Stochastic analysis of the past is an exploratory exercise in counterfactual history, considering alternative possible scenarios. In particular, the dynamic perturbations that might have transitioned a volcano from an unrest to an eruptive state need to be considered. The stochastic modeling of past volcanic crises leads to estimates of eruption probability that can illuminate historical volcanic crisis decisions. It can also inform future economic risk management decisions in regions where there has been some volcanic unrest, but no actual eruption for at least hundreds of years. Furthermore, the availability of a library of past eruption probabilities would provide benchmark support for estimates of eruption probability in future volcanic crises.

A stochastic SIRS epidemic model with infectious force under intervention strategies

Science.gov (United States)

Cai, Yongli; Kang, Yun; Banerjee, Malay; Wang, Weiming

2015-12-01

In this paper, we extend a classical SIRS epidemic model with the infectious forces under intervention strategies from a deterministic framework to a stochastic differential equation (SDE) one through introducing random fluctuations. The value of our study lies in two aspects. Mathematically, by using the Markov semigroups theory, we prove that the reproduction number R0S can be used to govern the stochastic dynamics of SDE model. If R0S 1, under mild extra conditions, it has an endemic stationary distribution which leads to the stochastical persistence of the disease. Epidemiologically, we find that random fluctuations can suppress disease outbreak, which can provide us some useful control strategies to regulate disease dynamics.

A Stochastic Operational Planning Model for Smart Power Systems

Directory of Open Access Journals (Sweden)

Sh. Jadid

2014-12-01

Full Text Available Smart Grids are result of utilizing novel technologies such as distributed energy resources, and communication technologies in power system to compensate some of its defects. Various power resources provide some benefits for operation domain however, power system operator should use a powerful methodology to manage them. Renewable resources and load add uncertainty to the problem. So, independent system operator should use a stochastic method to manage them. A Stochastic unit commitment is presented in this paper to schedule various power resources such as distributed generation units, conventional thermal generation units, wind and PV farms, and demand response resources. Demand response resources, interruptible loads, distributed generation units, and conventional thermal generation units are used to provide required reserve for compensating stochastic nature of various resources and loads. In the presented model, resources connected to distribution network can participate in wholesale market through aggregators. Moreover, a novel three-program model which can be used by aggregators is presented in this article. Loads and distributed generation can contract with aggregators by these programs. A three-bus test system and the IEEE RTS are used to illustrate usefulness of the presented model. The results show that ISO can manage the system effectively by using this model

Mapping of the stochastic Lotka-Volterra model to models of population genetics and game theory

Science.gov (United States)

Constable, George W. A.; McKane, Alan J.

2017-08-01

The relationship between the M -species stochastic Lotka-Volterra competition (SLVC) model and the M -allele Moran model of population genetics is explored via timescale separation arguments. When selection for species is weak and the population size is large but finite, precise conditions are determined for the stochastic dynamics of the SLVC model to be mappable to the neutral Moran model, the Moran model with frequency-independent selection, and the Moran model with frequency-dependent selection (equivalently a game-theoretic formulation of the Moran model). We demonstrate how these mappings can be used to calculate extinction probabilities and the times until a species' extinction in the SLVC model.

Characterizing economic trends by Bayesian stochastic model specification search

DEFF Research Database (Denmark)

Grassi, Stefano; Proietti, Tommaso

We extend a recently proposed Bayesian model selection technique, known as stochastic model specification search, for characterising the nature of the trend in macroeconomic time series. In particular, we focus on autoregressive models with possibly time-varying intercept and slope and decide on ...

Letter: Modeling reactive shock waves in heterogeneous solids at the continuum level with stochastic differential equations

Science.gov (United States)

Kittell, D. E.; Yarrington, C. D.; Lechman, J. B.; Baer, M. R.

2018-05-01

A new paradigm is introduced for modeling reactive shock waves in heterogeneous solids at the continuum level. Inspired by the probability density function methods from turbulent reactive flows, it is hypothesized that the unreacted material microstructures lead to a distribution of heat release rates from chemical reaction. Fluctuations in heat release, rather than velocity, are coupled to the reactive Euler equations which are then solved via the Riemann problem. A numerically efficient, one-dimensional hydrocode is used to demonstrate this new approach, and simulation results of a representative impact calculation (inert flyer into explosive target) are discussed.

Intimate Partner Violence: A Stochastic Model.

Science.gov (United States)

Guidi, Elisa; Meringolo, Patrizia; Guazzini, Andrea; Bagnoli, Franco

2017-01-01

Intimate partner violence (IPV) has been a well-studied problem in the past psychological literature, especially through its classical methodology such as qualitative, quantitative and mixed methods. This article introduces two basic stochastic models as an alternative approach to simulate the short and long-term dynamics of a couple at risk of IPV. In both models, the members of the couple may assume a finite number of states, updating them in a probabilistic way at discrete time steps. After defining the transition probabilities, we first analyze the evolution of the couple in isolation and then we consider the case in which the individuals modify their behavior depending on the perceived violence from other couples in their environment or based on the perceived informal social support. While high perceived violence in other couples may converge toward the own presence of IPV by means a gender-specific transmission, the gender differences fade-out in the case of received informal social support. Despite the simplicity of the two stochastic models, they generate results which compare well with past experimental studies about IPV and they give important practical implications for prevention intervention in this field. Copyright: © 2016 by Fabrizio Serra editore, Pisa · Roma.

Anomaly detection in random heterogeneous media Feynman-Kac formulae, stochastic homogenization and statistical inversion

CERN Document Server

Simon, Martin

2015-01-01

This monograph is concerned with the analysis and numerical solution of a stochastic inverse anomaly detection problem in electrical impedance tomography (EIT). Martin Simon studies the problem of detecting a parameterized anomaly in an isotropic, stationary and ergodic conductivity random field whose realizations are rapidly oscillating. For this purpose, he derives Feynman-Kac formulae to rigorously justify stochastic homogenization in the case of the underlying stochastic boundary value problem. The author combines techniques from the theory of partial differential equations and functional analysis with probabilistic ideas, paving the way to new mathematical theorems which may be fruitfully used in the treatment of the problem at hand. Moreover, the author proposes an efficient numerical method in the framework of Bayesian inversion for the practical solution of the stochastic inverse anomaly detection problem.   Contents Feynman-Kac formulae Stochastic homogenization Statistical inverse problems  Targe...

Realizations of highly heterogeneous collagen networks via stochastic reconstruction for micromechanical analysis of tumor cell invasion

Science.gov (United States)

Nan, Hanqing; Liang, Long; Chen, Guo; Liu, Liyu; Liu, Ruchuan; Jiao, Yang

2018-03-01

Three-dimensional (3D) collective cell migration in a collagen-based extracellular matrix (ECM) is among one of the most significant topics in developmental biology, cancer progression, tissue regeneration, and immune response. Recent studies have suggested that collagen-fiber mediated force transmission in cellularized ECM plays an important role in stress homeostasis and regulation of collective cellular behaviors. Motivated by the recent in vitro observation that oriented collagen can significantly enhance the penetration of migrating breast cancer cells into dense Matrigel which mimics the intravasation process in vivo [Han et al. Proc. Natl. Acad. Sci. USA 113, 11208 (2016), 10.1073/pnas.1610347113], we devise a procedure for generating realizations of highly heterogeneous 3D collagen networks with prescribed microstructural statistics via stochastic optimization. Specifically, a collagen network is represented via the graph (node-bond) model and the microstructural statistics considered include the cross-link (node) density, valence distribution, fiber (bond) length distribution, as well as fiber orientation distribution. An optimization problem is formulated in which the objective function is defined as the squared difference between a set of target microstructural statistics and the corresponding statistics for the simulated network. Simulated annealing is employed to solve the optimization problem by evolving an initial network via random perturbations to generate realizations of homogeneous networks with randomly oriented fibers, homogeneous networks with aligned fibers, heterogeneous networks with a continuous variation of fiber orientation along a prescribed direction, as well as a binary system containing a collagen region with aligned fibers and a dense Matrigel region with randomly oriented fibers. The generation and propagation of active forces in the simulated networks due to polarized contraction of an embedded ellipsoidal cell and a small group

A stochastic-bayesian model for the fracture probability of PWR pressure vessels

Energy Technology Data Exchange (ETDEWEB)

Francisco, Alexandre S.; Duran, Jorge Alberto R., E-mail: afrancisco@metal.eeimvr.uff.br, E-mail: duran@metal.eeimvr.uff.br [Universidade Federal Fluminense (UFF), Volta Redonda, RJ (Brazil). Dept. de Engenharia Mecanica

2013-07-01

Fracture probability of pressure vessels containing cracks can be obtained by methodologies of easy understanding, which require a deterministic treatment, complemented by statistical methods. However, more accurate results are required, methodologies need to be better formulated. This paper presents a new methodology to address this problem. First, a more rigorous methodology is obtained by means of the relationship of probability distributions that model crack incidence and nondestructive inspection efficiency using the Bayes' theorem. The result is an updated crack incidence distribution. Further, the accuracy of the methodology is improved by using a stochastic model for the crack growth. The stochastic model incorporates the statistical variability of the crack growth process, combining the stochastic theory with experimental data. Stochastic differential equations are derived by the randomization of empirical equations. From the solution of this equation, a distribution function related to the crack growth is derived. The fracture probability using both probability distribution functions is in agreement with theory, and presents realistic value for pressure vessels. (author)

A stochastic-bayesian model for the fracture probability of PWR pressure vessels

International Nuclear Information System (INIS)

Francisco, Alexandre S.; Duran, Jorge Alberto R.

2013-01-01

Fracture probability of pressure vessels containing cracks can be obtained by methodologies of easy understanding, which require a deterministic treatment, complemented by statistical methods. However, more accurate results are required, methodologies need to be better formulated. This paper presents a new methodology to address this problem. First, a more rigorous methodology is obtained by means of the relationship of probability distributions that model crack incidence and nondestructive inspection efficiency using the Bayes' theorem. The result is an updated crack incidence distribution. Further, the accuracy of the methodology is improved by using a stochastic model for the crack growth. The stochastic model incorporates the statistical variability of the crack growth process, combining the stochastic theory with experimental data. Stochastic differential equations are derived by the randomization of empirical equations. From the solution of this equation, a distribution function related to the crack growth is derived. The fracture probability using both probability distribution functions is in agreement with theory, and presents realistic value for pressure vessels. (author)

Estimation of Dynamic Panel Data Models with Stochastic Volatility Using Particle Filters

Directory of Open Access Journals (Sweden)

Wen Xu

2016-10-01

Full Text Available Time-varying volatility is common in macroeconomic data and has been incorporated into macroeconomic models in recent work. Dynamic panel data models have become increasingly popular in macroeconomics to study common relationships across countries or regions. This paper estimates dynamic panel data models with stochastic volatility by maximizing an approximate likelihood obtained via Rao-Blackwellized particle filters. Monte Carlo studies reveal the good and stable performance of our particle filter-based estimator. When the volatility of volatility is high, or when regressors are absent but stochastic volatility exists, our approach can be better than the maximum likelihood estimator which neglects stochastic volatility and generalized method of moments (GMM estimators.

Stochastic processes

CERN Document Server

Parzen, Emanuel

1962-01-01

Well-written and accessible, this classic introduction to stochastic processes and related mathematics is appropriate for advanced undergraduate students of mathematics with a knowledge of calculus and continuous probability theory. The treatment offers examples of the wide variety of empirical phenomena for which stochastic processes provide mathematical models, and it develops the methods of probability model-building.Chapter 1 presents precise definitions of the notions of a random variable and a stochastic process and introduces the Wiener and Poisson processes. Subsequent chapters examine

Development of stochastic indicator models of lithology, Yucca Mountain, Nevada

International Nuclear Information System (INIS)

Rautman, C.A.; Robey, T.H.

1994-01-01

Indicator geostatistical techniques have been used to produce a number of fully three-dimensional stochastic simulations of large-scale lithologic categories at the Yucca Mountain site. Each realization reproduces the available drill hole data used to condition the simulation. Information is propagated away from each point of observation in accordance with a mathematical model of spatial continuity inferred through soft data taken from published geologic cross sections. Variations among the simulated models collectively represent uncertainty in the lithology at unsampled locations. These stochastic models succeed in capturing many major features of welded-nonwelded lithologic framework of Yucca Mountain. However, contacts between welded and nonwelded rock types for individual simulations appear more complex than suggested by field observation, and a number of probable numerical artifacts exist in these models. Many of the apparent discrepancies between the simulated models and the general geology of Yucca Mountain represent characterization uncertainty, and can be traced to the sparse site data used to condition the simulations. Several vertical stratigraphic columns have been extracted from the three-dimensional stochastic models for use in simplified total-system performance assessment exercises. Simple, manual adjustments are required to eliminate the more obvious simulation artifacts and to impose a secondary set of deterministic geologic features on the overall stratigraphic framework provided by the indictor models

Modelling the stochastic behaviour of primary nucleation.

Science.gov (United States)

Maggioni, Giovanni Maria; Mazzotti, Marco

2015-01-01

We study the stochastic nature of primary nucleation and how it manifests itself in a crystallisation process at different scales and under different operating conditions. Such characteristics of nucleation are evident in many experiments where detection times of crystals are not identical, despite identical experimental conditions, but instead are distributed around an average value. While abundant experimental evidence has been reported in the literature, a clear theoretical understanding and an appropriate modelling of this feature is still missing. In this contribution, we present two models describing a batch cooling crystallisation, where the interplay between stochastic nucleation and deterministic crystal growth is described differently in each. The nucleation and growth rates of the two models are estimated by a comprehensive set of measurements of paracetamol crystallisation from aqueous solution in a 1 mL vessel [Kadam et al., Chemical Engineering Science, 2012, 72, 10-19]. Both models are applied to the cooling crystallisation process above under different operating conditions, i.e. different volumes, initial concentrations, cooling rates. The advantages and disadvantages of the two approaches are illustrated and discussed, with particular reference to their use across scales of nucleation rate measured in very small crystallisers.

Persistence and extinction for stochastic logistic model with Levy noise and impulsive perturbation

Directory of Open Access Journals (Sweden)

Chun Lu

2015-09-01

Full Text Available This article investigates a stochastic logistic model with Levy noise and impulsive perturbation. In the model, the impulsive perturbation and Levy noise are taken into account simultaneously. This model is new and more feasible and more accordance with the actual. The definition of solution to a stochastic differential equation with Levy noise and impulsive perturbation is established. Based on this definition, we show that our model has a unique global positive solution and obtains its explicit expression. Sufficient conditions for extinction are established as well as nonpersistence in the mean, weak persistence and stochastic permanence. The threshold between weak persistence and extinction is obtained.

Discontinuous finite element and characteristics methods for neutrons transport equation solution in heterogeneous grids; Resolution de l'equation du transport des neutrons par les methodes des elements finis discontinus et des caracteristiques structurees appliquees a des maillages heterogenes

Energy Technology Data Exchange (ETDEWEB)

Masiello, E

2006-07-01

The principal goal of this manuscript is devoted to the investigation of a new type of heterogeneous mesh adapted to the shape of the fuel pins (fuel-clad-moderator). The new heterogeneous mesh guarantees the spatial modelling of the pin-cell with a minimum of regions. Two methods are investigated for the spatial discretization of the transport equation: the discontinuous finite element method and the method of characteristics for structured cells. These methods together with the new representation of the pin-cell result in an appreciable reduction of calculation points. They allow an exact modelling of the fuel pin-cell without spatial homogenization. A new synthetic acceleration technique based on an angular multigrid is also presented for the speed up of the inner iterations. These methods are good candidates for transport calculations for a nuclear reactor core. A second objective of this work is the application of method of characteristics for non-structured geometries to the study of double heterogeneity problem. The letters is characterized by fuel material with a stochastic dispersion of heterogeneous grains, and until now was solved with a model based on collision probabilities. We propose a new statistical model based on renewal-Markovian theory, which makes possible to take into account the stochastic nature of the problem and to avoid the approximations of the collision probability model. The numerical solution of this model is guaranteed by the method of characteristics. (author)

Stochastic modeling of central apnea events in preterm infants

International Nuclear Information System (INIS)

Clark, Matthew T; Lake, Douglas E; Randall Moorman, J; Delos, John B; Lee, Hoshik; Fairchild, Karen D; Kattwinkel, John

2016-01-01

A near-ubiquitous pathology in very low birth weight infants is neonatal apnea, breathing pauses with slowing of the heart and falling blood oxygen. Events of substantial duration occasionally occur after an infant is discharged from the neonatal intensive care unit (NICU). It is not known whether apneas result from a predictable process or from a stochastic process, but the observation that they occur in seemingly random clusters justifies the use of stochastic models. We use a hidden-Markov model to analyze the distribution of durations of apneas and the distribution of times between apneas. The model suggests the presence of four breathing states, ranging from very stable (with an average lifetime of 12 h) to very unstable (with an average lifetime of 10 s). Although the states themselves are not visible, the mathematical analysis gives estimates of the transition rates among these states. We have obtained these transition rates, and shown how they change with post-menstrual age; as expected, the residence time in the more stable breathing states increases with age. We also extrapolated the model to predict the frequency of very prolonged apnea during the first year of life. This paradigm—stochastic modeling of cardiorespiratory control in neonatal infants to estimate risk for severe clinical events—may be a first step toward personalized risk assessment for life threatening apnea events after NICU discharge. (paper)

Stochastic modeling of central apnea events in preterm infants.

Science.gov (United States)

Clark, Matthew T; Delos, John B; Lake, Douglas E; Lee, Hoshik; Fairchild, Karen D; Kattwinkel, John; Moorman, J Randall

2016-04-01

A near-ubiquitous pathology in very low birth weight infants is neonatal apnea, breathing pauses with slowing of the heart and falling blood oxygen. Events of substantial duration occasionally occur after an infant is discharged from the neonatal intensive care unit (NICU). It is not known whether apneas result from a predictable process or from a stochastic process, but the observation that they occur in seemingly random clusters justifies the use of stochastic models. We use a hidden-Markov model to analyze the distribution of durations of apneas and the distribution of times between apneas. The model suggests the presence of four breathing states, ranging from very stable (with an average lifetime of 12 h) to very unstable (with an average lifetime of 10 s). Although the states themselves are not visible, the mathematical analysis gives estimates of the transition rates among these states. We have obtained these transition rates, and shown how they change with post-menstrual age; as expected, the residence time in the more stable breathing states increases with age. We also extrapolated the model to predict the frequency of very prolonged apnea during the first year of life. This paradigm-stochastic modeling of cardiorespiratory control in neonatal infants to estimate risk for severe clinical events-may be a first step toward personalized risk assessment for life threatening apnea events after NICU discharge.

Analysis of dynamic regimes in stochastically forced Kaldor model

International Nuclear Information System (INIS)

Bashkirtseva, Irina; Ryazanova, Tatyana; Ryashko, Lev

2015-01-01

We consider the business cycle Kaldor model forced by random noise. Detailed parametric analysis of deterministic system is carried out and zones of coexisting stable equilibrium and stable limit cycle are found. Noise-induced transitions between these attractors are studied using stochastic sensitivity function technique and confidence domains method. Critical values of noise intensity corresponding to noise-induced transitions “equilibrium → cycle” and “cycle → equilibrium” are estimated. Dominants in combined stochastic regimes are discussed.

Live imaging-based model selection reveals periodic regulation of the stochastic G1/S phase transition in vertebrate axial development.

Directory of Open Access Journals (Sweden)

Mayu Sugiyama

2014-12-01

Full Text Available In multicellular organism development, a stochastic cellular response is observed, even when a population of cells is exposed to the same environmental conditions. Retrieving the spatiotemporal regulatory mode hidden in the heterogeneous cellular behavior is a challenging task. The G1/S transition observed in cell cycle progression is a highly stochastic process. By taking advantage of a fluorescence cell cycle indicator, Fucci technology, we aimed to unveil a hidden regulatory mode of cell cycle progression in developing zebrafish. Fluorescence live imaging of Cecyil, a zebrafish line genetically expressing Fucci, demonstrated that newly formed notochordal cells from the posterior tip of the embryonic mesoderm exhibited the red (G1 fluorescence signal in the developing notochord. Prior to their initial vacuolation, these cells showed a fluorescence color switch from red to green, indicating G1/S transitions. This G1/S transition did not occur in a synchronous manner, but rather exhibited a stochastic process, since a mixed population of red and green cells was always inserted between newly formed red (G1 notochordal cells and vacuolating green cells. We termed this mixed population of notochordal cells, the G1/S transition window. We first performed quantitative analyses of live imaging data and a numerical estimation of the probability of the G1/S transition, which demonstrated the existence of a posteriorly traveling regulatory wave of the G1/S transition window. To obtain a better understanding of this regulatory mode, we constructed a mathematical model and performed a model selection by comparing the results obtained from the models with those from the experimental data. Our analyses demonstrated that the stochastic G1/S transition window in the notochord travels posteriorly in a periodic fashion, with doubled the periodicity of the neighboring paraxial mesoderm segmentation. This approach may have implications for the characterization of

Modeling animal movements using stochastic differential equations

Science.gov (United States)

Haiganoush K. Preisler; Alan A. Ager; Bruce K. Johnson; John G. Kie

2004-01-01

We describe the use of bivariate stochastic differential equations (SDE) for modeling movements of 216 radiocollared female Rocky Mountain elk at the Starkey Experimental Forest and Range in northeastern Oregon. Spatially and temporally explicit vector fields were estimated using approximating difference equations and nonparametric regression techniques. Estimated...

Modeling reliability of power systems substations by using stochastic automata networks

International Nuclear Information System (INIS)

Šnipas, Mindaugas; Radziukynas, Virginijus; Valakevičius, Eimutis

2017-01-01

In this paper, stochastic automata networks (SANs) formalism to model reliability of power systems substations is applied. The proposed strategy allows reducing the size of state space of Markov chain model and simplifying system specification. Two case studies of standard configurations of substations are considered in detail. SAN models with different assumptions were created. SAN approach is compared with exact reliability calculation by using a minimal path set method. Modeling results showed that total independence of automata can be assumed for relatively small power systems substations with reliable equipment. In this case, the implementation of Markov chain model by a using SAN method is a relatively easy task. - Highlights: • We present the methodology to apply stochastic automata network formalism to create Markov chain models of power systems. • The stochastic automata network approach is combined with minimal path sets and structural functions. • Two models of substation configurations with different model assumptions are presented to illustrate the proposed methodology. • Modeling results of system with independent automata and functional transition rates are similar. • The conditions when total independence of automata can be assumed are addressed.

Development of a Stochastically-driven, Forward Predictive Performance Model for PEMFCs

Science.gov (United States)

Harvey, David Benjamin Paul

A one-dimensional multi-scale coupled, transient, and mechanistic performance model for a PEMFC membrane electrode assembly has been developed. The model explicitly includes each of the 5 layers within a membrane electrode assembly and solves for the transport of charge, heat, mass, species, dissolved water, and liquid water. Key features of the model include the use of a multi-step implementation of the HOR reaction on the anode, agglomerate catalyst sub-models for both the anode and cathode catalyst layers, a unique approach that links the composition of the catalyst layer to key properties within the agglomerate model and the implementation of a stochastic input-based approach for component material properties. The model employs a new methodology for validation using statistically varying input parameters and statistically-based experimental performance data; this model represents the first stochastic input driven unit cell performance model. The stochastic input driven performance model was used to identify optimal ionomer content within the cathode catalyst layer, demonstrate the role of material variation in potential low performing MEA materials, provide explanation for the performance of low-Pt loaded MEAs, and investigate the validity of transient-sweep experimental diagnostic methods.

Information-theoretic model selection for optimal prediction of stochastic dynamical systems from data

Science.gov (United States)

Darmon, David

2018-03-01

In the absence of mechanistic or phenomenological models of real-world systems, data-driven models become necessary. The discovery of various embedding theorems in the 1980s and 1990s motivated a powerful set of tools for analyzing deterministic dynamical systems via delay-coordinate embeddings of observations of their component states. However, in many branches of science, the condition of operational determinism is not satisfied, and stochastic models must be brought to bear. For such stochastic models, the tool set developed for delay-coordinate embedding is no longer appropriate, and a new toolkit must be developed. We present an information-theoretic criterion, the negative log-predictive likelihood, for selecting the embedding dimension for a predictively optimal data-driven model of a stochastic dynamical system. We develop a nonparametric estimator for the negative log-predictive likelihood and compare its performance to a recently proposed criterion based on active information storage. Finally, we show how the output of the model selection procedure can be used to compare candidate predictors for a stochastic system to an information-theoretic lower bound.

Toward Development of a Stochastic Wake Model: Validation Using LES and Turbine Loads

Directory of Open Access Journals (Sweden)

Jae Sang Moon

2017-12-01

Full Text Available Wind turbines within an array do not experience free-stream undisturbed flow fields. Rather, the flow fields on internal turbines are influenced by wakes generated by upwind unit and exhibit different dynamic characteristics relative to the free stream. The International Electrotechnical Commission (IEC standard 61400-1 for the design of wind turbines only considers a deterministic wake model for the design of a wind plant. This study is focused on the development of a stochastic model for waked wind fields. First, high-fidelity physics-based waked wind velocity fields are generated using Large-Eddy Simulation (LES. Stochastic characteristics of these LES waked wind velocity field, including mean and turbulence components, are analyzed. Wake-related mean and turbulence field-related parameters are then estimated for use with a stochastic model, using Multivariate Multiple Linear Regression (MMLR with the LES data. To validate the simulated wind fields based on the stochastic model, wind turbine tower and blade loads are generated using aeroelastic simulation for utility-scale wind turbine models and compared with those based directly on the LES inflow. The study’s overall objective is to offer efficient and validated stochastic approaches that are computationally tractable for assessing the performance and loads of turbines operating in wakes.

Stochastic modeling of thermal fatigue crack growth

CERN Document Server

Radu, Vasile

2015-01-01

The book describes a systematic stochastic modeling approach for assessing thermal-fatigue crack-growth in mixing tees, based on the power spectral density of temperature fluctuation at the inner pipe surface. It shows the development of a frequency-temperature response function in the framework of single-input, single-output (SISO) methodology from random noise/signal theory under sinusoidal input. The frequency response of stress intensity factor (SIF) is obtained by a polynomial fitting procedure of thermal stress profiles at various instants of time. The method, which takes into account the variability of material properties, and has been implemented in a real-world application, estimates the probabilities of failure by considering a limit state function and Monte Carlo analysis, which are based on the proposed stochastic model. Written in a comprehensive and accessible style, this book presents a new and effective method for assessing thermal fatigue crack, and it is intended as a concise and practice-or...

Fuzzy Stochastic Optimization Theory, Models and Applications

CERN Document Server

Wang, Shuming

2012-01-01

Covering in detail both theoretical and practical perspectives, this book is a self-contained and systematic depiction of current fuzzy stochastic optimization that deploys the fuzzy random variable as a core mathematical tool to model the integrated fuzzy random uncertainty. It proceeds in an orderly fashion from the requisite theoretical aspects of the fuzzy random variable to fuzzy stochastic optimization models and their real-life case studies.   The volume reflects the fact that randomness and fuzziness (or vagueness) are two major sources of uncertainty in the real world, with significant implications in a number of settings. In industrial engineering, management and economics, the chances are high that decision makers will be confronted with information that is simultaneously probabilistically uncertain and fuzzily imprecise, and optimization in the form of a decision must be made in an environment that is doubly uncertain, characterized by a co-occurrence of randomness and fuzziness. This book begins...

Powering stochastic reliability models by discrete event simulation

DEFF Research Database (Denmark)

Kozine, Igor; Wang, Xiaoyun

2012-01-01

it difficult to find a solution to the problem. The power of modern computers and recent developments in discrete-event simulation (DES) software enable to diminish some of the drawbacks of stochastic models. In this paper we describe the insights we have gained based on using both Markov and DES models...

Portfolio Optimization with Stochastic Dividends and Stochastic Volatility

Science.gov (United States)

Varga, Katherine Yvonne

2015-01-01

We consider an optimal investment-consumption portfolio optimization model in which an investor receives stochastic dividends. As a first problem, we allow the drift of stock price to be a bounded function. Next, we consider a stochastic volatility model. In each problem, we use the dynamic programming method to derive the Hamilton-Jacobi-Bellman…

Response spectrum analysis of a stochastic seismic model

International Nuclear Information System (INIS)

Kimura, Koji; Sakata, Masaru; Takemoto, Shinichiro.

1990-01-01

The stochastic response spectrum approach is presented for predicting the dynamic behavior of structures to earthquake excitation expressed by a random process, one of whose sample functions can be regarded as a recorded strong-motion earthquake accelerogram. The approach consists of modeling recorded ground motion by a random process and the root-mean-square response (rms) analysis of a single-degree-of-freedom system by using the moment equations method. The stochastic response spectrum is obtained as a plot of the maximum rms response versus the natural period of the system and is compared with the conventional response spectrum. (author)

Discrete stochastic analogs of Erlang epidemic models.

Science.gov (United States)

Getz, Wayne M; Dougherty, Eric R

2018-12-01

Erlang differential equation models of epidemic processes provide more realistic disease-class transition dynamics from susceptible (S) to exposed (E) to infectious (I) and removed (R) categories than the ubiquitous SEIR model. The latter is itself is at one end of the spectrum of Erlang SE[Formula: see text]I[Formula: see text]R models with [Formula: see text] concatenated E compartments and [Formula: see text] concatenated I compartments. Discrete-time models, however, are computationally much simpler to simulate and fit to epidemic outbreak data than continuous-time differential equations, and are also much more readily extended to include demographic and other types of stochasticity. Here we formulate discrete-time deterministic analogs of the Erlang models, and their stochastic extension, based on a time-to-go distributional principle. Depending on which distributions are used (e.g. discretized Erlang, Gamma, Beta, or Uniform distributions), we demonstrate that our formulation represents both a discretization of Erlang epidemic models and generalizations thereof. We consider the challenges of fitting SE[Formula: see text]I[Formula: see text]R models and our discrete-time analog to data (the recent outbreak of Ebola in Liberia). We demonstrate that the latter performs much better than the former; although confining fits to strict SEIR formulations reduces the numerical challenges, but sacrifices best-fit likelihood scores by at least 7%.

Discontinuous finite element and characteristics methods for neutrons transport equation solution in heterogeneous grids; Resolution de l'equation du transport des neutrons par les methodes des elements finis discontinus et des caracteristiques structurees appliquees a des maillages heterogenes

Energy Technology Data Exchange (ETDEWEB)

Masiello, E

2006-07-01

The principal goal of this manuscript is devoted to the investigation of a new type of heterogeneous mesh adapted to the shape of the fuel pins (fuel-clad-moderator). The new heterogeneous mesh guarantees the spatial modelling of the pin-cell with a minimum of regions. Two methods are investigated for the spatial discretization of the transport equation: the discontinuous finite element method and the method of characteristics for structured cells. These methods together with the new representation of the pin-cell result in an appreciable reduction of calculation points. They allow an exact modelling of the fuel pin-cell without spatial homogenization. A new synthetic acceleration technique based on an angular multigrid is also presented for the speed up of the inner iterations. These methods are good candidates for transport calculations for a nuclear reactor core. A second objective of this work is the application of method of characteristics for non-structured geometries to the study of double heterogeneity problem. The letters is characterized by fuel material with a stochastic dispersion of heterogeneous grains, and until now was solved with a model based on collision probabilities. We propose a new statistical model based on renewal-Markovian theory, which makes possible to take into account the stochastic nature of the problem and to avoid the approximations of the collision probability model. The numerical solution of this model is guaranteed by the method of characteristics. (author)

Stochastic modeling of oligodendrocyte generation in cell culture: model validation with time-lapse data

Directory of Open Access Journals (Sweden)

Noble Mark

2006-05-01

Full Text Available Abstract Background The purpose of this paper is two-fold. The first objective is to validate the assumptions behind a stochastic model developed earlier by these authors to describe oligodendrocyte generation in cell culture. The second is to generate time-lapse data that may help biomathematicians to build stochastic models of cell proliferation and differentiation under other experimental scenarios. Results Using time-lapse video recording it is possible to follow the individual evolutions of different cells within each clone. This experimental technique is very laborious and cannot replace model-based quantitative inference from clonal data. However, it is unrivalled in validating the structure of a stochastic model intended to describe cell proliferation and differentiation at the clonal level. In this paper, such data are reported and analyzed for oligodendrocyte precursor cells cultured in vitro. Conclusion The results strongly support the validity of the most basic assumptions underpinning the previously proposed model of oligodendrocyte development in cell culture. However, there are some discrepancies; the most important is that the contribution of progenitor cell death to cell kinetics in this experimental system has been underestimated.

Bond and CDS Pricing via the Stochastic Recovery Black-Cox Model

Directory of Open Access Journals (Sweden)

Albert Cohen

2017-04-01

Full Text Available Building on recent work incorporating recovery risk into structural models by Cohen & Costanzino (2015, we consider the Black-Cox model with an added recovery risk driver. The recovery risk driver arises naturally in the context of imperfect information implicit in the structural framework. This leads to a two-factor structural model we call the Stochastic Recovery Black-Cox model, whereby the asset risk driver At defines the default trigger and the recovery risk driver Rt defines the amount recovered in the event of default. We then price zero-coupon bonds and credit default swaps under the Stochastic Recovery Black-Cox model. Finally, we compare our results with the classic Black-Cox model, give explicit expressions for the recovery risk premium in the Stochastic Recovery Black-Cox model, and detail how the introduction of separate but correlated risk drivers leads to a decoupling of the default and recovery risk premiums in the credit spread. We conclude this work by computing the effect of adding coupons that are paid continuously until default, and price perpetual (consol bonds in our two-factor firm value model, extending calculations in the seminal paper by Leland (1994.

Stochastic modeling of pitting corrosion: A new model for initiation and growth of multiple corrosion pits

International Nuclear Information System (INIS)

Valor, A.; Caleyo, F.; Alfonso, L.; Rivas, D.; Hallen, J.M.

2007-01-01

In this work, a new stochastic model capable of simulating pitting corrosion is developed and validated. Pitting corrosion is modeled as the combination of two stochastic processes: pit initiation and pit growth. Pit generation is modeled as a nonhomogeneous Poisson process, in which induction time for pit initiation is simulated as the realization of a Weibull process. In this way, the exponential and Weibull distributions can be considered as the possible distributions for pit initiation time. Pit growth is simulated using a nonhomogeneous Markov process. Extreme value statistics is used to find the distribution of maximum pit depths resulting from the combination of the initiation and growth processes for multiple pits. The proposed model is validated using several published experiments on pitting corrosion. It is capable of reproducing the experimental observations with higher quality than the stochastic models available in the literature for pitting corrosion

Stochastic modeling of pitting corrosion: A new model for initiation and growth of multiple corrosion pits

Energy Technology Data Exchange (ETDEWEB)

Valor, A. [Facultad de Fisica, Universidad de La Habana, San Lazaro y L, Vedado, 10400 Havana (Cuba); Caleyo, F. [Departamento de Ingenieria, Metalurgica, IPN-ESIQIE, UPALM Edif. 7, Zacatenco, Mexico DF 07738 (Mexico)]. E-mail: fcaleyo@gmail.com; Alfonso, L. [Departamento de Ingenieria, Metalurgica, IPN-ESIQIE, UPALM Edif. 7, Zacatenco, Mexico DF 07738 (Mexico); Rivas, D. [Departamento de Ingenieria, Metalurgica, IPN-ESIQIE, UPALM Edif. 7, Zacatenco, Mexico DF 07738 (Mexico); Hallen, J.M. [Departamento de Ingenieria, Metalurgica, IPN-ESIQIE, UPALM Edif. 7, Zacatenco, Mexico DF 07738 (Mexico)

2007-02-15

In this work, a new stochastic model capable of simulating pitting corrosion is developed and validated. Pitting corrosion is modeled as the combination of two stochastic processes: pit initiation and pit growth. Pit generation is modeled as a nonhomogeneous Poisson process, in which induction time for pit initiation is simulated as the realization of a Weibull process. In this way, the exponential and Weibull distributions can be considered as the possible distributions for pit initiation time. Pit growth is simulated using a nonhomogeneous Markov process. Extreme value statistics is used to find the distribution of maximum pit depths resulting from the combination of the initiation and growth processes for multiple pits. The proposed model is validated using several published experiments on pitting corrosion. It is capable of reproducing the experimental observations with higher quality than the stochastic models available in the literature for pitting corrosion.

Fitting Nonlinear Ordinary Differential Equation Models with Random Effects and Unknown Initial Conditions Using the Stochastic Approximation Expectation-Maximization (SAEM) Algorithm.

Science.gov (United States)

Chow, Sy-Miin; Lu, Zhaohua; Sherwood, Andrew; Zhu, Hongtu

2016-03-01

The past decade has evidenced the increased prevalence of irregularly spaced longitudinal data in social sciences. Clearly lacking, however, are modeling tools that allow researchers to fit dynamic models to irregularly spaced data, particularly data that show nonlinearity and heterogeneity in dynamical structures. We consider the issue of fitting multivariate nonlinear differential equation models with random effects and unknown initial conditions to irregularly spaced data. A stochastic approximation expectation-maximization algorithm is proposed and its performance is evaluated using a benchmark nonlinear dynamical systems model, namely, the Van der Pol oscillator equations. The empirical utility of the proposed technique is illustrated using a set of 24-h ambulatory cardiovascular data from 168 men and women. Pertinent methodological challenges and unresolved issues are discussed.

GillesPy: A Python Package for Stochastic Model Building and Simulation.

Science.gov (United States)

Abel, John H; Drawert, Brian; Hellander, Andreas; Petzold, Linda R

2016-09-01

GillesPy is an open-source Python package for model construction and simulation of stochastic biochemical systems. GillesPy consists of a Python framework for model building and an interface to the StochKit2 suite of efficient simulation algorithms based on the Gillespie stochastic simulation algorithms (SSA). To enable intuitive model construction and seamless integration into the scientific Python stack, we present an easy to understand, action-oriented programming interface. Here, we describe the components of this package and provide a detailed example relevant to the computational biology community.

Modeling and Prediction Using Stochastic Differential Equations

DEFF Research Database (Denmark)

Juhl, Rune; Møller, Jan Kloppenborg; Jørgensen, John Bagterp

2016-01-01

Pharmacokinetic/pharmakodynamic (PK/PD) modeling for a single subject is most often performed using nonlinear models based on deterministic ordinary differential equations (ODEs), and the variation between subjects in a population of subjects is described using a population (mixed effects) setup...... deterministic and can predict the future perfectly. A more realistic approach would be to allow for randomness in the model due to e.g., the model be too simple or errors in input. We describe a modeling and prediction setup which better reflects reality and suggests stochastic differential equations (SDEs...

Levy-Student processes for a stochastic model of beam halos

Energy Technology Data Exchange (ETDEWEB)

Petroni, N. Cufaro [Department of Mathematics, University of Bari, and INFN Sezione di Bari, via E. Orabona 4, 70125 Bari (Italy)]. E-mail: cufaro@ba.infn.it; De Martino, S. [Department of Physics, University of Salerno, and INFN Sezione di Napoli (gruppo di Salerno), Via S. Allende, I-84081 Baronissi (SA) (Italy); De Siena, S. [Department of Physics, University of Salerno, and INFN Sezione di Napoli (gruppo di Salerno), Via S. Allende, I-84081 Baronissi (SA) (Italy); Illuminati, F. [Department of Physics, University of Salerno, and INFN Sezione di Napoli (gruppo di Salerno), Via S. Allende, I-84081 Baronissi (SA) (Italy)

2006-06-01

We describe the transverse beam distribution in particle accelerators within the controlled, stochastic dynamical scheme of the stochastic mechanics which produces time reversal invariant diffusion processes. In this paper we analyze the consequences of introducing the generalized Student laws, namely non-Gaussian, Levy infinitely divisible (but not stable) distributions. We will analyze this idea from two different standpoints: (a) first by supposing that the stationary distribution of our (Wiener powered) stochastic model is a Student distribution; (b) by supposing that our model is based on a (non-Gaussian) Levy process whose increments are Student distributed. In the case (a) the longer tails of the power decay of the Student laws, and in the case (b) the discontinuities of the Levy-Student process can well account for the rare escape of particles from the beam core, and hence for the formation of a halo in intense beams.

Levy-Student processes for a stochastic model of beam halos

International Nuclear Information System (INIS)

Petroni, N. Cufaro; De Martino, S.; De Siena, S.; Illuminati, F.

2006-01-01

We describe the transverse beam distribution in particle accelerators within the controlled, stochastic dynamical scheme of the stochastic mechanics which produces time reversal invariant diffusion processes. In this paper we analyze the consequences of introducing the generalized Student laws, namely non-Gaussian, Levy infinitely divisible (but not stable) distributions. We will analyze this idea from two different standpoints: (a) first by supposing that the stationary distribution of our (Wiener powered) stochastic model is a Student distribution; (b) by supposing that our model is based on a (non-Gaussian) Levy process whose increments are Student distributed. In the case (a) the longer tails of the power decay of the Student laws, and in the case (b) the discontinuities of the Levy-Student process can well account for the rare escape of particles from the beam core, and hence for the formation of a halo in intense beams

Assessing Exhaustiveness of Stochastic Sampling for Integrative Modeling of Macromolecular Structures.

Science.gov (United States)

Viswanath, Shruthi; Chemmama, Ilan E; Cimermancic, Peter; Sali, Andrej

2017-12-05

Modeling of macromolecular structures involves structural sampling guided by a scoring function, resulting in an ensemble of good-scoring models. By necessity, the sampling is often stochastic, and must be exhaustive at a precision sufficient for accurate modeling and assessment of model uncertainty. Therefore, the very first step in analyzing the ensemble is an estimation of the highest precision at which the sampling is exhaustive. Here, we present an objective and automated method for this task. As a proxy for sampling exhaustiveness, we evaluate whether two independently and stochastically generated sets of models are sufficiently similar. The protocol includes testing 1) convergence of the model score, 2) whether model scores for the two samples were drawn from the same parent distribution, 3) whether each structural cluster includes models from each sample proportionally to its size, and 4) whether there is sufficient structural similarity between the two model samples in each cluster. The evaluation also provides the sampling precision, defined as the smallest clustering threshold that satisfies the third, most stringent test. We validate the protocol with the aid of enumerated good-scoring models for five illustrative cases of binary protein complexes. Passing the proposed four tests is necessary, but not sufficient for thorough sampling. The protocol is general in nature and can be applied to the stochastic sampling of any set of models, not just structural models. In addition, the tests can be used to stop stochastic sampling as soon as exhaustiveness at desired precision is reached, thereby improving sampling efficiency; they may also help in selecting a model representation that is sufficiently detailed to be informative, yet also sufficiently coarse for sampling to be exhaustive. Copyright © 2017 Biophysical Society. Published by Elsevier Inc. All rights reserved.

AUTOMATIC CALIBRATION OF A STOCHASTIC-LAGRANGIAN TRANSPORT MODEL (SLAM)

Science.gov (United States)

Numerical models are a useful tool in evaluating and designing NAPL remediation systems. Traditional constitutive finite difference and finite element models are complex and expensive to apply. For this reason, this paper presents the application of a simplified stochastic-Lagran...

Stochastic programming framework for Lithuanian pension payout modelling

Directory of Open Access Journals (Sweden)

Audrius Kabašinskas

2014-12-01

Full Text Available The paper provides a scientific approach to the problem of selecting a pension fund by taking into account some specific characteristics of the Lithuanian Republic (LR pension accumulation system. The decision making model, which can be used to plan a long-term pension accrual of the Lithuanian Republic (LR citizens, in an optimal way is presented. This model focuses on factors that influence the sustainability of the pension system selection under macroeconomic, social and demographic uncertainty. The model is formalized as a single stage stochastic optimization problem where the long-term optimal strategy can be obtained based on the possible scenarios generated for a particular participant. Stochastic programming methods allow including the pension fund rebalancing moment and direction of investment, and taking into account possible changes of personal income, changes of society and the global financial market. The collection of methods used to generate scenario trees was found useful to solve strategic planning problems.

The water budget of heterogeneous areas : impact of soil and rainfall variability

NARCIS (Netherlands)

Kim, C.P.

1995-01-01

In this thesis the heterogeneity of the soil water budget components is investigated. Heterogeneity of soil hydraulic properties and rainfall rate are taken into account by using stochastic methods. The importance of lateral groundwater flow in causing heterogeneity of the water budget

The threshold of a stochastic avian-human influenza epidemic model with psychological effect

Science.gov (United States)

Zhang, Fengrong; Zhang, Xinhong

2018-02-01

In this paper, a stochastic avian-human influenza epidemic model with psychological effect in human population and saturation effect within avian population is investigated. This model describes the transmission of avian influenza among avian population and human population in random environments. For stochastic avian-only system, persistence in the mean and extinction of the infected avian population are studied. For the avian-human influenza epidemic system, sufficient conditions for the existence of an ergodic stationary distribution are obtained. Furthermore, a threshold of this stochastic model which determines the outcome of the disease is obtained. Finally, numerical simulations are given to support the theoretical results.

Deterministic and stochastic trends in the Lee-Carter mortality model

DEFF Research Database (Denmark)

Callot, Laurent; Haldrup, Niels; Kallestrup-Lamb, Malene

The Lee and Carter (1992) model assumes that the deterministic and stochastic time series dynamics loads with identical weights when describing the development of age specific mortality rates. Effectively this means that the main characteristics of the model simplifies to a random walk model...... that characterizes mortality data. We find empirical evidence that this feature of the Lee-Carter model overly restricts the system dynamics and we suggest to separate the deterministic and stochastic time series components at the benefit of improved fit and forecasting performance. In fact, we find...... that the classical Lee-Carter model will otherwise over estimate the reduction of mortality for the younger age groups and will under estimate the reduction of mortality for the older age groups. In practice, our recommendation means that the Lee-Carter model instead of a one-factor model should be formulated...

Modeling stochastic frontier based on vine copulas

Science.gov (United States)

Constantino, Michel; Candido, Osvaldo; Tabak, Benjamin M.; da Costa, Reginaldo Brito

2017-11-01

This article models a production function and analyzes the technical efficiency of listed companies in the United States, Germany and England between 2005 and 2012 based on the vine copula approach. Traditional estimates of the stochastic frontier assume that data is multivariate normally distributed and there is no source of asymmetry. The proposed method based on vine copulas allow us to explore different types of asymmetry and multivariate distribution. Using data on product, capital and labor, we measure the relative efficiency of the vine production function and estimate the coefficient used in the stochastic frontier literature for comparison purposes. This production vine copula predicts the value added by firms with given capital and labor in a probabilistic way. It thereby stands in sharp contrast to the production function, where the output of firms is completely deterministic. The results show that, on average, S&P500 companies are more efficient than companies listed in England and Germany, which presented similar average efficiency coefficients. For comparative purposes, the traditional stochastic frontier was estimated and the results showed discrepancies between the coefficients obtained by the application of the two methods, traditional and frontier-vine, opening new paths of non-linear research.

An Innovative Real-time Environment for Unified Deterministic and Stochastic Groundwater Modeling

Science.gov (United States)

Li, S.; Liu, Q.

2003-12-01

Despite an exponential growth of computational capability over the last two decades-one that has allowed computational science and engineering to become a unique, powerful tool for scientific discovery-the extreme cost of groundwater modeling continues to limit its use. This occurs primarily because the modeling paradigm that has been employed for decades limits our ability to take full advantage of recent developments in computer, communication, graphic, and visualization technologies. In this presentation we introduce an innovative and sophisticated computational environment for groundwater modeling that promises to eliminate the current bottleneck and greatly expand the utility of computational tools for scientific discovery related to groundwater. Based on a set of efficient and robust computational algorithms, the new software system, called Interactive Groundwater (IGW), allows simulating complex flow and transport in aquifers subject to both systematic and "randomly" varying stresses and geological and chemical heterogeneity. Adopting a new paradigm, IGW eliminates a major bottleneck inherent in the traditional fragmented modeling technologies and enables real-time modeling, real-time visualization, real-time analysis, and real-time presentation. IGW functions as a "numerical laboratory" in which a researcher can freely explore in real-time: creating visually an aquifer of desired configurations, interactively imposing desired stresses, and then immediately investigating and visualizing the geology and the processes of flow and contaminant transport and transformation. A modeler can pause to edit at any time and interact on-line with any aspects (e.g., conceptual and numerical representation, boundary conditions, model solvers, and ways of visualization and analysis) of the integrated modeling process; he/she can initiate or stop, whenever needed, particle tracking, plume modeling, subscale modeling, cross-sectional modeling, stochastic modeling, monitoring

Evapotranspiration Estimates for a Stochastic Soil-Moisture Model

Science.gov (United States)

Chaleeraktrakoon, Chavalit; Somsakun, Somrit

2009-03-01

Potential evapotranspiration is information that is necessary for applying a widely used stochastic model of soil moisture (I. Rodriguez Iturbe, A. Porporato, L. Ridolfi, V. Isham and D. R. Cox, Probabilistic modelling of water balance at a point: The role of climate, soil and vegetation, Proc. Roy. Soc. London A455 (1999) 3789-3805). An objective of the present paper is thus to find a proper estimate of the evapotranspiration for the stochastic model. This estimate is obtained by comparing the calculated soil-moisture distribution resulting from various techniques, such as Thornthwaite, Makkink, Jensen-Haise, FAO Modified Penman, and Blaney-Criddle, with an observed one. The comparison results using five sequences of daily soil-moisture for a dry season from November 2003 to April 2004 (Udornthani Province, Thailand) have indicated that all methods can be used if the weather information required is available. This is because their soil-moisture distributions are alike. In addition, the model is shown to have its ability in approximately describing the phenomenon at a weekly or biweekly time scale which is desirable for agricultural engineering applications.

Qualitative and Quantitative Integrated Modeling for Stochastic Simulation and Optimization

Directory of Open Access Journals (Sweden)

Xuefeng Yan

2013-01-01

Full Text Available The simulation and optimization of an actual physics system are usually constructed based on the stochastic models, which have both qualitative and quantitative characteristics inherently. Most modeling specifications and frameworks find it difficult to describe the qualitative model directly. In order to deal with the expert knowledge, uncertain reasoning, and other qualitative information, a qualitative and quantitative combined modeling specification was proposed based on a hierarchical model structure framework. The new modeling approach is based on a hierarchical model structure which includes the meta-meta model, the meta-model and the high-level model. A description logic system is defined for formal definition and verification of the new modeling specification. A stochastic defense simulation was developed to illustrate how to model the system and optimize the result. The result shows that the proposed method can describe the complex system more comprehensively, and the survival probability of the target is higher by introducing qualitative models into quantitative simulation.

News Impact Curve for Stochastic Volatility Models

OpenAIRE

Makoto Takahashi; Yasuhiro Omori; Toshiaki Watanabe

2012-01-01

This paper proposes a new method to compute the news impact curve for stochastic volatility (SV) models. The new method incorporates the joint movement of return and volatility, which has been ignored by the extant literature, by simply adding a couple of steps to the Bayesian MCMC estimation procedures for SV models. This simple procedure is versatile and applicable to various SV type models. Contrary to the monotonic news impact functions in the extant literature, the new method gives a U-s...

From medium heterogeneity to flow and transport: A time-domain random walk approach

Science.gov (United States)

Hakoun, V.; Comolli, A.; Dentz, M.

2017-12-01

The prediction of flow and transport processes in heterogeneous porous media is based on the qualitative and quantitative understanding of the interplay between 1) spatial variability of hydraulic conductivity, 2) groundwater flow and 3) solute transport. Using a stochastic modeling approach, we study this interplay through direct numerical simulations of Darcy flow and advective transport in heterogeneous media. First, we study flow in correlated hydraulic permeability fields and shed light on the relationship between the statistics of log-hydraulic conductivity, a medium attribute, and the flow statistics. Second, we determine relationships between Eulerian and Lagrangian velocity statistics, this means, between flow and transport attributes. We show how Lagrangian statistics and thus transport behaviors such as late particle arrival times are influenced by the medium heterogeneity on one hand and the initial particle velocities on the other. We find that equidistantly sampled Lagrangian velocities can be described by a Markov process that evolves on the characteristic heterogeneity length scale. We employ a stochastic relaxation model for the equidistantly sampled particle velocities, which is parametrized by the velocity correlation length. This description results in a time-domain random walk model for the particle motion, whose spatial transitions are characterized by the velocity correlation length and temporal transitions by the particle velocities. This approach relates the statistical medium and flow properties to large scale transport, and allows for conditioning on the initial particle velocities and thus to the medium properties in the injection region. The approach is tested against direct numerical simulations.

Averaging principle for second-order approximation of heterogeneous models with homogeneous models.

Science.gov (United States)

Fibich, Gadi; Gavious, Arieh; Solan, Eilon

2012-11-27

Typically, models with a heterogeneous property are considerably harder to analyze than the corresponding homogeneous models, in which the heterogeneous property is replaced by its average value. In this study we show that any outcome of a heterogeneous model that satisfies the two properties of differentiability and symmetry is O(ε(2)) equivalent to the outcome of the corresponding homogeneous model, where ε is the level of heterogeneity. We then use this averaging principle to obtain new results in queuing theory, game theory (auctions), and social networks (marketing).

Averaging principle for second-order approximation of heterogeneous models with homogeneous models

Science.gov (United States)

Fibich, Gadi; Gavious, Arieh; Solan, Eilon

2012-01-01

Typically, models with a heterogeneous property are considerably harder to analyze than the corresponding homogeneous models, in which the heterogeneous property is replaced by its average value. In this study we show that any outcome of a heterogeneous model that satisfies the two properties of differentiability and symmetry is O(ɛ2) equivalent to the outcome of the corresponding homogeneous model, where ɛ is the level of heterogeneity. We then use this averaging principle to obtain new results in queuing theory, game theory (auctions), and social networks (marketing). PMID:23150569

Stochastic model of energetic nuclear reactor

International Nuclear Information System (INIS)

Bojko, R.V.; Ryazanov, V.V.

2002-01-01

Behaviour of nuclear reactor was treated using the theory of branching processes. As mathematical model descriptive the neutron number in time the Markov occasional process is proposed. Application of branching occasional processes with variable regime to the description of neutron behaviour in the reactor makes possible conducting strong description of critical operation regime and demonstrates the severity of the process. Three regimes of the critical behaviour depending on the sign of manipulated variables and feedbacks were discovered. Probability regularities peculiar to the behaviour of the reactor are embodied to the suggested stochastic model [ru

A coupled stochastic inverse-management framework for dealing with nonpoint agriculture pollution under groundwater parameter uncertainty

Science.gov (United States)