Structural estimation of jump-diffusion processes in macroeconomics

This paper shows how to solve and estimate a continuous-time dynamic stochastic general equilibrium (DSGE) model with jumps. It also shows that a continuous-time formulation can make it simpler (relative to its discrete-time version) to compute and estimate the deep parameters using the likelihood function when non-linearities and/or non-normalities are considered. We illustrate our approach by solving and estimating the stochastic AK and the neoclassical growth models. Our Monte Carlo experiments demonstrate that non-normalities can be detected for this class of models. Moreover, we provide strong empirical evidence for jumps in aggregate US data.     

Stochastic differential equation modeling the prevalence of rare chronic diseases with an application to systemic lupus erythematosus

In this article we describe a system of stochastic differential equations to model the age‐specific prevalence of rare chronic diseases from incidence and mortality rates. As an application, the age profile of the prevalence of systemic lupus erythematosus in England and Wales in1995 is calculated. The results are in good agreement with the observed epidemiological measures.     

Optimal policy in Markov-switching rational expectations models

In this paper we consider the optimal quadratic control problem of Markov-switching linear rational expectation models. These models are general and flexible tools for modelling not only regime but also model or parameter uncertainty. We show, first, how to find the solution of a Markov-switching linear rational expectation model. Based on this solution we then show how to apply dynamic programming to find the optimal time-consistent policy and the resulting Nash-Stackelberg equilibrium. Suitable modifications of the algorithm allow to deal with the (non-RE) case in which the policymaker and the private sector hold different beliefs or probabilities over regime change. We also show how the optimisation procedure can be employed to obtain the optimal policy under commitment. As an illustration we compute the optimal policy in a small open economy subject to stochastic structural breaks in some of its key parameters.     

Comparing solution methods for dynamic equilibrium economies

This paper compares solution methods for dynamic equilibrium economies. We compute and simulate the stochastic neoclassical growth model with leisure choice using first, second, and fifth order perturbations in levels and in logs, the finite elements method, Chebyshev polynomials, and value function iteration for several calibrations. We document the performance of the methods in terms of computing time, implementation complexity, and accuracy, and we present some conclusions based on the reported evidence.     

Estimation of the Option Prime: Microsimulation of Backward Stochastic Differential Equations

A mathematical statistical model is needed to obtain an option prime and create a hedging strategy. With formulas derived from stochastic differential equations, the primes for US Dollar/Chilean Pesos currency options using a prime calculator are obtained. Furthermore, a backward simulation of the option prime trajectory is used with a numerical method created for backward stochastic differential equations. The use of statistics in finance is highly important in order to develop complex products.     

Nowcasting GDP with a pool of factor models and a fast estimation algorithm

We propose a novel mixed-frequency dynamic factor model with time-varying parameters and stochastic volatility for macroeconomic nowcasting and develop a fast estimation algorithm. This enables us to generate forecast densities based on a large space of factor models. We apply our framework to nowcast US GDP growth in real time. Our results reveal that stochastic volatility seems to improve the accuracy of point forecasts the most, compared to the constant-parameter factor model. These gains are most prominent during unstable periods such as the Covid-19 pandemic. Finally, we highlight indicators driving the US GDP growth forecasts and associated downside risks in real time.     

Dynamic Modeling with Structural Equations and Stochastic Differential Equations: Applications with the German Socio-economic Panel

The paper discusses an application of linear dynamic models to multi-wave longitudinal data. Starting from three-wave and four-wave simplex models using standard structural equations, linear dynamic state space models with stochastic differential equations are presented. The main differences between longitudinal structural equations (static view) and stochastic differential equations (dynamic view) are emphasized. Substantively, the models prove the relation, stability and change of two concepts in a period of 10 years: National Identity and Intention to stay in Germany. Data from a sample of migrant workers in Germany included in the German Socio-economic Panel (GSOEP) are used for the analyses. Results and further developments of dynamic models are discussed in the final section.The authors thank Hermann Singer for his comments and discussions on applications of dynamic models.     

The stochastic lake game: A numerical solution

In this paper, we numerically solve a stochastic dynamic programming problem for the solution of a stochastic dynamic game for which there is a potential function. The players select a mean level of control. The state transition dynamics is a function of the current state of the system and a multiplicative noise factor on the control variables of the players. The particular application is for lake water usage. The control variables are the levels of phosphorus discharged (typically by farmers) into the watershed of the lake, and the random shock is the rainfall that washes the phosphorus into the lake. The state of the system is the accumulated level of phosphorus in the lake. The system dynamics are sufficiently nonlinear so that there can be two Nash equilibria. A Skiba-like point can be present in the optimal control solution.We analyze (numerically) how the dynamics and the Skiba-like point change as the variance of the noise (the rain) increases. The numerical analysis uses a result of Dechert (1978. Optimal control problems from second order difference equations. Journal of Economic Theory 19, 50–63) to construct a potential function for the dynamic game. This greatly reduces the computational burden in finding Nash equilibria solutions for the dynamic game.     

