This topic is currently not treated at the WIAS

The valuation and risk assessment of financial products requires appropriate analysis of the characteristics involved in the modeling. Particular examples are volatility and characteristics of extreme events. Usually the models under consideration are multi-dimensional and have time dependent parameters.

The developed approaches for estimating time dependent volatility end extreme events have to be extended to the multivariate case. In this context suitable dimension reduction methods will play an important role. The results will be needed for modeling complex financial structures.


Statistical procedures for inference of non-stationary time series are developed in the univariate case. In particular, new procedures for estimating time dependent volatilities based on local adaptive smoothing techniques are proposed.


  Articles in Refereed Journals

  • M. Ladkau, J.G.M. Schoenmakers, J. Zhang, Libor model with expiry-wise stochastic volatility and displacement, International Journal of Portfolio Analysis and Management, 1 (2013) pp. 224--249.
    We develop a multi-factor stochastic volatility Libor model with displacement, where each individual forward Libor is driven by its own square-root stochastic volatility process. The main advantage of this approach is that, maturity-wise, each square-root process can be calibrated to the corresponding cap(let)vola-strike panel at the market. However, since even after freezing the Libors in the drift of this model, the Libor dynamics are not affine, new affine approximations have to be developed in order to obtain Fourier based (approximate) pricing procedures for caps and swaptions. As a result, we end up with a Libor modeling package that allows for efficient calibration to a complete system of cap/swaption market quotes that performs well even in crises times, where structural breaks in vola-strike-maturity panels are typically observed.

  • G. Dos Reis, A. Réveillac, J. Zhang, FBSDEs with time delayed generators: $L^p$-solutions, differentiability, representation formulas and path regularity, Stochastic Processes and their Applications, 121 (2011) pp. 2114--2150.

  • V. Krätschmer, H. Zähle, Sensitivity of risk measures with respect to the normal approximation of total claim distributions, Insurance: Mathematics & Economics, 49 (2011) pp. 335--344.
    A simple and commonly used method to approximate the total claim distribution of a (possible weakly dependent) insurance collective is the normal approximation. In this article, we investigate the error made when the normal approximation is plugged in a fairly general distribution-invariant risk measure. We focus on the rate of the convergence of the error relative to the number of clients, we specify the relative error's asymptotic distribution, and we illustrate our results by means of a numerical example. Regarding the risk measure, we take into account distortion risk measures as well as distribution-invariant coherent risk measures.

  • Y. Chen, W. Härdle, V. Spokoiny, GHICA --- Risk analysis with GH distributions and independent components, Journal of Empirical Finance, 17 (2010) pp. 255--269.

  • A.E. Kyprianou, R.L. Loeffen, Refracted Lévy processes, Annales de l'Institut Henri Poincare. Probabilites et Statistiques, 46 (2010) pp. 24--44.

  • P. Čížek, W. Härdle, V. Spokoiny, Adaptive pointwise estimation in time-inhomogeneous conditional heteroscedasticity models, The Econometrics Journal, 12 (2009) pp. 248--271.

  • P. Friz, S. Benaim, Regular variation and smile asymptotics, Mathematical Finance. An International Journal of Mathematics, Statistics and Financial Economics, 19 (2009) pp. 1--12.

  • E. Giacomini, W. Härdle, V. Spokoiny, Inhomogeneous dependency modelling with time varying copulae, Journal of Business & Economic Statistics, 27 (2009) pp. 224--234.
    Measuring dependence in a multivariate time series is tantamount to modelling its dynamic structure in space and time. In the context of a multivariate normally distributed time series, the evolution of the covariance (or correlation) matrix over time describes this dynamic. A wide variety of applications, though, requires a modelling framework different from the multivariate normal. In risk management the non-normal behaviour of most financial time series calls for nonlinear dependency. The correct modelling of non-gaussian dependencies is therefore a key issue in the analysis of multivariate time series. In this paper we use copulae functions with adaptively estimated time varying parameters for modelling the distribution of returns, free from the usual normality assumptions. Further, we apply copulae to estimation of Value-at-Risk (VaR) of a portfolio and show its better performance over the RiskMetrics approach, a widely used methodology for VaR estimation

  • V. Spokoiny, Multiscale local change point detection with applications to Value-at-Risk, The Annals of Statistics, 37 (2009) pp. 1405--1436.
    This paper offers a new procedure for nonparametric estimation and forecasting of time series with applications to volatility modeling for financial data. The approach is based on the assumption of local homogeneity: for every time point there exists a historical emphinterval of homogeneity, in which the volatility parameter can be well approximated by a constant. The procedure recovers this interval from the data using the local change point (LCP) analysis. Afterwards the estimate of the volatility can be simply obtained by local averaging. The approach carefully addresses the question of choosing the tuning parameters of the procedure using the so called ``propagation'' condition. The main result claims a new ``oracle'' inequality in terms of the modeling bias which measures the quality of the local constant approximation. This result yields the optimal rate of estimation for smooth and piecewise constant volatility functions. Then the new procedure is applied to some data sets and a comparison with a standard GARCH model is also provided. Finally we discuss applications of the new method to the Value at Risk problem. The numerical results demonstrate a very reasonable performance of the new method.

  • I.G. Grama, V. Spokoiny, Statistics of extremes by oracle estimation, The Annals of Statistics, 36 (2008) pp. 1619--1648.
    We use the fitted Pareto law to construct an accompanying approximation of the excess distribution function. A selection rule of the location of the excess distribution function is proposed based on a stagewise lack-of-fit testing procedure. Our main result is an oracle type inequality for the Kullback-Leibler loss of the obtained adaptive estimator.

  • Y. Chen, W. Härdle, V. Spokoiny, Portfolio value at risk based on independent components analysis, Journal of Computational and Applied Mathematics, 205 (2007) pp. 594--607.
    Risk management technology applied to high-dimensional portfolios needs simple and fast methods for calculation of value at risk (VaR). The multivariate normal framework provides a simple off-the-shelf methodology but lacks the heavy-tailed distributional properties that are observed in data. A principle component-based method (tied closely to the elliptical structure of the distribution) is therefore expected to be unsatisfactory. Here, we propose and analyze a technology that is based on independent component analysis (ICA). We study the proposed ICVaR methodology in an extensive simulation study and apply it to a high-dimensional portfolio situation. Our analysis yields very accurate VaRs.

  • G.N. Milstein, J.G.M. Schoenmakers, V. Spokoiny, Forward and reverse representations for Markov chains, Stochastic Processes and their Applications, 117 (2007) pp. 1052--1075.
    In this paper we carry over the concept of reverse probabilistic representations developed in Milstein, Schoenmakers, Spokoiny (2004) for diffusion processes, to discrete time Markov chains. We outline the construction of reverse chains in several situations and apply this to processes which are connected with jump-diffusion models and finite state Markov chains. By combining forward and reverse representations we then construct transition density estimators for chains which have root-N accuracy in any dimension and consider some applications.

  • D. Spivakovskaya, A.W. Heemink, J.G.M. Schoenmakers, Two-particle models for the estimation of the mean and standard deviation of concentrations in coastal waters, Sochastic Environmental Research and Risk Assessment (SERRA), 21 (2007) pp. 235--251.

  • D. Belomestny, V. Spokoiny, Spatial aggregation of local likelihood estimates with applications to classification, The Annals of Statistics, 35 (2007) pp. 2287--2311.
    This paper presents a new method for spatially adaptive local (constant) likelihood estimation which applies to a broad class of nonparametric models, including the Gaussian, Poisson and binary response models. The main idea of the method is given a sequence of local likelihood estimates (''weak'' estimates), to construct a new aggregated estimate whose pointwise risk is of order of the smallest risk among all ``weak'' estimates. We also propose a new approach towards selecting the parameters of the procedure by providing the prescribed behavior of the resulting estimate in the simple parametric situation. We establish a number of important theoretical results concerning the optimality of the aggregated estimate. In particular, our ``oracle'' results claims that its risk is up to some logarithmic multiplier equal to the smallest risk for the given family of estimates. The performance of the procedure is illustrated by application to the classification problem. A numerical study demonstrates its nice performance in simulated and real life examples.

  • E. VAN DEN Berg, A.W. Heemink, H.X. Lin, J.G.M. Schoenmakers, Probability density estimation in stochastic environmental models using reverse representations, Stochastic Environmental Research & Risk Assessment, 20 (2006) pp. 126--139.
    The estimation of probability densities of variables described by stochastic differential equations has long been done using forward time estimators, which rely on the generation of forward in time realizations of the model. Recently, an estimator based on the combination of forward and reverse time estimators has been developed. This estimator has a higher order of convergence than the classical one. In this article, we explore the new estimator and compare the forward and forward? reverse estimators by applying them to a biochemical oxygen demand model. Finally, we show that the computational efficiency of the forward?reverse estimator is superior to the classical one, and discuss the algorithmic aspects of the estimator.

  • D. Spivakovskaya, A.W. Heemink, G.N. Milstein, J.G.M. Schoenmakers, Simulation of the transport of particles in coastal waters using forward and reverse time diffusion, Advances in Water Resources, 28 (2005) pp. 927--938.
    Particle models are often used to simulate the spreading of a pollutant in coastal waters in case of a calamity at sea. Here many different particle tracks starting at the point of release are generated to determine the particle concentration at some critical locations after the release. This Monte Carlo method however consumes a large CPU time. Recently, Milstein, Schoenmakers and Spokoiny (2003) introduced the concept of reverse-time diffusion. They derived a reverse system from the original forward simulation model and showed that the Monte Carlo estimator can also be based on realizations of this reverse system. In this paper we apply this concept to estimate particle concentrations in coastal waters. The results for the experiments considered show that the CPU time compared with the classical method is reduced orders of magnitude.

  • M. Giurcanu, V. Spokoiny, Confidence estimation of the covariance function of stationary and locally stationary processes, Statistics & Decisions. International Journal for Statistical Theory and Related Fields, 22 (2004) pp. 283--300.

  • H. Haaf, O. Reiss, J.G.M. Schoenmakers, Numerically stable computation of CreditRisk+, Phys. Rev. E (3), 6 (2004) pp. 1--10.
    The CreditRisk+ model launched by Credit Suisse First Boston in 1997 is widely used by practitioners in the banking sector as a simple means for the quantification of credit risk, primarily of the loan book. We present an alternative numerical recursion scheme for CreditRisk+, equivalent to an algorithm recently proposed by Giese, that is based on well-known expansions of the logarithm and the exponential of a power series. We show that it is advantageous for the Panjer recursion advocated in the original CreditRisk+ document, in that it is numerically stable. The crucial stability arguments are explained in detail. We explain how to apply the suggested recursion scheme to incorporate stochastic exposures into the CreditRisk+ model as introduced by Tasche (2004). Finally, the computational complexity of the resulting algorithm is stated and compared with other methods for computing the CreditRisk+ loss distribution.

  • G.N. Milstein, J.G.M. Schoenmakers, V. Spokoiny, Transition density estimation for stochastic differential equations via forward-reverse representations, Bernoulli. Official Journal of the Bernoulli Society for Mathematical Statistics and Probability, 10 (2004) pp. 281--312.
    The general reverse diffusion equations are derived and applied to the problem of transition density estimation of diffusion processes between two fixed states. For this problem we propose density estimation based on forward?reverse representations and show that this method allows essentially better results to be achieved than the usual kernel or projection estimation based on forward representations only.

  • V. Spokoiny, D. Mercurio, Statistical inference for time-inhomogeneous volatility models, The Annals of Statistics, 32 (2004) pp. 577--602.

  • M.-Y. Cheng, J. Fan, V. Spokoiny, Dynamic nonparametric filtering with application to volatility estimation, , (2003) pp. 315-333.

  • W. Härdle, H. Herwatz, V. Spokoiny, Time inhomogeneous multiple volatility modelling, , 1 (2003) pp. 55-95.

  • J.G.M. Schoenmakers, A.W. Heemink, K. Ponnambalm, P.E. Kloeden, Variance reduction for Monte Carlo simulation of stochastic environmental models, Applied Mathematical Modelling. Simulation and Computation for Engineering and Environmental Systems. Elsevier Science Inc., New York, NY. English, English abstracts., 26 (2002) pp. 787--795.
    To determine the probability of exceedence Monte Carlo simulation of stochastic models is often used. Mathematically this requires the evaluation of an expectation of some function of a solution of a stochastic model. This can be reformulated as a Kolmogorov final value problem. It can thus be calculated numerically by either solving a deterministic partial differential equation (Kolmogorov's Backwards equations) or by simulating a large number of trajectories of the stochastic differential equation. Here we discuss a composite method of variance reduced Monte Carlo simulation. The variance reduction is obtained by the Girsanov transformation to modify the stochastic model by a correction term that is obtained from an approximate solution of the partial differential equation computed by a classical numerical method. The composite method is more efficient than either the standard Monte Carlo or the classical numerical method. The approach is applied to estimate the probability of exceedence in a model for biochemical-oxygen demand.

  Contributions to Collected Editions

  • M. Elagin, V. Spokoiny, Locally time homogeneous time series modelling, in: Applied Quantitative Finance (2nd edition), W. Härdle, N. Hautsch, L. Overbeck, eds., Springer, Berlin, 2008, pp. 345--362.

  • H. Haaf, O. Reiss, J.G.M. Schoenmakers, Numerically stable computation of Credit Risk+, in: CreditRisk+ in the Banking Industry, M. Gundlach, F. Lehrbass, eds., XII, Springer, Berlin Heidelberg, 2004, pp. 67-76.

  • D. Spivakovskaya, A.W. Heemink, J.G.M. Schoenmakers, G.N. Milstein, Stochastic modeling of transport in coastal waters using forward and reverse time diffusion, in: Computational Methods in Water Resources, Proceedings of the 15th International Conference on Computational Methods in Water Resources (CMWR XV), June 13--17 2004, Chapel Hill, North Carolina, C.T. Miller, M.W. Farthing, W.G. Gray, G.F. Pinder, eds., 55 of Developments in Water Science, Elsevier Science, Amsterdam, 2004, pp. 1813--1824.

  Talks, Poster

  • M. Ladkau, A new multi-factor stochastic volatility model with displacement, PreMoLab Workshop on: Advances in predictive modeling and optimization, May 16 - 17, 2013, WIAS-Berlin, May 16, 2013.

  • V. Panov, Affine stochastic volatility models: asymptotic behavior of the characteristic function and estimation of the Blumenthal - Getoor index, Tagung des SFB 649 "Ökonomisches Risiko" in Motzen, June 30 - July 2, 2011, Humboldt-Universität zu Berlin, Wirtschaftswissenschaftliche Fakultät, July 1, 2011.

  • J. Zhang, Lp-solutions of BSDEs with time delayed generators, Mathematikkolloquium, Universität Innsbruck, Fachbereich Mathematik, Austria, April 28, 2011.

  • J. Zhang, Solvability and numerical simulation of BSDEs related to BSPDEs with applications to utility maximization, Universität Innsbruck, Fachbereich Mathematik, Austria, April 26, 2011.

  • V. Krätschmer, A uniform central limit theorem for distortions of empirical distributions with applications to nonparametric estimation of distribution-invariant risk measures, Leipziger Stochastik Tage, March 1 - 5, 2010, Universität Leipzig, Fakultät für Mathematik und Informatik, March 2, 2010.

  • V. Krätschmer, Error caused by normal approximation of the total claim distribution when plugged in risk measures, 4th International Conference Mathematical and Statistical Methods for Actuarial Sciences and Finance, April 7 - 9, 2010, Università degli Studi di Salerno, Dipartimento di Scienze Economiche e Statistiche, Ravello, Italy, April 8, 2010.

  • V. Krätschmer, Error caused by normal approximation of the total claim distribution when plugged in risk measures, DAGStat2010, March 23 - 26, 2010, Technische Universität Dortmund, Fakultät Statistik, March 23, 2010.

  • V. Krätschmer, Nichtparametrische Schätzung verteilungsinvarianter Risikomaße, Universität Bayreuth, Fachbereich Mathematik, July 27, 2010.

  • V. Krätschmer, Nichtparametrische Schätzung verteilungsinvarianter Risikomaße, Universität Mannheim, Institut für Mathematik und Informatik, May 4, 2010.

  • V. Krätschmer, Nichtparametrische Schätzung verteilungsinvarianter Risikomaße, Heinrich Heine Universität Düsseldorf, Institut für Mathematik, April 23, 2010.

  • V. Krätschmer, Nichtparametrische Schätzung von verteilungsinvarianten Risikomaßen, Mathematisches Kolloquium, Universität Duisburg-Essen, Fakultät für Mathematik, October 19, 2010.

  • V. Krätschmer, Nonparametric estimation of low invariant risk measures for time series, Statistische Woche Nürnberg 2010, September 14 - 17, 2010, Deutschen Statistischen Gesellschaft (DStatG), September 16, 2010.

  • P. Friz, Stochastic analysis and quantitative finance, Microsoft Research Cambridge, UK, January 5, 2010.

  • V. Spokoiny, Local parametric estimation, October 18 - 22, 2010, École Nationale de la Statistique et de l'Analyse de l'Information (ENSAI), Rennes, France.

  • J.G.M. Schoenmakers, Statistical and numerical methods for evaluation for financial derivates and risk, Center Days 2009 (DFG Research Center scshape Matheon), March 30 - April 1, 2009, Technische Universität Berlin, March 31, 2009.

  • V. Spokoiny, Adaptive local parametric estimation, Université Joseph Fourier Grenoble I, Équipe de Statistique et Modélisation Stochastique, Laboratoire Jean Kuntzmann, France, February 26, 2009.

  • V. Spokoiny, Modern nonparametric statistics (block lecture), October 2 - 13, 2009, École Nationale de la Statistique et de l'Analyse de l'Information (ENSAI), Rennes, France.

  • V. Spokoiny, Modern nonparametric statistics (block lecture), October 18 - 29, 2009, Yale University, New Haven, USA.

  • V. Spokoiny, Modern nonparametric statistics (block lecture), January 13 - 16, 2009, École Nationale de la Statistique et de l'Analyse de l'Information (ENSAI), Rennes, France.

  • J. Polzehl, On a nonstationary structural adaptive approach to volatility estimation, University of Gothenburg, Centre for Finance, Sweden, May 5, 2004.

  • J.G.M. Schoenmakers, Numerically stable computation of CreditRisk+, Karlsruher Stochastik-Tage 2004, March 23 - 26, 2004, Universität Karlsruhe, March 25, 2004.

  • J.G.M. Schoenmakers, Transition density estimation for stochastic differential equations via forward-reverse representations, Tandem-Workshop Stochastik-Numerik, June 11 - March 26, 2004, Humboldt-Universität zu Berlin, June 11, 2004.

  • J. Polzehl, Local modeling by structural adaptation, The Art of Semiparametrics, October 19 - 21, 2003, Berlin, October 20, 2003.

  • J. Polzehl, Structural adaptive smoothing methods and applications in imaging, Magnetic Resonance Seminar, Physikalisch-Technische Bundesanstalt, March 13, 2003.

  • J.G.M. Schoenmakers, Kreditrisiko Portfolio-Modelle, Kreditanstalt für Wiederaufbau, Frankfurt, July 18, 2003.

  • J.G.M. Schoenmakers, Transition density estimation for stochastic differential equations via forward reverse representations, IV IMACS Seminar on Monte Carlo Methods (MCM 2003), September 15 - 19, 2003, Berlin, September 16, 2003.

  • J. Polzehl, Structural adaptive smoothing and its applications in imaging and time series, Uppsala University, Department of Mathematics, Sweden, May 2, 2002.

  External Preprints

  • R.L. Loeffen, I. Czarna, Z. Palmowski, Parisian ruin probability for spectrally negative Lévy processes, Preprint no. arXiv:1102.4055, Cornell University Library,, 2011.

  • A. Kyprianou, R.L. Loeffen, J.-L. Perez, Optimal control with absolutely continuous strategies for spectrally negative Levy processes, Preprint no. arXiv:1008.2363, Cornell University Library,, 2010.
    In the last few years there has been renewed interest in the classical control problem of de Finetti for the case that underlying source of randomness is a spectrally negative Levy process. In particular a significant step forward is made in an article of Loeffen where it is shown that a natural and very general condition on the underlying Levy process which allows one to proceed with the analysis of the associated Hamilton-Jacobi-Bellman equation is that its Levy measure is absolutely continuous, having completely monotone density. In this paper we consider de Finetti's control problem but now with the restriction that control strategies are absolutely continuous with respect to Lebesgue measure. This problem has been considered by Asmussen and Taksar, Jeanblanc and Shiryaev and Boguslavskaya in the diffusive case and Gerber and Shiu for the case of a Cramer-Lundberg process with exponentially distributed jumps. We show the robustness of the condition that the underlying Levy measure has a completely monotone density and establish an explicit optimal strategy for this case that envelopes the aforementioned existing results. The explicit optimal strategy in question is the so-called refraction strategy.

  • E. VAN  den Berg, A.W. Heemink, H.X. Lin, J.G.M. Schoenmakers, Probability density estimation in stochastic environmental models using reverse representations, Report no. 03-06, Delft University of Technology, 2003.

  • E. VAN DEN Berg, A.W. Heemink, H.X. Lin, J.G.M. Schoenmakers, Probability density estimation in stochastic environmental models using reverse representations, Report no. 6, TU Delft, The Netherlands, Applied Mathematical Analysis, 2003.