Effects of Survey Geometry on Power Spectrum Covariance for Cosmic Shear Survey

Size: px
Start display at page:

Download "Effects of Survey Geometry on Power Spectrum Covariance for Cosmic Shear Survey"

Transcription

1 Effects of Survey Geometry on Power Spectrum Covariance for Cosmic Shear Survey Ryuichi Takahashi (Hirosaki U) with Shunji Soma, Masahiro Takada and Issha Kayo RT+, in preparation

2 Abstract What survey shape is the BEST (or WORST) to extract the cosmological information in weak-lensing cosmic shear survey?

3 Quiz Which observational survey geometry is better to extract the cosmological information in weak-lensing cosmic shear survey? Square VS Belt Elongated Rectangular the surface areas are same

4 Quiz Which observational survey geometry is better to extract the cosmological information in weak-lensing cosmic shear survey? Square VS Belt Elongated Rectangular better the surface areas are same

5 Quiz2 Which observational survey geometry is better to extract the cosmological information in weak-lensing cosmic shear survey? Square Sparse Sampling VS the surface areas are same

6 Quiz2 Which observational survey geometry is better to extract the cosmological information in weak-lensing cosmic shear survey? Square Sparse Sampling VS better the surface areas are same

7 Summary Sparse Sampling Belt Elongated Rectangular The above observational shapes have advantages at both large and small scales: larger-scale fluctuations can be scales (e.g. Kaiser 1986, 1998; Kilbinger & Schneider 2004; Chiang+ 2013) the non-gaussian error can be scales

8 Summary Sparse Sampling Belt Elongated Rectangular The above observational shapes have advantages at both large and small scales: larger-scale fluctuations can be scales (e.g. Kaiser 1986, 1998; Kilbinger & Schneider 2004; Chiang+ 2013) the non-gaussian error can be scales

9 Summary Sparse Sampling Belt Elongated Rectangular Signal-to-Noise in power spectrum measurement can be 3 times higher than the S/N in the square observational shape ( 9 times larger survey area)

10 Introduction Weak-lensing cosmic shear is a powerful observational tool to constrain the cosmological models (including dark energy). previous surveys : SDSS, CHFTLS, COSMOS, on-going & future surveys : Subaru HSC, DES, LSST, from LSST homepage

11 Convergence field κ(θ) : projected surface mass density at angular position θ Convergence power spectrum C(l) (=shear power spectrum) C l = 1 N l κ (l) 2 l l κ(l) N l : multipole : convergence in Fourier domain : # of mode

12 Measurement large scales Gaussian error C(l) C(l) = 2 N l N l : # of mode observational area Gaussian error 1 (observational area) 1/2 Measurement errors of cosmological parameters 1 (observational area) 1/2

13 Measurement small scales Small-scale fluctuations are correlated due to non-linear mode coupling (Meiksin & White 1999; Scoccimarro+ 1999) Measurement error = Gaussian error + non-gaussian error depends on non-linearity of fluctuations (scale, redshift) & survey shape What observational shape is the best to reduce the non-gaussian error?

14 Log-normal Convergence field Convergence field is assumed to follow log-normal distribution which is supported by gravitational-lensing ray-tracing simulations (e.g. Jain+ 2000; Taruya+ 2002; Das & Ostriker 2006; RT+ 2011) κ 0 2 σ κ : minimum convergence empty beam : variance the mean is zero dp dκ κ 0 0 κ

15 203deg Log-normal Convergence Maps input power spectrum at redshift 0.9 halo-fit nonlinear model each map is a square shape 203x203 deg 2 (= 4π stradian) angular resolution = 1arcmin 1000 maps prepared θ 2 cut a region & calculate the power spectrum 0 203deg θ 1

16 Log-normal Convergence Maps We calculate the power spectra & covariance from the 1000 maps SN(Signal-to-Noise ratio) in power spectrum measurement Cumulative SN up to l max S N 2 = C ij 1 l i,l j <l max C W l i C W (l i ) C W (l) : power spectrum C ij : covariance

17 Results rectangular observational shape with an area of 100deg^2 θ A θ B θ A θ B = 100 deg 2

18 Cumulative SN for rectangular observational shapes Survey area = 100deg^2 Cumulative SN up to l_max z s = 0.9 Maximum multipole

19 Cumulative SN for rectangular observational shapes Survey area = 100deg^2 Cumulative SN up to l_max 2 times better 3 times better z s = 0.9 Maximum multipole

20 What observational survey shape is the best (or worst) to suppress the non-gaussian error in the power spectrum measurement Problem setting : 1. prepare 100 small patches, each patch has an area of 1x1 deg^2 2. place these 100 patches on the all sky (203x203 deg^2) survey area is 100deg^2 in total

21 100 deg^2 area in total

22 Cumulative SN up to l_max 3 times better 3 times better Maximum multipole

23 The non-gaussian error depends on variance of convergence in the survey region: (Hilbert+ 2011; Takada & Hu 2013) σ W 2 = d 2 q 2π 2 W q 2 C(q) = 1 S W 2 d2 θ 1 d 2 θ 2 W θ 1 W(θ 2 )ξ( θ 1 θ 2 ) Two-point correlation function Small non-gaussian error Small variance σ W 2

