Statistical analysis of the inverse problem
|
|
- Michael Gibson
- 5 years ago
- Views:
Transcription
1 Statistical analysis of the inverse problem Hongyuan Cao and Jayanta Kumar Pal SAMSI SAMSI/CRSC Undergraduate Workshop at NCSU May 24, 2007
2 Outline 1 Introduction 2 Preliminary analysis Estimation of parameters Residual analysis
3 Introduction Outline 1 Introduction 2 Preliminary analysis Estimation of parameters Residual analysis
4 Introduction Estimation problem We discuss the statistical procedure to solve the inverse problem. We use the spring model : d 2 y(t) dt 2 + C dy(t) dt + Ky(t) = 0. We have observations : (t 1, y 1 ),..., (t m, y m ). The goal is to estimate the unknown parameters C and K.
5 Introduction Estimation problem We discuss the statistical procedure to solve the inverse problem. We use the spring model : d 2 y(t) dt 2 + C dy(t) dt + Ky(t) = 0. We have observations : (t 1, y 1 ),..., (t m, y m ). The goal is to estimate the unknown parameters C and K.
6 Introduction Estimation problem We discuss the statistical procedure to solve the inverse problem. We use the spring model : d 2 y(t) dt 2 + C dy(t) dt + Ky(t) = 0. We have observations : (t 1, y 1 ),..., (t m, y m ). The goal is to estimate the unknown parameters C and K.
7 Introduction Differential equation For any C, K, the differential equation has a solution given initial value y(0) = y 0. We call that particular solution y(t; C, K). Notice the dependence on C and K! For fixed C, K it is just a function of t, as we would have liked.
8 Introduction Differential equation For any C, K, the differential equation has a solution given initial value y(0) = y 0. We call that particular solution y(t; C, K). Notice the dependence on C and K! For fixed C, K it is just a function of t, as we would have liked.
9 Introduction Cost function : least squares Fix C and K. We define the cost function L(C, K) = m (y i y(t i ; C, K)) 2 i=1 We seek to minimize this cost function L(C, K) over all possible values of C and K. The least squares estimator : Ĉ, ˆK.
10 Introduction Cost function : least squares Fix C and K. We define the cost function L(C, K) = m (y i y(t i ; C, K)) 2 i=1 We seek to minimize this cost function L(C, K) over all possible values of C and K. The least squares estimator : Ĉ, ˆK.
11 Introduction Non-linear regression This is a non-linear regression model. Assumption : we have equal variance measurement errors. Model y i = y(t i ; C, K) + ɛ i ɛ i are independent, identically distributed Normal errors with mean zero and variance σ 2. This is non-linear, since the solution to the ODE is a function involving exponentials and trigonometric forms.
12 Introduction Non-linear regression This is a non-linear regression model. Assumption : we have equal variance measurement errors. Model y i = y(t i ; C, K) + ɛ i ɛ i are independent, identically distributed Normal errors with mean zero and variance σ 2. This is non-linear, since the solution to the ODE is a function involving exponentials and trigonometric forms.
13 Introduction Linear and non-linear regression Examples : Linear model y = β 0 + β 1 x + β 2 x 2 y = β 0 + β 1 cos(x) + β 2 sin(x) + β 3 exp( x) Example : non-linear model y = exp( β 1 x) + α cos(β 2 x) The non-linearity is in terms of the parameters, not the other variables.
14 Introduction Linear and non-linear regression Examples : Linear model y = β 0 + β 1 x + β 2 x 2 y = β 0 + β 1 cos(x) + β 2 sin(x) + β 3 exp( x) Example : non-linear model y = exp( β 1 x) + α cos(β 2 x) The non-linearity is in terms of the parameters, not the other variables.
15 Outline 1 Introduction 2 Preliminary analysis Estimation of parameters Residual analysis
16 Preliminary analysis MATLAB commands I saved all my MATLAB commands in the file costminanalysis.m
17 Preliminary analysis Preliminary analysis Extracting the data from the.mat file U = load( data1 ) time = U.trace_x ; disp = U.trace_y(4,:); N = length(time); Plotting the data plot(time, disp); xlabel( t ); ylabel( y(t) );
18 Preliminary analysis Preliminary analysis Extracting the data from the.mat file U = load( data1 ) time = U.trace_x ; disp = U.trace_y(4,:); N = length(time); Plotting the data plot(time, disp); xlabel( t ); ylabel( y(t) );
19 Preliminary analysis Plot of y(t) vs t 8 x y(t) t
20 Preliminary analysis Cleaning the data for our purpose Observe that the data has some pattern at the beginning, that we will truncate. We will start at the point where we have maximum amplitude and the data dampens out gradually. Also, the variable disp is in row vector form. We change it to column vector by transposing. disp = disp ; [maxd,index] = max(disp); time = time(index:n)-time(index); disp = disp(index:n);
21 Preliminary analysis Cleaning the data for our purpose Observe that the data has some pattern at the beginning, that we will truncate. We will start at the point where we have maximum amplitude and the data dampens out gradually. Also, the variable disp is in row vector form. We change it to column vector by transposing. disp = disp ; [maxd,index] = max(disp); time = time(index:n)-time(index); disp = disp(index:n);
22 Preliminary analysis Centering the data We also just center the data so that the mean is 0, which is necessary for the analysis. disp = disp - mean(disp); We plot the new data again.
23 Preliminary analysis Plot of y(t) 8 x y(t) t
24 Estimation of parameters Solving a differential equation We write ẏ = dy dt, and ÿ = d2 y. So, the ODE of dt 2 spring-mass-dashpot becomes ÿ + Cẏ + Ky = 0 Or alternately, we write in a vector form, i.e. ( ) ( ) d y1 y2 = dt y 2 Ky 1 Cy 2
25 Estimation of parameters Solving a differential equation We write ẏ = dy dt, and ÿ = d2 y. So, the ODE of dt 2 spring-mass-dashpot becomes ÿ + Cẏ + Ky = 0 Or alternately, we write in a vector form, i.e. ( ) ( ) d y1 y2 = dt y 2 Ky 1 Cy 2
26 Estimation of parameters MATLAB ODE solver Here we find the function y(t; C, K) which solves the equation ÿ + Cẏ + Ky = 0 under the initial condition y(0) = y 0. The routine ode23 solves the equation system z = ode_model(t,z) when we call it. First, we code the ODE system in MATLAB function dy = ode_model(t,y,c,k) dy = zeros(2,1); dy(1) = y(2); dy(2) = -K*y(1) - C*y(2) Save this file in the working directory as ode_model.m.
27 Estimation of parameters MATLAB ODE solver Here we find the function y(t; C, K) which solves the equation ÿ + Cẏ + Ky = 0 under the initial condition y(0) = y 0. The routine ode23 solves the equation system z = ode_model(t,z) when we call it. First, we code the ODE system in MATLAB function dy = ode_model(t,y,c,k) dy = zeros(2,1); dy(1) = y(2); dy(2) = -K*y(1) - C*y(2) Save this file in the working directory as ode_model.m.
28 Estimation of parameters Cost function of the vibrating beam problem Next, we solve the ODE using the MATLAB routine, using the initial values, and also calculate the cost L(C, K) as defined before. The following function takes as an input the parameter values C and K with the data, and computes the cost L(C, K) = m (y i y(t i ; C, K)) 2 i=1 as defined before.
29 Estimation of parameters The MATLAB function to compute L(C, K) First we define a function to compute L(C, K). function [cost,d_model]=cost_beam(ck,time,disp) x0 = [disp(1) 0]; C = CK(1); K = CK(2); [t,x] = ode23(@ode_model, time, x0, [], C, K); d_model = x(:,1); cost = sum((d_model-disp).^2); Save this file as cost_beam.m
30 Estimation of parameters Plot of the function L(C, K) Let us plot the function over a wide range of values for C and K. C0 = [0:.05:2]; K0 = [1000:2:2000]; LCK = zeros(length(c0), length(k0)); for i = 1:length(C0) for j = 1:length(K0) CK0 = [C0(i);K0(j)]; [c,d] = cost_beam(ck0,time,disp); LCK(i,j) = c; end end surf(c0,k0, LCK ) xlabel( C ) ylabel( K )
31 Estimation of parameters Plot of L(C, K) 1.2 x K C
32 Estimation of parameters Optimization routine We use the MATLAB routine to minimize the cost function L(C, K) over a wide range of values of C and K. We have to give some initialization values. C0 = 1; K0 = 1500; CK0 = [C0;K0]; [CK,cost]=fminsearch(@cost_beam,CK0,[],time,disp C = CK(1) K = CK(2) Finally, we get the estimators Ĉ and ˆK. Ĉ =.9052, ˆK =
33 Estimation of parameters Optimization routine We use the MATLAB routine to minimize the cost function L(C, K) over a wide range of values of C and K. We have to give some initialization values. C0 = 1; K0 = 1500; CK0 = [C0;K0]; [CK,cost]=fminsearch(@cost_beam,CK0,[],time,disp C = CK(1) K = CK(2) Finally, we get the estimators Ĉ and ˆK. Ĉ =.9052, ˆK =
34 Estimation of parameters Optimization routine We use the MATLAB routine to minimize the cost function L(C, K) over a wide range of values of C and K. We have to give some initialization values. C0 = 1; K0 = 1500; CK0 = [C0;K0]; [CK,cost]=fminsearch(@cost_beam,CK0,[],time,disp C = CK(1) K = CK(2) Finally, we get the estimators Ĉ and ˆK. Ĉ =.9052, ˆK =
35 Estimation of parameters Estimation of σ 2 Recall that we have the model y i = y(t i ; C, K) + ɛ i ɛ i are independent, identically distributed Normal errors with mean zero and variance σ 2. An estimate of σ 2 is 1 m 2 m (y i y(t i, Ĉ, ˆK)) 2 i=1
36 Estimation of parameters Estimation of σ 2 Recall that we have the model y i = y(t i ; C, K) + ɛ i ɛ i are independent, identically distributed Normal errors with mean zero and variance σ 2. An estimate of σ 2 is 1 m 2 m (y i y(t i, Ĉ, ˆK)) 2 i=1
37 Estimation of parameters To get the residuals The MATLAB code is [c, d_model] =cost_beam(ck, time,disp); d_res = disp - d_model; m = length(disp); sigma2 = sum(d_res.^2)/m; We save the residuals for exploratory analysis. We estimate ˆσ = from the data
38 Residual analysis Plotting the residuals We start plotting the residuals to see whether they have any pattern. Residual vs Time plot(time,d_res) xlabel( time ) ylabel( residuals )
39 Residual analysis Residuals against the time 8 x residuals time The residuals show some pattern, and the independence assumption seems tenuous!!
40 Residual analysis Fitted vs residuals Next, we check for the homoscedasticity assumption : that is, whether the residuals have the same variance or not. Residual vs Fitted plot(d_model,d_res, k. ) xlabel( fitted ) ylabel( residuals )
41 Residual analysis Checking homoscedasticity 8 x residuals fitted x 10 5 Here also, some pattern, though less evident, can be observed.
42 Residual analysis Quantile-quantile plot Finally, we check the normality assumptions. We use quantile-quantile plot for that. We plot the sorted residuals against the quantiles for the normal distribution. We expect the points stay close to a straight line. qqplot(d_res) ylabel( residuals )
43 Residual analysis Q-Q plot 8 x 10 5 QQ Plot of Sample Data versus Standard Normal residuals Standard Normal Quantiles The assumption of normality is not that strong too!!!
44 Residual analysis Improvement for the model? The spring model does not seem so appropriate for the data. Residual shows dependence structure. Variances may not be equal. Normality assumption may not hold. We may have the same underlying model, but with different assumptions for the error structure. Or we may have an altogether different model. Is there a possible improvement on the physical spring model we use here?
45 Residual analysis Improvement for the model? The spring model does not seem so appropriate for the data. Residual shows dependence structure. Variances may not be equal. Normality assumption may not hold. We may have the same underlying model, but with different assumptions for the error structure. Or we may have an altogether different model. Is there a possible improvement on the physical spring model we use here?
46 Residual analysis Improvement for the model? The spring model does not seem so appropriate for the data. Residual shows dependence structure. Variances may not be equal. Normality assumption may not hold. We may have the same underlying model, but with different assumptions for the error structure. Or we may have an altogether different model. Is there a possible improvement on the physical spring model we use here?
February 2 Math 2335 sec 51 Spring 2016
February 2 Math 2335 sec 51 Spring 2016 Section 3.1: Root Finding, Bisection Method Many problems in the sciences, business, manufacturing, etc. can be framed in the form: Given a function f (x), find
More informationFeb. 4 Math 2335 sec 001 Spring 2014
Feb. 4 Math 2335 sec 001 Spring 2014 Propagated Error in Function Evaluation Let f (x) be some differentiable function. Suppose x A is an approximation to x T, and we wish to determine the function value
More informationProblem Set 3 Due by Sat 12:00, week 3
Business 35150 John H. Cochrane Problem Set 3 Due by Sat 12:00, week 3 Part I. Reading questions: These refer to the reading assignment in the syllabus. Please hand in short answers. Where appropriate,
More informationErrata, Mahler Study Aids for Exam 3/M, Spring 2010 HCM, 1/26/13 Page 1
Errata, Mahler Study Aids for Exam 3/M, Spring 2010 HCM, 1/26/13 Page 1 1B, p. 72: (60%)(0.39) + (40%)(0.75) = 0.534. 1D, page 131, solution to the first Exercise: 2.5 2.5 λ(t) dt = 3t 2 dt 2 2 = t 3 ]
More informationTHE UNIVERSITY OF TEXAS AT AUSTIN McCombs School of Business
THE UNIVERSIT OF TEXAS AT AUSTIN McCombs School of Business STA 37.5 Tom Shively SIMPLE EXPONENTIAL SMOOTHING MODELS The statistical model for simple exponential smoothing is t = M t- + ε t ε t iid N(0,
More informationECON 6022B Problem Set 2 Suggested Solutions Fall 2011
ECON 60B Problem Set Suggested Solutions Fall 0 September 7, 0 Optimal Consumption with A Linear Utility Function (Optional) Similar to the example in Lecture 3, the household lives for two periods and
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 informationMultiple regression - a brief introduction
Multiple regression - a brief introduction Multiple regression is an extension to regular (simple) regression. Instead of one X, we now have several. Suppose, for example, that you are trying to predict
More informationPortfolio Optimization. Prof. Daniel P. Palomar
Portfolio Optimization Prof. Daniel P. Palomar The Hong Kong University of Science and Technology (HKUST) MAFS6010R- Portfolio Optimization with R MSc in Financial Mathematics Fall 2018-19, HKUST, Hong
More informationLecture 3: Factor models in modern portfolio choice
Lecture 3: Factor models in modern portfolio choice Prof. Massimo Guidolin Portfolio Management Spring 2016 Overview The inputs of portfolio problems Using the single index model Multi-index models Portfolio
More informationECE 567 STATISTICAL SIGNAL PROCESSING FALL Homework Assignment #4 Solutions
ECE 567 STATISTICAL SIGNAL PROCESSING FALL 017 Homework Assignment #4 Solutions 1. a As noted in the solution to problem a of Homework Assignment #3, the test is snxn b if A 0 but is if A 0, so that there
More informationBooth School of Business, University of Chicago Business 41202, Spring Quarter 2014, Mr. Ruey S. Tsay. Solutions to Midterm
Booth School of Business, University of Chicago Business 41202, Spring Quarter 2014, Mr. Ruey S. Tsay Solutions to Midterm Problem A: (30 pts) Answer briefly the following questions. Each question has
More information(5) Multi-parameter models - Summarizing the posterior
(5) Multi-parameter models - Summarizing the posterior Spring, 2017 Models with more than one parameter Thus far we have studied single-parameter models, but most analyses have several parameters For example,
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2009, Mr. Ruey S. Tsay. Solutions to Final Exam
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2009, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (42 pts) Answer briefly the following questions. 1. Questions
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 informationIntroduction to Population Modeling
Introduction to Population Modeling In addition to estimating the size of a population, it is often beneficial to estimate how the population size changes over time. Ecologists often uses models to create
More informationLecture 5: Fundamentals of Statistical Analysis and Distributions Derived from Normal Distributions
Lecture 5: Fundamentals of Statistical Analysis and Distributions Derived from Normal Distributions ELE 525: Random Processes in Information Systems Hisashi Kobayashi Department of Electrical Engineering
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay. Solutions to Final Exam.
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (32 pts) Answer briefly the following questions. 1. Suppose
More informationEC316a: Advanced Scientific Computation, Fall Discrete time, continuous state dynamic models: solution methods
EC316a: Advanced Scientific Computation, Fall 2003 Notes Section 4 Discrete time, continuous state dynamic models: solution methods We consider now solution methods for discrete time models in which decisions
More informationEconometric Methods for Valuation Analysis
Econometric Methods for Valuation Analysis Margarita Genius Dept of Economics M. Genius (Univ. of Crete) Econometric Methods for Valuation Analysis Cagliari, 2017 1 / 26 Correlation Analysis Simple Regression
More informationOptimal Portfolio Selection
Optimal Portfolio Selection We have geometrically described characteristics of the optimal portfolio. Now we turn our attention to a methodology for exactly identifying the optimal portfolio given a set
More informationLecture 3: Review of Probability, MATLAB, Histograms
CS 4980/6980: Introduction to Data Science c Spring 2018 Lecture 3: Review of Probability, MATLAB, Histograms Instructor: Daniel L. Pimentel-Alarcón Scribed and Ken Varghese This is preliminary work and
More informationChapter 6. Transformation of Variables
6.1 Chapter 6. Transformation of Variables 1. Need for transformation 2. Power transformations: Transformation to achieve linearity Transformation to stabilize variance Logarithmic transformation MACT
More informationOptimal Trading Strategy With Optimal Horizon
Optimal Trading Strategy With Optimal Horizon Financial Math Festival Florida State University March 1, 2008 Edward Qian PanAgora Asset Management Trading An Integral Part of Investment Process Return
More informationLecture 12: The Bootstrap
Lecture 12: The Bootstrap Reading: Chapter 5 STATS 202: Data mining and analysis October 20, 2017 1 / 16 Announcements Midterm is on Monday, Oct 30 Topics: chapters 1-5 and 10 of the book everything until
More informationLimit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies
Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies George Tauchen Duke University Viktor Todorov Northwestern University 2013 Motivation
More informationREGIONAL WORKSHOP ON TRAFFIC FORECASTING AND ECONOMIC PLANNING
International Civil Aviation Organization 27/8/10 WORKING PAPER REGIONAL WORKSHOP ON TRAFFIC FORECASTING AND ECONOMIC PLANNING Cairo 2 to 4 November 2010 Agenda Item 3 a): Forecasting Methodology (Presented
More informationA Monte Carlo Based Analysis of Optimal Design Criteria
A Monte Carlo Based Analysis of Optimal Design Criteria H. T. Banks, Kathleen J. Holm and Franz Kappel Center for Quantitative Sciences in Biomedicine Center for Research in Scientific Computation North
More informationMultiple Regression. Review of Regression with One Predictor
Fall Semester, 2001 Statistics 621 Lecture 4 Robert Stine 1 Preliminaries Multiple Regression Grading on this and other assignments Assignment will get placed in folder of first member of Learning Team.
More informationAnalysis of Variance in Matrix form
Analysis of Variance in Matrix form The ANOVA table sums of squares, SSTO, SSR and SSE can all be expressed in matrix form as follows. week 9 Multiple Regression A multiple regression model is a model
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Simulation Efficiency and an Introduction to Variance Reduction Methods Martin Haugh Department of Industrial Engineering and Operations Research Columbia University
More informationAgricultural and Applied Economics 637 Applied Econometrics II
Agricultural and Applied Economics 637 Applied Econometrics II Assignment I Using Search Algorithms to Determine Optimal Parameter Values in Nonlinear Regression Models (Due: February 3, 2015) (Note: Make
More informationSystems of Ordinary Differential Equations. Lectures INF2320 p. 1/48
Systems of Ordinary Differential Equations Lectures INF2320 p. 1/48 Lectures INF2320 p. 2/48 ystems of ordinary differential equations Last two lectures we have studied models of the form y (t) = F(y),
More informationThe Use of Importance Sampling to Speed Up Stochastic Volatility Simulations
The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations Stan Stilger June 6, 1 Fouque and Tullie use importance sampling for variance reduction in stochastic volatility simulations.
More informationThe Two Sample T-test with One Variance Unknown
The Two Sample T-test with One Variance Unknown Arnab Maity Department of Statistics, Texas A&M University, College Station TX 77843-343, U.S.A. amaity@stat.tamu.edu Michael Sherman Department of Statistics,
More informationQuantile Regression due to Skewness. and Outliers
Applied Mathematical Sciences, Vol. 5, 2011, no. 39, 1947-1951 Quantile Regression due to Skewness and Outliers Neda Jalali and Manoochehr Babanezhad Department of Statistics Faculty of Sciences Golestan
More informationWindow Width Selection for L 2 Adjusted Quantile Regression
Window Width Selection for L 2 Adjusted Quantile Regression Yoonsuh Jung, The Ohio State University Steven N. MacEachern, The Ohio State University Yoonkyung Lee, The Ohio State University Technical Report
More informationMATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS
MATH307/37 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS School of Mathematics and Statistics Semester, 04 Tutorial problems should be used to test your mathematical skills and understanding of the lecture material.
More informationPASS Sample Size Software
Chapter 850 Introduction Cox proportional hazards regression models the relationship between the hazard function λ( t X ) time and k covariates using the following formula λ log λ ( t X ) ( t) 0 = β1 X1
More informationHomework Assignments for BusAdm 713: Business Forecasting Methods. Assignment 1: Introduction to forecasting, Review of regression
Homework Assignments for BusAdm 713: Business Forecasting Methods Note: Problem points are in parentheses. Assignment 1: Introduction to forecasting, Review of regression 1. (3) Complete the exercises
More informationA RIDGE REGRESSION ESTIMATION APPROACH WHEN MULTICOLLINEARITY IS PRESENT
Fundamental Journal of Applied Sciences Vol. 1, Issue 1, 016, Pages 19-3 This paper is available online at http://www.frdint.com/ Published online February 18, 016 A RIDGE REGRESSION ESTIMATION APPROACH
More informationIntroduction to Computational Finance and Financial Econometrics Descriptive Statistics
You can t see this text! Introduction to Computational Finance and Financial Econometrics Descriptive Statistics Eric Zivot Summer 2015 Eric Zivot (Copyright 2015) Descriptive Statistics 1 / 28 Outline
More informationCase Study: Heavy-Tailed Distribution and Reinsurance Rate-making
Case Study: Heavy-Tailed Distribution and Reinsurance Rate-making May 30, 2016 The purpose of this case study is to give a brief introduction to a heavy-tailed distribution and its distinct behaviors in
More informationLinear regression model
Regression Model Assumptions (Solutions) STAT-UB.0003: Regression and Forecasting Models Linear regression model 1. Here is the least squares regression fit to the Zagat restaurant data: 10 15 20 25 10
More informationModelling volatility - ARCH and GARCH models
Modelling volatility - ARCH and GARCH models Beáta Stehlíková Time series analysis Modelling volatility- ARCH and GARCH models p.1/33 Stock prices Weekly stock prices (library quantmod) Continuous returns:
More informationLESSON 7 INTERVAL ESTIMATION SAMIE L.S. LY
LESSON 7 INTERVAL ESTIMATION SAMIE L.S. LY 1 THIS WEEK S PLAN Part I: Theory + Practice ( Interval Estimation ) Part II: Theory + Practice ( Interval Estimation ) z-based Confidence Intervals for a Population
More informationMeasuring Financial Risk using Extreme Value Theory: evidence from Pakistan
Measuring Financial Risk using Extreme Value Theory: evidence from Pakistan Dr. Abdul Qayyum and Faisal Nawaz Abstract The purpose of the paper is to show some methods of extreme value theory through analysis
More informationAP Stats: 3B ~ Least Squares Regression and Residuals. Objectives:
Objectives: INTERPRET the slope and y intercept of a least-squares regression line USE the least-squares regression line to predict y for a given x CALCULATE and INTERPRET residuals and their standard
More informationF UNCTIONAL R ELATIONSHIPS BETWEEN S TOCK P RICES AND CDS S PREADS
F UNCTIONAL R ELATIONSHIPS BETWEEN S TOCK P RICES AND CDS S PREADS Amelie Hüttner XAIA Investment GmbH Sonnenstraße 19, 80331 München, Germany amelie.huettner@xaia.com March 19, 014 Abstract We aim to
More informationThe 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 informationEfficient Parameter Estimation for ODE Models from
MASAMB 2016 Efficient Parameter Estimation for ODE Models from Relative Data Using Hierarchical Optimization Sabrina Krause, Carolin Loos, Jan Hasenauer Helmholtz Zentrum München Institute of Computational
More informationEE/AA 578 Univ. of Washington, Fall Homework 8
EE/AA 578 Univ. of Washington, Fall 2016 Homework 8 1. Multi-label SVM. The basic Support Vector Machine (SVM) described in the lecture (and textbook) is used for classification of data with two labels.
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 informationQuestion 1 Consider an economy populated by a continuum of measure one of consumers whose preferences are defined by the utility function:
Question 1 Consider an economy populated by a continuum of measure one of consumers whose preferences are defined by the utility function: β t log(c t ), where C t is consumption and the parameter β satisfies
More informationA Monte Carlo Based Analysis of Optimal Design Criteria
A Monte Carlo Based Analysis of Optimal Design Criteria H. T. Banks, Kathleen J. Holm and Franz Kappel Center for Quantitative Sciences in Biomedicine Center for Research in Scientific Computation North
More informationChapter 8: CAPM. 1. Single Index Model. 2. Adding a Riskless Asset. 3. The Capital Market Line 4. CAPM. 5. The One-Fund Theorem
Chapter 8: CAPM 1. Single Index Model 2. Adding a Riskless Asset 3. The Capital Market Line 4. CAPM 5. The One-Fund Theorem 6. The Characteristic Line 7. The Pricing Model Single Index Model 1 1. Covariance
More informationJacob: What data do we use? Do we compile paid loss triangles for a line of business?
PROJECT TEMPLATES FOR REGRESSION ANALYSIS APPLIED TO LOSS RESERVING BACKGROUND ON PAID LOSS TRIANGLES (The attached PDF file has better formatting.) {The paid loss triangle helps you! distinguish between
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2010, Mr. Ruey S. Tsay Solutions to Final Exam
The University of Chicago, Booth School of Business Business 410, Spring Quarter 010, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (4 pts) Answer briefly the following questions. 1. Questions 1
More informationLecture 6: Non Normal Distributions
Lecture 6: Non Normal Distributions and their Uses in GARCH Modelling Prof. Massimo Guidolin 20192 Financial Econometrics Spring 2015 Overview Non-normalities in (standardized) residuals from asset return
More informationStock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models
Stock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models David Prager 1 1 Associate Professor of Mathematics Anderson University (SC) Based on joint work with Professor Qing Zhang,
More informationRisk Reduction Potential
Risk Reduction Potential Research Paper 006 February, 015 015 Northstar Risk Corp. All rights reserved. info@northstarrisk.com Risk Reduction Potential In this paper we introduce the concept of risk reduction
More informationThis is a open-book exam. Assigned: Friday November 27th 2009 at 16:00. Due: Monday November 30th 2009 before 10:00.
University of Iceland School of Engineering and Sciences Department of Industrial Engineering, Mechanical Engineering and Computer Science IÐN106F Industrial Statistics II - Bayesian Data Analysis Fall
More informationFinite Element Method
In Finite Difference Methods: the solution domain is divided into a grid of discrete points or nodes the PDE is then written for each node and its derivatives replaced by finite-divided differences In
More informationMatlab Based Stochastic Processes in Stochastic Asset Portfolio Optimization
Matlab Based Stochastic Processes in Stochastic Asset Portfolio Optimization Workshop OR und Statistische Analyse mit Mathematischen Tools, May 13th, 2014 Dr. Georg Ostermaier, Ömer Kuzugüden Agenda Introduction
More informationMatlab Workshop MFE 2006 Lecture 4
Matlab Workshop MFE 2006 Lecture 4 Stefano Corradin Peng Liu http://faculty.haas.berkeley.edu/peliu/computing Haas School of Business, Berkeley, MFE 2006 Applications in Finance II 4.1 Optimum toolbox.
More informationLab 2 - Decision theory
Lab 2 - Decision theory Edvin Listo Zec 920625-2976 edvinli@student.chalmers.se September 29, 2014 Co-worker: Jessica Fredby Introduction The goal of this computer assignment is to analyse a given set
More informationBivariate Birnbaum-Saunders Distribution
Department of Mathematics & Statistics Indian Institute of Technology Kanpur January 2nd. 2013 Outline 1 Collaborators 2 3 Birnbaum-Saunders Distribution: Introduction & Properties 4 5 Outline 1 Collaborators
More informationHomework Assignment Section 3
Homework Assignment Section 3 Tengyuan Liang Business Statistics Booth School of Business Problem 1 A company sets different prices for a particular stereo system in eight different regions of the country.
More informationThe 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 informationApproximating a life table by linear combinations of exponential distributions and valuing life-contingent options
Approximating a life table by linear combinations of exponential distributions and valuing life-contingent options Zhenhao Zhou Department of Statistics and Actuarial Science The University of Iowa Iowa
More informationDynamic Replication of Non-Maturing Assets and Liabilities
Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland
More informationTests for the Difference Between Two Linear Regression Intercepts
Chapter 853 Tests for the Difference Between Two Linear Regression Intercepts Introduction Linear regression is a commonly used procedure in statistical analysis. One of the main objectives in linear regression
More informationarxiv: v1 [math.st] 18 Sep 2018
Gram Charlier and Edgeworth expansion for sample variance arxiv:809.06668v [math.st] 8 Sep 08 Eric Benhamou,* A.I. SQUARE CONNECT, 35 Boulevard d Inkermann 900 Neuilly sur Seine, France and LAMSADE, Universit
More informationSYSM 6304 Risk and Decision Analysis Lecture 2: Fitting Distributions to Data
SYSM 6304 Risk and Decision Analysis Lecture 2: Fitting Distributions to Data M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu September 5, 2015
More informationSupplementary Appendix to The Risk Premia Embedded in Index Options
Supplementary Appendix to The Risk Premia Embedded in Index Options Torben G. Andersen Nicola Fusari Viktor Todorov December 214 Contents A The Non-Linear Factor Structure of Option Surfaces 2 B Additional
More informationFinancial Risk Management
Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #4 1 Correlation and copulas 1. The bivariate Gaussian copula is given
More informationExtend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty
Extend the ideas of Kan and Zhou paper on Optimal Portfolio Construction under parameter uncertainty George Photiou Lincoln College University of Oxford A dissertation submitted in partial fulfilment for
More informationForecasting Chapter 14
Forecasting Chapter 14 14-01 Forecasting Forecast: A prediction of future events used for planning purposes. It is a critical inputs to business plans, annual plans, and budgets Finance, human resources,
More informationWHERE HAS ALL THE BIG DATA GONE?
WHERE HAS ALL THE BIG DATA GONE? Maryam Farboodi Princeton Adrien Matray Princeton Laura Veldkamp NYU Stern School of Business 2018 MOTIVATION Increase in big data in financial sector 1. data processing
More informationTangent 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 informationSupplementary Material for Combinatorial Partial Monitoring Game with Linear Feedback and Its Application. A. Full proof for Theorems 4.1 and 4.
Supplementary Material for Combinatorial Partial Monitoring Game with Linear Feedback and Its Application. A. Full proof for Theorems 4.1 and 4. If the reader will recall, we have the following problem-specific
More informationApproximation of functions and American options
Approximation of functions and American options Responsible teacher: Anatoliy Malyarenko December 8, 2003 Abstract Contents of the lecture: Approximation of functions. Gorner s scheme. American options.
More informationAn analysis of momentum and contrarian strategies using an optimal orthogonal portfolio approach
An analysis of momentum and contrarian strategies using an optimal orthogonal portfolio approach Hossein Asgharian and Björn Hansson Department of Economics, Lund University Box 7082 S-22007 Lund, Sweden
More informationRegression and Simulation
Regression and Simulation This is an introductory R session, so it may go slowly if you have never used R before. Do not be discouraged. A great way to learn a new language like this is to plunge right
More informationFinancial Econometrics (FinMetrics04) Time-series Statistics Concepts Exploratory Data Analysis Testing for Normality Empirical VaR
Financial Econometrics (FinMetrics04) Time-series Statistics Concepts Exploratory Data Analysis Testing for Normality Empirical VaR Nelson Mark University of Notre Dame Fall 2017 September 11, 2017 Introduction
More informationGARCH Models. Instructor: G. William Schwert
APS 425 Fall 2015 GARCH Models Instructor: G. William Schwert 585-275-2470 schwert@schwert.ssb.rochester.edu Autocorrelated Heteroskedasticity Suppose you have regression residuals Mean = 0, not autocorrelated
More informationNumerical Simulation of Stochastic Differential Equations: Lecture 1, Part 1. Overview of Lecture 1, Part 1: Background Mater.
Numerical Simulation of Stochastic Differential Equations: Lecture, Part Des Higham Department of Mathematics University of Strathclyde Course Aim: Give an accessible intro. to SDEs and their numerical
More informationCE 513: STATISTICAL METHODS
/CE 608 CE 513: STATISTICAL METHODS IN CIVIL ENGINEERING Lecture-1: Introduction & Overview Dr. Budhaditya Hazra Room: N-307 Department of Civil Engineering 1 Schedule of Lectures Last class before puja
More informationTHE CHINESE UNIVERSITY OF HONG KONG Department of Mathematics MMAT5250 Financial Mathematics Homework 2 Due Date: March 24, 2018
THE CHINESE UNIVERSITY OF HONG KONG Department of Mathematics MMAT5250 Financial Mathematics Homework 2 Due Date: March 24, 2018 Name: Student ID.: I declare that the assignment here submitted is original
More informationVaR Estimation under Stochastic Volatility Models
VaR Estimation under Stochastic Volatility Models Chuan-Hsiang Han Dept. of Quantitative Finance Natl. Tsing-Hua University TMS Meeting, Chia-Yi (Joint work with Wei-Han Liu) December 5, 2009 Outline Risk
More informationA potentially useful approach to model nonlinearities in time series is to assume different behavior (structural break) in different subsamples
1.3 Regime switching models A potentially useful approach to model nonlinearities in time series is to assume different behavior (structural break) in different subsamples (or regimes). If the dates, the
More informationRegression Review and Robust Regression. Slides prepared by Elizabeth Newton (MIT)
Regression Review and Robust Regression Slides prepared by Elizabeth Newton (MIT) S-Plus Oil City Data Frame Monthly Excess Returns of Oil City Petroleum, Inc. Stocks and the Market SUMMARY: The oilcity
More informationMean-Variance Analysis
Mean-Variance Analysis Mean-variance analysis 1/ 51 Introduction How does one optimally choose among multiple risky assets? Due to diversi cation, which depends on assets return covariances, the attractiveness
More informationV(0.1) V( 0.5) 0.6 V(0.5) V( 0.5)
In-class exams are closed book, no calculators, except for one 8.5"x11" page, written in any density (student may bring a magnifier). Students are bound by the University of Florida honor code. Exam papers
More informationFinal Exam Suggested Solutions
University of Washington Fall 003 Department of Economics Eric Zivot Economics 483 Final Exam Suggested Solutions This is a closed book and closed note exam. However, you are allowed one page of handwritten
More informationStat 328, Summer 2005
Stat 328, Summer 2005 Exam #2, 6/18/05 Name (print) UnivID I have neither given nor received any unauthorized aid in completing this exam. Signed Answer each question completely showing your work where
More informationBusiness Statistics: A First Course
Business Statistics: A First Course Fifth Edition Chapter 12 Correlation and Simple Linear Regression Business Statistics: A First Course, 5e 2009 Prentice-Hall, Inc. Chap 12-1 Learning Objectives In this
More informationDescription Quick start Menu Syntax Options Remarks and examples Stored results Methods and formulas Acknowledgment References Also see
Title stata.com tssmooth shwinters Holt Winters seasonal smoothing Description Quick start Menu Syntax Options Remarks and examples Stored results Methods and formulas Acknowledgment References Also see
More informationu (x) < 0. and if you believe in diminishing return of the wealth, then you would require
Chapter 8 Markowitz Portfolio Theory 8.7 Investor Utility Functions People are always asked the question: would more money make you happier? The answer is usually yes. The next question is how much more
More informationStat 101 Exam 1 - Embers Important Formulas and Concepts 1
1 Chapter 1 1.1 Definitions Stat 101 Exam 1 - Embers Important Formulas and Concepts 1 1. Data Any collection of numbers, characters, images, or other items that provide information about something. 2.
More information