THE NEGATIVE BINOMIAL-ERLANG DISTRIBUTION WITH APPLICATIONS. Kasetsart University Chatuchak, Bangkok, 10900, THAILAND
|
|
- Basil Phelps
- 6 years ago
- Views:
Transcription
1 International Journal of Pure and Applied Mathematis Volume 92 No , ISSN: (printed version; ISSN: (on-line version url: doi: PAipam.eu THE NEGATIVE BINOMIAL-ERLANG DISTRIBUTION WITH APPLICATIONS Siriporn Kongrod 1, Winai Bodhisuwan 2, Prasit Payakkapong 3 1,2,3 Department of Statistis Kasetsart University Chatuhak, Bangkok, 10900, THAILAND Abstrat: This paper introdues a new three-parameter of the mied negative binomial distribution whih is alled the negative binomial-erlang distribution. This distribution obtained by miing the negative binomial distribution with the Erlang distribution. The negative binomial-erlang distribution an be used to desribe ount data with a large number of zeros. The negative binomialeponential is presented as speial ases of the negative binomial-erlang distribution. In addition, we present some properties of the negative binomial-erlang distribution inluding fatorial moments, mean, variane, skewness and kurtosis. The parameter estimation for negative binomial-erlang distribution by the maimum likelihood estimation are provided. Appliations of the the negative binomial-erlang distribution are arried out two real ount data sets. The result shown that the negative binomial-erlang distribution is better than fit when ompared the Poisson and negative binomial distributions. AMS Subet Classifiation: 60E05 Key Words: mied distribution, ount data, Erlang distribution, the negative binomial-erlang distribution, overdispersion Reeived: Deember 4, 2013 Correspondene author 2014 Aademi Publiations, Ltd. url:
2 390 S. Kongrod, W. Bodhisuwan, P. Payakkapong 1. Introdution The Poisson distribution is a disrete probability distribution for the ounts of events that our randomly in a given interval of time, sine data in multiple researh fields often ahieve the Poisson distribution. However, ount data often ehibit overdispersion, with a variane larger than the mean. The negative binomial (NB distribution is employed as a funtional from that whih relaes the overdispersion restrition of the Poisson distribution. The NB distribution is a miture of Poisson distribution by miing the Poisson distribution and gamma distribution. The NB distribution has beome inreasingly popular as a more fleible alternative to ount data with overdispersion. The NB distribution was first introdued by Pasal [2]. Montmort [13] applied the NB distribution to represent the number of times that a oin should be flipped to obtain a ertain number of heads. Student [16] used the NB distribution as an alternative to the Poisson distribution. Eggenberger and Polya [5] has been shown to be the limiting form of the NB distribution. Ever sine, there has been an inreasing number of appliations of the NB distribution have been developed in a parallel way. For instant, the Poisson-inverse Gaussian distribution [15], [9], the negative binomial-inverse Gaussian distributions [4], [15] and the negative binomial-beta Eponential distribution [3]. In this paper we introdue a new mied negative binomial distribution, whih is alled the negative binomial-erlang (NB-EL distribution. The NB- EL distribution an be used to desribe ount data distribution with a large number of zeros in Poisson distribution. The Erlang (EL distribution was introdued by Agner Erlang [1] was the first author to etend the eponential distribution with his method of stages. He defined a non-negative random variable as the time taken to move through a fied number of stages, spending an eponential amount of time with a fied rate in eah one. The NB-EL distribution is obtained by miing the NB distribution whose probability of suess parameter p follows the EL distribution. In other words, if the probability of suess parameter p of the NB distribution has the EL distribution, then the resulting distribution is referred to as the NB-EL with parameters r, k and. The negative binomial-eponential has three-parameter is speial ase of the NB-EL distribution. In partiular, various strutural properties of the new distribution are derived, inluding epansions for its fatorial moments, mean, variane, skewness and kurtosis.
3 THE NEGATIVE BINOMIAL-ERLANG The Negative Binomial Distribution The negative binomial distribution is a natural and more fleible etension of Poisson distribution and allows for overdispersion relative to the Poisson distribution. If X is the NB distribution with parameter r > 0 and 0 < p < 1 then we an write the pmf as ( r + 1 f(;r,p = p r (1 p ; = 0,1,2,... (1 Note that p is a probability of suess, the eperiment is repeated as many times as required to obtain r suesses. The first two moments about zero and the fatorial moment of order m of the NB distribution are given respetively E(X = r(1 p, p Var = r(1 p p 2, µ [m] (X = Γ(r +m] Γ(r (1 pm p m, m = 1,2,... (2 The log likelihood funtion of the NB distribution is given by L(r,p = n [( r +i 1 i p r (1 p i ], then, we an write the log likelihood funtion of the NB(r,p as logl(r,p = n ( r +i 1 log +nrlog(p i ( n log(1 p i. (3
4 392 S. Kongrod, W. Bodhisuwan, P. Payakkapong 3. The Erlang Distribution The EL distribution was introdued by Agner Erlang [1]. It is a ontinuous probability distribution with wide appliability mainly due to its relative to the eponential distributions. The probability density funtion (pdf of the EL distribution is given by g( = k k 1 e, (4 (k 1! for > 0 and k, > 0. Note that parameter k is alled the shape parameter and the parameter is alled the rate parameter. The mean of the EL distribution is given as The variane of the EL distribution is then given as E(X = k. (5 Var(X = k 2. (6 The moment generating funtion of the EL distribution is defined by M X (t = ( 1 t k for t > 0. (7 4. The Negative Binomial-Erlang Distribution Definition 1. Let X be a random variable of the NB-EL (r,k, distribution, when the NB distribution have parameters r > 0 and p = ep( λ, where λ is distributed as the EL distribution with positive parameters k and, i.e., X λ NB(r,p = ep( λ and λ EL(k,. Theorem 2. Let X NB-EL(r,k,. The pmf of X is given by ( r+ 1 f(;r,k, = for > 0 and k, > 0. ( ( ( 1 +(r + k, (8
5 THE NEGATIVE BINOMIAL-ERLANG Proof. If X λ NB(r,p = ep( λ in Eq. 1 and λ EL(k, in Eq. 4, then the pmf of X an be obtained by f(;r,k, = where, f 1 ( λ is defined by ( r + 1 f 1 ( λ = ( r + 1 = 0 f 1 ( λg(λ;k,dλ (9 e λ (1 e λ By substituting Eq. 10 into Eq. 9, we obtain ( ( ( r + 1 f( λ = ( 1 0 ( ( r + 1 = ( ( 1 e λ(r+. (10 e λ(r+ f 2 (λ;α,βdλ ( 1 (M λ ( (r +. (11 Substituting the moment generating funtion of EL distribution Eq. 7 into Eq. 11, the pmf of NB-EL (r,k, is finally given as ( ( ( r+ 1 f(;r,k, = ( 1 k. +(r + Many well known distributions are subsumed by the NB-EL distribution. Net, we display graphs of the pmf of the NB distribution with various values of parameters are shown in Figure 1. Corollary 3. Ifk = 1thentheNB-EL distributionreduestothenegative binomial-eponential (NB-E distribution with pmf given by ( ( ( r + 1 +(r + f(;r, = ( 1, (12 where = 0,1,2,... for r and > 0. Proof. If X λ NB(r,p = ep λ and λ EL(k = 1,, then the pmf of X is ( ( ( r + 1 +(r + f(;r, = ( 1
6 394 S. Kongrod, W. Bodhisuwan, P. Payakkapong (a (b Probability Probability ( (d Probability Probability (e (f Probability Probability Figure 1: (a-f: The probability mass funtion of the NB-EL distribution of some values of parameters: (a r=0.5, k=0.5, =0.5, (b r=0.5, k=1, =0.5, ( r=0.5, k=2, =0.5, (dr=5, k=1, =1.5, (e r=5, k=1, =2 and (f r=5, k=2, =5
7 THE NEGATIVE BINOMIAL-ERLANG ( ( ( r + 1 +(r + = ( 1. From Corollary 3, we find the negative binomial-eponential distribution displayed in Eq. 12, whih introdued by Panger and Willmot [17]. 5. Properties of the NB-EL Distribution In this setion, we have studied the properties of the NB-EL distribution, whih inludes the fatorial moments, mean, variane, skewness and kurtosis are given as follow Theorem 4. If X NB-EL(r, k,, then the fatorial moment of order m of X is given by µ [m] (X = Γ(r +m Γ(r m where = 0,1,2,... for r,k and > 0. ( ( m ( 1 k, (13 (m Proof. If X λ NB(r,p = ep( λ in Eq. 1 and λ EL(k, in Eq. 4, then the fatorial moment of order m of X an be obtained by µ [k] (X = E λ [µ k (X λ]. The fatorial moment of order m of the NB distribution in Eq. 2, µ [m] (X beomes ( Γ(r +m (1 e λ m Γ(r +m µ [m] (X = E λ Γ(r e λm = E λ (e λ 1 m. Γ(r Using a binomial epansion of (e λ 1 m, then shows that µ [m] (X an be written as µ [m] (X = = Γ(r +m Γ(r Γ(r +m Γ(r m ( m ( 1 E λ (e λ(m m ( m ( 1 M λ (m.
8 396 S. Kongrod, W. Bodhisuwan, P. Payakkapong From the moment generating funtion of the EL distribution in Eq. 7, we have finally that µ [m] (X an be written as µ [m] (X = Γ(r +m Γ(r m ( ( m ( 1 k. (m From the fatorial moments of the NB-EL distribution, it is straightforward to redue the first four moments given in Eq Eq. 17, variane in Eq. 18, skewness in Eq. 19 and kurtosis in Eq. 20. E(X = r(ϕ 1 1, (14 E(X 2 = (r 2 +rϕ 2 (2r 2 +rϕ 1 +r 2, (15 E(X 3 = (r 3 +3r 2 +2rϕ 3 (3r 3 +6r 2 +3rϕ 2 +(3r 3 +3r 2 +rϕ 1 r 3, (16 E(X 4 = (r 4 +6r 3 +11r 2 +6rϕ 4 (4r 4 +18r 3 +26r 2 +12rϕ 3 +(6r 4 +18r 3 +19r 2 +7rϕ 2 (4r 4 +6r 3 +4r 2 +rϕ 1 +r 4, (17 Var(X = E(X 2 (E(X 2 = (r 2 +rϕ 2 rϕ 1 (1+rϕ 1, (18 Skewness(X = [ E(X 3 3E(X 2 E(X+2[E(X] 3 / ] σ 3 [ = (r 3 +3r 2 +2rϕ 3 (3r 2 +3rϕ 2 +rϕ 1 +3r 2 (ϕ 2 2 (3r 3 +3r 2 ϕ 1 ϕ 2 +2r 3 (ϕ 2 3] /σ 3, (19 Kurtosis(X = = [ E(X 4 4E(X 3 E(X+6E(X 2 (E(X 2 3[E(X] 4 / ] σ 4 [ (r 4 +6r 3 +11r 2 +6rϕ 4 (6r 3 +18r 2 +12rϕ 3 3r 4 (ϕ 1 4 +(7r 2 +7rϕ 2 4rϕ 1 6r 3 (ϕ 1 3 +(12r 3 +12r 2 ϕ 1 ϕ 2 (4r 4 +12r 3 +8r 2 ϕ 1 ϕ 3 +(6r 4 +6r 3 (ϕ 1 2 ϕ 2 ]/σ 4, (20
9 THE NEGATIVE BINOMIAL-ERLANG where ( i k ϕ i =. 6. Parameters Estimation The estimation of parameters for NB-EL distribution via the maimum-likelihood estimation (MLE method proedure will be disussed. The likelihood funtion of the NB-EL(r,k, is given by n i L(r,k, = ( r+i 1 i ( ( i ( 1 with orresponding log-likelihood funtion n ( r+i 1 L = logl(r,k, = log + n i ( i log ( 1 ( i +(r + k. (21 k. (22 +(r + The first order onditions for finding the optimal values of the parameters obtained by differentiating Eq. 22 with respet to r,k and give rise to the following differential equations and L r = + L n k = L r n i 1 k=0 n i = log(r +k k i i ( i [ ( 1 i n i ( i ( 1 k [ ( 1 ( i +(r+ ( i ( 1 i [ (+(r+ ] k+1 k +(r+ ] klog [ +(r+ +(r+ ] k ( i [ ( 1 k k 1 (r+ [ ] k ( 1 ( i (+r+ k+1 ] +(r+, (23 ], (24. (25
10 398 S. Kongrod, W. Bodhisuwan, P. Payakkapong We estimated parameters of the NB-EL distribution with equating Eq Eq. 25 to zero, the MLE solutions of ˆr,ˆk and ĉ an be obtained by solving the resulting equations simultaneously using a numerial proedure, the Newton- Raphson method. 7. Illustrative Eamples The NB-EL distribution is applied on two real ount data sets. The details are presented in the following eamples. Eample 1. We use a real data set is numbers of inured from the aident on maor road in Bangkok in 2007 [3]. The data was olleted by Department of Highways, Ministry of Transport, Thailand. We use the real data are fitted by the Poisson, NB and NB-EL distributions (Table 1, whih show the observed and epeted frequenies, grouped in lasses of epeted frequeny greater than five for the hi-square goodness of fit test. Based on the p-value, the maimum likelihood method provides very poor fit for the Poisson and the NB distribution. We an see that NB-EL distribution provides the highest p-value of fitting for this data set. Table 1: Observed and epeted frequenies for the aident data Number of Observed Fitting distribution inured frequeny Poisson NB NB-EL } Estimated parameters ˆλ= ˆr = ˆr = ˆp = ˆk= ĉ= Chi-squares Degrees of freedom p-value < Eample 2. We used a data set whih was obtained from Klugman et
11 THE NEGATIVE BINOMIAL-ERLANG al. [15], provides the data set of 9,461 automobile insurane poliies where by the number of aidents of eah poliy has been reorded. We apply this data set to fit with the Poisson, NB and NB-EL distributions. Obtained results (Table 2, shows that the observed and epeted frequenies, grouped in lasses of epeted frequeny greater than five for the hi-square goodness of fit test. The maimum likelihood method provides very poor fit for the Poisson and the NB distribution. Based on p-value we found that the NB-EL distribution seen to be best fit among there three distributions, with p-value Table 2: Observed and epeted frequenies for the aident data Number of Observed Fitting distribution laims frequeny Poisson NB NB-EL Estimated parameters ˆλ= ˆr = ˆr = ˆp = ˆk= ĉ= Chi-squares Degrees of freedom p-value < Conlusion In this paper, we introdue a new three-parameter negative binomial-erlang distribution, NB-EL(r, k,. This distribution obtained by miing the NB with the EL distribution (when the NB distribution have parameters r > 0 and p = ep( λ, where λ is distributed as the EL distribution with positive parameters k and. We showed that the negative binomial-eponential distributions isaspeial ase ofthenb-el distribution. Inaddition, themoments of thenb- EL distribution whih inludes the fatorial moments, mean, variane, skewness
12 400 S. Kongrod, W. Bodhisuwan, P. Payakkapong and kurtosis are derived. Moreover, the parameter estimation of the NB-EL using the maimum likelihood estimation are developed. We inlude an appliation of the NB-EL distribution to fit two real data sets. We found that this data set is best fit with NB-EL distribution among there three distributions. We hope that NB-EL distribution may attrat wider appliations in analyzing ount data. Aknowledgments We would like to thank Dr.Chookait Pudprommarat, who has given many useful omments and valuable suggestions. This work was supported by researh funding from the Researh Professional Development Proet under the Siene Ahievement Sholarship of Thailand (SAST. Referenes [1] A.K. Erlang, Solution of some problems in the theory of probabilities of signifiane in automati telephone ehanges, Elektrotkeknikeren, 13(1917, 513. [2] Blaise Pasal, Varia opera mathematia D.Petri de Fermat, Tolossae, (1679. [3] C. Pudprommarat, W. Bodhisuwan, A new mied negative binomial distribution, Journal of Applied Sienes, 12(2012, , doi: /as [4] E. Gomez-Deniz, J.M. Sarabia and E. Calderin-Oeda, Univariate and multivariate versions of the negative binomial-inverse Gaussian distributions with appliations, Insurane: Mathematis and Eonomis, 4(2008, doi: /J.insmatheo [5] F. Eggenberger and G.Polya, Uber die Statistik verketetter Votgange, Zeitshrift fur Angewandte Mathematik and Mehanik, 1 (1923, , doi: /zamm [6] G.E. Willmot, The Poisson-inverse Gaussian distribution as an alternative to the negative binomial, Sandinavian Atuar, 2 (1987, , doi: /
13 THE NEGATIVE BINOMIAL-ERLANG [7] H. Ranger and E. Willmot, Finite sum evaluation of the negative binomial eponential model, Astin Bulletin, 12(1981, [8] J. Lemaire, How to define a bonus-malus system with an eponential utility funtion. ASTIN Bull. 10(1979, [9] L. Tremblay, Using the Poisson inverse Gaussian in bonus-malus systems, ASTIN Bull, 22(1992, [10] L.J. Simon, Fitting negative binomial distributions by the method of maimum likelihood, Pro. Casual. Atuar. So., 48(1961, [11] M. Greenwood and G.U. Yule, An inquiry into the Nature of frequeny distributions representative of multiple happenings with partiular referene to the ourrene of multiple attaks of disease or of repeated aidents, Journal of the Royal Statistial Soiety, 83(1920, [12] M.D. Bressoud, A radial approah to real analysis, Cambridge University Press, Cambridge (2005. [13] P.R. de Montmort, Essay d Analyse sur les Jeu de Hazard, Seonde Edition, Revue Augmentee de plusieurs Lettre, Quillau, Paris (1980. [14] S. Meng,Y. Wei and G.A. Whitmore, Aounting for individual overdispersion in a bonus-malus system, ASTIN Bull.(1999. [15] S.A. Klugman. H.H. Paner and G.E.Willmot, Loss Models: From Data to Deision, John Wiley and Sons USA, (2008. [16] W.S. Gosset Student, On the error of ounting with a Haemaytometer, Biometrika, 5(1907, [17] H.H. Paner and G.E. Willmot, Finite sum eqvluation of the negative binomail-eponential model, Astin Bulletin, 12(1981,
14 402
Bonus-Malus System with the Claim Frequency Distribution is Geometric and the Severity Distribution is Truncated Weibull
IOP Conferene Series: Earth and Environmental Siene PAPER OPEN ACCESS Bonus-Malus System with the Claim Frequeny Distribution is Geometri and the Severity Distribution is Trunated Weibull To ite this artile:
More informationUsing the Average of the Extreme Values of a Triangular Distribution for a Transformation, and Its Approximant via the Continuous Uniform Distribution
British Journal of Mathematis & Computer Siene 4(4):., 014 ISSN: 31-0851 SCIENCEDOMAIN international www.sienedomain.org Using the Average of the Extreme Values of a Triangular Distribution for a Transformation,
More informationBinomial Mixture of Erlang Distribution
International Research Journal of Scientific Findings. Vol. 1 (6), Pp. 209-215, September, 2014. www.ladderpublishers.org Research Paper Binomial Mixture of Erlang Distribution * Kareema Abed Al-Kadim
More informationLimiting Limited Liability
Limiting Limited Liability Giuseppe Dari-Mattiai Amsterdam Center for Law & Eonomis Working Paper No. 2005-05 This paper an be downloaded without harge from the Soial Siene Researh Network Eletroni Paper
More informationEconomics 325 Intermediate Macroeconomic Analysis Practice Problem Set 1 Suggested Solutions Professor Sanjay Chugh Spring 2011
Department of Eonomis Universit of Marland Eonomis 35 Intermediate Maroeonomi Analsis Pratie Problem Set Suggested Solutions Professor Sanja Chugh Spring 0. Partial Derivatives. For eah of the following
More informationEconomics 2202 (Section 05) Macroeconomic Theory Practice Problem Set 3 Suggested Solutions Professor Sanjay Chugh Fall 2014
Department of Eonomis Boston College Eonomis 2202 (Setion 05) Maroeonomi Theory Pratie Problem Set 3 Suggested Solutions Professor Sanjay Chugh Fall 2014 1. Interation of Consumption Tax and Wage Tax.
More informationOptional Section: Continuous Probability Distributions
6 Optional Setion: Continuous Probability Distributions 6.5 The Normal Approximation to the Binomial Distribution For eah retangle, the width is 1 and the height is equal to the probability assoiated with
More informationRetirement Benefits Schemes (Miscellaneous Amendments) RETIREMENT BENEFITS SCHEMES (MISCELLANEOUS AMENDMENTS) REGULATIONS 2014
Retirement Benefits Shemes (Misellaneous Amendments) Index RETIREMENT BENEFITS SCHEMES (MISCELLANEOUS AMENDMENTS) REGULATIONS 2014 Index Regulation Page 1 Title... 3 2 Commenement... 3 3 Amendment of the
More informationClass Notes: Week 6. Multinomial Outcomes
Ronald Hek Class Notes: Week 6 1 Class Notes: Week 6 Multinomial Outomes For the next ouple of weeks or so, we will look at models where there are more than two ategories of outomes. Multinomial logisti
More informationPage 80. where C) refers to estimation cell (defined by industry and, for selected industries, region)
Nonresponse Adjustment in the Current Statistis Survey 1 Kennon R. Copeland U.S. Bureau of Labor Statistis 2 Massahusetts Avenue, N.E. Washington, DC 20212 (Copeland.Kennon@bls.gov) I. Introdution The
More informationContinuous Distributions
Quantitative Methods 2013 Continuous Distributions 1 The most important probability distribution in statistics is the normal distribution. Carl Friedrich Gauss (1777 1855) Normal curve A normal distribution
More informationKurtosis Statistics with Reference to Power Function Distribution
ISSN 68-80 Journal o Statistis Volume, 06. pp. -0 Abstrat Kurtosis Statistis with Reerene to Power Funtion Distribution Azaz Ahmad and Ahmed Saeed Akhter Pearson statistis o skewness and kurtosis gave
More informationDecision-making Method for Low-rent Housing Construction Investment. Wei Zhang*, Liwen You
5th International Conferene on Civil Enineerin and Transportation (ICCET 5) Deision-makin Method for Low-rent Housin Constrution Investment Wei Zhan*, Liwen You University of Siene and Tehnoloy Liaonin,
More informationVersion A. Problem 1. Let X be the continuous random variable defined by the following pdf: 1 x/2 when 0 x 2, f(x) = 0 otherwise.
Math 224 Q Exam 3A Fall 217 Tues Dec 12 Version A Problem 1. Let X be the continuous random variable defined by the following pdf: { 1 x/2 when x 2, f(x) otherwise. (a) Compute the mean µ E[X]. E[X] x
More informationSAMPLE CHAPTERS UNESCO EOLSS INVESTMENT MODELS. Ulrich Rieder University of Ulm, Germany
INVESMEN MODELS Ulrih Rieder University of Ulm, Germany Keywords: meanvariane ortfolio seletion, Markowitz model, minimum variane ortfolio, twofund searation, HARAutility, BlakSholes model, stohasti dynami
More informationSimulation Wrap-up, Statistics COS 323
Simulation Wrap-up, Statistics COS 323 Today Simulation Re-cap Statistics Variance and confidence intervals for simulations Simulation wrap-up FYI: No class or office hours Thursday Simulation wrap-up
More informationThe Impact of Capacity Costs on Bidding Strategies in Procurement Auctions
Review of Aounting Studies, 4, 5 13 (1999) 1999 Kluwer Aademi Publishers, Boston. Manufatured in The Netherlands. The Impat of Capaity Costs on Bidding Strategies in Prourement Autions JÖRG BUDDE University
More informationAn EOQ Model with Parabolic Demand Rate and Time Varying Selling Price
Annals of Pure and Applied Mathematis Vol.,.,, 3-43 ISSN: 79-87X (P),79-888(online) Published on 5 September www.researhmathsi.org Annals of An EOQ Model with Paraboli Demand Rate and ime Varying Selling
More informationmathematics Prey Predator Models with Variable Carrying Capacity Article Mariam K. A. Al-Moqbali, Nasser S. Al-Salti and Ibrahim M.
mathematis Artile Prey Predator Models with Variable Carrying Capaity Mariam K. A. Al-Moqbali, Nasser S. Al-Salti and Ibrahim M. Elmojtaba * Department of Mathematis and Statistis, College of Siene, Sultan
More informationFree Space Path Loss of UWB Communications
Free Spae Path Loss of UWB Communiations Pihaya Supanakoon 1, Sathit Aroonpraparat 1, Sathaporn Promwong 1 and Jun-ihi Takada 1 Department of Information Engineering, Faulty of Engineering, King Mongkut
More informationAvailability Analysis with Opportunistic Maintenance of a Two Component Deteriorating System
Analysis with Maintenane of a Two Component Deteriorating System Joel P. Varghese and Girish Kumar Abstrat This paper desribes the opportunisti maintenane model for availability analysis of two omponent
More informationValuation of Bermudan-DB-Underpin Option
Valuation of Bermudan-DB-Underpin Option Mary, Hardy 1, David, Saunders 1 and Xiaobai, Zhu 1 1 Department of Statistis and Atuarial Siene, University of Waterloo Marh 31, 2017 Abstrat The study of embedded
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 informationLecture 10: Point Estimation
Lecture 10: Point Estimation MSU-STT-351-Sum-17B (P. Vellaisamy: MSU-STT-351-Sum-17B) Probability & Statistics for Engineers 1 / 31 Basic Concepts of Point Estimation A point estimate of a parameter θ,
More informationi e AT 16 of 2008 INSURANCE ACT 2008
i e AT 16 of 2008 INSURANCE ACT 2008 Insurane At 2008 Index i e INSURANCE ACT 2008 Index Setion Page PART 1 REGULATORY OBJECTIVES 9 1 Regulatory objetives... 9 2 [Repealed]... 9 PART 2 ADMINISTRATION
More informationRandom Variables Handout. Xavier Vilà
Random Variables Handout Xavier Vilà Course 2004-2005 1 Discrete Random Variables. 1.1 Introduction 1.1.1 Definition of Random Variable A random variable X is a function that maps each possible outcome
More informationEDWARD W. (JED) FREES 1, GLENN MEYERS 2 AND A. DAVID CUMMINGS 2 ABSTRACT
DEPENDENT MULTI-PERIL RATEMAKING MODELS BY EDWARD W. (JED) FREES, GLENN MEYERS 2 AND A. DAVID CUMMINGS 2 ABSTRACT This paper onsiders insurane laims that are available by ause of loss, or peril. Using
More informationCHAPTER 9 BUDGETARY PLANNING SUMMARY OF QUESTIONS BY STUDY OBJECTIVES AND BLOOM S TAXONOMY. True-False Statements. Multiple Choice Questions
HTER 9 BUDGETARY PLANNING SUMMARY OF QUESTIONS BY STUDY OBJETIVES AND BLOOM S TAXONOMY Item SO BT Item SO BT Item SO BT Item SO BT 4 6 6 6 1 11. 11. 114. 11. 116. 117. 118. 119. 10. 11. 1. 1. 14. 1. 16.
More informationEfficient Pricing of European Options with Stochastic Interest Rate Using Fourier Transform Method
Amerian Journal of Applied Mathematis 2016; 4(4): 181-185 http://www.sienepublishinggroup.om/j/ajam doi: 10.11648/j.ajam.20160404.13 ISSN: 2330-0043 (Print); ISSN: 2330-006X (Online) Effiient Priing of
More informationPoint Estimation. Copyright Cengage Learning. All rights reserved.
6 Point Estimation Copyright Cengage Learning. All rights reserved. 6.2 Methods of Point Estimation Copyright Cengage Learning. All rights reserved. Methods of Point Estimation The definition of unbiasedness
More informationi e SD No.2017/0343 PAYMENT SERVICES (AMENDMENT) REGULATIONS 2017
i e SD No.2017/0343 PAYMENT SERVICES (AMENDMENT) REGULATIONS 2017 Payment Servies (Amendment) Regulations 2017 Index PAYMENT SERVICES (AMENDMENT) REGULATIONS 2017 Index Regulation Page 1 Title... 5 2
More informationChapter 4: Commonly Used Distributions. Statistics for Engineers and Scientists Fourth Edition William Navidi
Chapter 4: Commonly Used Distributions Statistics for Engineers and Scientists Fourth Edition William Navidi 2014 by Education. This is proprietary material solely for authorized instructor use. Not authorized
More informationProbability Theory and Simulation Methods. April 9th, Lecture 20: Special distributions
April 9th, 2018 Lecture 20: Special distributions Week 1 Chapter 1: Axioms of probability Week 2 Chapter 3: Conditional probability and independence Week 4 Chapters 4, 6: Random variables Week 9 Chapter
More informationME3620. Theory of Engineering Experimentation. Spring Chapter III. Random Variables and Probability Distributions.
ME3620 Theory of Engineering Experimentation Chapter III. Random Variables and Probability Distributions Chapter III 1 3.2 Random Variables In an experiment, a measurement is usually denoted by a variable
More informationGOVERNMENT GAZETTE REPUBLIC OF NAMIBIA
.. GOVERNMENT GAZETTE OF THE REPUBLIC OF NAMIBIA N$2.25 WINDHOEK -22 Deember 2001 CONTENTS Page GOVERNMENT NOTICE No. 249 Promulgation of Environment Investment Fund of Namibia At, 2001 (At No. 13 of 200
More information6. Genetics examples: Hardy-Weinberg Equilibrium
PBCB 206 (Fall 2006) Instructor: Fei Zou email: fzou@bios.unc.edu office: 3107D McGavran-Greenberg Hall Lecture 4 Topics for Lecture 4 1. Parametric models and estimating parameters from data 2. Method
More informationLicensing and Patent Protection
Kennesaw State University DigitalCommons@Kennesaw State University Faulty Publiations 00 Liensing and Patent Protetion Arijit Mukherjee University of Nottingham Aniruddha Baghi Kennesaw State University,
More informationCAS Course 3 - Actuarial Models
CAS Course 3 - Actuarial Models Before commencing study for this four-hour, multiple-choice examination, candidates should read the introduction to Materials for Study. Items marked with a bold W are available
More informationTransport tax reforms, two-part tariffs, and revenue recycling. - A theoretical result
Transport tax reforms, to-part tariffs, and revenue reyling - A theoretial result Abstrat Jens Erik Nielsen Danish Transport Researh Institute We explore the interation beteen taxes on onership and on
More information4.3 Normal distribution
43 Normal distribution Prof Tesler Math 186 Winter 216 Prof Tesler 43 Normal distribution Math 186 / Winter 216 1 / 4 Normal distribution aka Bell curve and Gaussian distribution The normal distribution
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 informationPractice Exam 1. Loss Amount Number of Losses
Practice Exam 1 1. You are given the following data on loss sizes: An ogive is used as a model for loss sizes. Determine the fitted median. Loss Amount Number of Losses 0 1000 5 1000 5000 4 5000 10000
More informationChapter 5. Statistical inference for Parametric Models
Chapter 5. Statistical inference for Parametric Models Outline Overview Parameter estimation Method of moments How good are method of moments estimates? Interval estimation Statistical Inference for Parametric
More informationMarket Power Rents and Climate Change Mitigation. A Rationale for Export Taxes on Coal? Philipp M. Richter, Roman Mendelevitch, Frank Jotzo
Market Power Rents and Climate Change Mitigation A Rationale for Export Taxes on Coal? Philipp M. Rihter, Roman Mendelevith, Frank Jotzo Roman Mendelevith 9 th Trans-Atlanti Infraday, FERC, Washington
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 informationNormal distribution Approximating binomial distribution by normal 2.10 Central Limit Theorem
1.1.2 Normal distribution 1.1.3 Approimating binomial distribution by normal 2.1 Central Limit Theorem Prof. Tesler Math 283 Fall 216 Prof. Tesler 1.1.2-3, 2.1 Normal distribution Math 283 / Fall 216 1
More informationClipping Coupons: Redemption of Offers with Forward-Looking Consumers
Clipping Coupons: Redemption of Offers with Forward-Looking Consumers Kissan Joseph Oksana Loginova Marh 6, 2019 Abstrat Consumer redemption behavior pertaining to oupons, gift ertifiates, produt sampling,
More informationThe diversification delta: A different perspective. Author. Published. Journal Title. Version DOI. Copyright Statement.
The diversifiation delta: A different perspetive Author Salazar Flores, Yuri, Bianhi, Robert, Drew, Mihael, Truk, Stefan Published 07 Journal Title Journal of Portfolio Management Version Post-print DOI
More informationRanking dynamics and volatility. Ronald Rousseau KU Leuven & Antwerp University, Belgium
Ranking dynamis and volatility Ronald Rousseau KU Leuven & Antwerp University, Belgium ronald.rousseau@kuleuven.be Joint work with Carlos Garía-Zorita, Sergio Marugan Lazaro and Elias Sanz-Casado Department
More informationA 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 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 informationCalculus VCT plc. For investors looking for regular, tax-free income. Please send completed application packs to:
Calulus VCT pl For investors looking for regular, tax-free inome Please send ompleted appliation paks to: Calulus EIS Fund, 104 Park Street, London, W1K 6NF A portfolio of entrepreneurial, growing UK ompanies
More informationExplanatory Memorandum
IN THE KEYS HEAVILY INDEBTED POOR COUNTRIES (LIMITATION ON DEBT RECOVERY) BILL 202 Explanatory Memorandum. This Bill is promoted by the Counil of Ministers. 2. Clause provides for the short title of the
More informationPROSPECTUS May 1, Agency Shares
Dreyfus Institutional Reserves Funds Dreyfus Institutional Reserves Money Fund Class/Tiker Ageny shares DRGXX Dreyfus Institutional Reserves Treasury Fund Class/Tiker Ageny shares DGYXX Dreyfus Institutional
More informationChapter 3 Common Families of Distributions. Definition 3.4.1: A family of pmfs or pdfs is called exponential family if it can be expressed as
Lecture 0 on BST 63: Statistical Theory I Kui Zhang, 09/9/008 Review for the previous lecture Definition: Several continuous distributions, including uniform, gamma, normal, Beta, Cauchy, double exponential
More information2.1 Random variable, density function, enumerative density function and distribution function
Risk Theory I Prof. Dr. Christian Hipp Chair for Science of Insurance, University of Karlsruhe (TH Karlsruhe) Contents 1 Introduction 1.1 Overview on the insurance industry 1.1.1 Insurance in Benin 1.1.2
More informationLecture 5: Fundamentals of Statistical Analysis and Distributions Derived from Normal Distributions
Lecture 5: Fundamentals of Statistical Analysis and Distributions Derived from Normal Distributions ELE 525: Random Processes in Information Systems Hisashi Kobayashi Department of Electrical Engineering
More informationSampling & populations
Sampling & populations Sample proportions Sampling distribution - small populations Sampling distribution - large populations Sampling distribution - normal distribution approximation Mean & variance of
More informationدرس هفتم یادگیري ماشین. (Machine Learning) دانشگاه فردوسی مشهد دانشکده مهندسی رضا منصفی
یادگیري ماشین توزیع هاي نمونه و تخمین نقطه اي پارامترها Sampling Distributions and Point Estimation of Parameter (Machine Learning) دانشگاه فردوسی مشهد دانشکده مهندسی رضا منصفی درس هفتم 1 Outline Introduction
More informationON TRANSACTION COSTS IN STOCK TRADING
QUANTITATIVE METHODS IN ECONOMICS Volume XVIII, No., 07, pp. 58 67 ON TRANSACTION COSTS IN STOCK TRADING Marek Andrzej Koiński Faulty of Applied Informatis and Mathematis Warsaw University of Life Sienes
More informationFINANCIAL SERVICES (FEES) ORDER 2018
Finanial Servies (Fees) Order 2018 Artile 1 Statutory Doument No. 2018/0060 Finanial Servies At 2008 FINANCIAL SERVICES (FEES) ORDER 2018 Approved by Tynwald: 20 Marh 2018 Coming into Operation: 1 April
More informationStatistics 6 th Edition
Statistics 6 th Edition Chapter 5 Discrete Probability Distributions Chap 5-1 Definitions Random Variables Random Variables Discrete Random Variable Continuous Random Variable Ch. 5 Ch. 6 Chap 5-2 Discrete
More informationi e SD No.2015/0206 PAYMENT SERVICES REGULATIONS 2015
i e SD No.2015/0206 PAYMENT SERVICES REGULATIONS 2015 Payment Servies Regulations 2015 Index PAYMENT SERVICES REGULATIONS 2015 Index Regulation Page PART 1 INTRODUCTION 7 1 Title... 7 2 Commenement...
More informationSTAT/MATH 395 PROBABILITY II
STAT/MATH 395 PROBABILITY II Distribution of Random Samples & Limit Theorems Néhémy Lim University of Washington Winter 2017 Outline Distribution of i.i.d. Samples Convergence of random variables The Laws
More informationHomework Problems Stat 479
Chapter 2 1. Model 1 is a uniform distribution from 0 to 100. Determine the table entries for a generalized uniform distribution covering the range from a to b where a < b. 2. Let X be a discrete random
More informationAUTHOR COPY. The co-production approach to service: a theoretical background
Journal of the Operational Researh Soiety (213), 1 8 213 Operational Researh Soiety td. All rights reserved. 16-5682/13 www.palgrave-journals.om/jors/ The o-prodution approah to servie: a theoretial bakground
More informationCIVL Learning Objectives. Definitions. Discrete Distributions
CIVL 3103 Discrete Distributions Learning Objectives Define discrete distributions, and identify common distributions applicable to engineering problems. Identify the appropriate distribution (i.e. binomial,
More informationEngineering Statistics ECIV 2305
Engineering Statistics ECIV 2305 Section 5.3 Approximating Distributions with the Normal Distribution Introduction A very useful property of the normal distribution is that it provides good approximations
More informationProblems from 9th edition of Probability and Statistical Inference by Hogg, Tanis and Zimmerman:
Math 224 Fall 207 Homework 5 Drew Armstrong Problems from 9th edition of Probability and Statistical Inference by Hogg, Tanis and Zimmerman: Section 3., Exercises 3, 0. Section 3.3, Exercises 2, 3, 0,.
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 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 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 informationA Comparison Between Skew-logistic and Skew-normal Distributions
MATEMATIKA, 2015, Volume 31, Number 1, 15 24 c UTM Centre for Industrial and Applied Mathematics A Comparison Between Skew-logistic and Skew-normal Distributions 1 Ramin Kazemi and 2 Monireh Noorizadeh
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 informationAUDITING COST OVERRUN CLAIMS *
AUDITING COST OVERRUN CLAIMS * David Pérez-Castrillo # University of Copenhagen & Universitat Autònoma de Barelona Niolas Riedinger ENSAE, Paris Abstrat: We onsider a ost-reimbursement or a ost-sharing
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 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 information5. In fact, any function of a random variable is also a random variable
Random Variables - Class 11 October 14, 2012 Debdeep Pati 1 Random variables 1.1 Expectation of a function of a random variable 1. Expectation of a function of a random variable 2. We know E(X) = x xp(x)
More information4-1. Chapter 4. Commonly Used Distributions by The McGraw-Hill Companies, Inc. All rights reserved.
4-1 Chapter 4 Commonly Used Distributions 2014 by The Companies, Inc. All rights reserved. Section 4.1: The Bernoulli Distribution 4-2 We use the Bernoulli distribution when we have an experiment which
More informationProbability Theory. Mohamed I. Riffi. Islamic University of Gaza
Probability Theory Mohamed I. Riffi Islamic University of Gaza Table of contents 1. Chapter 2 Discrete Distributions The binomial distribution 1 Chapter 2 Discrete Distributions Bernoulli trials and the
More informationAggregation and capital allocation for portfolios of dependent risks
Aggregation and capital allocation for portfolios of dependent risks... with bivariate compound distributions Etienne Marceau, Ph.D. A.S.A. (Joint work with Hélène Cossette and Mélina Mailhot) Luminy,
More informationSTATISTICS and PROBABILITY
Introduction to Statistics Atatürk University STATISTICS and PROBABILITY LECTURE: PROBABILITY DISTRIBUTIONS Prof. Dr. İrfan KAYMAZ Atatürk University Engineering Faculty Department of Mechanical Engineering
More informationContinuous probability distribution
Microarray Center BIOSTATISTICS Lecture 6 Continuous Probability Distributions 16-4-1 Lecture 6. Continuous probability distributions Dr. Petr Nazarov petr.nazarov@crp-sante.lu OUTLINE Lecture 1 Continuous
More informationA simple Consumption-based Capital Asset Pricing Model
simle Consumtion-based Caital sset Priing Model Integrated wit MCandless and Wallae Kjetil Storesletten Setember 3, 200 6 Introdution Purose of leture: understand te onsumtion-based aital asset riing model
More informationDefinition 9.1 A point estimate is any function T (X 1,..., X n ) of a random sample. We often write an estimator of the parameter θ as ˆθ.
9 Point estimation 9.1 Rationale behind point estimation When sampling from a population described by a pdf f(x θ) or probability function P [X = x θ] knowledge of θ gives knowledge of the entire population.
More informationRoad Transport Regulations 2018 ROAD TRANSPORT REGULATIONS Title Commencement Interpretation... 5
Road Transport Regulations 2018 Index ROAD TRANSPORT REGULATIONS 2018 Index Regulation Page 1 Title... 5 2 Commenement... 5 3 Interpretation... 5 PART 2 ROAD TRANSPORT LICENSING COMMITTEE 6 MEETINGS OF
More informationLecture Notes 6. Assume F belongs to a family of distributions, (e.g. F is Normal), indexed by some parameter θ.
Sufficient Statistics Lecture Notes 6 Sufficiency Data reduction in terms of a particular statistic can be thought of as a partition of the sample space X. Definition T is sufficient for θ if the conditional
More informationThe Bernoulli distribution
This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike License. Your use of this material constitutes acceptance of that license and the conditions of use of materials on this
More informationSOCIETY OF ACTUARIES EXAM STAM SHORT-TERM ACTUARIAL MATHEMATICS EXAM STAM SAMPLE QUESTIONS
SOCIETY OF ACTUARIES EXAM STAM SHORT-TERM ACTUARIAL MATHEMATICS EXAM STAM SAMPLE QUESTIONS Questions 1-307 have been taken from the previous set of Exam C sample questions. Questions no longer relevant
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 informationAnti-Money Laundering and Countering the Financing of Terrorism ANTI-MONEY LAUNDERING AND COUNTERING THE FINANCING OF TERRORISM CODE 2015
Anti-Money Laundering and Countering the Finaning of Terrorism Index ANTI-MONEY LAUNDERING AND COUNTERING THE FINANCING OF TERRORISM CODE 2015 Index Paragraph Page PART 1 INTRODUCTORY 3 1 Title... 3 2
More informationGeneral Registry (Miscellaneous Fees) Order 2016 GENERAL REGISTRY (MISCELLANEOUS FEES) ORDER 2016
General Registry (Misellaneous Fees) Order 2016 Index GENERAL REGISTRY (MISCELLANEOUS FEES) ORDER 2016 Index Artile Page 1 Title... 3 2 Commenement... 3 3 Misellaneous Fees in the General Registry... 3
More informationD Accounting for management
QU EST IONS ounting for management 17 mins 1.1 Whih of the following statements about qualities of good information is false? It should be relevant for its purposes It should be ommuniated to the right
More informationThe Central Exposed to Risk as a Hedging Strategy: A Case Study of a Kenyan Pension Scheme
The Central Eposed to Risk as a Hedging Strategy: A Case Study of a Kenyan Pension Sheme Abstrat Titus Kipkoeh Rotih Shool of Siene, University of Eldoret, P.O BOX 1125-30100, Eldoret, Kenya. This study
More informationCambridge University Press Risk Modelling in General Insurance: From Principles to Practice Roger J. Gray and Susan M.
adjustment coefficient, 272 and Cramér Lundberg approximation, 302 existence, 279 and Lundberg s inequality, 272 numerical methods for, 303 properties, 272 and reinsurance (case study), 348 statistical
More informationMATH 3200 Exam 3 Dr. Syring
. Suppose n eligible voters are polled (randomly sampled) from a population of size N. The poll asks voters whether they support or do not support increasing local taxes to fund public parks. Let M be
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 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 informationMODELLING OF INCOME AND WAGE DISTRIBUTION USING THE METHOD OF L-MOMENTS OF PARAMETER ESTIMATION
International Days of Statistics and Economics, Prague, September -3, MODELLING OF INCOME AND WAGE DISTRIBUTION USING THE METHOD OF L-MOMENTS OF PARAMETER ESTIMATION Diana Bílková Abstract Using L-moments
More informationก ก ก ก ก ก ก. ก (Food Safety Risk Assessment Workshop) 1 : Fundamental ( ก ( NAC 2010)) 2 3 : Excel and Statistics Simulation Software\
ก ก ก ก (Food Safety Risk Assessment Workshop) ก ก ก ก ก ก ก ก 5 1 : Fundamental ( ก 29-30.. 53 ( NAC 2010)) 2 3 : Excel and Statistics Simulation Software\ 1 4 2553 4 5 : Quantitative Risk Modeling Microbial
More information