24 Two-point correlation function of convergence ξ has a minimum at θ=15deg Best separation among the patches is ~15deg Separation angle

25 100 deg^2 area in total Each patch is separated by ~15deg separation

26 Sparse sampling 1x1 deg^2 patches are placed in a configuration 10x10 with a separation θsep

27 Cumulative SN up to l_max Maximum multipole

28 Cumulative SN up to l_max more sparse distribution is better Maximum multipole

29 Summary Sparse Sampling Belt Elongated Rectangular The above observational shapes have advantages at both large and small scales: larger-scale fluctuations can be scales (e.g. Kaiser 1986, 1998; Kilbinger & Schneider 2004; Chiang+ 2013) the non-gaussian error can be scales

30 Power Spectrum Covariance for Log-normal Field Survey window function Total survey area W θ = 1 0 inside the survey region otherwise S W = d 2 θ W(θ) Power spectrum convolved with the window function C W (l) = 1 S W d 2 q W l q 2 C(q) Circular-averaged band power spectrum CW l = 1 S W l l d 2 l N l d 2 q 2π 2 W l q 2 C(q) N l = 2πl l : # of modes in (l l 2, l + l 2)

31 Power spectrum covariance (Takada & Hu 2013) C ij = CW l i, C W l j = 1 S W 2 N l C W(l i ) 2 δ ij + T W (l i, l j ) Gaussian non-gaussian terms = C G ij+ C T0 ij + C σ W ij 1 / (survey area) C G ij = 1 S W 2 N l C W(l i ) 2 δ ij depends on the survey geometry (beat-coupling term) C T0 ij = 1 S W l l i d 2 l N li d2 l T(l, l, l, l ) l l j N lj C σ W ij = 4 κ 0 2 σ W 2 C l i C(l j ) σ W 2 = d 2 q 2π 2 W q 2 C(q) Variance of convergence in the survey region

32 The trispectrum can be written in term of power spectra for the log-normal field Hilbert+ (2011) showed that a four point correlation function of log-normal field can be written in terms of its two point function Power spectrum covariance can be written in terms of its power spectrum

33

34

35

36 Gaussian three independent data Non-Gaussian Cls are fully correlated at different k single data Non-Gaussian error depends on survey geometry Optimal survey geometry can reduce the non-gaussian errors

37 シミュレーションマップを使った解析 Cylindrically-averaged power spectrum estimator C W,r l = 1 S W 1 N l κ W,r(l ) 2 l l Mean power spectrum C W l = 1 N r C N W,r (l) r r=1 N r = 1000 : # of realizations Covariance matrix C ij = CW l i, C W l j = 1 N r 1 N r r=1 C W,r l i C W l i C W,r (l i ) C W l j Cumulative SN(Signal-to-Noise ratio) up to l max S N 2 = Cij 1 l<l max C W l i C W (l i )

38 Ray-tracing simulation (Taruya+ 2002)

39 Covariance matrix (diagonal elements) Variance Multipole

40 Covariance matrix (off-diagonal elements) C ij C ii C jj Multipole

Project: IEEE P Working Group for Wireless Personal Area Networks (WPANs)

Project: IEEE P Working Group for Wireless Personal Area Networks (WPANs) Project: IEEE P802.15 Working Group for Wireless Personal Area Networks (WPANs) Title: Channel Model for Intra-Device Communications Date Submitted: 15 January 2016 Source: Alexander Fricke, Thomas Kürner,

More information

Geometry Factors pt. II, Error Analysis, and Optimizations

Geometry Factors pt. II, Error Analysis, and Optimizations Geometry Factors pt. II, Error Analysis, and Optimizations Christopher Coppola November 11, 2014 Simulation Method Error Analysis Optimizations Simulation Goals -Calculate geometry factors -Optimize pressure

More information

IDENTIFYING BROAD AND NARROW FINANCIAL RISK FACTORS VIA CONVEX OPTIMIZATION: PART II

IDENTIFYING 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 information

Dynamic Response of Jackup Units Re-evaluation of SNAME 5-5A Four Methods

Dynamic Response of Jackup Units Re-evaluation of SNAME 5-5A Four Methods ISOPE 2010 Conference Beijing, China 24 June 2010 Dynamic Response of Jackup Units Re-evaluation of SNAME 5-5A Four Methods Xi Ying Zhang, Zhi Ping Cheng, Jer-Fang Wu and Chee Chow Kei ABS 1 Main Contents

More information

High Dimensional Bayesian Optimisation and Bandits via Additive Models

High Dimensional Bayesian Optimisation and Bandits via Additive Models 1/20 High Dimensional Bayesian Optimisation and Bandits via Additive Models Kirthevasan Kandasamy, Jeff Schneider, Barnabás Póczos ICML 15 July 8 2015 2/20 Bandits & Optimisation Maximum Likelihood inference

More information

Computer Exercise 2 Simulation

