There are two sources of information available in empirical research in finance: one corresponding to historical data and the other to prices currently observed in the markets. When proposing a model, it is desirable to use information from both sources. However in modern finance, where stochastic differential equations have been one of the main modeling tools, the common models are typically different for historical data and for current market data. The former are usually assumed to be time homogeneous, while the latter are typically time in-homogeneous. This practice can be explained by the fact that a time-homogeneous model is stationary and easier to estimate, while time-inhomogeneous model are required in order to replicate market data sufficiently well without creating arbitrage opportunities. In this thesis, we study methods of statistical inference, both parametric and non-parametric, for stochastic differential equations with time-dependent parameters. In the first part, we propose a new class of stochastic differential equation with time-dependent drift and diffusion terms, where some of the parameters change according to a hidden Markov process. We show that under some technical conditions this innovative way of modeling switching times renders the resulting model stationary. We also explore different approaches to estimate parameters in our proposed model. Our simulation studies demonstrate that the parameters of the model can be efficiently estimated by using a version of the filtering method proposed in the literature. We illustrate our model and the proposed estimation method by applying them to interest rate data, and we detect significant time variations in early 1980s, when targets of the monetary policy in the United States were changed. One of the known drawbacks of parametric models is the risk of model misspecification. In the second part of the thesis, we allow the drift to be time-dependent and nonparametric, and our objective is to estimate it using a single trajectory of the process. The main idea underlying this method is to approximate the time-dependent function with a sequence of polynomials. Since we can estimate efficiently only a finite number of parameters for any finite length of data, in our method we propose to relate the number of parameters to the length of the observed trajectory. This idea is similar to the method of sieves proposed by Grenander (Abstract Inference, 1981). The asymptotic analysis that we present is based on the assumption that the length of available data $T$ increases to infinity. We investigate two cases, one is a Brownian motion with time-dependent drift and the other corresponds to a class of mean-reverting stochastic differential equations with time-dependent mean-reversion level. In both cases we prove asymptotic consistency and normality of a modified maximum likelihood estimator of the projected time-dependent component. The main challenge in proving our results in the second case stems from two features of the problem: one is due to the fact that coefficients of projections change with $T$ and the other is related to the confounding effect between the mean-reversion speed and the level function. By applying our method to the same interest rate data we use in the first part, we find another evidence of time-variation in the drift term.

There are two sources of information available in empirical research in finance: one corresponding to historical data and the other to prices currently observed in the markets. When proposing a model, it is desirable to use information from both sources. However in modern finance, where stochastic differential equations have been one of the main modeling tools, the common models are typically different for historical data and for current market data. The former are usually assumed to be time homogeneous, while the latter are typically time in-homogeneous. This practice can be explained by the fact that a time-homogeneous model is stationary and easier to estimate, while time-inhomogeneous model are required in order to replicate market data sufficiently well without creating arbitrage opportunities. In this thesis, we study methods of statistical inference, both parametric and non-parametric, for stochastic differential equations with time-dependent parameters. In the first part, we propose a new class of stochastic differential equation with time-dependent drift and diffusion terms, where some of the parameters change according to a hidden Markov process. We show that under some technical conditions this innovative way of modeling switching times renders the resulting model stationary. We also explore different approaches to estimate parameters in our proposed model. Our simulation studies demonstrate that the parameters of the model can be efficiently estimated by using a version of the filtering method proposed in the literature. We illustrate our model and the proposed estimation method by applying them to interest rate data, and we detect significant time variations in early 1980s, when targets of the monetary policy in the United States were changed. One of the known drawbacks of parametric models is the risk of model misspecification. In the second part of the thesis, we allow the drift to be time-dependent and nonparametric, and our objective is to estimate it using a single trajectory of the process. The main idea underlying this method is to approximate the time-dependent function with a sequence of polynomials. Since we can estimate efficiently only a finite number of parameters for any finite length of data, in our method we propose to relate the number of parameters to the length of the observed trajectory. This idea is similar to the method of sieves proposed by Grenander (Abstract Inference, 1981). The asymptotic analysis that we present is based on the assumption that the length of available data $T$ increases to infinity. We investigate two cases, one is a Brownian motion with time-dependent drift and the other corresponds to a class of mean-reverting stochastic differential equations with time-dependent mean-reversion level. In both cases we prove asymptotic consistency and normality of a modified maximum likelihood estimator of the projected time-dependent component. The main challenge in proving our results in the second case stems from two features of the problem: one is due to the fact that coefficients of projections change with $T$ and the other is related to the confounding effect between the mean-reversion speed and the level function. By applying our method to the same interest rate data we use in the first part, we find another evidence of time-variation in the drift term.

