Conjugate Bayesian Models for Massive Spatial Data
|
|
- Pearl Baker
- 5 years ago
- Views:
Transcription
1 Conjugate Bayesian Models for Massive Spatial Data Abhi Datta 1, Sudipto Banerjee 2 and Andrew O. Finley 3 July 31, Department of Biostatistics, Bloomberg School of Public Health, Johns Hopkins University, Baltimore, Maryland. 2 Department of Biostatistics, Fielding School of Public Health, University of California, Los Angeles. 3 Departments of Forestry and Geography, Michigan State University, East Lansing, Michigan.
2 Case Study: Alaska Tanana Valley Forest Height Dataset Forest height and tree cover Forest fire history Forest height (red lines) data from LiDAR at locations Knowledge of forest height is important for biomass assessment, carbon management etc 1
3 Case Study: Alaska Tanana Valley Forest Height Dataset Forest height and tree cover Forest fire history Goal: High-resolution domainwide prediction maps of forest height Covariates: Domainwide tree cover (grey) and forest fire history (red patches) in the last 20 years 1
4 Analyzing the data Models used: Non-spatial regression: y FH (s) = β 0 + β tree x tree + β fire x fire + ɛ(s) Figure: Variogram of the residuals from non-spatial regression indicates strong spatial pattern 2
5 NNGP models Collapsed NNGP: y FH (s) = β 0 + β tree x tree + β fire x fire + w(s) + ɛ(s) w(s) NNGP(0, C(, σ 2, φ)) y FH N(Xβ, C + τ 2 I) where C is the NNGP covariance matrix derived from C Response NNGP: y FH (s) NNGP(β 0 + β tree x tree + β fire x fire, Σ(, σ 2, φ, τ 2 )) y FH N(Xβ, Σ) where Σ is the NNGP covariance matrix derived from Σ = C + τ 2 I 3
6 NNGP models Non-spatial regression Collapsed NNGP Response NNGP CRPS RMSPE CP 93% 94% 94% CIW Table: Model comparison metrics for the Tanana valley dataset NNGP models perform significantly better than the non-spatial model MCMC run time for the NNGP models: Collapsed model: 319 hours Response model: 38 hours For massive spatial data, full Bayesian output for even NNGP models require substantial time 4
7 Another look at the response model Original full GP model: y(s) ind N(x(s) β + w(s), τ 2 ) w(s) GP with a stationary covariance function C(, σ 2, φ) Cov(w) = σ 2 R(φ) Full GP model: y N(Xβ, Σ) where Σ = σ 2 M M = R(φ) + αi α = τ 2 /σ 2 is the ratio of the noise to signal variance Response NNGP model: y N(Xβ, Σ) Σ = σ 2 M where M is the NNGP approximation for M 5
8 Conjugate NNGP y N(Xβ, σ 2 M) If φ and α are known, M, and hence M, are known matrices The model becomes a standard Bayesian linear model Assume a Normal Inverse Gamma (NIG) prior for (β, σ 2 ) (β, σ 2 ) NIG(µ β, V β, a σ, b σ ), i.e., β σ 2 N(µ β, σ 2 V β ) and σ 2 IG(a σ, b σ ) 6
9 Conjugate NNGP y N(Xβ, σ 2 M), M is known Joint likelihood: N(y Xβ, σ 2 M) N(β µ β, σ 2 V β ) IG(σ 2 a σ, b σ ) 7
10 Conjugate NNGP y N(Xβ, σ 2 M), M is known Joint likelihood: N(y Xβ, σ 2 M) N(β µ β, σ 2 V β ) IG(σ 2 a σ, b σ ) Conjugate posterior distribution (β, σ 2 ) y NIG(µ β, V β, a σ, b σ) Expressions for µ β, V β, a σ and b σ can be calculated in O(n) time 7
11 Conjugate NNGP (β, σ 2 ) y NIG(µ β, V β, a σ, b σ) Marginal posterior: β y MVt 2a σ (µ β, b σ a σ V β ) MVt k (m, V ) is the multivariate t distribution with degrees of k, mean m and scale matrix V E(β y) = µ β, Var(β y) = b σ aσ 1V β Marginal posterior: σ 2 y IG(aσ, bσ) E(σ 2 y) = b σ aσ 1, Var(σ2 y) = b 2 σ (a σ 1)2 (a σ 2) Exact posterior distributions of β and σ 2 are available 8
12 Predictive distributions y(s) y t 2a σ (m(s), b σ a σ v(s)) E(y(s) y) = m(s), Var(y(s) y) = b σ a σ 1v(s) m(s) and v(s) can be computed using O(m) flops Exact posterior predictive distributions of y(s) y for any s No MCMC required for parameter estimation or prediction 9
13 Choosing α and φ φ and α are chosen using K-fold cross validation over a grid of possible values Unlike MCMC, cross-validation can be completely parallelized Resolution of the grid for φ and α can be decided based on computing resources available In practice, a reasonably coarse grid often suffices 10
14 Choosing α and φ α φ RMSPE RMSPE Figure: Simulation experiment: True value (+) of (α, φ) and estimated value ( ) using 5-fold cross validation 11
15 Scalability Computation and storage requirements are O(n) One evaluation time similar to the response NNGP model Unlike response NNGP, does not involve any serial MCMC iterations For K fold cross validation and G combinations of φ and α, total number of evaluations is KG Embarassingly parallel: Each of the KG evaluations can proceed in parallel 12
16 Scalability Figure: Run times of different NNGP models with increasing sample size 13
17 Alaska Tanana Valley dataset Conjugate NNGP Collapsed NNGP Response NNGP β (2.35, 2.47) 2.37 (2.31,2.42) β TC (0.02, 0.02) 0.02 (0.02, 0.02) β Fire (0.34, 0.43) 0.43 (0.39, 0.48) σ (18.50, 18.81) (17.13, 17.41) τ (1.55, 1.56) 1.55 (1.54, 1.55) φ (3.70, 3.77) 4.15 (4.13, 4.19) CRPS RMSPE time (hrs.) Table: Parameter estimates and model comparison metrics for the Tanana valley dataset Conjugate model produces estimates and model comparison numbers very similar to the MCMC based NNGP models For locations, conjugate model takes 7 seconds 14
18 Summary MCMC free exact Bayesian approach by fixing some covariance parameters Conjugate posterior distributions of the parameters and posterior predictive distributions available in closed form Embarassingly parallel cross validation to identify best choices for fixed parameters Runs in seconds for massive spatial dataset with millions of locations Available in the spnngp package in R 15
Bayesian Linear Model: Gory Details
Bayesian Linear Model: Gory Details Pubh7440 Notes By Sudipto Banerjee Let y y i ] n i be an n vector of independent observations on a dependent variable (or response) from n experimental units. Associated
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 informationOutline. Review Continuation of exercises from last time
Bayesian Models II Outline Review Continuation of exercises from last time 2 Review of terms from last time Probability density function aka pdf or density Likelihood function aka likelihood Conditional
More informationST440/550: Applied Bayesian Analysis. (5) Multi-parameter models - Summarizing the posterior
(5) Multi-parameter models - Summarizing the posterior Models with more than one parameter Thus far we have studied single-parameter models, but most analyses have several parameters For example, consider
More informationBayesian Normal Stuff
Bayesian Normal Stuff - Set-up of the basic model of a normally distributed random variable with unknown mean and variance (a two-parameter model). - Discuss philosophies of prior selection - Implementation
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 informationConjugate Models. Patrick Lam
Conjugate Models Patrick Lam Outline Conjugate Models What is Conjugacy? The Beta-Binomial Model The Normal Model Normal Model with Unknown Mean, Known Variance Normal Model with Known Mean, Unknown Variance
More informationNon-informative Priors Multiparameter Models
Non-informative Priors Multiparameter Models Statistics 220 Spring 2005 Copyright c 2005 by Mark E. Irwin Prior Types Informative vs Non-informative There has been a desire for a prior distributions that
More informationMachine 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 informationBeating the market, using linear regression to outperform the market average
Radboud University Bachelor Thesis Artificial Intelligence department Beating the market, using linear regression to outperform the market average Author: Jelle Verstegen Supervisors: Marcel van Gerven
More informationWeight Smoothing with Laplace Prior and Its Application in GLM Model
Weight Smoothing with Laplace Prior and Its Application in GLM Model Xi Xia 1 Michael Elliott 1,2 1 Department of Biostatistics, 2 Survey Methodology Program, University of Michigan National Cancer Institute
More informationCS340 Machine learning Bayesian statistics 3
CS340 Machine learning Bayesian statistics 3 1 Outline Conjugate analysis of µ and σ 2 Bayesian model selection Summarizing the posterior 2 Unknown mean and precision The likelihood function is p(d µ,λ)
More information1 Bayesian Bias Correction Model
1 Bayesian Bias Correction Model Assuming that n iid samples {X 1,...,X n }, were collected from a normal population with mean µ and variance σ 2. The model likelihood has the form, P( X µ, σ 2, T n >
More informationPosterior Inference. , where should we start? Consider the following computational procedure: 1. draw samples. 2. convert. 3. compute properties
Posterior Inference Example. Consider a binomial model where we have a posterior distribution for the probability term, θ. Suppose we want to make inferences about the log-odds γ = log ( θ 1 θ), where
More informationSTAT 425: Introduction to Bayesian Analysis
STAT 45: Introduction to Bayesian Analysis Marina Vannucci Rice University, USA Fall 018 Marina Vannucci (Rice University, USA) Bayesian Analysis (Part 1) Fall 018 1 / 37 Lectures 9-11: Multi-parameter
More informationMissing Data. EM Algorithm and Multiple Imputation. Aaron Molstad, Dootika Vats, Li Zhong. University of Minnesota School of Statistics
Missing Data EM Algorithm and Multiple Imputation Aaron Molstad, Dootika Vats, Li Zhong University of Minnesota School of Statistics December 4, 2013 Overview 1 EM Algorithm 2 Multiple Imputation Incomplete
More informationModel 0: We start with a linear regression model: log Y t = β 0 + β 1 (t 1980) + ε, with ε N(0,
Stat 534: Fall 2017. Introduction to the BUGS language and rjags Installation: download and install JAGS. You will find the executables on Sourceforge. You must have JAGS installed prior to installing
More informationChapter 7: Estimation Sections
1 / 31 : Estimation Sections 7.1 Statistical Inference Bayesian Methods: 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions 7.4 Bayes Estimators Frequentist Methods: 7.5 Maximum Likelihood
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 informationExtended Model: Posterior Distributions
APPENDIX A Extended Model: Posterior Distributions A. Homoskedastic errors Consider the basic contingent claim model b extended by the vector of observables x : log C i = β log b σ, x i + β x i + i, i
More informationσ e, which will be large when prediction errors are Linear regression model
Linear regression model we assume that two quantitative variables, x and y, are linearly related; that is, the population of (x, y) pairs are related by an ideal population regression line y = α + βx +
More informationEstimation 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 informationKernel Conditional Quantile Estimation via Reduction Revisited
Kernel Conditional Quantile Estimation via Reduction Revisited Novi Quadrianto Novi.Quad@gmail.com The Australian National University, Australia NICTA, Statistical Machine Learning Program, Australia Joint
More informationModelling Returns: the CER and the CAPM
Modelling Returns: the CER and the CAPM Carlo Favero Favero () Modelling Returns: the CER and the CAPM 1 / 20 Econometric Modelling of Financial Returns Financial data are mostly observational data: they
More information1. You are given the following information about a stationary AR(2) model:
Fall 2003 Society of Actuaries **BEGINNING OF EXAMINATION** 1. You are given the following information about a stationary AR(2) model: (i) ρ 1 = 05. (ii) ρ 2 = 01. Determine φ 2. (A) 0.2 (B) 0.1 (C) 0.4
More informationRobust Regression for Capital Asset Pricing Model Using Bayesian Approach
Thai Journal of Mathematics : 016) 71 8 Special Issue on Applied Mathematics : Bayesian Econometrics http://thaijmath.in.cmu.ac.th ISSN 1686-009 Robust Regression for Capital Asset Pricing Model Using
More informationGeostatistical Inference under Preferential Sampling
Geostatistical Inference under Preferential Sampling Marie Ozanne and Justin Strait Diggle, Menezes, and Su, 2010 October 12, 2015 Marie Ozanne and Justin Strait Preferential Sampling October 12, 2015
More informationChapter 7: Estimation Sections
1 / 40 Chapter 7: Estimation Sections 7.1 Statistical Inference Bayesian Methods: Chapter 7 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions 7.4 Bayes Estimators Frequentist Methods:
More informationAnalysis of extreme values with random location Abstract Keywords: 1. Introduction and Model
Analysis of extreme values with random location Ali Reza Fotouhi Department of Mathematics and Statistics University of the Fraser Valley Abbotsford, BC, Canada, V2S 7M8 Ali.fotouhi@ufv.ca Abstract Analysis
More informationChapter 8: Sampling distributions of estimators Sections
Chapter 8 continued Chapter 8: Sampling distributions of estimators Sections 8.1 Sampling distribution of a statistic 8.2 The Chi-square distributions 8.3 Joint Distribution of the sample mean and sample
More informationSession 8: The Markowitz problem p. 1
Session 8: The Markowitz problem Susan Thomas http://www.igidr.ac.in/ susant susant@mayin.org IGIDR Bombay Session 8: The Markowitz problem p. 1 Portfolio optimisation Session 8: The Markowitz problem
More informationSTART HERE: Instructions. 1 Exponential Family [Zhou, Manzil]
START HERE: Instructions Thanks a lot to John A.W.B. Constanzo and Shi Zong for providing and allowing to use the latex source files for quick preparation of the HW solution. The homework was due at 9:00am
More informationThe Time-Varying Effects of Monetary Aggregates on Inflation and Unemployment
経営情報学論集第 23 号 2017.3 The Time-Varying Effects of Monetary Aggregates on Inflation and Unemployment An Application of the Bayesian Vector Autoregression with Time-Varying Parameters and Stochastic Volatility
More informationCommon Misconceptions about "Beta" Hedging, Estimation and Horizon Effects 1
QuantNugget3 Common Misconceptions about "Beta" Hedging, Estimation and Horizon Effects 1 Attilio Meucci 2 attilio_meucci@symmys.com this version: eptember 27 2010 last version available at: http://ssrn.com/abstract=1619923
More informationLog-linear Modeling Under Generalized Inverse Sampling Scheme
Log-linear Modeling Under Generalized Inverse Sampling Scheme Soumi Lahiri (1) and Sunil Dhar (2) (1) Department of Mathematical Sciences New Jersey Institute of Technology University Heights, Newark,
More informationBayesian course - problem set 3 (lecture 4)
Bayesian course - problem set 3 (lecture 4) Ben Lambert November 14, 2016 1 Ticked off Imagine once again that you are investigating the occurrence of Lyme disease in the UK. This is a vector-borne disease
More informationThe Monte Carlo Method in High Performance Computing
The Monte Carlo Method in High Performance Computing Dieter W. Heermann Monte Carlo Methods 2015 Dieter W. Heermann (Monte Carlo Methods)The Monte Carlo Method in High Performance Computing 2015 1 / 1
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 informationTRΛNSPΛRΣNCY ΛNΛLYTICS
TRΛNSPΛRΣNCY ΛNΛLYTICS RISK-AI, LLC PRESENTATION INTRODUCTION I. Transparency Analytics is a state-of-the-art risk management analysis and research platform for Investment Advisors, Funds of Funds, Family
More informationTechnical Appendix: Policy Uncertainty and Aggregate Fluctuations.
Technical Appendix: Policy Uncertainty and Aggregate Fluctuations. Haroon Mumtaz Paolo Surico July 18, 2017 1 The Gibbs sampling algorithm Prior Distributions and starting values Consider the model to
More informationINSTITUTE AND FACULTY OF ACTUARIES. Curriculum 2019 SPECIMEN EXAMINATION
INSTITUTE AND FACULTY OF ACTUARIES Curriculum 2019 SPECIMEN EXAMINATION Subject CS1A Actuarial Statistics Time allowed: Three hours and fifteen minutes INSTRUCTIONS TO THE CANDIDATE 1. Enter all the candidate
More informationTop-down particle filtering for Bayesian decision trees
Top-down particle filtering for Bayesian decision trees Balaji Lakshminarayanan 1, Daniel M. Roy 2 and Yee Whye Teh 3 1. Gatsby Unit, UCL, 2. University of Cambridge and 3. University of Oxford Outline
More informationApplication of MCMC Algorithm in Interest Rate Modeling
Application of MCMC Algorithm in Interest Rate Modeling Xiaoxia Feng and Dejun Xie Abstract Interest rate modeling is a challenging but important problem in financial econometrics. This work is concerned
More informationBayesian Forecasting & Scalable Multivariate Volatility Analysis Using Simultaneous Graphical Dynamic Models
Bayesian Forecasting & Scalable Multivariate Volatility Analysis Using Simultaneous Graphical Dynamic Models Lutz F. Gruber 1, Mike West 2 Duke University Abstract The recently introduced class of simultaneous
More informationOptimal weights for the MSCI North America index. Optimal weights for the MSCI Europe index
Portfolio construction with Bayesian GARCH forecasts Wolfgang Polasek and Momtchil Pojarliev Institute of Statistics and Econometrics University of Basel Holbeinstrasse 12 CH-4051 Basel email: Momtchil.Pojarliev@unibas.ch
More informationEstimating and Accounting for the Output Gap with Large Bayesian Vector Autoregressions
Estimating and Accounting for the Output Gap with Large Bayesian Vector Autoregressions James Morley 1 Benjamin Wong 2 1 University of Sydney 2 Reserve Bank of New Zealand The view do not necessarily represent
More informationarxiv: v1 [stat.ap] 27 Jun 2016 Abstract
Bayesian Forecasting & Scalable Multivariate Volatility Analysis Using Simultaneous Graphical Dynamic Models Lutz F. Gruber 1, Mike West 2 Duke University arxiv:1606.08291v1 [stat.ap] 27 Jun 2016 Abstract
More informationStochastic Volatility Models. Hedibert Freitas Lopes
Stochastic Volatility Models Hedibert Freitas Lopes SV-AR(1) model Nonlinear dynamic model Normal approximation R package stochvol Other SV models STAR-SVAR(1) model MSSV-SVAR(1) model Volume-volatility
More informationHigh 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 informationAppendix for Beazer, Quintin H. & Byungwon Woo IMF Conditionality, Government Partisanship, and the Progress of Economic Reforms
Appendix for Beazer, Quintin H. & Byungwon Woo. 2015. IMF Conditionality, Government Partisanship, and the Progress of Economic Reforms This appendix contains the additional analyses that space considerations
More informationTwo hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER
Two hours MATH20802 To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER STATISTICAL METHODS Answer any FOUR of the SIX questions.
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 informationSubject CS1 Actuarial Statistics 1 Core Principles. Syllabus. for the 2019 exams. 1 June 2018
` Subject CS1 Actuarial Statistics 1 Core Principles Syllabus for the 2019 exams 1 June 2018 Copyright in this Core Reading is the property of the Institute and Faculty of Actuaries who are the sole distributors.
More informationBayesian Multinomial Model for Ordinal Data
Bayesian Multinomial Model for Ordinal Data Overview This example illustrates how to fit a Bayesian multinomial model by using the built-in mutinomial density function (MULTINOM) in the MCMC procedure
More informationDuangporn Jearkpaporn, Connie M. Borror Douglas C. Montgomery and George C. Runger Arizona State University Tempe, AZ
Process Monitoring for Correlated Gamma Distributed Data Using Generalized Linear Model Based Control Charts Duangporn Jearkpaporn, Connie M. Borror Douglas C. Montgomery and George C. Runger Arizona State
More informationCommon one-parameter models
Common one-parameter models In this section we will explore common one-parameter models, including: 1. Binomial data with beta prior on the probability 2. Poisson data with gamma prior on the rate 3. Gaussian
More informationHierarchical Generalized Linear Models. Measurement Incorporated Hierarchical Linear Models Workshop
Hierarchical Generalized Linear Models Measurement Incorporated Hierarchical Linear Models Workshop Hierarchical Generalized Linear Models So now we are moving on to the more advanced type topics. To begin
More information(11) Case Studies: Adaptive clinical trials. ST440/540: Applied Bayesian Analysis
Use of Bayesian methods in clinical trials Bayesian methods are becoming more common in clinical trials analysis We will study how to compute the sample size for a Bayesian clinical trial We will then
More informationMoral Hazard: Dynamic Models. Preliminary Lecture Notes
Moral Hazard: Dynamic Models Preliminary Lecture Notes Hongbin Cai and Xi Weng Department of Applied Economics, Guanghua School of Management Peking University November 2014 Contents 1 Static Moral Hazard
More informationThis chapter introduces Markov chain Monte Carlo (MCMC) methods for empirical corporate
Damien (Typeset by SPi) 516 of 696 26 Markov chain Monte Carlo methods in corporate finance arthur korteweg This chapter introduces Markov chain Monte Carlo (MCMC) methods for empirical corporate finance.
More informationComputational Statistics Handbook with MATLAB
«H Computer Science and Data Analysis Series Computational Statistics Handbook with MATLAB Second Edition Wendy L. Martinez The Office of Naval Research Arlington, Virginia, U.S.A. Angel R. Martinez Naval
More informationcontinuous rv Note for a legitimate pdf, we have f (x) 0 and f (x)dx = 1. For a continuous rv, P(X = c) = c f (x)dx = 0, hence
continuous rv Let X be a continuous rv. Then a probability distribution or probability density function (pdf) of X is a function f(x) such that for any two numbers a and b with a b, P(a X b) = b a f (x)dx.
More informationModeling 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 informationInstitute of Actuaries of India Subject CT6 Statistical Methods
Institute of Actuaries of India Subject CT6 Statistical Methods For 2014 Examinations Aim The aim of the Statistical Methods subject is to provide a further grounding in mathematical and statistical techniques
More informationPractice Exercises for Midterm Exam ST Statistical Theory - II The ACTUAL exam will consists of less number of problems.
Practice Exercises for Midterm Exam ST 522 - Statistical Theory - II The ACTUAL exam will consists of less number of problems. 1. Suppose X i F ( ) for i = 1,..., n, where F ( ) is a strictly increasing
More informationObjective Bayesian Analysis for Heteroscedastic Regression
Analysis for Heteroscedastic Regression & Esther Salazar Universidade Federal do Rio de Janeiro Colóquio Inter-institucional: Modelos Estocásticos e Aplicações 2009 Collaborators: Marco Ferreira and Thais
More informationA Multivariate Analysis of Intercompany Loss Triangles
A Multivariate Analysis of Intercompany Loss Triangles Peng Shi School of Business University of Wisconsin-Madison ASTIN Colloquium May 21-24, 2013 Peng Shi (Wisconsin School of Business) Intercompany
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 informationSELECTION OF VARIABLES INFLUENCING IRAQI BANKS DEPOSITS BY USING NEW BAYESIAN LASSO QUANTILE REGRESSION
Vol. 6, No. 1, Summer 2017 2012 Published by JSES. SELECTION OF VARIABLES INFLUENCING IRAQI BANKS DEPOSITS BY USING NEW BAYESIAN Fadel Hamid Hadi ALHUSSEINI a Abstract The main focus of the paper is modelling
More informationUniversity of California, Los Angeles Department of Statistics. Final exam 07 June 2013
University of California, Los Angeles Department of Statistics Statistics C183/C283 Instructor: Nicolas Christou Final exam 07 June 2013 Name: Problem 1 (20 points) a. Suppose the variable X follows the
More informationTherefore, statistical modelling tools are required which make thorough space-time analyses of insurance regression data possible and allow to explore
Bayesian space time analysis of health insurance data Stefan Lang, Petra Kragler, Gerhard Haybach and Ludwig Fahrmeir University of Munich, Ludwigstr. 33, 80539 Munich email: lang@stat.uni-muenchen.de
More informationChapter 5 Univariate time-series analysis. () Chapter 5 Univariate time-series analysis 1 / 29
Chapter 5 Univariate time-series analysis () Chapter 5 Univariate time-series analysis 1 / 29 Time-Series Time-series is a sequence fx 1, x 2,..., x T g or fx t g, t = 1,..., T, where t is an index denoting
More informationSTAT758. Final Project. Time series analysis of daily exchange rate between the British Pound and the. US dollar (GBP/USD)
STAT758 Final Project Time series analysis of daily exchange rate between the British Pound and the US dollar (GBP/USD) Theophilus Djanie and Harry Dick Thompson UNR May 14, 2012 INTRODUCTION Time Series
More informationThe Usefulness of Bayesian Optimal Designs for Discrete Choice Experiments
The Usefulness of Bayesian Optimal Designs for Discrete Choice Experiments Roselinde Kessels Joint work with Bradley Jones, Peter Goos and Martina Vandebroek Outline 1. Motivating example from healthcare
More informationBlack-Litterman Model
Institute of Financial and Actuarial Mathematics at Vienna University of Technology Seminar paper Black-Litterman Model by: Tetyana Polovenko Supervisor: Associate Prof. Dipl.-Ing. Dr.techn. Stefan Gerhold
More informationThe Normal Distribution
The Normal Distribution The normal distribution plays a central role in probability theory and in statistics. It is often used as a model for the distribution of continuous random variables. Like all models,
More informationModelling the Sharpe ratio for investment strategies
Modelling the Sharpe ratio for investment strategies Group 6 Sako Arts 0776148 Rik Coenders 0777004 Stefan Luijten 0783116 Ivo van Heck 0775551 Rik Hagelaars 0789883 Stephan van Driel 0858182 Ellen Cardinaels
More information**BEGINNING OF EXAMINATION** A random sample of five observations from a population is:
**BEGINNING OF EXAMINATION** 1. You are given: (i) A random sample of five observations from a population is: 0.2 0.7 0.9 1.1 1.3 (ii) You use the Kolmogorov-Smirnov test for testing the null hypothesis,
More informationSemiparametric Modeling, Penalized Splines, and Mixed Models
Semi 1 Semiparametric Modeling, Penalized Splines, and Mixed Models David Ruppert Cornell University http://wwworiecornelledu/~davidr January 24 Joint work with Babette Brumback, Ray Carroll, Brent Coull,
More informationIdiosyncratic risk, insurance, and aggregate consumption dynamics: a likelihood perspective
Idiosyncratic risk, insurance, and aggregate consumption dynamics: a likelihood perspective Alisdair McKay Boston University June 2013 Microeconomic evidence on insurance - Consumption responds to idiosyncratic
More informationBayesian Dynamic Linear Models for Strategic Asset Allocation
Bayesian Dynamic Linear Models for Strategic Asset Allocation Jared Fisher Carlos Carvalho, The University of Texas Davide Pettenuzzo, Brandeis University April 18, 2016 Fisher (UT) Bayesian Risk Prediction
More informationAustralian Journal of Basic and Applied Sciences. Conditional Maximum Likelihood Estimation For Survival Function Using Cox Model
AENSI Journals Australian Journal of Basic and Applied Sciences Journal home page: wwwajbaswebcom Conditional Maximum Likelihood Estimation For Survival Function Using Cox Model Khawla Mustafa Sadiq University
More informationHeterogeneous 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 informationSemiparametric Modeling, Penalized Splines, and Mixed Models David Ruppert Cornell University
Semiparametric Modeling, Penalized Splines, and Mixed Models David Ruppert Cornell University Possible Model SBMD i,j is spinal bone mineral density on ith subject at age equal to age i,j lide http://wwworiecornelledu/~davidr
More informationStochastic Volatility (SV) Models Lecture 9. Morettin & Toloi, 2006, Section 14.6 Tsay, 2010, Section 3.12 Tsay, 2013, Section 4.
Stochastic Volatility (SV) Models Lecture 9 Morettin & Toloi, 2006, Section 14.6 Tsay, 2010, Section 3.12 Tsay, 2013, Section 4.13 Stochastic volatility model The canonical stochastic volatility model
More informationEquilibrium 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 informationFinal 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 informationA 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 informationConsistent estimators for multilevel generalised linear models using an iterated bootstrap
Multilevel Models Project Working Paper December, 98 Consistent estimators for multilevel generalised linear models using an iterated bootstrap by Harvey Goldstein hgoldstn@ioe.ac.uk Introduction Several
More informationVolatility Models and Their Applications
HANDBOOK OF Volatility Models and Their Applications Edited by Luc BAUWENS CHRISTIAN HAFNER SEBASTIEN LAURENT WILEY A John Wiley & Sons, Inc., Publication PREFACE CONTRIBUTORS XVII XIX [JQ VOLATILITY MODELS
More informationAnalysis of the Bitcoin Exchange Using Particle MCMC Methods
Analysis of the Bitcoin Exchange Using Particle MCMC Methods by Michael Johnson M.Sc., University of British Columbia, 2013 B.Sc., University of Winnipeg, 2011 Project Submitted in Partial Fulfillment
More informationChapter 7: Estimation Sections
Chapter 7: Estimation Sections 7.1 Statistical Inference Bayesian Methods: 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions Frequentist Methods: 7.5 Maximum Likelihood Estimators
More informationAnalyzing Oil Futures with a Dynamic Nelson-Siegel Model
Analyzing Oil Futures with a Dynamic Nelson-Siegel Model NIELS STRANGE HANSEN & ASGER LUNDE DEPARTMENT OF ECONOMICS AND BUSINESS, BUSINESS AND SOCIAL SCIENCES, AARHUS UNIVERSITY AND CENTER FOR RESEARCH
More informationOnline Appendix to ESTIMATING MUTUAL FUND SKILL: A NEW APPROACH. August 2016
Online Appendix to ESTIMATING MUTUAL FUND SKILL: A NEW APPROACH Angie Andrikogiannopoulou London School of Economics Filippos Papakonstantinou Imperial College London August 26 C. Hierarchical mixture
More informationBayesian Hierarchical/ Multilevel and Latent-Variable (Random-Effects) Modeling
Bayesian Hierarchical/ Multilevel and Latent-Variable (Random-Effects) Modeling 1: Formulation of Bayesian models and fitting them with MCMC in WinBUGS David Draper Department of Applied Mathematics and
More informationChapter 5 Macroeconomics and Finance
Macro II Chapter 5 Macro and Finance 1 Chapter 5 Macroeconomics and Finance Main references : - L. Ljundqvist and T. Sargent, Chapter 7 - Mehra and Prescott 1985 JME paper - Jerman 1998 JME paper - J.
More informationARCH and GARCH models
ARCH and GARCH models Fulvio Corsi SNS Pisa 5 Dic 2011 Fulvio Corsi ARCH and () GARCH models SNS Pisa 5 Dic 2011 1 / 21 Asset prices S&P 500 index from 1982 to 2009 1600 1400 1200 1000 800 600 400 200
More informationDiscrete & continuous characters: The threshold model. Liam J. Revell Anthrotree, 2014
Discrete & continuous characters: The threshold model Liam J. Revell Anthrotree, 2014 Discrete & continuous characters: the threshold model So far we have discussed continuous & discrete character models
More informationMCMC Estimation of Multiscale Stochastic Volatility Models
MCMC Estimation of Multiscale Stochastic Volatility Models German Molina, Chuan-Hsiang Han and Jean-Pierre Fouque Technical Report #23-6 June 3, 23 This material was based upon work supported by the National
More informationActuarial Society of India EXAMINATIONS
Actuarial Society of India EXAMINATIONS 7 th June 005 Subject CT6 Statistical Models Time allowed: Three Hours (0.30 am 3.30 pm) INSTRUCTIONS TO THE CANDIDATES. Do not write your name anywhere on the answer
More information