Computer Exercise 2 Simulation Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Fall 2017 Computer Exercise 2 Simulation This lab deals with pricing

More information

arxiv: v1 [astro-ph.im] 22 Jan 2009

arxiv: v1 [astro-ph.im] 22 Jan 2009 Astronomy & Astrophysics manuscript no. sparsity c ESO 009 January, 009 Optimal PSF modelling for weak lensing : complexity and sparsity arxiv:0901.3557v1 [astro-ph.im] Jan 009 S. Paulin-Henriksson 1,

More information

Log-Robust Portfolio Management

Log-Robust Portfolio Management Log-Robust Portfolio Management Dr. Aurélie Thiele Lehigh University Joint work with Elcin Cetinkaya and Ban Kawas Research partially supported by the National Science Foundation Grant CMMI-0757983 Dr.

More information

dt+ ρσ 2 1 ρ2 σ 2 κ i and that A is a rather lengthy expression that we may or may not need. (Brigo & Mercurio Lemma Thm , p. 135.

dt+ ρσ 2 1 ρ2 σ 2 κ i and that A is a rather lengthy expression that we may or may not need. (Brigo & Mercurio Lemma Thm , p. 135. A 2D Gaussian model (akin to Brigo & Mercurio Section 4.2) Suppose where ( κ1 0 dx(t) = 0 κ 2 r(t) = δ 0 +X 1 (t)+x 2 (t) )( X1 (t) X 2 (t) ) ( σ1 0 dt+ ρσ 2 1 ρ2 σ 2 )( dw Q 1 (t) dw Q 2 (t) ) In this

More information

Modelling strategies for bivariate circular data

Modelling strategies for bivariate circular data Modelling strategies for bivariate circular data John T. Kent*, Kanti V. Mardia, & Charles C. Taylor Department of Statistics, University of Leeds 1 Introduction On the torus there are two common approaches

More information

dt + ρσ 2 1 ρ2 σ 2 B i (τ) = 1 e κ iτ κ i

dt + ρσ 2 1 ρ2 σ 2 B i (τ) = 1 e κ iτ κ i A 2D Gaussian model (akin to Brigo & Mercurio Section 4.2) Suppose where dx(t) = ( κ1 0 0 κ 2 ) ( X1 (t) X 2 (t) In this case we find (BLACKBOARD) that r(t) = δ 0 + X 1 (t) + X 2 (t) ) ( σ1 0 dt + ρσ 2

More information

Multi-Period Trading via Convex Optimization

Multi-Period Trading via Convex Optimization Multi-Period Trading via Convex Optimization Stephen Boyd Enzo Busseti Steven Diamond Ronald Kahn Kwangmoo Koh Peter Nystrup Jan Speth Stanford University & Blackrock City University of Hong Kong September

More information

Copulas? What copulas? R. Chicheportiche & J.P. Bouchaud, CFM

Copulas? What copulas? R. Chicheportiche & J.P. Bouchaud, CFM Copulas? What copulas? R. Chicheportiche & J.P. Bouchaud, CFM Multivariate linear correlations Standard tool in risk management/portfolio optimisation: the covariance matrix R ij = r i r j Find the portfolio

More information

Nonlinear Manifold Learning for Financial Markets Integration

Nonlinear Manifold Learning for Financial Markets Integration Nonlinear Manifold Learning for Financial Markets Integration George Tzagkarakis 1 & Thomas Dionysopoulos 1,2 1 EONOS Investment Technologies, Paris (FR) 2 Dalton Strategic Partnership, London (UK) Nice,

More information

Structural credit risk models and systemic capital

Structural credit risk models and systemic capital Structural credit risk models and systemic capital Somnath Chatterjee CCBS, Bank of England November 7, 2013 Structural credit risk model Structural credit risk models are based on the notion that both

More information

9.1 Principal Component Analysis for Portfolios

9.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 information

Project 1: Double Pendulum

Project 1: Double Pendulum Final Projects Introduction to Numerical Analysis II http://www.math.ucsb.edu/ atzberg/winter2009numericalanalysis/index.html Professor: Paul J. Atzberger Due: Friday, March 20th Turn in to TA s Mailbox:

More information

Tangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford.

Tangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford. Tangent Lévy Models Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford June 24, 2010 6th World Congress of the Bachelier Finance Society Sergey

More information

A Stochastic Reserving Today (Beyond Bootstrap)

A Stochastic Reserving Today (Beyond Bootstrap) A Stochastic Reserving Today (Beyond Bootstrap) Presented by Roger M. Hayne, PhD., FCAS, MAAA Casualty Loss Reserve Seminar 6-7 September 2012 Denver, CO CAS Antitrust Notice The Casualty Actuarial Society

More information

An Online Appendix of Technical Trading: A Trend Factor

An Online Appendix of Technical Trading: A Trend Factor An Online Appendix of Technical Trading: A Trend Factor In this online appendix, we provide a comparative static analysis of the theoretical model as well as further robustness checks on the trend factor.

More information

Accelerated Stochastic Gradient Descent Praneeth Netrapalli MSR India

Accelerated Stochastic Gradient Descent Praneeth Netrapalli MSR India Accelerated Stochastic Gradient Descent Praneeth Netrapalli MSR India Presented at OSL workshop, Les Houches, France. Joint work with Prateek Jain, Sham M. Kakade, Rahul Kidambi and Aaron Sidford Linear

More information

Resource Allocation within Firms and Financial Market Dislocation: Evidence from Diversified Conglomerates

Resource Allocation within Firms and Financial Market Dislocation: Evidence from Diversified Conglomerates Resource Allocation within Firms and Financial Market Dislocation: Evidence from Diversified Conglomerates Gregor Matvos and Amit Seru (RFS, 2014) Corporate Finance - PhD Course 2017 Stefan Greppmair,

More information

Statistical and Computational Inverse Problems with Applications Part 5B: Electrical impedance tomography

Statistical 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 information

Machine Learning for Quantitative Finance

Machine Learning for Quantitative Finance Machine Learning for Quantitative Finance Fast derivative pricing Sofie Reyners Joint work with Jan De Spiegeleer, Dilip Madan and Wim Schoutens Derivative pricing is time-consuming... Vanilla option pricing

More information

Robust Longevity Risk Management

Robust Longevity Risk Management Robust Longevity Risk Management Hong Li a,, Anja De Waegenaere a,b, Bertrand Melenberg a,b a Department of Econometrics and Operations Research, Tilburg University b Netspar Longevity 10 3-4, September,

More information

SAQ KONTROLL AB Box 49306, STOCKHOLM, Sweden Tel: ; Fax:

SAQ KONTROLL AB Box 49306, STOCKHOLM, Sweden Tel: ; Fax: ProSINTAP - A Probabilistic Program for Safety Evaluation Peter Dillström SAQ / SINTAP / 09 SAQ KONTROLL AB Box 49306, 100 29 STOCKHOLM, Sweden Tel: +46 8 617 40 00; Fax: +46 8 651 70 43 June 1999 Page

More information

THE DISTRIBUTION OF LOAN PORTFOLIO VALUE * Oldrich Alfons Vasicek

THE DISTRIBUTION OF LOAN PORTFOLIO VALUE * Oldrich Alfons Vasicek HE DISRIBUION OF LOAN PORFOLIO VALUE * Oldrich Alfons Vasicek he amount of capital necessary to support a portfolio of debt securities depends on the probability distribution of the portfolio loss. Consider

More information

Robust Optimization Applied to a Currency Portfolio

Robust Optimization Applied to a Currency Portfolio Robust Optimization Applied to a Currency Portfolio R. Fonseca, S. Zymler, W. Wiesemann, B. Rustem Workshop on Numerical Methods and Optimization in Finance June, 2009 OUTLINE Introduction Motivation &

More information

MAKING OPTIMISATION TECHNIQUES ROBUST WITH AGNOSTIC RISK PARITY

MAKING OPTIMISATION TECHNIQUES ROBUST WITH AGNOSTIC RISK PARITY Technical Note May 2017 MAKING OPTIMISATION TECHNIQUES ROBUST WITH AGNOSTIC RISK PARITY Introduction The alternative investment industry is becoming ever more accessible to those wishing to diversify away

More information

IDENTIFYING BROAD AND NARROW FINANCIAL RISK FACTORS VIA CONVEX OPTIMIZATION: PART I

IDENTIFYING 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

Probabilistic Meshless Methods for Bayesian Inverse Problems. Jon Cockayne July 8, 2016

Probabilistic Meshless Methods for Bayesian Inverse Problems. Jon Cockayne July 8, 2016 Probabilistic Meshless Methods for Bayesian Inverse Problems Jon Cockayne July 8, 2016 1 Co-Authors Chris Oates Tim Sullivan Mark Girolami 2 What is PN? Many problems in mathematics have no analytical

More information

Amath 546/Econ 589 Univariate GARCH Models

Amath 546/Econ 589 Univariate GARCH Models Amath 546/Econ 589 Univariate GARCH Models Eric Zivot April 24, 2013 Lecture Outline Conditional vs. Unconditional Risk Measures Empirical regularities of asset returns Engle s ARCH model Testing for ARCH

More information

Calibration Lecture 4: LSV and Model Uncertainty

Calibration Lecture 4: LSV and Model Uncertainty Calibration Lecture 4: LSV and Model Uncertainty March 2017 Recap: Heston model Recall the Heston stochastic volatility model ds t = rs t dt + Y t S t dw 1 t, dy t = κ(θ Y t ) dt + ξ Y t dw 2 t, where

More information

Automated PSF measurement and homogenization in DESDM

Automated PSF measurement and homogenization in DESDM Automated PSF measurement and homogenization in DESDM E.Bertin (IAP) E. Bertin DES Munich meeting 05/2010 1 PSF homogenization History Science requirements PSFEX internals Point source selection PSF modeling

More information

Choice Probabilities. Logit Choice Probabilities Derivation. Choice Probabilities. Basic Econometrics in Transportation.

Choice Probabilities. Logit Choice Probabilities Derivation. Choice Probabilities. Basic Econometrics in Transportation. 1/31 Choice Probabilities Basic Econometrics in Transportation Logit Models Amir Samimi Civil Engineering Department Sharif University of Technology Primary Source: Discrete Choice Methods with Simulation

More information

12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006.

12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006. 12. Conditional heteroscedastic models (ARCH) MA6622, Ernesto Mordecki, CityU, HK, 2006. References for this Lecture: Robert F. Engle. Autoregressive Conditional Heteroscedasticity with Estimates of Variance

More information

Earnings Inequality and the Minimum Wage: Evidence from Brazil

Earnings Inequality and the Minimum Wage: Evidence from Brazil Earnings Inequality and the Minimum Wage: Evidence from Brazil Niklas Engbom June 16, 2016 Christian Moser World Bank-Bank of Spain Conference This project Shed light on drivers of earnings inequality

More information

A Structural Model of Continuous Workout Mortgages (Preliminary Do not cite)

A Structural Model of Continuous Workout Mortgages (Preliminary Do not cite) A Structural Model of Continuous Workout Mortgages (Preliminary Do not cite) Edward Kung UCLA March 1, 2013 OBJECTIVES The goal of this paper is to assess the potential impact of introducing alternative

More information

Quantitative Introduction ro Risk and Uncertainty in Business Module 5: Hypothesis Testing Examples

Quantitative Introduction ro Risk and Uncertainty in Business Module 5: Hypothesis Testing Examples Quantitative Introduction ro Risk and Uncertainty in Business Module 5: Hypothesis Testing Examples M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu

More information

Riemannian Geometry, Key to Homework #1

Riemannian Geometry, Key to Homework #1 Riemannian Geometry Key to Homework # Let σu v sin u cos v sin u sin v cos u < u < π < v < π be a parametrization of the unit sphere S {x y z R 3 x + y + z } Fix an angle < θ < π and consider the parallel

More information

Monte Carlo Methods for Uncertainty Quantification

Monte Carlo Methods for Uncertainty Quantification Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)

More information

A Correlated Sampling Method for Multivariate Normal and Log-normal Distributions

A 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 information

Portfolio Management and Optimal Execution via Convex Optimization

Portfolio Management and Optimal Execution via Convex Optimization Portfolio Management and Optimal Execution via Convex Optimization Enzo Busseti Stanford University April 9th, 2018 Problems portfolio management choose trades with optimization minimize risk, maximize

More information

Overnight Index Rate: Model, calibration and simulation

Overnight Index Rate: Model, calibration and simulation Research Article Overnight Index Rate: Model, calibration and simulation Olga Yashkir and Yuri Yashkir Cogent Economics & Finance (2014), 2: 936955 Page 1 of 11 Research Article Overnight Index Rate: Model,

More information

Practical Section 02 May 02, Part 1: Analytic transforms versus FFT algorithm. AnalBoxCar = 2*AB*BW*sin(2*pi*BW*f).

Practical Section 02 May 02, Part 1: Analytic transforms versus FFT algorithm. AnalBoxCar = 2*AB*BW*sin(2*pi*BW*f). 12.714 Practical Section 02 May 02, 2012 Part 1: Analytic transforms versus FFT algorithm (a) For a box car time domain signal with width 2*BW and amplitude BA, compare the analytic version of the Fourier

More information

The method of Maximum Likelihood.

The method of Maximum Likelihood. Maximum Likelihood The method of Maximum Likelihood. In developing the least squares estimator - no mention of probabilities. Minimize the distance between the predicted linear regression and the observed

More information

1 Explaining Labor Market Volatility

1 Explaining Labor Market Volatility Christiano Economics 416 Advanced Macroeconomics Take home midterm exam. 1 Explaining Labor Market Volatility The purpose of this question is to explore a labor market puzzle that has bedeviled business

More information

Chapter 7: Portfolio Theory

Chapter 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 information

Estimation Appendix to Dynamics of Fiscal Financing in the United States

Estimation Appendix to Dynamics of Fiscal Financing in the United States Estimation Appendix to Dynamics of Fiscal Financing in the United States Eric M. Leeper, Michael Plante, and Nora Traum July 9, 9. Indiana University. This appendix includes tables and graphs of additional

More information

ISI, RIN and MPN Modeling: Some Clarifications

ISI, RIN and MPN Modeling: Some Clarifications ISI, RIN and MPN Modeling: Some Clarifications Kasyapa Balemarthy Robert Lingle Jr. July 5, 01 ISI, MPN and RIN: the Big Picture : ISI with a single-mode VCSEL Worst-case received waveform due to single-mode

More information

Dynamic Portfolio Execution Detailed Proofs

Dynamic 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 information

A local RBF method based on a finite collocation approach

A local RBF method based on a finite collocation approach Boundary Elements and Other Mesh Reduction Methods XXXVIII 73 A local RBF method based on a finite collocation approach D. Stevens & H. Power Department of Mechanical Materials and Manufacturing Engineering,

More information

Equilibrium Asset Pricing: With Non-Gaussian Factors and Exponential Utilities

Equilibrium Asset Pricing: With Non-Gaussian Factors and Exponential Utilities Equilibrium Asset Pricing: With Non-Gaussian Factors and Exponential Utilities Dilip Madan Robert H. Smith School of Business University of Maryland Madan Birthday Conference September 29 2006 1 Motivation

More information

Regret-based Selection

Regret-based Selection Regret-based Selection David Puelz (UT Austin) Carlos M. Carvalho (UT Austin) P. Richard Hahn (Chicago Booth) May 27, 2017 Two problems 1. Asset pricing: What are the fundamental dimensions (risk factors)

More information

A model reduction approach to numerical inversion for parabolic partial differential equations

A 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 information

Parameter estimation in SDE:s

Parameter estimation in SDE:s Lund University Faculty of Engineering Statistics in Finance Centre for Mathematical Sciences, Mathematical Statistics HT 2011 Parameter estimation in SDE:s This computer exercise concerns some estimation

More information

ARCH Models and Financial Applications

ARCH Models and Financial Applications Christian Gourieroux ARCH Models and Financial Applications With 26 Figures Springer Contents 1 Introduction 1 1.1 The Development of ARCH Models 1 1.2 Book Content 4 2 Linear and Nonlinear Processes 5

More information

LECTURE 3: FREE CENTRAL LIMIT THEOREM AND FREE CUMULANTS

LECTURE 3: FREE CENTRAL LIMIT THEOREM AND FREE CUMULANTS LECTURE 3: FREE CENTRAL LIMIT THEOREM AND FREE CUMULANTS Recall from Lecture 2 that if (A, φ) is a non-commutative probability space and A 1,..., A n are subalgebras of A which are free with respect to

More information

Extended Libor Models and Their Calibration

Extended Libor Models and Their Calibration Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Vienna, 16 November 2007 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Vienna, 16 November

More information

Impact of Calendar Effects. in the Volatility of Vale Shares

Impact of Calendar Effects. in the Volatility of Vale Shares Journal of Finance and Investment Analysis, vol.2, no.3, 2013, 1-16 ISSN: 2241-0988 (print version), 2241-0996 (online) Scienpress Ltd, 2013 Impact of Calendar Effects in the Volatility of Vale Shares

More information

EX-POST VERIFICATION OF PREDICTION MODELS OF WAGE DISTRIBUTIONS

EX-POST VERIFICATION OF PREDICTION MODELS OF WAGE DISTRIBUTIONS EX-POST VERIFICATION OF PREDICTION MODELS OF WAGE DISTRIBUTIONS LUBOŠ MAREK, MICHAL VRABEC University of Economics, Prague, Faculty of Informatics and Statistics, Department of Statistics and Probability,

More information

An Introduction to Statistical Extreme Value Theory

An Introduction to Statistical Extreme Value Theory An Introduction to Statistical Extreme Value Theory Uli Schneider Geophysical Statistics Project, NCAR January 26, 2004 NCAR Outline Part I - Two basic approaches to extreme value theory block maxima,

More information

Application of Moment Expansion Method to Option Square Root Model

Application of Moment Expansion Method to Option Square Root Model Application of Moment Expansion Method to Option Square Root Model Yun Zhou Advisor: Professor Steve Heston University of Maryland May 5, 2009 1 / 19 Motivation Black-Scholes Model successfully explain

More information

Modeling Volatility Risk in Equity Options: a Cross-sectional approach

Modeling Volatility Risk in Equity Options: a Cross-sectional approach ICBI Global Derivatives, Amsterdam, 2014 Modeling Volatility Risk in Equity Options: a Cross-sectional approach Marco Avellaneda NYU & Finance Concepts Doris Dobi* NYU * This work is part of Doris Dobi

More information

Modeling Co-movements and Tail Dependency in the International Stock Market via Copulae

Modeling Co-movements and Tail Dependency in the International Stock Market via Copulae Modeling Co-movements and Tail Dependency in the International Stock Market via Copulae Katja Ignatieva, Eckhard Platen Bachelier Finance Society World Congress 22-26 June 2010, Toronto K. Ignatieva, E.

More information

Multilevel quasi-monte Carlo path simulation

Multilevel quasi-monte Carlo path simulation Multilevel quasi-monte Carlo path simulation Michael B. Giles and Ben J. Waterhouse Lluís Antoni Jiménez Rugama January 22, 2014 Index 1 Introduction to MLMC Stochastic model Multilevel Monte Carlo Milstein

More information

In terms of covariance the Markowitz portfolio optimisation problem is:

In terms of covariance the Markowitz portfolio optimisation problem is: Markowitz portfolio optimisation Solver To use Solver to solve the quadratic program associated with tracing out the efficient frontier (unconstrained efficient frontier UEF) in Markowitz portfolio optimisation

More information

ASTIN Colloquium Understanding Split Credibility. Ira Robbin, PhD AVP and Senior Pricing Actuary Endurance US Insurance Operations

ASTIN Colloquium Understanding Split Credibility. Ira Robbin, PhD AVP and Senior Pricing Actuary Endurance US Insurance Operations ASTIN Colloquium Understanding Split Credibility Ira Robbin, PhD AVP and Senior Pricing Actuary Endurance US Insurance Operations Ground Rules Follow US Anti-trust Laws, s il vous plait! Violators will

More information

Imperfect Information and Market Segmentation Walsh Chapter 5

Imperfect Information and Market Segmentation Walsh Chapter 5 Imperfect Information and Market Segmentation Walsh Chapter 5 1 Why Does Money Have Real Effects? Add market imperfections to eliminate short-run neutrality of money Imperfect information keeps price from

More information

Final exam solutions

Final exam solutions EE365 Stochastic Control / MS&E251 Stochastic Decision Models Profs. S. Lall, S. Boyd June 5 6 or June 6 7, 2013 Final exam solutions This is a 24 hour take-home final. Please turn it in to one of the

More information

CHAPTER II LITERATURE STUDY

CHAPTER II LITERATURE STUDY CHAPTER II LITERATURE STUDY 2.1. Risk Management Monetary crisis that strike Indonesia during 1998 and 1999 has caused bad impact to numerous government s and commercial s bank. Most of those banks eventually

More information

I. Return Calculations (20 pts, 4 points each)

I. Return Calculations (20 pts, 4 points each) University of Washington Winter 015 Department of Economics Eric Zivot Econ 44 Midterm Exam Solutions This is a closed book and closed note exam. However, you are allowed one page of notes (8.5 by 11 or

More information

Mathematics of Finance Final Preparation December 19. To be thoroughly prepared for the final exam, you should

Mathematics of Finance Final Preparation December 19. To be thoroughly prepared for the final exam, you should Mathematics of Finance Final Preparation December 19 To be thoroughly prepared for the final exam, you should 1. know how to do the homework problems. 2. be able to provide (correct and complete!) definitions

More information

Direct Methods for linear systems Ax = b basic point: easy to solve triangular systems

Direct Methods for linear systems Ax = b basic point: easy to solve triangular systems NLA p.1/13 Direct Methods for linear systems Ax = b basic point: easy to solve triangular systems... 0 0 0 etc. a n 1,n 1 x n 1 = b n 1 a n 1,n x n solve a n,n x n = b n then back substitution: takes n

More information

Nobel Symposium Money and Banking

Nobel Symposium Money and Banking Nobel Symposium Money and Banking https://www.houseoffinance.se/nobel-symposium May 26-28, 2018 Clarion Hotel Sign, Stockholm Money and Banking: Some DSGE Challenges Nobel Symposium on Money and Banking

More information

Graduate 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 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 information

IEOR E4703: Monte-Carlo Simulation

IEOR E4703: Monte-Carlo Simulation IEOR E4703: Monte-Carlo Simulation Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com

More information

Operational Risk Aggregation

Operational 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 information

Lecture IV Portfolio management: Efficient portfolios. Introduction to Finance Mathematics Fall Financial mathematics

Lecture IV Portfolio management: Efficient portfolios. Introduction to Finance Mathematics Fall Financial mathematics Lecture IV Portfolio management: Efficient portfolios. Introduction to Finance Mathematics Fall 2014 Reduce the risk, one asset Let us warm up by doing an exercise. We consider an investment with σ 1 =

More information

Genetics and/of basket options

Genetics and/of basket options Genetics and/of basket options Wolfgang Karl Härdle Elena Silyakova Ladislaus von Bortkiewicz Chair of Statistics Humboldt-Universität zu Berlin http://lvb.wiwi.hu-berlin.de Motivation 1-1 Basket derivatives

More information

The rth moment of a real-valued random variable X with density f(x) is. x r f(x) dx

The rth moment of a real-valued random variable X with density f(x) is. x r f(x) dx 1 Cumulants 1.1 Definition The rth moment of a real-valued random variable X with density f(x) is µ r = E(X r ) = x r f(x) dx for integer r = 0, 1,.... The value is assumed to be finite. Provided that

More information

Option-Implied Information in Asset Allocation Decisions

Option-Implied Information in Asset Allocation Decisions Option-Implied Information in Asset Allocation Decisions Grigory Vilkov Goethe University Frankfurt 12 December 2012 Grigory Vilkov Option-Implied Information in Asset Allocation 12 December 2012 1 / 32

More information

Learning from cycles

Learning from cycles Learning from cycles Mehmet Dalkır April 30, 2007 preliminary copy Abstract In this paper, Bayesian learning rules with distributions on non-euclidean spaces are discussed. The findings of this study is

More information

Heterogeneous Firm, Financial Market Integration and International Risk Sharing

Heterogeneous Firm, Financial Market Integration and International Risk Sharing Heterogeneous Firm, Financial Market Integration and International Risk Sharing Ming-Jen Chang, Shikuan Chen and Yen-Chen Wu National DongHwa University Thursday 22 nd November 2018 Department of Economics,

More information

Multi-dimensional Term Structure Models

Multi-dimensional Term Structure Models Multi-dimensional Term Structure Models We will focus on the affine class. But first some motivation. A generic one-dimensional model for zero-coupon yields, y(t; τ), looks like this dy(t; τ) =... dt +

More information

Estimating Market Power in Differentiated Product Markets

Estimating Market Power in Differentiated Product Markets Estimating Market Power in Differentiated Product Markets Metin Cakir Purdue University December 6, 2010 Metin Cakir (Purdue) Market Equilibrium Models December 6, 2010 1 / 28 Outline Outline Estimating

More information

Course information FN3142 Quantitative finance

Course information FN3142 Quantitative finance Course information 015 16 FN314 Quantitative finance This course is aimed at students interested in obtaining a thorough grounding in market finance and related empirical methods. Prerequisite If taken

More information

A model reduction approach to numerical inversion for parabolic partial differential equations

A 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 information

Determining source cumulants in femtoscopy with Gram-Charlier and Edgeworth series

Determining source cumulants in femtoscopy with Gram-Charlier and Edgeworth series Determining source cumulants in femtoscopy with Gram-Charlier and Edgeworth series M.B. de Kock a H.C. Eggers a J. Schmiegel b a University of Stellenbosch, South Africa b Aarhus University, Denmark VI

More information

Simulation Wrap-up, Statistics COS 323

Simulation Wrap-up, Statistics COS 323 Simulation Wrap-up, Statistics COS 323 Today Simulation Re-cap Statistics Variance and confidence intervals for simulations Simulation wrap-up FYI: No class or office hours Thursday Simulation wrap-up

More information

Premia 14 HESTON MODEL CALIBRATION USING VARIANCE SWAPS PRICES

Premia 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 information

Phys 731: String Theory - Assignment 2

Phys 731: String Theory - Assignment 2 Phys 73: String Theory - Assignment Andrzej Pokraka October 8, 07 Problem Part a) We are given the action S S P + S J, S P 4πα d σδ ab a Xσ) b Xσ) S J d σjσ)xσ). Here, σ a, ) like in a regular QFT. First,