<p>Noise-induced stabilization occurs when an unstable deterministic system is stabilized by the addition of white noise. Proving that this phenomenon occurs for a particular system is often manifested through the construction of a global Lyapunov function. However, the procedure for constructing a Lyapunov function is often quite ad hoc, involving much time and tedium. In this thesis, a systematic algorithm for the construction of a global Lyapunov function for planar systems is presented. The general methodology is to construct a sequence of local Lyapunov functions in different regions of the plane, where the regions are delineated by different behaviors of the deterministic dynamics. A priming region, where the deterministic drift is directed inward, is first identified where there is an obvious choice for a local Lyapunov function. This priming Lyapunov function is then propagated to the other regions through a series of Poisson equations. The local Lyapunov functions are lastly patched together to form one smooth global Lyapunov function.</p><p>The algorithm is applied to a model problem which displays finite time blow up in the deterministic setting in order to prove that the system exhibits noise-induced stabilization. Moreover, the Lyapunov function constructed is in fact what we define to be a super Lyapunov function. We prove that the existence of a super Lyapunov function, along with a minorization condition, implies that the corresponding system converges to a unique invariant probability measure at an exponential rate that is independent of the initial condition.</p> / Dissertation

In this dissertation, we have studied diffusion models and their applications in risk theory and insurance. Let Xt be a d-dimensional diffusion process satisfying a system of Stochastic Differential Equations defined on an open set G Rd, and let Ut be a utility function of Xt with U0 = u0. Let T be the first time that Ut reaches a level u^*. We study the Laplace transform of the distribution of T, as well as the probability of ruin, psileft(u_{0}right)=Prleft{ T&lt;inftyright} , and other important probabilities. A class of exponential martingales is constructed to analyze the asymptotic properties of all probabilities. In addition, we prove that the expected discounted penalty function, a generalization of the probability of ultimate ruin, satisfies an elliptic partial differential equation, subject to some initial boundary conditions. Two examples from areas of actuarial work to which martingales have been applied are given to illustrate our methods and results: 1. Insurer's insolvency. 2. Terrorism risk. In particular, we study insurer's insolvency for the Cram'{e}r-Lundberg model with investments whose price follows a geometric Brownian motion. We prove the conjecture proposed by Constantinescu and Thommann. / Mathematics

This thesis investigates the possibility of approximating stationary solutions of stochastic differential equations using numerical methods. We consider a particular class of stochastic differential equations, which are known to generate random dynamical systems. The existence of stochastic stationary solution is proved using global attractor approach. Euler's numerical method, applied to the stochastic differential equation, is proved to generate a discrete random dynamical system. The existence of stationary solution is proved again using global attractor approach. At last we prove that the approximate stationary point converges in mean-square sense to the exact one as the time step of the numerical scheme diminishes.

博士（経営） / 85 p. / 一橋大学

Diese Diplomarbeit untersucht die Optimierung eines Mean-Variance Portfolios auf einem vollständigen Markt unter der Bedingung, dass die Insolvenz des Investors ausgeschlossen ist. Hierbei wird die duale Methode (auch Martingalmethode genannt)

There exist many numerical methods for numerical solutions of the systems of stochastic differential equations. We choose the method of deterministic quadrature formulae proposed by Müller–Gronbach, and Yaroslavtseva in 2016. The idea is to apply a simplified version of the cubature in Wiener space. We explain the method and check how good it works in the simplest case of the classical Black–Scholes model.

In this paper, we develop a model for predation movements of a subsocial spider species, Anelosimus studiosus. We expand on a previous model to include multiple spider interaction on the web as well as a latency period during predation. We then use the model to test different spatial configurations to determine the optimal spacing of spiders within a colony for successful capture during predation. The model simulations indicate that spiders uniformly spacing out along the edge of the web results in the most successful predation strategy. This is similar to the behavior observed by Ross (2013) in which it was determined to be statistically significant that during certain times of the day, spiders were positioned along the edge more than expected under complete spatial randomness.

We consider two distinct techniques for estimating random parameters in random differential equation (RDE) models. In one approach, the solution to a RDE is represented by a collection of solution trajectories in the form of sample deterministic equations. In a second approach we employ pointwise equivalent stochastic differential equation (SDE) representations for certain RDEs. Each of the approaches is tested using deterministic model comparison techniques for a logistic growth model which is viewed as a special case of a more general Bernoulli growth model. We demonstrate efficacy of the preferred method with experimental data using algae growth model comparisons.