Physics-informed Gaussian process regression for states estimation and forecasting in power grids

Real-time state estimation and forecasting are critical for the efficient operation of power grids. In this paper, a physics-informed Gaussian process regression (PhI-GPR) method is presented and used for forecasting and estimating the phase angle, angular speed, and wind mechanical power of a three-generator power grid system using sparse measurements. In standard data-driven Gaussian process regression (GPR), parameterized models for the prior statistics are fit by maximizing the marginal likelihood of observed data. In the PhI-GPR method, we propose to compute the prior statistics offline by solving stochastic differential equations (SDEs) governing the power grid dynamics. The short-term forecast of a power grid system dominated by wind generation is complicated by the stochastic nature of the wind and the resulting uncertainty in wind mechanical power. Here, we assume that the power grid dynamics are governed by swing equations, with the wind mechanical power fluctuating randomly in time. We solve these equations for the mean and covariances of the power grid states using the Monte Carlo simulation method.We demonstrate that the proposed PhI-GPR method can accurately forecast and estimate observed and unobserved states. For the considered problem, PhI-GPR has computational advantages over the ensemble Kalman filter (EnKF) method: In PhI-GPR, ensembles are computed offline and independently of the data acquisition process, whereas for EnFK, ensembles are computed online with data acquisition, rendering real-time forecast more challenging. We also demonstrate that the PhI-GPR forecast is more accurate than the EnKF forecast when the random mechanical wind power is non-Markovian. In contrast, the two methods produce similar forecasts for the Markovian mechanical wind power.For observed states, we show that PhI-GPR provides a forecast comparable to the standard data-driven GPR; both forecasts are significantly more accurate than the autoregressive integrated moving average (ARIMA) forecast. We also show that the ARIMA forecast is more sensitive to observation frequency and measurement errors than the PhI-GPR forecast.     

A class of statistical models to weaken independence in two-way contingency tables

In this paper we study a new class of statistical models for contingency tables. We define this class of models through a subset of the binomial equations of the classical independence model. We prove that they are log-linear and we use some notions from Algebraic Statistics to compute their sufficient statistic and their parametric representation. Moreover, we show how to compute maximum likelihood estimates and to perform exact inference through the Diaconis-Sturmfels algorithm. Examples show that these models can be useful in a wide range of applications.     

Stochastic models of resonating markets

This paper describes a way to model a seasonally and irregularly peaking price dynamics, as that originated in commodity and energy markets, using a system of coupled nonlinear stochastic differential equations. The specific case of an electric power market is used to show which microeconomic features this approach is able to model. Critical point analysis is used in a simple way to show how the interaction between dynamic criticality and stochasticity can be used to develop further models, useful to explore more deeply other types of peaking price dynamics.     

Markovian lifts of positive semidefinite affine Volterra-type processes

We consider stochastic partial differential equations appearing as Markovian lifts of matrix-valued (affine) Volterra-type processes from the point of view of the generalized Feller property (see, e.g., Dörsek and Teichmann in A semigroup point of view on splitting schemes for stochastic (partial) differential equations, 2010. arXiv:1011.2651). We introduce in particular Volterra Wishart processes with fractional kernels and values in the cone of positive semidefinite matrices. They are constructed from matrix products of infinite dimensional Ornstein–Uhlenbeck processes whose state space is the set of matrix-valued measures. Parallel to that we also consider positive definite Volterra pure jump processes, giving rise to multivariate Hawkes-type processes. We apply these affine covariance processes for multivariate (rough) volatility modeling and introduce a (rough) multivariate Volterra Heston-type model.

     

Analysis of optimal dynamic withdrawal policies in withdrawal guarantee products