More information

HIGH ORDER DISCONTINUOUS GALERKIN METHODS FOR 1D PARABOLIC EQUATIONS. Ahmet İzmirlioğlu. BS, University of Pittsburgh, 2004

HIGH ORDER DISCONTINUOUS GALERKIN METHODS FOR 1D PARABOLIC EQUATIONS. Ahmet İzmirlioğlu. BS, University of Pittsburgh, 2004 HIGH ORDER DISCONTINUOUS GALERKIN METHODS FOR D PARABOLIC EQUATIONS by Ahmet İzmirlioğlu BS, University of Pittsburgh, 24 Submitted to the Graduate Faculty of Art and Sciences in partial fulfillment of

More information

Synthetic CDO Pricing Using the Student t Factor Model with Random Recovery

Synthetic CDO Pricing Using the Student t Factor Model with Random Recovery Synthetic CDO Pricing Using the Student t Factor Model with Random Recovery UNSW Actuarial Studies Research Symposium 2006 University of New South Wales Tom Hoedemakers Yuri Goegebeur Jurgen Tistaert Tom

More information

Chapter 5 The Production Process and Costs

Chapter 5 The Production Process and Costs Managerial Economics & Business Strategy Chapter 5 The Production Process and Costs McGraw-Hill/Irwin Copyright 2010 by the McGraw-Hill Companies, Inc. All rights reserved. I. Production Analysis Overview

