Estimation of a Ramsay-Curve IRT Model using the Metropolis-Hastings Robbins-Monro Algorithm
|
|
- Blaze Williams
- 5 years ago
- Views:
Transcription
1 1 / 34 Estimation of a Ramsay-Curve IRT Model using the Metropolis-Hastings Robbins-Monro Algorithm Scott Monroe & Li Cai IMPS 2012, Lincoln, Nebraska
2 Outline 2 / 34 1 Introduction and Motivation 2 Review of Ramsay Curves 3 Review of the MH-RM Algorithm 4 Simulation Study 5 Empirical Application 6 Conclusion
3 Outline 3 / 34 1 Introduction and Motivation 2 Review of Ramsay Curves 3 Review of the MH-RM Algorithm 4 Simulation Study 5 Empirical Application 6 Conclusion
4 Ramsay-Curve IRT 4 / 34 The need to estimate the latent trait distribution has received increasing attention in recent years. One approach is Ramsay Curve IRT, (RC-IRT, Woods & Thissen, 2006), which uses Bock and Aitkin s (1981) EM algorithm for estimation. However, the EM algorithm does not yield the observed information matrix upon convergence. Hence, no standard errors are available.
5 Ramsay-Curve IRT 5 / 34 The need to estimate the latent trait distribution has received increasing attention in recent years. One approach is Ramsay Curve IRT, (RC-IRT, Woods & Thissen, 2006), which uses Bock and Aitkin s (1981) EM algorithm for estimation. However, the EM algorithm does not yield the observed information matrix upon convergence. Hence, no standard errors are available. A solution: use the Metropolis-Hastings Robbins-Monro (MH-RM, Cai, 2010) algorithm.
6 Outline 6 / 34 1 Introduction and Motivation 2 Review of Ramsay Curves 3 Review of the MH-RM Algorithm 4 Simulation Study 5 Empirical Application 6 Conclusion
7 Review of Ramsay Curves For full details, see Woods and Thissen (2006). Given observed latent trait scores, RCs can estimate the distribution s shape. Essentially, the method uses B-spline regression, with a small number of parameters Densities Ramsay Curve Emp. Hist / 34
8 Outline 8 / 34 1 Introduction and Motivation 2 Review of Ramsay Curves 3 Review of the MH-RM Algorithm 4 Simulation Study 5 Empirical Application 6 Conclusion
9 9 / 34 A Review of MH-RM Metropolis-Hastings Robbins-Monro (Cai, 2010) is motivated by: Fisher (1925) The gradient of the observed log likelihood = the expectation of the gradient of the complete log likelihood.
10 A Review of MH-RM 10 / 34 Metropolis-Hastings Robbins-Monro (Cai, 2010) is motivated by: Fisher (1925) The gradient of the observed log likelihood = the expectation of the gradient of the complete log likelihood. Iterative approach: Use MH sampler to find an MC approximation of the expectation. Use this approximation to update the parameter estimates.
11 A Review of MH-RM 11 / 34 Metropolis-Hastings Robbins-Monro (Cai, 2010) is motivated by: Fisher (1925) The gradient of the observed log likelihood = the expectation of the gradient of the complete log likelihood. Iterative approach: Use MH sampler to find an MC approximation of the expectation. Use this approximation to update the parameter estimates. Since the approximations are noisy, use the RM method to filter. The filtering operates through decreasing gain constants.
12 MH-RM Estimation in a Picture 12 / 34 Stage 1 (iterations 1-800): move to neighborhood of MLE. Parameter Estimate Iteration Number
13 13 / 34 MH-RM Estimation in a Picture Stage 1 (iterations 1-800): move to neighborhood of MLE. Stage 2 (iterations ): find starting values for Stage 3 by averaging estimates across Stage 2 iterations. Parameter Estimate Iteration Number
14 14 / 34 MH-RM Estimation in a Picture Stage 1 (iterations 1-800): move to neighborhood of MLE. Stage 2 (iterations ): find starting values for Stage 3 by averaging estimates across Stage 2 iterations. Stage 3 (iterations 1001-convergence): obtain MLE by using decreasing gain constants to average out the noise. Parameter Estimate Iteration Number
15 15 / 34 Standard Errors in MH-RM Louis (1982) The observed information can be expressed as a function of the complete log likelihood, involving its gradient and hessian. With MH-RM, the pieces used in the approximation can be calculated by: recursive approximation, as in Cai (2010); or Monte Carlo approximation, as in Diebolt and Ip (1996).
16 Outline 16 / 34 1 Introduction and Motivation 2 Review of Ramsay Curves 3 Review of the MH-RM Algorithm 4 Simulation Study 5 Empirical Application 6 Conclusion
17 Simulation Study 17 / 34 Simulation Study Conditions Sample Size 1000 Number of Items 25 Slope Dist. N (1.7, 0.3) Location Dist. N (0, 1), truncated at ±2 True g(θ η) Normal, Skewed, and Bimodal Replications 100 Estimation Methods MH-RM and BA-EM
18 18 / 34 True g(θ η) used in Simulation Study Simulation Densities Normal Skewed Bimodal The Skewed and Bimodal densities were found as mixtures of normals. Latent trait scores were drawn by rejection sampling. Item responses were simulated in accordance with the 2PL model.
19 Parameter Recovery Log likelihood values for MH-RM (x-axis) and BA-EM (y-axis) Normal MH RM EM Skewed MH RM EM Bimodal MH RM EM Both MH-RM and EM perform maximum marginal likelihood estimation, so we shouldn t expect substantial differences in estimation performance. 19 / 34
20 Parameter Recovery Average RMSE for item parameters within a replication for MH-RM (x-axis) and BA-EM (y-axis) Normal MH RM EM Skewed MH RM EM Bimodal MH RM EM Takeaway: MH-RM and BA-EM are producing comparable results. 20 / 34
21 21 / 34 Ramsay Curve Recovery: Normal g(θ η) Density θ estimated RCs by replication shown in different colors true g(θ η) shown in black
22 22 / 34 Ramsay Curve Recovery: Skewed g(θ η) Density θ estimated RCs by replication shown in different colors true g(θ η) shown in black
23 23 / 34 Ramsay Curve Recovery: Bimodal g(θ η) Density θ estimated RCs by replication shown in different colors true g(θ η) shown in black
24 Standard Error Estimation 24 / 34 For item parameters, average standard errors (x-axis) and Monte Carlo standard deviations (y-axis). Normal Skewed Bimodal Monte Carlo SD Monte Carlo SD Monte Carlo SD Standard Error Standard Error Standard Error Note: For all three g(θ η) shapes, the MC standard deviations are slightly larger.
25 25 / 34 Standard Error Estimation Examining coverage probabilities (based on 100 replications): 68% coverage 95% coverage g(θ η) a c all a c all Normal Skewed Bimodal a = slope c = intercept all = slopes and intercepts
26 Outline 26 / 34 1 Introduction and Motivation 2 Review of Ramsay Curves 3 Review of the MH-RM Algorithm 4 Simulation Study 5 Empirical Application 6 Conclusion
27 Example: National Comorbidity Survey 27 / 34 Background: The National Comorbidity Survey (NCS, 1994) was a nationwide household survey of the United States. Ten Likert items (4-categories) measured current emotional distress. Example Items: During the past 30 days how often did you... worry too much about things? feel exhausted for no good reason? feel hopeless about the future? Response options: never (0), rarely (1), sometimes (2), often (3).
28 28 / 34 Distress Scale Results: Ramsay Curve Densities Distress Normal Note: there appears to be a sizeable group of respondents at θ = 3.
29 Distress Scale Results: Slope Estimates and Standard Errors 29 / 34 Item Slope Estimates Standard Errors Blue values correspond to RC estimation. Black values correspond to IRT estimation with a normal distribution.
30 Outline 30 / 34 1 Introduction and Motivation 2 Review of Ramsay Curves 3 Review of the MH-RM Algorithm 4 Simulation Study 5 Empirical Application 6 Conclusion
31 31 / 34 Discussion Overall, results support the validity and utility of the MH-RM implementation. Embedding RC-IRT in MH-RM (as opposed to EM) enables the computation of the observed information matrix. This, in turn, facilitates: test assembly limited-information goodness-of-fit testing differential item functioning
32 References I 32 / 34 Bock, R. D., & Aitkin, M. (1981). Marginal maximum likelihood estimation of item parameters: Application of an EM algorithm. Psychometrika, 46, Cai, L. (2010). High-dimensional exploratory item factor analysis by a Metropolis-Hastings Robbins-Monro algorithm. Psychometrika, 75, Cai, L. (2012). flexmirt: Flexible multilevel item factor analysis and test scoring [Computer software]. Seattle, WA: Vector Psychometric Group, LLC. Diebolt, J., & Ip, E. H. S. (1996). Stochastic EM: method and application. In W. Gilks, S. Richardson, & D. Spiegelhalter (Eds.), Markov chain Monte Carlo in practice (p ). London: Chapman and Hall. Fisher, R. A. (1925). Theory of statistical estimation. Proceedings of the Cambridge Philosophical Society, 22,
33 References II 33 / 34 Louis, T. A. (1982). Finding the observed information matrix when using the EM algorithm. Journal of the Royal Statistical Society - Series B, 44(2), Mislevy, R. J. (1984). Estimating latent distributions. Psychometrika, 49(3), R Development Core Team. (2010). R: A language and environment for statistical computing [Computer software manual]. Vienna, Austria. Available from (ISBN ) Woods, C. M., & Lin, N. (2008). Item response theory with estimation of the latent density using Davidian curves. Applied Psychological Measurement, 33(2), Woods, C. M., & Thissen, D. (2006). Item response theory with estimation of the latent population distribution using spline-based densities. Psychometrika, 71,
34 Acknowledgements 34 / 34 This research is supported by grants from the Institute of Education Sciences (R305B and R305D100039) and the National Institute on Drug Abuse (R01DA and R01DA030466). Thanks to Larry Thomas and Mark Hansen for providing datasets used in the empirical applications.
On Implementation of the Markov Chain Monte Carlo Stochastic Approximation Algorithm
On Implementation of the Markov Chain Monte Carlo Stochastic Approximation Algorithm Yihua Jiang, Peter Karcher and Yuedong Wang Abstract The Markov Chain Monte Carlo Stochastic Approximation Algorithm
More informationEstimation of the Markov-switching GARCH model by a Monte Carlo EM algorithm
Estimation of the Markov-switching GARCH model by a Monte Carlo EM algorithm Maciej Augustyniak Fields Institute February 3, 0 Stylized facts of financial data GARCH Regime-switching MS-GARCH Agenda Available
More informationCalibration of Interest Rates
WDS'12 Proceedings of Contributed Papers, Part I, 25 30, 2012. ISBN 978-80-7378-224-5 MATFYZPRESS Calibration of Interest Rates J. Černý Charles University, Faculty of Mathematics and Physics, Prague,
More informationWC-5 Just How Credible Is That Employer? Exploring GLMs and Multilevel Modeling for NCCI s Excess Loss Factor Methodology
Antitrust Notice The Casualty Actuarial Society is committed to adhering strictly to the letter and spirit of the antitrust laws. Seminars conducted under the auspices of the CAS are designed solely to
More informationCOS 513: Gibbs Sampling
COS 513: Gibbs Sampling Matthew Salesi December 6, 2010 1 Overview Concluding the coverage of Markov chain Monte Carlo (MCMC) sampling methods, we look today at Gibbs sampling. Gibbs sampling is a simple
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 informationStochastic Approximation Algorithms and Applications
Harold J. Kushner G. George Yin Stochastic Approximation Algorithms and Applications With 24 Figures Springer Contents Preface and Introduction xiii 1 Introduction: Applications and Issues 1 1.0 Outline
More informationMaximum Likelihood Estimation
Maximum Likelihood Estimation EPSY 905: Fundamentals of Multivariate Modeling Online Lecture #6 EPSY 905: Maximum Likelihood In This Lecture The basics of maximum likelihood estimation Ø The engine that
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 informationA NEW POINT ESTIMATOR FOR THE MEDIAN OF GAMMA DISTRIBUTION
Banneheka, B.M.S.G., Ekanayake, G.E.M.U.P.D. Viyodaya Journal of Science, 009. Vol 4. pp. 95-03 A NEW POINT ESTIMATOR FOR THE MEDIAN OF GAMMA DISTRIBUTION B.M.S.G. Banneheka Department of Statistics and
More informationThe histogram should resemble the uniform density, the mean should be close to 0.5, and the standard deviation should be close to 1/ 12 =
Chapter 19 Monte Carlo Valuation Question 19.1 The histogram should resemble the uniform density, the mean should be close to.5, and the standard deviation should be close to 1/ 1 =.887. Question 19. The
More informationRelevant parameter changes in structural break models
Relevant parameter changes in structural break models A. Dufays J. Rombouts Forecasting from Complexity April 27 th, 2018 1 Outline Sparse Change-Point models 1. Motivation 2. Model specification Shrinkage
More informationEstimating Mixed Logit Models with Large Choice Sets. Roger H. von Haefen, NC State & NBER Adam Domanski, NOAA July 2013
Estimating Mixed Logit Models with Large Choice Sets Roger H. von Haefen, NC State & NBER Adam Domanski, NOAA July 2013 Motivation Bayer et al. (JPE, 2007) Sorting modeling / housing choice 250,000 individuals
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 informationA Test of the Normality Assumption in the Ordered Probit Model *
A Test of the Normality Assumption in the Ordered Probit Model * Paul A. Johnson Working Paper No. 34 March 1996 * Assistant Professor, Vassar College. I thank Jahyeong Koo, Jim Ziliak and an anonymous
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 informationA comment on Christoffersen, Jacobs and Ornthanalai (2012), Dynamic jump intensities and risk premiums: Evidence from S&P500 returns and options
A comment on Christoffersen, Jacobs and Ornthanalai (2012), Dynamic jump intensities and risk premiums: Evidence from S&P500 returns and options Garland Durham 1 John Geweke 2 Pulak Ghosh 3 February 25,
More informationUsing MCMC and particle filters to forecast stochastic volatility and jumps in financial time series
Using MCMC and particle filters to forecast stochastic volatility and jumps in financial time series Ing. Milan Fičura DYME (Dynamical Methods in Economics) University of Economics, Prague 15.6.2016 Outline
More informationPARAMETRIC AND NON-PARAMETRIC BOOTSTRAP: A SIMULATION STUDY FOR A LINEAR REGRESSION WITH RESIDUALS FROM A MIXTURE OF LAPLACE DISTRIBUTIONS
PARAMETRIC AND NON-PARAMETRIC BOOTSTRAP: A SIMULATION STUDY FOR A LINEAR REGRESSION WITH RESIDUALS FROM A MIXTURE OF LAPLACE DISTRIBUTIONS Melfi Alrasheedi School of Business, King Faisal University, Saudi
More informationProbits. Catalina Stefanescu, Vance W. Berger Scott Hershberger. Abstract
Probits Catalina Stefanescu, Vance W. Berger Scott Hershberger Abstract Probit models belong to the class of latent variable threshold models for analyzing binary data. They arise by assuming that the
More informationMonitoring Processes with Highly Censored Data
Monitoring Processes with Highly Censored Data Stefan H. Steiner and R. Jock MacKay Dept. of Statistics and Actuarial Sciences University of Waterloo Waterloo, N2L 3G1 Canada The need for process monitoring
More informationAn Improved Saddlepoint Approximation Based on the Negative Binomial Distribution for the General Birth Process
Computational Statistics 17 (March 2002), 17 28. An Improved Saddlepoint Approximation Based on the Negative Binomial Distribution for the General Birth Process Gordon K. Smyth and Heather M. Podlich Department
More informationMarket Risk Analysis Volume I
Market Risk Analysis Volume I Quantitative Methods in Finance Carol Alexander John Wiley & Sons, Ltd List of Figures List of Tables List of Examples Foreword Preface to Volume I xiii xvi xvii xix xxiii
More informationWeek 1 Quantitative Analysis of Financial Markets Distributions B
Week 1 Quantitative Analysis of Financial Markets Distributions B Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 October
More informationLikelihood-based Optimization of Threat Operation Timeline Estimation
12th International Conference on Information Fusion Seattle, WA, USA, July 6-9, 2009 Likelihood-based Optimization of Threat Operation Timeline Estimation Gregory A. Godfrey Advanced Mathematics Applications
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 informationA Two-Step Estimator for Missing Values in Probit Model Covariates
WORKING PAPER 3/2015 A Two-Step Estimator for Missing Values in Probit Model Covariates Lisha Wang and Thomas Laitila Statistics ISSN 1403-0586 http://www.oru.se/institutioner/handelshogskolan-vid-orebro-universitet/forskning/publikationer/working-papers/
More informationLecture 2: Forecasting stock returns
Lecture 2: Forecasting stock returns Prof. Massimo Guidolin Advanced Financial Econometrics III Winter/Spring 2018 Overview The objective of the predictability exercise on stock index returns Predictability
More informationLecture 2: Forecasting stock returns
Lecture 2: Forecasting stock returns Prof. Massimo Guidolin Advanced Financial Econometrics III Winter/Spring 2016 Overview The objective of the predictability exercise on stock index returns Predictability
More informationChapter 2 Uncertainty Analysis and Sampling Techniques
Chapter 2 Uncertainty Analysis and Sampling Techniques The probabilistic or stochastic modeling (Fig. 2.) iterative loop in the stochastic optimization procedure (Fig..4 in Chap. ) involves:. Specifying
More informationInverse reinforcement learning from summary data
Inverse reinforcement learning from summary data Antti Kangasrääsiö, Samuel Kaski Aalto University, Finland ECML PKDD 2018 journal track Published in Machine Learning (2018), 107:1517 1535 September 12,
More informationIranian Journal of Economic Studies. Inflation Behavior in Top Sukuk Issuing Countries : Using a Bayesian Log-linear Model
Iranian Journal of Economic Studies, 6(1) 017, 9-46 Iranian Journal of Economic Studies Journal homepage: ijes.shirazu.ac.ir Inflation Behavior in Top Sukuk Issuing Countries : Using a Bayesian Log-linear
More informationImproving the accuracy of estimates for complex sampling in auditing 1.
Improving the accuracy of estimates for complex sampling in auditing 1. Y. G. Berger 1 P. M. Chiodini 2 M. Zenga 2 1 University of Southampton (UK) 2 University of Milano-Bicocca (Italy) 14-06-2017 1 The
More informationA Practical Implementation of the Gibbs Sampler for Mixture of Distributions: Application to the Determination of Specifications in Food Industry
A Practical Implementation of the for Mixture of Distributions: Application to the Determination of Specifications in Food Industry Julien Cornebise 1 Myriam Maumy 2 Philippe Girard 3 1 Ecole Supérieure
More informationEmpirical Distribution Testing of Economic Scenario Generators
1/27 Empirical Distribution Testing of Economic Scenario Generators Gary Venter University of New South Wales 2/27 STATISTICAL CONCEPTUAL BACKGROUND "All models are wrong but some are useful"; George Box
More informationA Markov Chain Monte Carlo Approach to Estimate the Risks of Extremely Large Insurance Claims
International Journal of Business and Economics, 007, Vol. 6, No. 3, 5-36 A Markov Chain Monte Carlo Approach to Estimate the Risks of Extremely Large Insurance Claims Wan-Kai Pang * Department of Applied
More informationExtracting Information from the Markets: A Bayesian Approach
Extracting Information from the Markets: A Bayesian Approach Daniel Waggoner The Federal Reserve Bank of Atlanta Florida State University, February 29, 2008 Disclaimer: The views expressed are the author
More informationAsymmetric Type II Compound Laplace Distributions and its Properties
CHAPTER 4 Asymmetric Type II Compound Laplace Distributions and its Properties 4. Introduction Recently there is a growing trend in the literature on parametric families of asymmetric distributions which
More informationGamma Distribution Fitting
Chapter 552 Gamma Distribution Fitting Introduction This module fits the gamma probability distributions to a complete or censored set of individual or grouped data values. It outputs various statistics
More informationInvesting through Economic Cycles with Ensemble Machine Learning Algorithms
Investing through Economic Cycles with Ensemble Machine Learning Algorithms Thomas Raffinot Silex Investment Partners Big Data in Finance Conference Thomas Raffinot (Silex-IP) Economic Cycles-Machine Learning
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 information1 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 informationA Hidden Markov Model Approach to Information-Based Trading: Theory and Applications
A Hidden Markov Model Approach to Information-Based Trading: Theory and Applications Online Supplementary Appendix Xiangkang Yin and Jing Zhao La Trobe University Corresponding author, Department of Finance,
More informationOn the Distributional Assumptions in the StoNED model
INSTITUTT FOR FORETAKSØKONOMI DEPARTMENT OF BUSINESS AND MANAGEMENT SCIENCE FOR 24 2015 ISSN: 1500-4066 September 2015 Discussion paper On the Distributional Assumptions in the StoNED model BY Xiaomei
More informationModelling component reliability using warranty data
ANZIAM J. 53 (EMAC2011) pp.c437 C450, 2012 C437 Modelling component reliability using warranty data Raymond Summit 1 (Received 10 January 2012; revised 10 July 2012) Abstract Accelerated testing is often
More informationUnobserved Heterogeneity Revisited
Unobserved Heterogeneity Revisited Robert A. Miller Dynamic Discrete Choice March 2018 Miller (Dynamic Discrete Choice) cemmap 7 March 2018 1 / 24 Distributional Assumptions about the Unobserved Variables
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 informationMarket Correlations in the Euro Changeover Period With a View to Portfolio Management
Preprint, April 2010 Market Correlations in the Euro Changeover Period With a View to Portfolio Management Gernot Müller Keywords: European Monetary Union European Currencies Markov Chain Monte Carlo Minimum
More informationLecture 9: Markov and Regime
Lecture 9: Markov and Regime Switching Models Prof. Massimo Guidolin 20192 Financial Econometrics Spring 2017 Overview Motivation Deterministic vs. Endogeneous, Stochastic Switching Dummy Regressiom Switching
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 informationGetting started with WinBUGS
1 Getting started with WinBUGS James B. Elsner and Thomas H. Jagger Department of Geography, Florida State University Some material for this tutorial was taken from http://www.unt.edu/rss/class/rich/5840/session1.doc
More informationOption Pricing Using Bayesian Neural Networks
Option Pricing Using Bayesian Neural Networks Michael Maio Pires, Tshilidzi Marwala School of Electrical and Information Engineering, University of the Witwatersrand, 2050, South Africa m.pires@ee.wits.ac.za,
More informationFast Convergence of Regress-later Series Estimators
Fast Convergence of Regress-later Series Estimators New Thinking in Finance, London Eric Beutner, Antoon Pelsser, Janina Schweizer Maastricht University & Kleynen Consultants 12 February 2014 Beutner Pelsser
More informationAdaptive Interest Rate Modelling
Modelling Mengmeng Guo Wolfgang Karl Härdle Ladislaus von Bortkiewicz Chair of Statistics C.A.S.E. - Center for Applied Statistics and Economics Humboldt-Universität zu Berlin http://lvb.wiwi.hu-berlin.de
More informationFitting parametric distributions using R: the fitdistrplus package
Fitting parametric distributions using R: the fitdistrplus package M. L. Delignette-Muller - CNRS UMR 5558 R. Pouillot J.-B. Denis - INRA MIAJ user! 2009,10/07/2009 Background Specifying the probability
More informationAppendix. A.1 Independent Random Effects (Baseline)
A Appendix A.1 Independent Random Effects (Baseline) 36 Table 2: Detailed Monte Carlo Results Logit Fixed Effects Clustered Random Effects Random Coefficients c Coeff. SE SD Coeff. SE SD Coeff. SE SD Coeff.
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 informationLecture 8: Markov and Regime
Lecture 8: Markov and Regime Switching Models Prof. Massimo Guidolin 20192 Financial Econometrics Spring 2016 Overview Motivation Deterministic vs. Endogeneous, Stochastic Switching Dummy Regressiom Switching
More informationAdaptive Metropolis-Hastings samplers for the Bayesian analysis of large linear Gaussian systems
Adaptive Metropolis-Hastings samplers for the Bayesian analysis of large linear Gaussian systems Stephen KH Yeung stephen.yeung@ncl.ac.uk Darren J Wilkinson d.j.wilkinson@ncl.ac.uk Department of Statistics,
More informationEmpirical Analysis of the US Swap Curve Gough, O., Juneja, J.A., Nowman, K.B. and Van Dellen, S.
WestminsterResearch http://www.westminster.ac.uk/westminsterresearch Empirical Analysis of the US Swap Curve Gough, O., Juneja, J.A., Nowman, K.B. and Van Dellen, S. This is a copy of the final version
More informationAn 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 informationInternational Journal of Computer Engineering and Applications, Volume XII, Issue II, Feb. 18, ISSN
Volume XII, Issue II, Feb. 18, www.ijcea.com ISSN 31-3469 AN INVESTIGATION OF FINANCIAL TIME SERIES PREDICTION USING BACK PROPAGATION NEURAL NETWORKS K. Jayanthi, Dr. K. Suresh 1 Department of Computer
More informationGENERATION OF STANDARD NORMAL RANDOM NUMBERS. Naveen Kumar Boiroju and M. Krishna Reddy
GENERATION OF STANDARD NORMAL RANDOM NUMBERS Naveen Kumar Boiroju and M. Krishna Reddy Department of Statistics, Osmania University, Hyderabad- 500 007, INDIA Email: nanibyrozu@gmail.com, reddymk54@gmail.com
More informationLOSS SEVERITY DISTRIBUTION ESTIMATION OF OPERATIONAL RISK USING GAUSSIAN MIXTURE MODEL FOR LOSS DISTRIBUTION APPROACH
LOSS SEVERITY DISTRIBUTION ESTIMATION OF OPERATIONAL RISK USING GAUSSIAN MIXTURE MODEL FOR LOSS DISTRIBUTION APPROACH Seli Siti Sholihat 1 Hendri Murfi 2 1 Department of Accounting, Faculty of Economics,
More informationUsing New SAS 9.4 Features for Cumulative Logit Models with Partial Proportional Odds Paul J. Hilliard, Educational Testing Service (ETS)
Using New SAS 9.4 Features for Cumulative Logit Models with Partial Proportional Odds Using New SAS 9.4 Features for Cumulative Logit Models with Partial Proportional Odds INTRODUCTION Multicategory Logit
More informationMonte 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 informationJournal of Economics and Financial Analysis, Vol:1, No:1 (2017) 1-13
Journal of Economics and Financial Analysis, Vol:1, No:1 (2017) 1-13 Journal of Economics and Financial Analysis Type: Double Blind Peer Reviewed Scientific Journal Printed ISSN: 2521-6627 Online ISSN:
More informationFinancial Econometrics Notes. Kevin Sheppard University of Oxford
Financial Econometrics Notes Kevin Sheppard University of Oxford Monday 15 th January, 2018 2 This version: 22:52, Monday 15 th January, 2018 2018 Kevin Sheppard ii Contents 1 Probability, Random Variables
More informationLecture 17: More on Markov Decision Processes. Reinforcement learning
Lecture 17: More on Markov Decision Processes. Reinforcement learning Learning a model: maximum likelihood Learning a value function directly Monte Carlo Temporal-difference (TD) learning COMP-424, Lecture
More informationInternational Journal of Computer Engineering and Applications, Volume XII, Issue II, Feb. 18, ISSN
International Journal of Computer Engineering and Applications, Volume XII, Issue II, Feb. 18, www.ijcea.com ISSN 31-3469 AN INVESTIGATION OF FINANCIAL TIME SERIES PREDICTION USING BACK PROPAGATION NEURAL
More informationMCMC Package Example
MCMC Package Example Charles J. Geyer April 4, 2005 This is an example of using the mcmc package in R. The problem comes from a take-home question on a (take-home) PhD qualifying exam (School of Statistics,
More informationAdjusted Priors for Bayes Factors Involving Reparameterized Order Constraints
Adjusted Priors for Bayes Factors Involving Reparameterized Order Constraints Supplementary Material Daniel W. Heck & Eric-Jan Wagenmakers April 29, 2016 Contents 1 The Product-Binomial Model 2 1.1 Parameter
More informationADVANCED OPERATIONAL RISK MODELLING IN BANKS AND INSURANCE COMPANIES
Small business banking and financing: a global perspective Cagliari, 25-26 May 2007 ADVANCED OPERATIONAL RISK MODELLING IN BANKS AND INSURANCE COMPANIES C. Angela, R. Bisignani, G. Masala, M. Micocci 1
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 informationELEMENTS OF MONTE CARLO SIMULATION
APPENDIX B ELEMENTS OF MONTE CARLO SIMULATION B. GENERAL CONCEPT The basic idea of Monte Carlo simulation is to create a series of experimental samples using a random number sequence. According to the
More informationImportance sampling and Monte Carlo-based calibration for time-changed Lévy processes
Importance sampling and Monte Carlo-based calibration for time-changed Lévy processes Stefan Kassberger Thomas Liebmann BFS 2010 1 Motivation 2 Time-changed Lévy-models and Esscher transforms 3 Applications
More informationEfficiency Measurement with the Weibull Stochastic Frontier*
OXFORD BULLETIN OF ECONOMICS AND STATISTICS, 69, 5 (2007) 0305-9049 doi: 10.1111/j.1468-0084.2007.00475.x Efficiency Measurement with the Weibull Stochastic Frontier* Efthymios G. Tsionas Department of
More informationIs the Ex ante Premium Always Positive? Evidence and Analysis from Australia
Is the Ex ante Premium Always Positive? Evidence and Analysis from Australia Kathleen D Walsh * School of Banking and Finance University of New South Wales This Draft: Oct 004 Abstract: An implicit assumption
More informationSupplementary material for the paper Identifiability and bias reduction in the skew-probit model for a binary response
Supplementary material for the paper Identifiability and bias reduction in the skew-probit model for a binary response DongHyuk Lee and Samiran Sinha Department of Statistics, Texas A&M University, College
More informationBayesian Estimation of the Markov-Switching GARCH(1,1) Model with Student-t Innovations
Bayesian Estimation of the Markov-Switching GARCH(1,1) Model with Student-t Innovations Department of Quantitative Economics, Switzerland david.ardia@unifr.ch R/Rmetrics User and Developer Workshop, Meielisalp,
More informationM.S. in Quantitative Finance & Risk Analytics (QFRA) Fall 2017 & Spring 2018
M.S. in Quantitative Finance & Risk Analytics (QFRA) Fall 2017 & Spring 2018 2 - Required Professional Development &Career Workshops MGMT 7770 Prof. Development Workshop 1/Career Workshops (Fall) Wed.
More informationThe Leveled Chain Ladder Model. for Stochastic Loss Reserving
The Leveled Chain Ladder Model for Stochastic Loss Reserving Glenn Meyers, FCAS, MAAA, CERA, Ph.D. Abstract The popular chain ladder model forms its estimate by applying age-to-age factors to the latest
More informationDiscrete Choice Methods with Simulation
Discrete Choice Methods with Simulation Kenneth E. Train University of California, Berkeley and National Economic Research Associates, Inc. iii To Daniel McFadden and in memory of Kenneth Train, Sr. ii
More informationEstimating a Dynamic Oligopolistic Game with Serially Correlated Unobserved Production Costs. SS223B-Empirical IO
Estimating a Dynamic Oligopolistic Game with Serially Correlated Unobserved Production Costs SS223B-Empirical IO Motivation There have been substantial recent developments in the empirical literature on
More informationMonte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50)
Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50) Magnus Wiktorsson Centre for Mathematical Sciences Lund University, Sweden Lecture 5 Sequential Monte Carlo methods I January
More informationOmitted Variables Bias in Regime-Switching Models with Slope-Constrained Estimators: Evidence from Monte Carlo Simulations
Journal of Statistical and Econometric Methods, vol. 2, no.3, 2013, 49-55 ISSN: 2051-5057 (print version), 2051-5065(online) Scienpress Ltd, 2013 Omitted Variables Bias in Regime-Switching Models with
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 informationSTARRY GOLD ACADEMY , , Page 1
ICAN KNOWLEDGE LEVEL QUANTITATIVE TECHNIQUE IN BUSINESS MOCK EXAMINATION QUESTIONS FOR NOVEMBER 2016 DIET. INSTRUCTION: ATTEMPT ALL QUESTIONS IN THIS SECTION OBJECTIVE QUESTIONS Given the following sample
More informationAlgorithmic Trading using Reinforcement Learning augmented with Hidden Markov Model
Algorithmic Trading using Reinforcement Learning augmented with Hidden Markov Model Simerjot Kaur (sk3391) Stanford University Abstract This work presents a novel algorithmic trading system based on reinforcement
More informationCalibration of PD term structures: to be Markov or not to be
CUTTING EDGE. CREDIT RISK Calibration of PD term structures: to be Markov or not to be A common discussion in credit risk modelling is the question of whether term structures of default probabilities can
More informationFitting financial time series returns distributions: a mixture normality approach
Fitting financial time series returns distributions: a mixture normality approach Riccardo Bramante and Diego Zappa * Abstract Value at Risk has emerged as a useful tool to risk management. A relevant
More information[D7] PROBABILITY DISTRIBUTION OF OUTSTANDING LIABILITY FROM INDIVIDUAL PAYMENTS DATA Contributed by T S Wright
Faculty and Institute of Actuaries Claims Reserving Manual v.2 (09/1997) Section D7 [D7] PROBABILITY DISTRIBUTION OF OUTSTANDING LIABILITY FROM INDIVIDUAL PAYMENTS DATA Contributed by T S Wright 1. Introduction
More informationModeling skewness and kurtosis in Stochastic Volatility Models
Modeling skewness and kurtosis in Stochastic Volatility Models Georgios Tsiotas University of Crete, Department of Economics, GR December 19, 2006 Abstract Stochastic volatility models have been seen as
More informationAccelerated Option Pricing Multiple Scenarios
Accelerated Option Pricing in Multiple Scenarios 04.07.2008 Stefan Dirnstorfer (stefan@thetaris.com) Andreas J. Grau (grau@thetaris.com) 1 Abstract This paper covers a massive acceleration of Monte-Carlo
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 informationRisk Measurement in Credit Portfolio Models
9 th DGVFM Scientific Day 30 April 2010 1 Risk Measurement in Credit Portfolio Models 9 th DGVFM Scientific Day 30 April 2010 9 th DGVFM Scientific Day 30 April 2010 2 Quantitative Risk Management Profit
More information1 Introduction. Term Paper: The Hall and Taylor Model in Duali 1. Yumin Li 5/8/2012
Term Paper: The Hall and Taylor Model in Duali 1 Yumin Li 5/8/2012 1 Introduction In macroeconomics and policy making arena, it is extremely important to have the ability to manipulate a set of control
More information1 Residual life for gamma and Weibull distributions
Supplement to Tail Estimation for Window Censored Processes Residual life for gamma and Weibull distributions. Gamma distribution Let Γ(k, x = x yk e y dy be the upper incomplete gamma function, and let
More informationCalibrating Financial Models Using Consistent Bayesian Estimators
Calibrating Financial Models Using Consistent Bayesian Estimators Christoph Reisinger Joint work with Alok Gupta June 25, 2010 Example model uncertainty A local volatility model, jump diffusion model,
More informationNo Predictable Components in G7 Stock Returns
No Predictable Components in G7 Stock Returns Prasad V. Bidarkota* Khurshid M. Kiani Abstract: We search for time-varying predictable components in monthly excess stock index returns over the risk free
More information