GENERATION OF STANDARD NORMAL RANDOM NUMBERS. Naveen Kumar Boiroju and M. Krishna Reddy
|
|
- Elizabeth Lindsey
- 6 years ago
- Views:
Transcription
1 GENERATION OF STANDARD NORMAL RANDOM NUMBERS Naveen Kumar Boiroju and M. Krishna Reddy Department of Statistics, Osmania University, Hyderabad , INDIA Abstract: In this paper, we propose an easy method to generate standard normal random variables and this method is compared with the other existing methods. It is observed that the proposed method is quite competitive with most of the existing methods. Randomness of the generated sample is tested using the autocorrelation function and Ljung-Box test. Normality of the samples is checked with the help of relation between testing of hypothesis and confidence intervals. Keywords: Autocorrelation function, Bootstrap confidence intervals, Central limit theorem, Ljung-Box Statistic, Standard normal distribution, Random numbers.. Introduction The normal distribution is also called as the Gaussian distribution and it is the most important distribution in statistics. The distribution is given by f ( x) ( x µ ) = exp ; <, < > 0 x µ and σ σ π σ
2 where µ is a location parameter, equal to the mean, and σ the standard deviation. For µ = 0 and σ =, we refer to this distribution as the standard normal distribution. Some characteristics of a standard normal distribution include the following:. Symmetric about mean and Mean=Median=Mode.. Mean=0, Standard deviation =, Skewness=0 and Kurtosis relative to normal ( β 3) =0. 3. First and third quartiles are and respectively. 4. Normal random variables X s with mean µ and variance σ are generated by the relationship X = µ + σ Z, where Z is the standard normal random variable. Generation or simulation of random numbers, using deterministic algorithms, is widely used by statisticians for several purposes. These include evaluating performance of estimation and test procedures, and also use of simulation as a part of statistical inference. Generating normal random numbers is an old and very important problem in the statistical literature. Standard normal random variables are frequently used in computer science, computational statistics, and in particular, in applications of the Monte Carlo method. Several algorithms are available in the literature to generate normal random numbers like method of sums of uniforms, Box-Muller method, Polar method, acceptance-rejection method, inversion method, etc. Howell and Rheinfurth (98), Devroye (986), Johnson et.al. (995), Thomas et.al. (007) and Rao et.al. (0) provides an extensive list of different algorithms for the generation of standard normal random numbers. In view of the enormous recent development in statistics that relies on the use of simulation, it is essential for every statistician to have some knowledge of the corresponding techniques. It is also useful to note that because of its extensive use, search for an appropriate simulation technique is itself an area of research. In Section, some of the existing methods to generate standard normal
3 random numbers discussed and a very simple and convenient method is proposed. Testing of randomness and normality of the generated samples is presented in Section 3. Final conclusion is presented in Section 4.. Methods for Generating Standard Normal Random Numbers Some existing methods for generating standard normal random numbers discussed in this section. A new algorithm to generate standard normal random numbers is also proposed and is named as method-9 in this article. Method-: Sum of Uniform Random Variables The simplest way of generating normal variables is an application of the central limit theorem. The central limit theorem is a weak convergence result that expresses the fact that any sum of many small independent random variables is approximately normally distributed. Use of the central limit theorem on U(0,) random variables provide a simple method for closely approximating normal random variates. The following algorithm is used to generate the standard normal variables.. Generate independent uniform random numbers,. Return Z = U i i= 6. U, U, K, U ( ) ~ iid U 0,. This method requires uniform random variables to generate a single standard normal random number. 3
4 Method-: Box-Muller Method This method is due to Box and Muller (958) and generates a pair of independent standard normal random variables using a pair of uniform random numbers as follows:. Generate two independent random numbers Uand U from U(0,) distribution.. Return Z = lnu Cos( πu ) and Z lnu Sin( πu ) =. Method-3: Polar Method Another form of the Box-Muller method is called the polar technique. This improves over the previous technique in being quicker as well as numerically more robust. The algorithm can be summarized as follows:. Generate two random numbers Uand U from U(0,) distribution.. Set V = U, V = U and S = V +. Note that V and V are U(-,). V 3. If S>, go to step, otherwise go to step Return two independent standard normal variables Z ln S S ln S =. S = V and Z V This algorithm also requires two uniform variables to generate a single standard normal random number. Method-4: Inversion Method follows: This method is due to Hastings (955) and the algorithm can be summarized as 4
5 . Generate U ~ U(0,).. Return Z = Φ ( U ) 3.if U 0.5, use Z = -Φ = -w + i= 0 3 j= 0 ( U ). i a w i i b w j ; U 0.5. where w = a a a 0 ln ( U ) =.5557; = ; = ; ; b b b b 3 0 = ; =.43788; = ; = This method requires single uniform random number but involves seven constants. Method-5: Acceptance-Rejection Method. If The following method is borrowed from the paper of Sigman (007).. Generate two independent exponentials at rate; Y ( Y ) /, set 3. Generate U.Set Z = 3 Z = Y ; otherwise go to back to step. Z if U 3 0.5, set Z = - Z if U Y = ln 3 ( U ) and Y = ln( U ). > 0.5. This algorithm is an iterative algorithm and requires at least three uniform random numbers to generate a single standard random number. Method-6: Using Generalized Exponential Distribution The following algorithm is given by Kundu et.al. (006).. Generate an uniform (0,) random variable, U ~. Generate X = -ln ( U ). 3. Compute the standard normal random variable as Z U(0,). ln = ( X)
6 This algorithm requires only one uniform random number for the generation of single standard normal random number. Method-7: Bol shev Formula A formula was suggested by Bol shev (959) is explained in the following steps.. Let the five uniform random numbers, U, U, U and U ~ ( 0,) [ U j ] 5. Compute = 3( ) X. 5 j= 3 3. Z = X 0.0( 3X X ) U., U This formula requires atleast five uniform random numbers to generate a single standard normal random number. Method-8: Inversion Method The following algorithm proposed by Rao et.al. (0) and it is explained in the following steps.. Generate U from the U(0,) distribution.. Return ln U Z =..70 This algorithm requires only a single uniform random variable to generate the standard normal variables and it is very much easier to apply and compute. 6
7 Method-9: Proposed Method The proposed algorithm is developed using feedforward neural networks and summarized as follows:. Generate uniform random number U from U(0,) distribution.. Set X = Tanh( U), X =Tanh( U), X 3 =Tanh( U) and X 4 =Tanh( U). 3. Return Z= X X X X 4. This method requires only one uniform random number to generate the standard normal random numbers. 3. Testing Randomness and Normality A comparative study of the above methods discussed in this section. 000 random numbers generated using each of the methods discussed in section and tested for the randomness and normality of the random numbers generated using the proposed method. Testing Randomness: For a random series, we expect to find that autocorrelation function (ACF) is zero for all non-zero values of lag k and none of the autocorrelations is significantly different from zero at the level Hence, the ACF is used as a graphical test procedure to test the 7
8 randomness of the generated data (Chatfield, 99). The sample ACF of the standard normal random numbers generated from the proposed method (Method-9) is presented in the following figure. Figure : Autocorrelation function of the sample generated from the proposed method As the results indicate, none of these autocorrelations of the standard normal random numbers generated from the proposed method is significantly different from zero at the level This proves that the generated sample is a random sample. Randomness of the sample generated from the proposed method is tested using Ljung- Box test (Chatfield, 99; Ljung and Box, 978). The following table presents the values of Ljung-Box statistic and its significant probability values for different lags of the autocorrelations. 8
9 Table : Ljung-Box test Lag Ljung-Box Statistic Value df P-value From the above table, it is observed that the P-values of the Ljung-Box statistic at different lags is greater than the level 0.05, therefore we accept the null hypothesis of randomness of the sample and we may conclude that the sample generated using the proposed method is a random sample. Testing of Normality: The histogram of the sample generated using the proposed method (Method-9) is presented in the following figure. 9
10 Figure : Histogram of the sample generated using the proposed method From the above histogram, it is observed that the sample generated using the proposed method is symmetric and normal. The following table presents the descriptive statistics of the samples generated using each of the methods. Table : Descriptive Statistics Descriptive Method Statistics M M M3 M4 M5 M6 M7 M8 M9 N Valid Mean Median Mode Std. Deviation Skewness Kurtosis Range Minimum Maximum Quartiles Q Q Q
11 From the above table, it is clear that the mean, median and mode are almost equal to zero and same for the samples generated from the each method. Standard deviation, skewness and kurtosis of the random numbers generated using each of the methods is near to one, zero and zero respectively. This indicates that the methods producing standard normal random numbers. The sample quartiles are also agreeing with that of the standard normal distribution. Method-4 generates a wide range of the standard random numbers whereas method-8 generates a low range of the standard normal random numbers. The relation between testing of hypothesis and confidence intervals is used to test the normality of the samples. The characteristics of the standard normal distribution is checked and tested with the help of bootstrap confidence intervals (Becher et.al.; 993 and Efron and Tibishirani; 993). If all the parameters (mean=0, standard deviation=, skewness=0 and kurtosis=0 ) of the standard normal distribution lies in between the lower and upper bounds of bootstrap confidence interval, then the normality of the sample is satisfied. Bootstrap confidence intervals for the mean, standard deviation, skewness and kurtosis of the samples generated from each of the method is presented in the following table.
12 Table 3: Bootstrap confidence intervals Bootstrap method (0,000 iterations) Method Descriptive Statistics Statistic 95% Confidence Interval Bias Std. Error Lower Upper Mean Method- Std. Deviation Skewness Kurtosis Mean Method- Std. Deviation Skewness Kurtosis Mean Method-3 Std. Deviation Skewness Kurtosis Mean Method-4 Std. Deviation Skewness Kurtosis Mean Method-5 Std. Deviation Skewness Kurtosis Mean Method-6 Std. Deviation Skewness Kurtosis Mean Method-7 Std. Deviation Skewness Kurtosis Mean Method-8 Std. Deviation Skewness Kurtosis Mean Method-9 Std. Deviation Skewness Kurtosis
13 From the above table, it is observed that the parameters of the standard normal distribution belong to the respective bootstrap confidence intervals for all the methods. This implies that all the methods considered in this article can be used to generate the standard normal random numbers. 4. Conclusion From the above study, it is observed that the proposed method possesses the randomness as well as the normality. The proposed method is an easy and convenient method to generate the standard normal random numbers. Hence, the proposed method can be viewed as a competitive method to generate the standard normal random numbers. References: [] Becher H., Hall P., Wilson S.R. (993), Bootstrap hypothesis testing procedures, Biometrics, 49(4): [] Bol shev, L.N. (959), On transformations of random variables, Teoriega Veroyatitnosteti i ee primeneniya, 4, 9-4. [3] Box, G. E. P. and Muller, M. E. (958), A note on the generation of random normal deviates, Annals of Mathematical Statistics, 9, [4] Chatfield, C. (99), The Analysis of Time Series: An Introduction, 5 th ed., Chapman and Hall, London. [5] Devroye L. (986), Non-uniform random variable distribution, Springer-Verlag: New York. [6] Efron B., Tibshirani R.J. (993), An introduction to the Bootstrap, Chapman & Hall, New York. 3
14 [7] Hastings, C.(955), Approximations for Digital Computers. Princeton University Press, Princeton, N.J. [8] Howell, L.W. and Rheinfurth, M.H. (98), Generation of Pseudo-Random Numbers, Technical Paper 05, NASA. [9] Johnson, N., Kotz, S. and Balakrishnan, N. (995), Continuous Univariate Distribution, Vol., John Wiley and Sons, New York. [0] Kundu, D., Gupta, R.D. and Manglick, A. (006), A convenient way of generating normal random variables using generalized exponential distribution J. Modern Appl. Statist. Methods, 5, (home.iitk.ac.in/~kundu/paper04.pdf.) [] Ljung, G. M. and Box, G. E. P. (978), On a measure of lack of fit in time series models, Biometrika, 65, [] Rao, K.R., Boiroju, N.K. and Reddy, M.K. (0), Generation of standard normal random variables, Indian Journal of Scientific Research, (4), [3] Sigman, K. (007), Simulating normal (Gaussian) rvs with applications to simulating Brownian motion and geometric Brownian motion in one and two dimensions, [4] Thomas, D.B., Luk, W., Leong, P.H.W., Villasenor, J.D. (007), Gaussian Random Number Generators, ACM Computing Surveys, Vol. 39, No. 4, Article. 4
A Convenient Way of Generating Normal Random Variables Using Generalized Exponential Distribution
A Convenient Way of Generating Normal Random Variables Using Generalized Exponential Distribution Debasis Kundu 1, Rameshwar D. Gupta 2 & Anubhav Manglick 1 Abstract In this paper we propose a very convenient
More informationOn the Distribution and Its Properties of the Sum of a Normal and a Doubly Truncated Normal
The Korean Communications in Statistics Vol. 13 No. 2, 2006, pp. 255-266 On the Distribution and Its Properties of the Sum of a Normal and a Doubly Truncated Normal Hea-Jung Kim 1) Abstract This paper
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 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 informationRobust Critical Values for the Jarque-bera Test for Normality
Robust Critical Values for the Jarque-bera Test for Normality PANAGIOTIS MANTALOS Jönköping International Business School Jönköping University JIBS Working Papers No. 00-8 ROBUST CRITICAL VALUES FOR THE
More informationA Skewed Truncated Cauchy Logistic. Distribution and its Moments
International Mathematical Forum, Vol. 11, 2016, no. 20, 975-988 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/imf.2016.6791 A Skewed Truncated Cauchy Logistic Distribution and its Moments Zahra
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 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 informationGENERATION OF APPROXIMATE GAMMA SAMPLES BY PARTIAL REJECTION
IASC8: December 5-8, 8, Yokohama, Japan GEERATIO OF APPROXIMATE GAMMA SAMPLES BY PARTIAL REJECTIO S.H. Ong 1 Wen Jau Lee 1 Institute of Mathematical Sciences, University of Malaya, 563 Kuala Lumpur, MALAYSIA
More informationECE 295: Lecture 03 Estimation and Confidence Interval
ECE 295: Lecture 03 Estimation and Confidence Interval Spring 2018 Prof Stanley Chan School of Electrical and Computer Engineering Purdue University 1 / 23 Theme of this Lecture What is Estimation? You
More informationKARACHI UNIVERSITY BUSINESS SCHOOL UNIVERSITY OF KARACHI BS (BBA) VI
88 P a g e B S ( B B A ) S y l l a b u s KARACHI UNIVERSITY BUSINESS SCHOOL UNIVERSITY OF KARACHI BS (BBA) VI Course Title : STATISTICS Course Number : BA(BS) 532 Credit Hours : 03 Course 1. Statistical
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 informationSample Size for Assessing Agreement between Two Methods of Measurement by Bland Altman Method
Meng-Jie Lu 1 / Wei-Hua Zhong 1 / Yu-Xiu Liu 1 / Hua-Zhang Miao 1 / Yong-Chang Li 1 / Mu-Huo Ji 2 Sample Size for Assessing Agreement between Two Methods of Measurement by Bland Altman Method Abstract:
More informationESTIMATION OF MODIFIED MEASURE OF SKEWNESS. Elsayed Ali Habib *
Electronic Journal of Applied Statistical Analysis EJASA, Electron. J. App. Stat. Anal. (2011), Vol. 4, Issue 1, 56 70 e-issn 2070-5948, DOI 10.1285/i20705948v4n1p56 2008 Università del Salento http://siba-ese.unile.it/index.php/ejasa/index
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 informationModified ratio estimators of population mean using linear combination of co-efficient of skewness and quartile deviation
CSIRO PUBLISHING The South Pacific Journal of Natural and Applied Sciences, 31, 39-44, 2013 www.publish.csiro.au/journals/spjnas 10.1071/SP13003 Modified ratio estimators of population mean using linear
More informationMean GMM. Standard error
Table 1 Simple Wavelet Analysis for stocks in the S&P 500 Index as of December 31 st 1998 ^ Shapiro- GMM Normality 6 0.9664 0.00281 11.36 4.14 55 7 0.9790 0.00300 56.58 31.69 45 8 0.9689 0.00319 403.49
More informationOn the Distribution of Kurtosis Test for Multivariate Normality
On the Distribution of Kurtosis Test for Multivariate Normality Takashi Seo and Mayumi Ariga Department of Mathematical Information Science Tokyo University of Science 1-3, Kagurazaka, Shinjuku-ku, Tokyo,
More informationSymmetricity of the Sampling Distribution of CV r for Exponential Samples
World Applied Sciences Journal 17 (Special Issue of Applied Math): 60-65, 2012 ISSN 1818-4952 IDOSI Publications, 2012 Symmetricity of the Sampling Distribution of CV r for Exponential Samples Fauziah
More informationMonte Carlo Methods in Option Pricing. UiO-STK4510 Autumn 2015
Monte Carlo Methods in Option Pricing UiO-STK4510 Autumn 015 The Basics of Monte Carlo Method Goal: Estimate the expectation θ = E[g(X)], where g is a measurable function and X is a random variable such
More informationGeneralized Modified Ratio Type Estimator for Estimation of Population Variance
Sri Lankan Journal of Applied Statistics, Vol (16-1) Generalized Modified Ratio Type Estimator for Estimation of Population Variance J. Subramani* Department of Statistics, Pondicherry University, Puducherry,
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 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 informationOptimal rebalancing of portfolios with transaction costs assuming constant risk aversion
Optimal rebalancing of portfolios with transaction costs assuming constant risk aversion Lars Holden PhD, Managing director t: +47 22852672 Norwegian Computing Center, P. O. Box 114 Blindern, NO 0314 Oslo,
More informationdiscussion Papers Some Flexible Parametric Models for Partially Adaptive Estimators of Econometric Models
discussion Papers Discussion Paper 2007-13 March 26, 2007 Some Flexible Parametric Models for Partially Adaptive Estimators of Econometric Models Christian B. Hansen Graduate School of Business at the
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 informationResults for option pricing
Results for option pricing [o,v,b]=optimal(rand(1,100000 Estimators = 0.4619 0.4617 0.4618 0.4613 0.4619 o = 0.46151 % best linear combination (true value=0.46150 v = 1.1183e-005 %variance per uniform
More informationCOMPARISON OF RATIO ESTIMATORS WITH TWO AUXILIARY VARIABLES K. RANGA RAO. College of Dairy Technology, SPVNR TSU VAFS, Kamareddy, Telangana, India
COMPARISON OF RATIO ESTIMATORS WITH TWO AUXILIARY VARIABLES K. RANGA RAO College of Dairy Technology, SPVNR TSU VAFS, Kamareddy, Telangana, India Email: rrkollu@yahoo.com Abstract: Many estimators of the
More informationM249 Diagnostic Quiz
THE OPEN UNIVERSITY Faculty of Mathematics and Computing M249 Diagnostic Quiz Prepared by the Course Team [Press to begin] c 2005, 2006 The Open University Last Revision Date: May 19, 2006 Version 4.2
More informationChapter 7. Inferences about Population Variances
Chapter 7. Inferences about Population Variances Introduction () The variability of a population s values is as important as the population mean. Hypothetical distribution of E. coli concentrations from
More informationA Robust Test for Normality
A Robust Test for Normality Liangjun Su Guanghua School of Management, Peking University Ye Chen Guanghua School of Management, Peking University Halbert White Department of Economics, UCSD March 11, 2006
More informationEFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS
Commun. Korean Math. Soc. 23 (2008), No. 2, pp. 285 294 EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS Kyoung-Sook Moon Reprinted from the Communications of the Korean Mathematical Society
More informationStatistical Analysis of Data from the Stock Markets. UiO-STK4510 Autumn 2015
Statistical Analysis of Data from the Stock Markets UiO-STK4510 Autumn 2015 Sampling Conventions We observe the price process S of some stock (or stock index) at times ft i g i=0,...,n, we denote it by
More informationدرس هفتم یادگیري ماشین. (Machine Learning) دانشگاه فردوسی مشهد دانشکده مهندسی رضا منصفی
یادگیري ماشین توزیع هاي نمونه و تخمین نقطه اي پارامترها Sampling Distributions and Point Estimation of Parameter (Machine Learning) دانشگاه فردوسی مشهد دانشکده مهندسی رضا منصفی درس هفتم 1 Outline Introduction
More informationChapter Seven. The Normal Distribution
Chapter Seven The Normal Distribution 7-1 Introduction Many continuous variables have distributions that are bellshaped and are called approximately normally distributed variables, such as the heights
More informationReview: Population, sample, and sampling distributions
Review: Population, sample, and sampling distributions A population with mean µ and standard deviation σ For instance, µ = 0, σ = 1 0 1 Sample 1, N=30 Sample 2, N=30 Sample 100000000000 InterquartileRange
More informationComputational Finance Improving Monte Carlo
Computational Finance Improving Monte Carlo School of Mathematics 2018 Monte Carlo so far... Simple to program and to understand Convergence is slow, extrapolation impossible. Forward looking method ideal
More informationOn Some Statistics for Testing the Skewness in a Population: An. Empirical Study
Available at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 12, Issue 2 (December 2017), pp. 726-752 Applications and Applied Mathematics: An International Journal (AAM) On Some Statistics
More informationA New Multivariate Kurtosis and Its Asymptotic Distribution
A ew Multivariate Kurtosis and Its Asymptotic Distribution Chiaki Miyagawa 1 and Takashi Seo 1 Department of Mathematical Information Science, Graduate School of Science, Tokyo University of Science, Tokyo,
More informationA Skewed Truncated Cauchy Uniform Distribution and Its Moments
Modern Applied Science; Vol. 0, No. 7; 206 ISSN 93-844 E-ISSN 93-852 Published by Canadian Center of Science and Education A Skewed Truncated Cauchy Uniform Distribution and Its Moments Zahra Nazemi Ashani,
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 information1) 3 points Which of the following is NOT a measure of central tendency? a) Median b) Mode c) Mean d) Range
February 19, 2004 EXAM 1 : Page 1 All sections : Geaghan Read Carefully. Give an answer in the form of a number or numeric expression where possible. Show all calculations. Use a value of 0.05 for any
More informationPower of t-test for Simple Linear Regression Model with Non-normal Error Distribution: A Quantile Function Distribution Approach
Available Online Publications J. Sci. Res. 4 (3), 609-622 (2012) JOURNAL OF SCIENTIFIC RESEARCH www.banglajol.info/index.php/jsr of t-test for Simple Linear Regression Model with Non-normal Error Distribution:
More informationJackknife Empirical Likelihood Inferences for the Skewness and Kurtosis
Georgia State University ScholarWorks @ Georgia State University Mathematics Theses Department of Mathematics and Statistics 5-10-2014 Jackknife Empirical Likelihood Inferences for the Skewness and Kurtosis
More informationClass 16. Daniel B. Rowe, Ph.D. Department of Mathematics, Statistics, and Computer Science. Marquette University MATH 1700
Class 16 Daniel B. Rowe, Ph.D. Department of Mathematics, Statistics, and Computer Science Copyright 013 by D.B. Rowe 1 Agenda: Recap Chapter 7. - 7.3 Lecture Chapter 8.1-8. Review Chapter 6. Problem Solving
More informationECE 340 Probabilistic Methods in Engineering M/W 3-4:15. Lecture 10: Continuous RV Families. Prof. Vince Calhoun
ECE 340 Probabilistic Methods in Engineering M/W 3-4:15 Lecture 10: Continuous RV Families Prof. Vince Calhoun 1 Reading This class: Section 4.4-4.5 Next class: Section 4.6-4.7 2 Homework 3.9, 3.49, 4.5,
More informationPoint Estimation. Some General Concepts of Point Estimation. Example. Estimator quality
Point Estimation Some General Concepts of Point Estimation Statistical inference = conclusions about parameters Parameters == population characteristics A point estimate of a parameter is a value (based
More informationABILITY OF VALUE AT RISK TO ESTIMATE THE RISK: HISTORICAL SIMULATION APPROACH
ABILITY OF VALUE AT RISK TO ESTIMATE THE RISK: HISTORICAL SIMULATION APPROACH Dumitru Cristian Oanea, PhD Candidate, Bucharest University of Economic Studies Abstract: Each time an investor is investing
More informationOptimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing
Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing Prof. Chuan-Ju Wang Department of Computer Science University of Taipei Joint work with Prof. Ming-Yang Kao March 28, 2014
More informationThe Economic and Social BOOTSTRAPPING Review, Vol. 31, No. THE 4, R/S October, STATISTIC 2000, pp
The Economic and Social BOOTSTRAPPING Review, Vol. 31, No. THE 4, R/S October, STATISTIC 2000, pp. 351-359 351 Bootstrapping the Small Sample Critical Values of the Rescaled Range Statistic* MARWAN IZZELDIN
More informationAN INFORMATION-BASED APPROACH TO CREDIT-RISK MODELLING. by Matteo L. Bedini Universitè de Bretagne Occidentale
AN INFORMATION-BASED APPROACH TO CREDIT-RISK MODELLING by Matteo L. Bedini Universitè de Bretagne Occidentale Matteo.Bedini@univ-brest.fr Agenda Credit Risk The Information-based Approach Defaultable Discount
More informationSimulation of probability distributions commonly used in hydrological frequency analysis
HYDROLOGICAL PROCESSES Hydrol. Process. 2, 5 6 (27) Published online May 26 in Wiley InterScience (www.interscience.wiley.com) DOI: 2/hyp.676 Simulation of probability distributions commonly used in hydrological
More informationEstimation of Value at Risk and ruin probability for diffusion processes with jumps
Estimation of Value at Risk and ruin probability for diffusion processes with jumps Begoña Fernández Universidad Nacional Autónoma de México joint work with Laurent Denis and Ana Meda PASI, May 21 Begoña
More informationHomework Assignments
Homework Assignments Week 1 (p. 57) #4.1, 4., 4.3 Week (pp 58 6) #4.5, 4.6, 4.8(a), 4.13, 4.0, 4.6(b), 4.8, 4.31, 4.34 Week 3 (pp 15 19) #1.9, 1.1, 1.13, 1.15, 1.18 (pp 9 31) #.,.6,.9 Week 4 (pp 36 37)
More informationProcess capability estimation for non normal quality characteristics: A comparison of Clements, Burr and Box Cox Methods
ANZIAM J. 49 (EMAC2007) pp.c642 C665, 2008 C642 Process capability estimation for non normal quality characteristics: A comparison of Clements, Burr and Box Cox Methods S. Ahmad 1 M. Abdollahian 2 P. Zeephongsekul
More informationBayesian Inference for Volatility of Stock Prices
Journal of Modern Applied Statistical Methods Volume 3 Issue Article 9-04 Bayesian Inference for Volatility of Stock Prices Juliet G. D'Cunha Mangalore University, Mangalagangorthri, Karnataka, India,
More informationFinancial Econometrics Jeffrey R. Russell. Midterm 2014 Suggested Solutions. TA: B. B. Deng
Financial Econometrics Jeffrey R. Russell Midterm 2014 Suggested Solutions TA: B. B. Deng Unless otherwise stated, e t is iid N(0,s 2 ) 1. (12 points) Consider the three series y1, y2, y3, and y4. Match
More informationModelling catastrophic risk in international equity markets: An extreme value approach. JOHN COTTER University College Dublin
Modelling catastrophic risk in international equity markets: An extreme value approach JOHN COTTER University College Dublin Abstract: This letter uses the Block Maxima Extreme Value approach to quantify
More informationContinous time models and realized variance: Simulations
Continous time models and realized variance: Simulations Asger Lunde Professor Department of Economics and Business Aarhus University September 26, 2016 Continuous-time Stochastic Process: SDEs Building
More informationCopyright 2005 Pearson Education, Inc. Slide 6-1
Copyright 2005 Pearson Education, Inc. Slide 6-1 Chapter 6 Copyright 2005 Pearson Education, Inc. Measures of Center in a Distribution 6-A The mean is what we most commonly call the average value. It is
More informationMEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL
MEASURING PORTFOLIO RISKS USING CONDITIONAL COPULA-AR-GARCH MODEL Isariya Suttakulpiboon MSc in Risk Management and Insurance Georgia State University, 30303 Atlanta, Georgia Email: suttakul.i@gmail.com,
More information2.1 Mathematical Basis: Risk-Neutral Pricing
Chapter Monte-Carlo Simulation.1 Mathematical Basis: Risk-Neutral Pricing Suppose that F T is the payoff at T for a European-type derivative f. Then the price at times t before T is given by f t = e r(t
More informationA New Test for Correlation on Bivariate Nonnormal Distributions
Journal of Modern Applied Statistical Methods Volume 5 Issue Article 8 --06 A New Test for Correlation on Bivariate Nonnormal Distributions Ping Wang Great Basin College, ping.wang@gbcnv.edu Ping Sa University
More informationOverview/Outline. Moving beyond raw data. PSY 464 Advanced Experimental Design. Describing and Exploring Data The Normal Distribution
PSY 464 Advanced Experimental Design Describing and Exploring Data The Normal Distribution 1 Overview/Outline Questions-problems? Exploring/Describing data Organizing/summarizing data Graphical presentations
More informationValue at Risk Ch.12. PAK Study Manual
Value at Risk Ch.12 Related Learning Objectives 3a) Apply and construct risk metrics to quantify major types of risk exposure such as market risk, credit risk, liquidity risk, regulatory risk etc., and
More informationOption Pricing Formula for Fuzzy Financial Market
Journal of Uncertain Systems Vol.2, No., pp.7-2, 28 Online at: www.jus.org.uk Option Pricing Formula for Fuzzy Financial Market Zhongfeng Qin, Xiang Li Department of Mathematical Sciences Tsinghua University,
More informationMAINTAINABILITY DATA DECISION METHODOLOGY (MDDM)
TECHNICAL REPORT NO. TR-2011-19 MAINTAINABILITY DATA DECISION METHODOLOGY (MDDM) June 2011 APPROVED FOR PUBLIC RELEASE; DISTRIBUTION IS UNLIMITED U.S. ARMY MATERIEL SYSTEMS ANALYSIS ACTIVITY ABERDEEN PROVING
More informationMONTE CARLO EXTENSIONS
MONTE CARLO EXTENSIONS School of Mathematics 2013 OUTLINE 1 REVIEW OUTLINE 1 REVIEW 2 EXTENSION TO MONTE CARLO OUTLINE 1 REVIEW 2 EXTENSION TO MONTE CARLO 3 SUMMARY MONTE CARLO SO FAR... Simple to program
More informationInferences on Correlation Coefficients of Bivariate Log-normal Distributions
Inferences on Correlation Coefficients of Bivariate Log-normal Distributions Guoyi Zhang 1 and Zhongxue Chen 2 Abstract This article considers inference on correlation coefficients of bivariate log-normal
More informationSAMPLE STANDARD DEVIATION(s) CHART UNDER THE ASSUMPTION OF MODERATENESS AND ITS PERFORMANCE ANALYSIS
Science SAMPLE STANDARD DEVIATION(s) CHART UNDER THE ASSUMPTION OF MODERATENESS AND ITS PERFORMANCE ANALYSIS Kalpesh S Tailor * * Assistant Professor, Department of Statistics, M K Bhavnagar University,
More informationRichardson Extrapolation Techniques for the Pricing of American-style Options
Richardson Extrapolation Techniques for the Pricing of American-style Options June 1, 2005 Abstract Richardson Extrapolation Techniques for the Pricing of American-style Options In this paper we re-examine
More informationChapter 7: Point Estimation and Sampling Distributions
Chapter 7: Point Estimation and Sampling Distributions Seungchul Baek Department of Statistics, University of South Carolina STAT 509: Statistics for Engineers 1 / 20 Motivation In chapter 3, we learned
More informationSome developments about a new nonparametric test based on Gini s mean difference
Some developments about a new nonparametric test based on Gini s mean difference Claudio Giovanni Borroni and Manuela Cazzaro Dipartimento di Metodi Quantitativi per le Scienze Economiche ed Aziendali
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 informationEX-POST VERIFICATION OF PREDICTION MODELS OF WAGE DISTRIBUTIONS
EX-POST VERIFICATION OF PREDICTION MODELS OF WAGE DISTRIBUTIONS LUBOŠ MAREK, MICHAL VRABEC University of Economics, Prague, Faculty of Informatics and Statistics, Department of Statistics and Probability,
More informationData Distributions and Normality
Data Distributions and Normality Definition (Non)Parametric Parametric statistics assume that data come from a normal distribution, and make inferences about parameters of that distribution. These statistical
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 informationSt. Xavier s College Autonomous Mumbai T.Y.B.A. Syllabus For 5 th Semester Courses in Statistics (June 2016 onwards)
St. Xavier s College Autonomous Mumbai T.Y.B.A. Syllabus For 5 th Semester Courses in Statistics (June 2016 onwards) Contents: Theory Syllabus for Courses: A.STA.5.01 Probability & Sampling Distributions
More informationMath 416/516: Stochastic Simulation
Math 416/516: Stochastic Simulation Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 13 Haijun Li Math 416/516: Stochastic Simulation Week 13 1 / 28 Outline 1 Simulation
More informationEdgeworth Binomial Trees
Mark Rubinstein Paul Stephens Professor of Applied Investment Analysis University of California, Berkeley a version published in the Journal of Derivatives (Spring 1998) Abstract This paper develops a
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 New Hybrid Estimation Method for the Generalized Pareto Distribution
A New Hybrid Estimation Method for the Generalized Pareto Distribution Chunlin Wang Department of Mathematics and Statistics University of Calgary May 18, 2011 A New Hybrid Estimation Method for the GPD
More informationA lower bound on seller revenue in single buyer monopoly auctions
A lower bound on seller revenue in single buyer monopoly auctions Omer Tamuz October 7, 213 Abstract We consider a monopoly seller who optimally auctions a single object to a single potential buyer, with
More informationNormal Probability Distributions
Normal Probability Distributions Properties of Normal Distributions The most important probability distribution in statistics is the normal distribution. Normal curve A normal distribution is a continuous
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 informationOn modelling of electricity spot price
, Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction
More informationProblem Set 4 Answer Key
Economics 31 Menzie D. Chinn Fall 4 Social Sciences 7418 University of Wisconsin-Madison Problem Set 4 Answer Key This problem set is due in lecture on Wednesday, December 1st. No late problem sets will
More informationReturn dynamics of index-linked bond portfolios
Return dynamics of index-linked bond portfolios Matti Koivu Teemu Pennanen June 19, 2013 Abstract Bond returns are known to exhibit mean reversion, autocorrelation and other dynamic properties that differentiate
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Generating Random Variables and Stochastic Processes Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationStrategies for Improving the Efficiency of Monte-Carlo Methods
Strategies for Improving the Efficiency of Monte-Carlo Methods Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction The Monte-Carlo method is a useful
More informationarxiv: v1 [q-fin.rm] 13 Dec 2016
arxiv:1612.04126v1 [q-fin.rm] 13 Dec 2016 The hierarchical generalized linear model and the bootstrap estimator of the error of prediction of loss reserves in a non-life insurance company Alicja Wolny-Dominiak
More informationMonte Carlo Simulation (Random Number Generation)
Monte Carlo Simulation (Random Number Generation) Revised: 10/11/2017 Summary... 1 Data Input... 1 Analysis Options... 6 Summary Statistics... 6 Box-and-Whisker Plots... 7 Percentiles... 9 Quantile Plots...
More informationTechnology Support Center Issue
United States Office of Office of Solid EPA/600/R-02/084 Environmental Protection Research and Waste and October 2002 Agency Development Emergency Response Technology Support Center Issue Estimation of
More informationEconomics 483. Midterm Exam. 1. Consider the following monthly data for Microsoft stock over the period December 1995 through December 1996:
University of Washington Summer Department of Economics Eric Zivot Economics 3 Midterm Exam This is a closed book and closed note exam. However, you are allowed one page of handwritten notes. Answer all
More informationMonte Carlo Simulation (General Simulation Models)
Monte Carlo Simulation (General Simulation Models) Revised: 10/11/2017 Summary... 1 Example #1... 1 Example #2... 10 Summary Monte Carlo simulation is used to estimate the distribution of variables when
More informationMonte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMS091)
Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMS091) Magnus Wiktorsson Centre for Mathematical Sciences Lund University, Sweden Lecture 3 Importance sampling January 27, 2015 M. Wiktorsson
More informationFractional Liu Process and Applications to Finance
Fractional Liu Process and Applications to Finance Zhongfeng Qin, Xin Gao Department of Mathematical Sciences, Tsinghua University, Beijing 84, China qzf5@mails.tsinghua.edu.cn, gao-xin@mails.tsinghua.edu.cn
More informationPoint Estimation. Stat 4570/5570 Material from Devore s book (Ed 8), and Cengage
6 Point Estimation Stat 4570/5570 Material from Devore s book (Ed 8), and Cengage Point Estimation Statistical inference: directed toward conclusions about one or more parameters. We will use the generic
More informationHypothesis Tests: One Sample Mean Cal State Northridge Ψ320 Andrew Ainsworth PhD
Hypothesis Tests: One Sample Mean Cal State Northridge Ψ320 Andrew Ainsworth PhD MAJOR POINTS Sampling distribution of the mean revisited Testing hypotheses: sigma known An example Testing hypotheses:
More informationPricing CDOs with the Fourier Transform Method. Chien-Han Tseng Department of Finance National Taiwan University
Pricing CDOs with the Fourier Transform Method Chien-Han Tseng Department of Finance National Taiwan University Contents Introduction. Introduction. Organization of This Thesis Literature Review. The Merton
More information