More information

Modeling the extremes of temperature time series. Debbie J. Dupuis Department of Decision Sciences HEC Montréal

Modeling the extremes of temperature time series. Debbie J. Dupuis Department of Decision Sciences HEC Montréal Modeling the extremes of temperature time series Debbie J. Dupuis Department of Decision Sciences HEC Montréal Outline Fig. 1: S&P 500. Daily negative returns (losses), Realized Variance (RV) and Jump

More information

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay. Solutions to Final Exam

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay. Solutions to Final Exam The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (40 points) Answer briefly the following questions. 1. Describe

More information

An Implementation of Markov Regime Switching GARCH Models in Matlab

An Implementation of Markov Regime Switching GARCH Models in Matlab An Implementation of Markov Regime Switching GARCH Models in Matlab Thomas Chuffart Aix-Marseille University (Aix-Marseille School of Economics), CNRS & EHESS Abstract MSGtool is a MATLAB toolbox which

More information

IMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS

IMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS IMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS DENIS BELOMESTNY AND MARKUS REISS 1. Introduction The aim of this report is to describe more precisely how the spectral calibration method

More information

Bayesian Finance. Christa Cuchiero, Irene Klein, Josef Teichmann. Obergurgl 2017

Bayesian Finance. Christa Cuchiero, Irene Klein, Josef Teichmann. Obergurgl 2017 Bayesian Finance Christa Cuchiero, Irene Klein, Josef Teichmann Obergurgl 2017 C. Cuchiero, I. Klein, and J. Teichmann Bayesian Finance Obergurgl 2017 1 / 23 1 Calibrating a Bayesian model: a first trial

More information