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 multidimensional 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.
Highlights
Statistical procedures for inference of nonstationary 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.Publications
Articles in Refereed Journals

S. Bürger, P. Mathé, Discretized Lavrent'ev regularization for the autoconvolution equation, Applicable Analysis. An International Journal, 96 (2017), pp. 16181637, DOI 10.1080/00036811.2016.1212336 .
Abstract
Lavrent?ev regularization for the autoconvolution equation was considered by Janno J. in Lavrent?ev regularization of illposed problems containing nonlinear neartomonotone operators with application to autoconvolution equation, Inverse Prob. 2000;16:333?348. Here this study is extended by considering discretization of the Lavrent?ev scheme by splines. It is shown how to maintain the known convergence rate by an appropriate choice of spline spaces and a proper choice of the discretization level. For piecewise constant splines the discretized equation allows for an explicit solver, in contrast to using higher order splines. This is used to design a fast implementation by means of postsmoothing, which provides results, which are indistinguishable from results obtained by direct discretization using cubic splines. 
M. Ladkau, J.G.M. Schoenmakers, J. Zhang, Libor model with expirywise stochastic volatility and displacement, International Journal of Portfolio Analysis and Management, 1 (2013), pp. 224249.
Abstract
We develop a multifactor stochastic volatility Libor model with displacement, where each individual forward Libor is driven by its own squareroot stochastic volatility process. The main advantage of this approach is that, maturitywise, each squareroot process can be calibrated to the corresponding cap(let)volastrike 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 volastrikematurity 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. 21142150.

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. 335344.
Abstract
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 distributioninvariant 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 distributioninvariant 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. 255269.

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

P. Čížek, W. Härdle, V. Spokoiny, Adaptive pointwise estimation in timeinhomogeneous conditional heteroscedasticity models, The Econometrics Journal, 12 (2009), pp. 248271.

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

E. Giacomini, W. Härdle, V. Spokoiny, Inhomogeneous dependency modelling with time varying copulae, Journal of Business & Economic Statistics, 27 (2009), pp. 224234.
Abstract
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 nonnormal behaviour of most financial time series calls for nonlinear dependency. The correct modelling of nongaussian 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 ValueatRisk (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 ValueatRisk, The Annals of Statistics, 37 (2009), pp. 14051436.
Abstract
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. 16191648.
Abstract
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 lackoffit testing procedure. Our main result is an oracle type inequality for the KullbackLeibler 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. 594607.
Abstract
Risk management technology applied to highdimensional portfolios needs simple and fast methods for calculation of value at risk (VaR). The multivariate normal framework provides a simple offtheshelf methodology but lacks the heavytailed distributional properties that are observed in data. A principle componentbased 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 highdimensional 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. 10521075.
Abstract
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 jumpdiffusion models and finite state Markov chains. By combining forward and reverse representations we then construct transition density estimators for chains which have rootN accuracy in any dimension and consider some applications. 
D. Spivakovskaya, A.W. Heemink, J.G.M. Schoenmakers, Twoparticle models for the estimation of the mean and standard deviation of concentrations in coastal waters, Stochastic Environmental Research and Risk Assessment (SERRA), 21 (2007), pp. 235251.

D. Belomestny, V. Spokoiny, Spatial aggregation of local likelihood estimates with applications to classification, The Annals of Statistics, 35 (2007), pp. 22872311.
Abstract
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. 126139.
Abstract
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. 927938.
Abstract
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 reversetime 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. 283300.

H. Haaf, O. Reiss, J.G.M. Schoenmakers, Numerically stable computation of CreditRisk+, Phys. Rev. E (3), 6 (2004), pp. 110.
Abstract
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 wellknown 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 forwardreverse representations, Bernoulli. Official Journal of the Bernoulli Society for Mathematical Statistics and Probability, 10 (2004), pp. 281312.
Abstract
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 timeinhomogeneous volatility models, The Annals of Statistics, 32 (2004), pp. 577602.

W. Härdle, H. Herwatz, V. Spokoiny, Time inhomogeneous multiple volatility modelling, Journal of Financial Econometrics, 1 (2003), pp. 5595.

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. 787795.
Abstract
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 biochemicaloxygen 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. 345362.

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. 6776.

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 1317 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. 18131824.
Preprints, Reports, Technical Reports

A. Lejay, P. Pigato, A threshold model for local volatility: Evidence of leverage and mean reversion effects on historical data, Preprint no. 2467, WIAS, Berlin, 2017, DOI 10.20347/WIAS.PREPRINT.2467 .
Abstract, PDF (534 kByte)
In financial markets, low prices are generally associated with high volatilities and viceversa, this well known stylized fact usually being referred to as leverage effect. We propose a local volatility model, given by a stochastic differential equation with piecewise constant coefficients, which accounts of leverage and meanreversion effects in the dynamics of the prices. This model exhibits a regime switch in the dynamics accordingly to a certain threshold. It can be seen as a continuous time version of the SelfExciting Threshold Autoregressive (SETAR) model. We propose an estimation procedure for the volatility and drift coefficients as well as for the threshold level. Tests are performed on the daily prices of 21 assets. They show empirical evidence for leverage and meanreversion effects, consistent with the results in the literature.
Talks, Poster

M. Ladkau, A new multifactor stochastic volatility model with displacement, PreMoLab Workshop on: Advances in predictive modeling and optimization, May 16  17, 2013, WIASBerlin, 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, HumboldtUniversität zu Berlin, Wirtschaftswissenschaftliche Fakultät, July 1, 2011.

J. Zhang, Lpsolutions 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 distributioninvariant 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 DuisburgEssen, 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 StochastikTage 2004, March 23  26, 2004, Universität Karlsruhe, March 25, 2004.

J.G.M. Schoenmakers, Transition density estimation for stochastic differential equations via forwardreverse representations, TandemWorkshop StochastikNumerik, June 11  March 26, 2004, HumboldtUniversitä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, PhysikalischTechnische Bundesanstalt, March 13, 2003.

J.G.M. Schoenmakers, Kreditrisiko PortfolioModelle, 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, arXiv.org, 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, arXiv.org, 2010.
Abstract
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 HamiltonJacobiBellman 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 CramerLundberg 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 socalled 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. 0306, 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.