Partial Lanczos SVD methods for R
|
|
- Meghan Patterson
- 6 years ago
- Views:
Transcription
1 Partial Lanczos SVD methods for R Bryan Lewis 1, adapted from the work of Jim Baglama 2 and Lothar Reichel 3 1 REvolution Computing 2 University of Rhode Island 3 Kent State University UseR 2009
2 Outline SVD and partial SVD Partial Lanczos bidiagonalization The irlba package
3 SVD Let A R l n, l n. A = n j=1 σ j u j v T j, { vj T v k = uj T 1 if j = k, u k = 0 o.w., u j R l, v j R n, j = 1, 2,..., n, and σ 1 σ 2 σ n 0.
4 Partial SVD Let k < n. Ã k := k j=1 σ j u j v T j
5 Partial Lanczos bi-diagonalization Start with a given vector p 1. Compute m steps of the Lanczos process: AP m = Q m B m A T Q m = P m Bm T + r m em, T B m R m m, P m R n m, Q m R l m, P T mp m = Q T mq m = I m, r m R n, P T mr m = 0, P m = [p 1, p 2,..., p m ].
6 Approximating partial SVD with partial Lanczos bi-diagonalization A T AP m = A T Q m B m = P m BmB T m + r m emb T m,
7 Approximating partial SVD with partial Lanczos bi-diagonalization A T AP m = A T Q m B m = P m BmB T m + r m emb T m, AA T Q m = AP m B T m + Ar m e T m, = Q m B m B T m + Ar m e T m.
8 Approximating partial SVD with partial Lanczos bi-diagonalization Let B m = m σj B uj B j=1 ( v B j ) T.
9 Approximating partial SVD with partial Lanczos bi-diagonalization Let B m = m σj B uj B j=1 ( v B j ) T. ( ) i.e., B m vj B = σj B uj B, and Bmu T j b = σj B vj B. Define:
10 Approximating partial SVD with partial Lanczos bi-diagonalization Let B m = m σj B uj B j=1 ( v B j ) T. ( ) i.e., B m vj B = σj B uj B, and Bmu T j b = σj B vj B. Define: σ j := σ B j, ũ j := Q m u B j, ṽ j := P m v B j.
11 Partial SVD approximation of A Aṽ j = AP m vj B = Q m B m vj B = σ B j Q m u B j = σ j ũ j,
12 Partial SVD approximation of A Aṽ j = AP m vj B = Q m B m vj B = σ B j Q m u B j = σ j ũ j, A T ũ j = A T Q m u B j = P m B T mu B j + r m e T mu B j = σ B j P m v B j + r m e T mu B j = σ j ṽ j + r m e T mu B j.
13 Augment and restart 1. Compute the Lanczos process up to step m. 2. Compute k < m approximate singular vectors. 3. Orthogonalize against the approximate singular vectors to get a new starting vector. 4. Continue the Lanczos process with the new starting vector for m more steps. 5. Check for convergence tolerance and exit if met. 6. GOTO 1.
14 artial Lanczos SVD methods for R Sketch of the augmented process... P k+1 := [ṽ 1, ṽ 2,..., ṽ k, p m+1 ],
15 artial Lanczos SVD methods for R Sketch of the augmented process... P k+1 := [ṽ 1, ṽ 2,..., ṽ k, p m+1 ], A P k+1 = [ σ 1 ũ 1, σ 1 ũ 2,..., σ k ũ k, Ap m+1 ]
16 artial Lanczos SVD methods for R Sketch of the augmented process... P k+1 := [ṽ 1, ṽ 2,..., ṽ k, p m+1 ], A P k+1 = [ σ 1 ũ 1, σ 1 ũ 2,..., σ k ũ k, Ap m+1 ] Orthogonalize Ap m+1 against {ũ j } k j=1 : Ap m+1 = k j=1 ρ jũj + r k.
17 artial Lanczos SVD methods for R Sketch of the augmented process... P k+1 := [ṽ 1, ṽ 2,..., ṽ k, p m+1 ], A P k+1 = [ σ 1 ũ 1, σ 1 ũ 2,..., σ k ũ k, Ap m+1 ] Orthogonalize Ap m+1 against {ũ j } k j=1 : Ap m+1 = k j=1 ρ jũj + r k. Q k+1 := [ũ 1, ũ 2,..., ũ k, r k / r k ], σ 1 ρ 1 σ 2 ρ 2 B k+1 :=... ρk. r k
18 artial Lanczos SVD methods for R Sketch of the augmented process... P k+1 := [ṽ 1, ṽ 2,..., ṽ k, p m+1 ], A P k+1 = [ σ 1 ũ 1, σ 1 ũ 2,..., σ k ũ k, Ap m+1 ] Orthogonalize Ap m+1 against {ũ j } k j=1 : Ap m+1 = k j=1 ρ jũj + r k. Q k+1 := [ũ 1, ũ 2,..., ũ k, r k / r k ], σ 1 ρ 1 σ 2 ρ 2 B k+1 :=... ρk. r k A P k+1 = Q k+1 Bk+1.
19 The irlba package Usage: irlba (A, nu = 5, nv = 5, adjust = 3, aug = "ritz", sigma = "ls", maxit = 1000, reorth = 1, tol = 1e-06, V = NULL)
20 Small example > A<-matrix (rnorm(5000*5000),5000,5000) > require (irlba) > system.time (L<-irlba (A,nu=5,nv=5)) user system elapsed > gc() used (Mb)... max used (Mb) Ncells Vcells
21 artial Lanczos SVD methods for R Small example (continued) > system.time (S<-svd(A,nu=5,nv=5)) user system elapsed > gc() used (Mb)... max used (Mb) Ncells Vcells
22 Small example (continued) > system.time (S<-svd(A,nu=5,nv=5)) user system elapsed > gc() used (Mb)... max used (Mb) Ncells Vcells > sqrt (crossprod(s$d[1:5]-l$d)/crossprod(s$d[1:5])) [,1] [1,] e-12
23 Large examples (live demo) The R implementation of IRLBA supports: Dense real/complex in-process matrices (normal R matrices) Sparse real in-process matrices (Matrix) Dense, real in- or out-of-process huge matrices with bigmemory + bigalgebra
24 References jbaglama 3. reichel 4. blewis 5.
Solving the Stochastic Steady-State Diffusion Problem Using Multigrid
Solving the Stochastic Steady-State Diffusion Problem Using Multigrid Tengfei Su Applied Mathematics and Scientific Computing Program Advisor: Howard Elman Department of Computer Science May 5, 2016 Tengfei
More informationModeling multiple runoff tables
Modeling multiple runoff tables Vincent Lous www.posthuma-partners.nl Motivation Why creating a "model" for all lines of business is important Motivation Why creating a "model" for all lines of business
More informationTowards efficient option pricing in incomplete markets
Towards efficient option pricing in incomplete markets GPU TECHNOLOGY CONFERENCE 2016 Shih-Hau Tan 1 2 1 Marie Curie Research Project STRIKE 2 University of Greenwich Apr. 6, 2016 (University of Greenwich)
More informationA way to improve incremental 2-norm condition estimation
A way to improve incremental 2-norm condition estimation Jurjen Duintjer Tebbens Institute of Computer Science Academy of Sciences of the Czech Republic duintjertebbens@cs.cas.cz Miroslav Tůma Institute
More informationELEMENTS OF MATRIX MATHEMATICS
QRMC07 9/7/0 4:45 PM Page 5 CHAPTER SEVEN ELEMENTS OF MATRIX MATHEMATICS 7. AN INTRODUCTION TO MATRICES Investors frequently encounter situations involving numerous potential outcomes, many discrete periods
More information1 Explicit Euler Scheme (or Euler Forward Scheme )
Numerical methods for PDE in Finance - M2MO - Paris Diderot American options January 2018 Files: https://ljll.math.upmc.fr/bokanowski/enseignement/2017/m2mo/m2mo.html We look for a numerical approximation
More informationPackage ELMSO. September 3, 2018
Type Package Package ELMSO September 3, 2018 Title Implementation of the Efficient Large-Scale Online Display Advertising Algorithm Version 1.0.0 Date 2018-8-31 Maintainer Courtney Paulson
More informationBROWNIAN MOTION II. D.Majumdar
BROWNIAN MOTION II D.Majumdar DEFINITION Let (Ω, F, P) be a probability space. For each ω Ω, suppose there is a continuous function W(t) of t 0 that satisfies W(0) = 0 and that depends on ω. Then W(t),
More informationGeometry of orthogonally invariant matrix varieties
Geometry of orthogonally invariant matrix varieties Dmitriy Drusvyatskiy Mathematics, University of Washington Joint work with H.-L. Lee (UW), G. Ottaviano (Florence), and R.R. Thomas (UW) UC Davis Algebra
More informationROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices
ROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices Bachelier Finance Society Meeting Toronto 2010 Henley Business School at Reading Contact Author : d.ledermann@icmacentre.ac.uk Alexander
More informationA model reduction approach to numerical inversion for parabolic partial differential equations
A model reduction approach to numerical inversion for parabolic partial differential equations Liliana Borcea Alexander V. Mamonov 2, Vladimir Druskin 3, Mikhail Zaslavsky 3 University of Michigan, Ann
More informationPackage samplesize4surveys
Type Package Package samplesize4surveys July 23, 2018 Title Sample Size Calculations for Complex Surveys Version 3.6.1.0 Date 2018-07-20 Author Hugo Andres Gutierrez Rojas Maintainer Hugo Andres Gutierrez
More informationLecture Note of Bus 41202, Spring 2010: Analysis of Multiple Series with Applications. x 1t x 2t. holdings (OIH) and energy select section SPDR (XLE).
Lecture Note of Bus 41202, Spring 2010: Analysis of Multiple Series with Applications Focus on two series (i.e., bivariate case) Time series: Data: x 1, x 2,, x T. X t = Some examples: (a) U.S. quarterly
More informationEfficient Portfolio and Introduction to Capital Market Line Benninga Chapter 9
Efficient Portfolio and Introduction to Capital Market Line Benninga Chapter 9 Optimal Investment with Risky Assets There are N risky assets, named 1, 2,, N, but no risk-free asset. With fixed total dollar
More informationBOUNDS FOR THE LEAST SQUARES RESIDUAL USING SCALED TOTAL LEAST SQUARES
BOUNDS FOR THE LEAST SQUARES RESIDUAL USING SCALED TOTAL LEAST SQUARES Christopher C. Paige School of Computer Science, McGill University Montreal, Quebec, Canada, H3A 2A7 paige@cs.mcgill.ca Zdeněk Strakoš
More informationPackage GenOrd. September 12, 2015
Package GenOrd September 12, 2015 Type Package Title Simulation of Discrete Random Variables with Given Correlation Matrix and Marginal Distributions Version 1.4.0 Date 2015-09-11 Author Alessandro Barbiero,
More informationExercise List: Proving convergence of the (Stochastic) Gradient Descent Method for the Least Squares Problem.
Exercise List: Proving convergence of the (Stochastic) Gradient Descent Method for the Least Squares Problem. Robert M. Gower. October 3, 07 Introduction This is an exercise in proving the convergence
More information9.1 Principal Component Analysis for Portfolios
Chapter 9 Alpha Trading By the name of the strategies, an alpha trading strategy is to select and trade portfolios so the alpha is maximized. Two important mathematical objects are factor analysis and
More informationStatistical and Computational Inverse Problems with Applications Part 5B: Electrical impedance tomography
Statistical and Computational Inverse Problems with Applications Part 5B: Electrical impedance tomography Aku Seppänen Inverse Problems Group Department of Applied Physics University of Eastern Finland
More informationReduced models for sparse grid discretizations of the multi-asset Black-Scholes equation
Reduced models for sparse grid discretizations of the multi-asset Black-Scholes equation The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters.
More informationA model reduction approach to numerical inversion for parabolic partial differential equations
A model reduction approach to numerical inversion for parabolic partial differential equations Liliana Borcea Alexander V. Mamonov 2, Vladimir Druskin 2, Mikhail Zaslavsky 2 University of Michigan, Ann
More informationPremia 14 HESTON MODEL CALIBRATION USING VARIANCE SWAPS PRICES
Premia 14 HESTON MODEL CALIBRATION USING VARIANCE SWAPS PRICES VADIM ZHERDER Premia Team INRIA E-mail: vzherder@mailru 1 Heston model Let the asset price process S t follows the Heston stochastic volatility
More informationA Correlated Sampling Method for Multivariate Normal and Log-normal Distributions
A Correlated Sampling Method for Multivariate Normal and Log-normal Distributions Gašper Žerovni, Andrej Trov, Ivan A. Kodeli Jožef Stefan Institute Jamova cesta 39, SI-000 Ljubljana, Slovenia gasper.zerovni@ijs.si,
More informationDesign of a Financial Application Driven Multivariate Gaussian Random Number Generator for an FPGA
Design of a Financial Application Driven Multivariate Gaussian Random Number Generator for an FPGA Chalermpol Saiprasert, Christos-Savvas Bouganis and George A. Constantinides Department of Electrical
More informationSupply Chain Outsourcing Under Exchange Rate Risk and Competition
Supply Chain Outsourcing Under Exchange Rate Risk and Competition Published in Omega 2011;39; 539-549 Zugang Liu and Anna Nagurney Department of Business and Economics The Pennsylvania State University
More informationShow that the column rank and the row rank of A are both equal to 3.
hapter Vectors and matrices.. Exercises. Let A 2 5 4 3 2 4 2 2 3 5 4 2 4 3 Show that the column rank and the row rank of A are both equal to 3. 2. Let x and y be column vectors of size n, andleti be the
More informationSteepest descent and conjugate gradient methods with variable preconditioning
Ilya Lashuk and Andrew Knyazev 1 Steepest descent and conjugate gradient methods with variable preconditioning Ilya Lashuk (the speaker) and Andrew Knyazev Department of Mathematics and Center for Computational
More informationModelling, Estimation and Hedging of Longevity Risk
IA BE Summer School 2016, K. Antonio, UvA 1 / 50 Modelling, Estimation and Hedging of Longevity Risk Katrien Antonio KU Leuven and University of Amsterdam IA BE Summer School 2016, Leuven Module II: Fitting
More informationAdjoint methods for option pricing, Greeks and calibration using PDEs and SDEs
Adjoint methods for option pricing, Greeks and calibration using PDEs and SDEs Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Oxford-Man Institute of Quantitative Finance
More informationIDENTIFYING BROAD AND NARROW FINANCIAL RISK FACTORS VIA CONVEX OPTIMIZATION: PART I
1 IDENTIFYING BROAD AND NARROW FINANCIAL RISK FACTORS VIA CONVEX OPTIMIZATION: PART I Lisa Goldberg lrg@berkeley.edu MMDS Workshop. June 22, 2016. joint with Alex Shkolnik and Jeff Bohn. Identifying Broad
More information(a) Describe the game in plain english and find its equivalent strategic form.
Risk and Decision Making (Part II - Game Theory) Mock Exam MIT/Portugal pages Professor João Soares 2007/08 1 Consider the game defined by the Kuhn tree of Figure 1 (a) Describe the game in plain english
More informationGraduate School of Business, University of Chicago Business 41202, Spring Quarter 2007, Mr. Ruey S. Tsay. Solutions to Final Exam
Graduate School of Business, University of Chicago Business 41202, Spring Quarter 2007, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (30 pts) Answer briefly the following questions. 1. Suppose that
More informationUser s Guide for the Matlab Library Implementing Closed Form MLE for Diffusions
User s Guide for the Matlab Library Implementing Closed Form MLE for Diffusions Yacine Aït-Sahalia Department of Economics and Bendheim Center for Finance Princeton University and NBER This Version: July
More informationA general approach to calculating VaR without volatilities and correlations
page 19 A general approach to calculating VaR without volatilities and correlations Peter Benson * Peter Zangari Morgan Guaranty rust Company Risk Management Research (1-212) 648-8641 zangari_peter@jpmorgan.com
More informationGraduate Macro Theory II: Notes on Value Function Iteration
Graduate Macro Theory II: Notes on Value Function Iteration Eric Sims University of Notre Dame Spring 07 Introduction These notes discuss how to solve dynamic economic models using value function iteration.
More informationPerturbation Bounds for Determinants and Characteristic Polynomials
Perturbation Bounds for Determinants and Characteristic Polynomials Ilse Ipsen North Carolina State University, Raleigh, USA Joint work with: Rizwana Rehman Characteristic Polynomials n n complex matrix
More informationAppendix to Dividend yields, dividend growth, and return predictability in the cross-section of. stocks
Appendix to Dividend yields, dividend growth, and return predictability in the cross-section of stocks Paulo Maio 1 Pedro Santa-Clara 2 This version: February 2015 1 Hanken School of Economics. E-mail:
More informationChapter 5 Finite Difference Methods. Math6911 W07, HM Zhu
Chapter 5 Finite Difference Methods Math69 W07, HM Zhu References. Chapters 5 and 9, Brandimarte. Section 7.8, Hull 3. Chapter 7, Numerical analysis, Burden and Faires Outline Finite difference (FD) approximation
More informationEc2723, Asset Pricing I Class Notes, Fall Complete Markets, Incomplete Markets, and the Stochastic Discount Factor
Ec2723, Asset Pricing I Class Notes, Fall 2005 Complete Markets, Incomplete Markets, and the Stochastic Discount Factor John Y. Campbell 1 First draft: July 30, 2003 This version: October 10, 2005 1 Department
More informationUnpublished Appendices to Market Reactions to Tangible and Intangible Information. Market Reactions to Different Types of Information
Unpublished Appendices to Market Reactions to Tangible and Intangible Information. This document contains the unpublished appendices for Daniel and Titman (006), Market Reactions to Tangible and Intangible
More informationDynamic Portfolio Execution Detailed Proofs
Dynamic Portfolio Execution Detailed Proofs Gerry Tsoukalas, Jiang Wang, Kay Giesecke March 16, 2014 1 Proofs Lemma 1 (Temporary Price Impact) A buy order of size x being executed against i s ask-side
More informationarxiv: v1 [math.st] 6 Jun 2014
Strong noise estimation in cubic splines A. Dermoune a, A. El Kaabouchi b arxiv:1406.1629v1 [math.st] 6 Jun 2014 a Laboratoire Paul Painlevé, USTL-UMR-CNRS 8524. UFR de Mathématiques, Bât. M2, 59655 Villeneuve
More informationFinancial Risk Forecasting Chapter 3 Multivariate volatility models
Financial Risk Forecasting Chapter 3 Multivariate volatility models Jon Danielsson 2017 London School of Economics To accompany Financial Risk Forecasting www.financialriskforecasting.com Published by
More informationRobust Loss Development Using MCMC: A Vignette
Robust Loss Development Using MCMC: A Vignette Christopher W. Laws Frank A. Schmid July 2, 2010 Abstract For many lines of insurance, the ultimate loss associated with a particular exposure (accident or
More informationAsymptotic methods in risk management. Advances in Financial Mathematics
Asymptotic methods in risk management Peter Tankov Based on joint work with A. Gulisashvili Advances in Financial Mathematics Paris, January 7 10, 2014 Peter Tankov (Université Paris Diderot) Asymptotic
More informationWeb-based Supplementary Materials for. A space-time conditional intensity model. for invasive meningococcal disease occurence
Web-based Supplementary Materials for A space-time conditional intensity model for invasive meningococcal disease occurence by Sebastian Meyer 1,2, Johannes Elias 3, and Michael Höhle 4,2 1 Department
More informationPackage ald. February 1, 2018
Type Package Title The Asymmetric Laplace Distribution Version 1.2 Date 2018-01-31 Package ald February 1, 2018 Author Christian E. Galarza and Victor H. Lachos
More informationConvex-Cardinality Problems Part II
l 1 -norm Methods for Convex-Cardinality Problems Part II total variation iterated weighted l 1 heuristic matrix rank constraints Prof. S. Boyd, EE364b, Stanford University Total variation reconstruction
More informationAitken Acceleration of the EM Algorithm
Aitken Acceleration for a Single Variable Suppose x = ξ (scalar for the moment is the solution of a problem, which is attained by successive approximations of the form [1] which converge tothe value ξ
More informationIDENTIFYING BROAD AND NARROW FINANCIAL RISK FACTORS VIA CONVEX OPTIMIZATION: PART II
1 IDENTIFYING BROAD AND NARROW FINANCIAL RISK FACTORS VIA CONVEX OPTIMIZATION: PART II Alexander D. Shkolnik ads2@berkeley.edu MMDS Workshop. June 22, 2016. joint with Jeffrey Bohn and Lisa Goldberg. Identifying
More informationOperational Risk Aggregation
Operational Risk Aggregation Professor Carol Alexander Chair of Risk Management and Director of Research, ISMA Centre, University of Reading, UK. Loss model approaches are currently a focus of operational
More informationQuasi-Monte Carlo Methods in Financial Engineering: An Equivalence Principle and Dimension Reduction
Quasi-Monte Carlo Methods in Financial Engineering: An Equivalence Principle and Dimension Reduction Xiaoqun Wang,2, and Ian H. Sloan 2,3 Department of Mathematical Sciences, Tsinghua University, Beijing
More informationECON 581. Introduction to Arrow-Debreu Pricing and Complete Markets. Instructor: Dmytro Hryshko
ECON 58. Introduction to Arrow-Debreu Pricing and Complete Markets Instructor: Dmytro Hryshko / 28 Arrow-Debreu economy General equilibrium, exchange economy Static (all trades done at period 0) but multi-period
More informationIdentification and Estimation of Dynamic Games when Players Belief Are Not in Equilibrium
Identification and Estimation of Dynamic Games when Players Belief Are Not in Equilibrium A Short Review of Aguirregabiria and Magesan (2010) January 25, 2012 1 / 18 Dynamics of the game Two players, {i,
More informationCS 234 Winter 2019 Assignment 1 Due: January 23 at 11:59 pm
CS 234 Winter 2019 Assignment 1 Due: January 23 at 11:59 pm For submission instructions please refer to website 1 Optimal Policy for Simple MDP [20 pts] Consider the simple n-state MDP shown in Figure
More informationLindner, Szimayer: A Limit Theorem for Copulas
Lindner, Szimayer: A Limit Theorem for Copulas Sonderforschungsbereich 386, Paper 433 (2005) Online unter: http://epub.ub.uni-muenchen.de/ Projektpartner A Limit Theorem for Copulas Alexander Lindner Alexander
More informationWrite legibly. Unreadable answers are worthless.
MMF 2021 Final Exam 1 December 2016. This is a closed-book exam: no books, no notes, no calculators, no phones, no tablets, no computers (of any kind) allowed. Do NOT turn this page over until you are
More informationProbability and Stochastics for finance-ii Prof. Joydeep Dutta Department of Humanities and Social Sciences Indian Institute of Technology, Kanpur
Probability and Stochastics for finance-ii Prof. Joydeep Dutta Department of Humanities and Social Sciences Indian Institute of Technology, Kanpur Lecture - 07 Mean-Variance Portfolio Optimization (Part-II)
More informationCurves on a Surface. (Com S 477/577 Notes) Yan-Bin Jia. Oct 24, 2017
Curves on a Surface (Com S 477/577 Notes) Yan-Bin Jia Oct 24, 2017 1 Normal and Geodesic Curvatures One way to examine how much a surface bends is to look at the curvature of curves on the surface. Let
More informationLecture 8: Linear Prediction: Lattice filters
1 Lecture 8: Linear Prediction: Lattice filters Overview New AR parametrization: Reflection coefficients; Fast computation of prediction errors; Direct and Inverse Lattice filters; Burg lattice parameter
More informationThe Use of Penultimate Approximations in Risk Management
The Use of Penultimate Approximations in Risk Management www.math.ethz.ch/ degen (joint work with P. Embrechts) 6th International Conference on Extreme Value Analysis Fort Collins CO, June 26, 2009 Penultimate
More informationCLASSIC TWO-STEP DURBIN-TYPE AND LEVINSON-TYPE ALGORITHMS FOR SKEW-SYMMETRIC TOEPLITZ MATRICES
CLASSIC TWO-STEP DURBIN-TYPE AND LEVINSON-TYPE ALGORITHMS FOR SKEW-SYMMETRIC TOEPLITZ MATRICES IYAD T ABU-JEIB Abstract We present ecient classic two-step Durbin-type and Levinsontype algorithms for even
More informationImprovement and Efficient Implementation of a Lattice-based Signature scheme
Improvement and Efficient Implementation of a Lattice-based Signature scheme, Johannes Buchmann Technische Universität Darmstadt TU Darmstadt August 2013 Lattice-based Signatures1 Outline Introduction
More informationFinancial Time Series Lecture 10: Analysis of Multiple Financial Time Series with Applications
Financial Time Series Lecture 10: Analysis of Multiple Financial Time Series with Applications Reference: Chapters 8 and 10 of the textbook. We shall focus on two series (i.e., the bivariate case) Time
More informationInfluence of Real Interest Rate Volatilities on Long-term Asset Allocation
200 2 Ó Ó 4 4 Dec., 200 OR Transactions Vol.4 No.4 Influence of Real Interest Rate Volatilities on Long-term Asset Allocation Xie Yao Liang Zhi An 2 Abstract For one-period investors, fixed income securities
More informationMelbourne Institute Working Paper Series Working Paper No. 22/07
Melbourne Institute Working Paper Series Working Paper No. 22/07 Permanent Structural Change in the US Short-Term and Long-Term Interest Rates Chew Lian Chua and Chin Nam Low Permanent Structural Change
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Other Miscellaneous Topics and Applications of Monte-Carlo Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More information1 Explicit Euler Scheme (or Euler Forward Scheme )
Numerical methods for PDE in Finance - M2MO - Paris Diderot American options January 2017 Files: https://ljll.math.upmc.fr/bokanowski/enseignement/2016/m2mo/m2mo.html We look for a numerical approximation
More informationInformation Asymmetry, Information Precision, and the Cost of Capital
University of Pennsylvania ScholarlyCommons Accounting Papers Wharton Faculty Research 2012 Information Asymmetry, Information Precision, and the Cost of Capital Richard A. Lambert University of Pennsylvania
More informationChapter 10: Mixed strategies Nash equilibria, reaction curves and the equality of payoffs theorem
Chapter 10: Mixed strategies Nash equilibria reaction curves and the equality of payoffs theorem Nash equilibrium: The concept of Nash equilibrium can be extended in a natural manner to the mixed strategies
More informationLinearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing
Linearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing Liuren Wu, Baruch College Joint work with Peter Carr and Xavier Gabaix at New York University Board of
More informationMSc in Financial Engineering
Department of Economics, Mathematics and Statistics MSc in Financial Engineering On Numerical Methods for the Pricing of Commodity Spread Options Damien Deville September 11, 2009 Supervisor: Dr. Steve
More informationPackage FADA. May 20, 2016
Type Package Package FADA May 20, 2016 Title Variable Selection for Supervised Classification in High Dimension Version 1.3.2 Date 2016-05-12 Author Emeline Perthame (INRIA, Grenoble, France), Chloe Friguet
More informationChapter 14. The Multi-Underlying Black-Scholes Model and Correlation
Chapter 4 The Multi-Underlying Black-Scholes Model and Correlation So far we have discussed single asset options, the payoff function depended only on one underlying. Now we want to allow multiple underlyings.
More informationAsian Option Pricing: Monte Carlo Control Variate. A discrete arithmetic Asian call option has the payoff. S T i N N + 1
Asian Option Pricing: Monte Carlo Control Variate A discrete arithmetic Asian call option has the payoff ( 1 N N + 1 i=0 S T i N K ) + A discrete geometric Asian call option has the payoff [ N i=0 S T
More informationSome Bounds for the Singular Values of Matrices
Applied Mathematical Sciences, Vol., 007, no. 49, 443-449 Some Bounds for the Singular Values of Matrices Ramazan Turkmen and Haci Civciv Department of Mathematics, Faculty of Art and Science Selcuk University,
More informationMS-E2114 Investment Science Exercise 4/2016, Solutions
Capital budgeting problems can be solved based on, for example, the benet-cost ratio (that is, present value of benets per present value of the costs) or the net present value (the present value of benets
More informationCSE325 Principles of Operating Systems. Deadlock. David P. Duggan. March 22, 2011
CSE325 Principles of Operating Systems Deadlock David P. Duggan dduggan@sandia.gov March 22, 2011 Outline What s deadlock? Necessary conditions for deadlock Deadlock prevention Deadlock avoidance Deadlock
More informationNo-Arbitrage, State Prices and Trade in Thin Financial Markets
No-Arbitrage, State Prices and Trade in Thin Financial Markets Andrés Carvajal Marek Weretka October 7, 2009 Abstract We examine how non-competitiveness in financial markets affects the choice of asset
More informationNOTICE TO MEMBERS No May 6, 2003 SYMBOL CONVERSION LONG-TERM EQUITY OPTIONS EXPIRING IN JANUARY 2004
NOTICE TO MEMBERS No. 2003-025 May 6, 2003 SYMBOL CONVERSION LONG-TERM EQUITY OPTIONS EXPIRING IN JANUARY Bourse de Montréal Inc. (the Bourse) and (CDCC) hereby inform you that at the opening on Monday,
More informationValidation of Credit Rating Models - A Preliminary Look at Methodology and Literature Review
JCIC JCIC Column 1 93 1-15 Validation of Credit Rating Models - A Preliminary Look at Methodology and Literature Review Ming-Yi Sun, Szu-Fang Wang JCIC Risk Research Team I. Introduction In preparing for
More information(-J. MQ Gateway Trust. Financial Report ARSN MACQUARIE. For the year ended 30 June 2011
(-J MACQUARIE MQ Gateway Trust Financial Report ARSN 123 784 930 The Responsible Entity of MQ Gateway Trust is MQ Portfolio Management Limited (ACN 092 552 611) The Responsible Entity's registered office
More informationInference of Several Log-normal Distributions
Inference of Several Log-normal Distributions Guoyi Zhang 1 and Bose Falk 2 Abstract This research considers several log-normal distributions when variances are heteroscedastic and group sizes are unequal.
More informationNumerical Methods For American Option Pricing. Peng Liu. June 2008
Numerical Methods For American Option Pricing Peng Liu June 2008 Abstract An analytic solution does not exist for evaluating the American put option. Usually, the value is obtained by applying numerical
More informationChapter 7: Portfolio Theory
Chapter 7: Portfolio Theory 1. Introduction 2. Portfolio Basics 3. The Feasible Set 4. Portfolio Selection Rules 5. The Efficient Frontier 6. Indifference Curves 7. The Two-Asset Portfolio 8. Unrestriceted
More informationMixedModR2 Erika Mudrak Thursday, August 30, 2018
MixedModR Erika Mudrak Thursday, August 3, 18 Generate the Data Generate data points from a population with one random effect: levels of Factor A, each sampled 5 times set.seed(39) siga
More information# generate data num.obs <- 100 y <- rnorm(num.obs,mean = theta.true, sd = sqrt(sigma.sq.true))
Posterior Sampling from Normal Now we seek to create draws from the joint posterior distribution and the marginal posterior distributions and Note the marginal posterior distributions would be used to
More information2D penalized spline (continuous-by-continuous interaction)
2D penalized spline (continuous-by-continuous interaction) Two examples (RWC, Section 13.1): Number of scallops caught off Long Island Counts are made at specific coordinates. Incidence of AIDS in Italian
More informationTopics in Computational Sustainability CS 325 Spring 2016
Topics in Computational Sustainability CS 325 Spring 2016 Note to other teachers and users of these slides. Andrew would be delighted if you found this source material useful in giving your own lectures.
More information2013 Proposed Non-Ad Valorem Recap
BL BELMONT LAKES CDD $ 4,200.00 42 42 $ 176,400.00 0 0 42 42 1 BROWARD COUNTY GARBAGE - 1 $ 270.00 182 190 $ 51,300.00 0 0 182 190 4 BROWARD COUNTY GARBAGE - 4 $ 270.00 653 858 $ 231,660.00 2 0 655 858
More informationUnderstanding Predictability (JPE, 2004)
Understanding Predictability (JPE, 2004) Lior Menzly, Tano Santos, and Pietro Veronesi Presented by Peter Gross NYU October 19, 2009 Presented by Peter Gross (NYU) Understanding Predictability October
More informationAdjusting Accounts When a Customer is also a Vendor
Adjusting Accounts When a Customer is also a Vendor INTRODUCTION The following example illustrates how to adjust accounts if you have a balance with a vendor, and the vendor wants to apply your balance
More informationKostas Kyriakoulis ECG 790: Topics in Advanced Econometrics Fall Matlab Handout # 5. Two step and iterative GMM Estimation
Kostas Kyriakoulis ECG 790: Topics in Advanced Econometrics Fall 2004 Matlab Handout # 5 Two step and iterative GMM Estimation The purpose of this handout is to describe the computation of the two-step
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationRobust portfolio optimization under multiperiod mean-standard deviation criterion
Robust portfolio optimization under multiperiod mean-standard deviation criterion Spiridon Penev 1 Pavel Shevchenko 2 Wei Wu 1 1 The University of New South Wales, Sydney Australia 2 Macquarie University
More informationImage analysis of malign melanoma: Waveles and svd
Image analysis of malign melanoma: Waveles and svd Dan Dolonius University of Gothenburg gusdolod@student.gu.se April 28, 2015 Dan Dolonius (Applied Mathematics) Image analysis of malign melanoma April
More informationLattice (Binomial Trees) Version 1.2
Lattice (Binomial Trees) Version 1. 1 Introduction This plug-in implements different binomial trees approximations for pricing contingent claims and allows Fairmat to use some of the most popular binomial
More informationLecture 10: Alignments with Affine Gaps. The Local Alignment Recurrence
Lecture 10: Alignments with Affine Gaps Study Chapter 6.9-6.10 1 The Local Alignment Recurrence The largest value of s i,j over the whole edit graph is the score of the best local alignment. Smith-Waterman
More informationIterative Refinement in Three Precisions
Iterative Refinement in Three Precisions for Research Fast and Accurate MattersSolution of Ill-Conditioned February Sparse 25, 2009 Linear Systems Nick Nick Higham Higham Director Schoolof ofresearch Mathematics
More informationFormulating SALCs with Projection Operators
Formulating SALCs with Projection Operators U The mathematical form of a SALC for a particular symmetry species cannot always be deduced by inspection (e.g., e 1g and e u pi-mos of benzene). U A projection
More information