Guaranteed Minimum Withdrawal Benefits (GMWB) are popular riders in variable annuities with withdrawal guarantees. With withdrawals spread over the life of the annuities contract, the benefit promises to return the entire initial annuitization amount irrespective of the market performance of the underlying fund portfolio. Treating the dynamic withdrawal rate as the control variable, the earlier works on GMWB have considered the construction of a continuous singular stochastic control model and the numerical solution of the resulting pricing model. This paper presents a more detailed characterization of the pricing properties of the GMWB and performs a full mathematical analysis of the optimal dynamic withdrawal policies under the competing factors of time value of fund, optionality value provided by the guarantee and penalty charge on excessive withdrawal. When a proportional penalty charge is applied on any withdrawal amount, we can reduce the pricing formulation to an optimal stopping problem with lower and upper obstacles. We then derive the integral equations for the determination of a pair of optimal withdrawal boundaries. When a proportional penalty charge is applied on the amount that is above the contractual withdrawal rate, we manage to characterize the behavior of the optimal withdrawal boundaries that separate the domain of the pricing models into three regions: no withdrawal, continuous withdrawal at the contractual rate and an immediate withdrawal of a finite amount. Under certain limiting scenarios such as a high policy fund value, the time close to expiry, or a low value of guarantee account, we manage to obtain analytical approximate solution to the singular stochastic control model of dynamic withdrawals.     

一族非李普希兹系数的随机微分方程

在非李普希兹条件下,对带跳随机微分方程数值方法的研究微乎其微。在非李普希兹系数下,对带跳随机微分方程数值方法进行了研究,并且得出近似解关于时间和开始点在Lp空间上一致收敛到解析解。     

Pricing VIX options with stochastic volatility and random jumps

This study presents an analytical exact solution for the price of VIX options under stochastic volatility model with simultaneous jumps in the asset price and volatility processes. We shall demonstrate that our new pricing formula can be used to efficiently compute the numerical values of a VIX option. While we also show that the numerical results obtained from our formula consistently match those obtained from Monte Carlo simulation perfectly as a verification of the correctness of our formula, numerical evidence is offered to illustrate that the correctness of the formula proposed in Lin and Chang (J Futur Markets 29(6), 523–543, 2009) is in serious doubt. Moreover, some important and distinct properties of VIX options (e.g., put-call parity, hedging ratios) are also examined and discussed.     

Using dynamic programming with adaptive grid scheme for optimal control problems in economics

The study of the solutions of dynamic models with optimizing agents has often been limited by a lack of available analytical techniques to explicitly find the global solution paths. On the other hand, the application of numerical techniques such as dynamic programming to find the solution in interesting regions of the state was restricted by the use of fixed grid size techniques. Following Grüne (Numer. Math. 75 (3) (1997) 319; University of Bayreuth, submitted, 2003), in this paper an adaptive grid scheme is used for finding the global solutions of discrete time Hamilton–Jacobi–Bellman equations. Local error estimates are established and an adapting iteration for the discretization of the state space is developed. The advantage of the use of adaptive grid scheme is demonstrated by computing the solutions of one- and two-dimensional economic models which exhibit steep curvature, complicated dynamics due to multiple equilibria, thresholds (Skiba sets) separating domains of attraction and periodic solutions. We consider deterministic and stochastic model variants. The studied examples are from economic growth, investment theory, environmental and resource economics.     

An analysis of the effect of noise in a heterogeneous agent financial market model

Heterogeneous agent models (HAMs) in finance and economics are often characterised by high dimensional nonlinear stochastic differential or difference systems. Because of the complexity of the interaction between the nonlinearities and noise, a commonly used, often called indirect, approach to the study of HAMs combines theoretical analysis of the underlying deterministic skeleton with numerical analysis of the stochastic model. However, it is well known that this indirect approach may not properly characterise the nature of the stochastic model. This paper aims to tackle this issue by developing a direct and analytical approach to the analysis of a stochastic model of speculative price dynamics involving two types of agents, fundamentalists and chartists, and the market price equilibria of which can be characterised by the stationary measures of a stochastic dynamical system. Using the stochastic method of averaging and stochastic bifurcation theory, we show that the stochastic model displays behaviour consistent with that of the underlying deterministic model when the time lag in the formation of price trends used by the chartists is far away from zero. However, when this lag approaches zero, such consistency breaks down.     

非高斯随机分布系统研究现状的分析

针对非高斯随机分布系统理论研究的发展问题,从非高斯随机分布系统的研究背景、系统静态模型的建立方法、系统动态模型的建立方法、非高斯随机分布系统的性能指标以及常用的控制算法等方面介绍了非高斯随机分布控制系统的研究现状,使读者对非高斯随机分布系统理论能进一步深入的了解。     

Rational expectations models: An approach using forward–backward stochastic differential equations

We show that a general class of continuous time rational expectations models can be reformulated as forward–backward stochastic differential equations (FBSDEs). Using this connection we obtain results on the conditions under which paths leading to, or keeping close to equilibrium exist, as well as their qualitative properties. We also provide a method for the construction of such paths through the connection of FBSDEs with quasilinear partial differential equations (PDEs). The theory is applied to specific macroeconomic models.