The Beta-Burr Type V Distribution: Its Properties and Application to Real Life Data

Total Page:16

File Type:pdf, Size:1020Kb

The Beta-Burr Type V Distribution: Its Properties and Application to Real Life Data International Journal of Advanced Statistics and Probability, 5 (2) (2017) 83-86 International Journal of Advanced Statistics and Probability Website: www.sciencepubco.com/index.php/IJASP doi: 10.14419/ijasp.v5i2.7862 Research paper The beta-burr type v distribution: its properties and application to real life data Husaini Garba Dikko *, Aliyu Yakubu, Alfa Aliyu Saidu Department of Statistics, Ahmadu Bello University Zaria-Nigeria *Corresponding author E-mail: [email protected] Abstract A new distribution called the beta-Burr type V distribution that extends the Burr type V distribution was defined, investigated and estab- lished. The properties examined provide a comprehensive mathematical treatment of the distribution. Additionally, various structural properties of the new distribution verified include probability density function verification, asymptotic behavior, Hazard Rate Function and the cumulative distribution. Subsequently, we used the maximum likelihood estimation procedure to estimate the parameters of the new distribution. Application of real data set indicates that this new distribution would serve as a good alternative distribution function to model real- life data in many areas. Keywords: Burr Type V; Beta Burr V; Hazard Rate; Maximum Likelihood Estimation. tan x 2 tan x 1 1. Introduction e sec x 1 e x 22 fx Burr's type V is one of the twelve distributions introduced in 1942 0 elsewhere by [5], which can be used to fit basically any given set of unimod- al data [7]. So many researchers have established beta –G distribu- tion; these include among others: Nadarajah and Gupta [20] de- (1) fined the beta-Frechet distribution; Famoye, et al. [5] defined the beta-Weibull; Nadaraja and Kotz [19] defined the beta- The cumulative distribution function of Burr type V distribution is: exponential distribution; Kong et al. [15] proposed the beta- gamma distribution. Fischer and Vaughan [12] introduced the tan x beta-hyperbolic secant distribution; beta-Gumbel (BGU) distribu- F x 1 (2) tion was introduced by Nadarajah and Kotz [18]; the beta-Pareto distribution was defined and studied by Akinsete et al. [4]; the Where α and β are the location parameters of the distribution. beta-Rayleigh distribution was proposed by Akinsete and Lowe The rest of the paper is organized as follows: in section two, we [3]; beta-Burr XII by Paranaiba et al. [21]; and very recently, briefly introduce the beta-burr type V (BBV) distribution, its Merovci et al. [17] developed the beta-Burr type X distribution properties is studied in section 3 while estimates of the parameters leaving the rest types of Burr distributions with little or no interest of the distribution using the maximum likelihood estimation tech- from researchers. This is the case with Burr type V distribution. It nique is obtained in section 4, in section 5 we look at the compari- is clear therefore, that this distribution is yet to catch the eyes of son of our new distribution with other known distributions and we researchers. Burr Type V distribution can be used to model real finally conclude in section 6. lifetime data. Ever since the introduction of the Burr distribution, it has some- what been neglected as an option in the analysis of lifetime data. 2. Beta- burr type V (BBV) Even though there are researchers that indicate that this distribu- tion possesses sufficient flexibility to make it a possible model for Beta distribution takes on several problems in reliability analysis various types of data; for example, it was used to examine the since it has been widely acclaim to be very powerful and applica- strengths of 1.5 cm glass fibers [17], Paranaiba et al. [21] studied ble probability distribution. In recent years, development focuses cancer recurrence by using the real data set. on new techniques for building meaningful distributions, which Burr's type V distribution extension with two parameters is include the use of the logit of beta (the link function of the beta through the beta-G generator pioneered by Eugene et al. [8]. We generalized distribution which was introduced by Jones [14] and investigate and discern the properties of the proposed distribution elsewhere [10]). The logit of beta introduced by Jones has proba- from G(x) which was used by Eugene et al [8]. It is known as the bility density function (pdf). beta generalized class of distribution, and it has two shape pa- rameters in the generator, and it is given in equation 3. f(x) g(x) = [F(x)]a−1[1 − F(x)]b−1 (3) The probability density function of Burr type V distribution is: B(a,b) Copyright © 2017HusainiGarbaDikko et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. 84 International Journal of Advanced Statistics and Probability Where 3.0 a=0.6, b=0.9, alpha =0.5, beta=0.2 F(x)is the cumulative distribution function of the baseline distri- a=3.0, b=2.5, alpha =0.5, beta=0.8 a=3.0, b=2.5, alpha =1.5, beta=0.8 bution and 2.5 a=3.0, b=2.5, alpha =1.5, beta=2.0 a=0.5, b=0.5, alpha =1.5, beta=2.0 f(x)is the probability density function of the baseline distribution a=1.5, b=1.5, alpha =1.5, beta=1.5 Many researchers have used this technique to come up with many 2.0 compound distributions which include: [1],[3],[4],[6],[8],[9],[12],[15],[16],[17],[18],[19],[20],[21],[23] 1.5 and many others. The aim of this paper is to propose a new model BBTV of Pdf called beta- Burr type V distribution. Using the logit defined in (3), 1.0 the pdf of the proposed distribution is given by equation (4): 0.5 − tan x 2 αβe sec x 0.0 f(x) = [1 + βe− tan x]−αa−1[1 − (1 + β(a,b) -50 0 50 βe− tan x)−α]b−1 (4) X Fig 1: Graph of Beta-Burr Type V Distribution for Various Parameter 3. Properties Values. 3.2. Asymptotic behaviour We will examine the statistical properties of the distribution we have put forth in the previous section. It will be scrutinized by In order to investigate the asymptotic behavior of the proposed verifying the asymptotic behaviour, true probability function, π model (Beta –Burr type V), we find limit as x → and limit as cumulative density and hazard rate function. 2 π x → − of the BBV distribution. 2 3.1. Investigation of proper PDF αβe− tan xsec2x ( ) ( − tan x)−αa−1 [ To verify that the proposed density function is a proper pdf, we limπ f x = limπ 1 + βe 1 x→− x→ − β (a, b) integrate (4) with respect to x and see whether the integral is equal 2 2 − tan x −α b−1 to unity or otherwise. If it is equal to unity, the density function is − (1 + βe ) ] a proper pdf otherwise it is not. This is equal to zero, since π π 2 2 αβ − tan x 2 − tan x −αa−1 ∫ π f(x)dx = ∫ π e sec x[1 + βe ] [1 − − 푡푎푛 푥 2 − 푡푎푛 푥 −훼−1 − − β(a,b) 푙푚휋 훼훽푒 푠푒푐 푥(1 + 훽푒 ) = 0 2 2 푥→ − − tan x −α b−1 2 (1 + βe ) ] dx (5) Similarly, 푙푚휋 푓(푥) = 0 To integrate (5), 푥→ 2 − tan x 2 dp 2 Letp = βe , then dp = −psec x dxand dx = 2 Since 푙푚휋 푠푒푐 푥 = 0 −psec x 푥→ 2 Hence, 3.3. Cumulative density function (CDF) of BBV distri- π ∞ bution 2 α ∫ f(x)dx = ∫ psec2x [1 + p]−αa−1 [1 π − 0 β(a, b) In this section, the cumulative density function of the BBV distri- 2 dp bution will be obtained − (1 + p)−α]b−1 psec2 x 푥 훼훽 − 푡푎푛 푦 2 −푡푎푛푦 –훼푎−1 퐹(푥) = ∫ 휋 푒 푠푒푐 푦[1 + 훽푒 ] [1 − − 훽(푎,푏) ∞ 2 α −푡푎푛푦 −훼 푏−1 = ∫ (1 + p)−αa−1 [1 − (1 + p)−α]b−1 dp (1 + 훽푒 ) ] 푑푦 (6) 0 β(a, b) To integrate (6), 푙푒푡 푝 = 훽푒푡푎푛 푦, then 1 −α − let R = (1 + p) ThenR α = 1 + p 2 푑푝 푑푝 = −푝푠푒푐 푦푑푦And푑푦 = −푝푠푒푐2푦 This implies Therefore, −α−1 dR dR = −α(1 + p) dpAnddp = 1 1+ −αR α ∞ 훼 퐹(푥) = ∫ 푝푠푒푐2푦[1 + 푝]−훼푎−1[1 Therefore, 훽푒− 푡푎푛 푥 훽(푎, 푏) −훼 푏−1 푑푝 π − (1 + 푝) ] 1 2 2 1 β(a, b) 푝푠푒푐 푦 ∫ f(x)dx = ∫ Ra−1(1 − R)b−1 dR = = 1 π − β(a, b) 0 β(a, b) ∞ 2 훼 −훼푎−1 −훼 푏−1 = ∫ [1 + 푝] [1 − (1 + 푝) ] 푑푝 − 푡푎푛 푥 훽(푎, 푏) Hence, f(x) is a proper pdf? 훽푒 The graph of the BBV distribution for various parameter values (a, Let푅 = (1 + 푝)−훼then, b, α and β) is given in fig. 1 훼+1 푑푅 푑푅 = −훼푅 훼 푑푝 And 푑푝 = 훼+1 −훼푅 훼 Hence, International Journal of Advanced Statistics and Probability 85 −훼 ∑푛 (1+훽푒− 푡푎푛 푥) 훼푛훽푛푒− 푖=1 푡푎푛 푥푖 훼 훼푎+1 푑푅 푛 2 −푡푎푛푥푖 −훼푎−1 푏−1 = 푛 ∏ 푠푒푐 푥 (1 + 훽푒 ) [1 − 퐹(푥) = ∫ 푅 훼 [1 − 푅] ( ) 푖=1 푖 훽(푎, 푏) 훼+1 (훽 푎,푏 ) 0 훼 − 푡푎푛 푥푖 −훼 푏−1 훼푅 (1 + 훽푒 ) ] (11) − 푡푎푛 푥 −훼 (1+훽푒 ) 훼푎+1−훼−1 훼 푏−1 푑푅 Taking log of (11), we have the log-likelihood given by: 퐹(푥) = ∫ 푅 훼 (1 − 푅) 0 훽(푎, 푏) 훼 푛 푙(푥; 푎, 푏, 훼, 훽) = 푛 푙표푔 훼 + 푛 푙표푔 훽 − ∑푖=1 푡푎푛푥푖 − 푛 푙표푔(⌈푎) − − 푡푎푛 푥 −훼 푛 2 1 (1+훽푒 ) 푛 푙표푔(⌈푏) + 푛 푙표푔(⌈푎 + 푏) + ∑푖=1 푙표푔(푠푒푐 푥푖) − (훼푎 + 퐹(푥) = ∫ 푅푎−1(1 − 푅)푏−1푑푅 1) ∑푛 푙표푔(1 + 훽푒푡푎푛푥푖) + (푏 − 1) ∑푛 푙표푔[1 − (1 + 훽(푎, 푏) 푖=1 푖=1 0 훽푒− 푡푎푛 푥푖)−훼] (12) Therefore, Differentiating (12) partially with respect to 푎, 푏, 훼푎푛푑훽 gives 훽[{1+훽푒− 푡푎푛 푥}−훼,푎,푏] (13), (14), (15) and (16): 퐹(푥) = (7) 훽(푎,푏) 휕푙(푥:푎,푏,훼,훽) = 푛휓(푎) + 푛휓(푎 + 푏) − 훼 (13) − 푡푎푛 푥 −훼 휕푎 Where훽[{1 + 훽푒 } , 푎, 푏] is an incomplete beta function? 휕푙(푥: 푎, 푏, 훼, 훽) 3.4.
Recommended publications
  • Severity Models - Special Families of Distributions
    Severity Models - Special Families of Distributions Sections 5.3-5.4 Stat 477 - Loss Models Sections 5.3-5.4 (Stat 477) Claim Severity Models Brian Hartman - BYU 1 / 1 Introduction Introduction Given that a claim occurs, the (individual) claim size X is typically referred to as claim severity. While typically this may be of continuous random variables, sometimes claim sizes can be considered discrete. When modeling claims severity, insurers are usually concerned with the tails of the distribution. There are certain types of insurance contracts with what are called long tails. Sections 5.3-5.4 (Stat 477) Claim Severity Models Brian Hartman - BYU 2 / 1 Introduction parametric class Parametric distributions A parametric distribution consists of a set of distribution functions with each member determined by specifying one or more values called \parameters". The parameter is fixed and finite, and it could be a one value or several in which case we call it a vector of parameters. Parameter vector could be denoted by θ. The Normal distribution is a clear example: 1 −(x−µ)2=2σ2 fX (x) = p e ; 2πσ where the parameter vector θ = (µ, σ). It is well known that µ is the mean and σ2 is the variance. Some important properties: Standard Normal when µ = 0 and σ = 1. X−µ If X ∼ N(µ, σ), then Z = σ ∼ N(0; 1). Sums of Normal random variables is again Normal. If X ∼ N(µ, σ), then cX ∼ N(cµ, cσ). Sections 5.3-5.4 (Stat 477) Claim Severity Models Brian Hartman - BYU 3 / 1 Introduction some claim size distributions Some parametric claim size distributions Normal - easy to work with, but careful with getting negative claims.
    [Show full text]
  • Package 'Distributional'
    Package ‘distributional’ February 2, 2021 Title Vectorised Probability Distributions Version 0.2.2 Description Vectorised distribution objects with tools for manipulating, visualising, and using probability distributions. Designed to allow model prediction outputs to return distributions rather than their parameters, allowing users to directly interact with predictive distributions in a data-oriented workflow. In addition to providing generic replacements for p/d/q/r functions, other useful statistics can be computed including means, variances, intervals, and highest density regions. License GPL-3 Imports vctrs (>= 0.3.0), rlang (>= 0.4.5), generics, ellipsis, stats, numDeriv, ggplot2, scales, farver, digest, utils, lifecycle Suggests testthat (>= 2.1.0), covr, mvtnorm, actuar, ggdist RdMacros lifecycle URL https://pkg.mitchelloharawild.com/distributional/, https: //github.com/mitchelloharawild/distributional BugReports https://github.com/mitchelloharawild/distributional/issues Encoding UTF-8 Language en-GB LazyData true Roxygen list(markdown = TRUE, roclets=c('rd', 'collate', 'namespace')) RoxygenNote 7.1.1 1 2 R topics documented: R topics documented: autoplot.distribution . .3 cdf..............................................4 density.distribution . .4 dist_bernoulli . .5 dist_beta . .6 dist_binomial . .7 dist_burr . .8 dist_cauchy . .9 dist_chisq . 10 dist_degenerate . 11 dist_exponential . 12 dist_f . 13 dist_gamma . 14 dist_geometric . 16 dist_gumbel . 17 dist_hypergeometric . 18 dist_inflated . 20 dist_inverse_exponential . 20 dist_inverse_gamma
    [Show full text]
  • Bivariate Burr Distributions by Extending the Defining Equation
    DECLARATION This thesis contains no material which has been accepted for the award of any other degree or diploma in any university and to the best ofmy knowledge and belief, it contains no material previously published by any other person, except where due references are made in the text ofthe thesis !~ ~ Cochin-22 BISMI.G.NADH June 2005 ACKNOWLEDGEMENTS I wish to express my profound gratitude to Dr. V.K. Ramachandran Nair, Professor, Department of Statistics, CUSAT for his guidance, encouragement and support through out the development ofthis thesis. I am deeply indebted to Dr. N. Unnikrishnan Nair, Fonner Vice Chancellor, CUSAT for the motivation, persistent encouragement and constructive criticism through out the development and writing of this thesis. The co-operation rendered by all teachers in the Department of Statistics, CUSAT is also thankfully acknowledged. I also wish to express great happiness at the steadfast encouragement rendered by my family. Finally I dedicate this thesis to the loving memory of my father Late. Sri. C.K. Gopinath. BISMI.G.NADH CONTENTS Page CHAPTER I INTRODUCTION 1 1.1 Families ofDistributions 1.2· Burr system 2 1.3 Present study 17 CHAPTER II BIVARIATE BURR SYSTEM OF DISTRIBUTIONS 18 2.1 Introduction 18 2.2 Bivariate Burr system 21 2.3 General Properties of Bivariate Burr system 23 2.4 Members ofbivariate Burr system 27 CHAPTER III CHARACTERIZATIONS OF BIVARIATE BURR 33 SYSTEM 3.1 Introduction 33 3.2 Characterizations ofbivariate Burr system using reliability 33 concepts 3.3 Characterization ofbivariate
    [Show full text]
  • Loss Distributions
    Munich Personal RePEc Archive Loss Distributions Burnecki, Krzysztof and Misiorek, Adam and Weron, Rafal 2010 Online at https://mpra.ub.uni-muenchen.de/22163/ MPRA Paper No. 22163, posted 20 Apr 2010 10:14 UTC Loss Distributions1 Krzysztof Burneckia, Adam Misiorekb,c, Rafał Weronb aHugo Steinhaus Center, Wrocław University of Technology, 50-370 Wrocław, Poland bInstitute of Organization and Management, Wrocław University of Technology, 50-370 Wrocław, Poland cSantander Consumer Bank S.A., Wrocław, Poland Abstract This paper is intended as a guide to statistical inference for loss distributions. There are three basic approaches to deriving the loss distribution in an insurance risk model: empirical, analytical, and moment based. The empirical method is based on a sufficiently smooth and accurate estimate of the cumulative distribution function (cdf) and can be used only when large data sets are available. The analytical approach is probably the most often used in practice and certainly the most frequently adopted in the actuarial literature. It reduces to finding a suitable analytical expression which fits the observed data well and which is easy to handle. In some applications the exact shape of the loss distribution is not required. We may then use the moment based approach, which consists of estimating only the lowest characteristics (moments) of the distribution, like the mean and variance. Having a large collection of distributions to choose from, we need to narrow our selection to a single model and a unique parameter estimate. The type of the objective loss distribution can be easily selected by comparing the shapes of the empirical and theoretical mean excess functions.
    [Show full text]
  • Normal Variance Mixtures: Distribution, Density and Parameter Estimation
    Normal variance mixtures: Distribution, density and parameter estimation Erik Hintz1, Marius Hofert2, Christiane Lemieux3 2020-06-16 Abstract Normal variance mixtures are a class of multivariate distributions that generalize the multivariate normal by randomizing (or mixing) the covariance matrix via multiplication by a non-negative random variable W . The multivariate t distribution is an example of such mixture, where W has an inverse-gamma distribution. Algorithms to compute the joint distribution function and perform parameter estimation for the multivariate normal and t (with integer degrees of freedom) can be found in the literature and are implemented in, e.g., the R package mvtnorm. In this paper, efficient algorithms to perform these tasks in the general case of a normal variance mixture are proposed. In addition to the above two tasks, the evaluation of the joint (logarithmic) density function of a general normal variance mixture is tackled as well, as it is needed for parameter estimation and does not always exist in closed form in this more general setup. For the evaluation of the joint distribution function, the proposed algorithms apply randomized quasi-Monte Carlo (RQMC) methods in a way that improves upon existing methods proposed for the multivariate normal and t distributions. An adaptive RQMC algorithm that similarly exploits the superior convergence properties of RQMC methods is presented for the task of evaluating the joint log-density function. In turn, this allows the parameter estimation task to be accomplished via an expectation-maximization- like algorithm where all weights and log-densities are numerically estimated. It is demonstrated through numerical examples that the suggested algorithms are quite fast; even for high dimensions around 1000 the distribution function can be estimated with moderate accuracy using only a few seconds of run time.
    [Show full text]
  • Handbook on Probability Distributions
    R powered R-forge project Handbook on probability distributions R-forge distributions Core Team University Year 2009-2010 LATEXpowered Mac OS' TeXShop edited Contents Introduction 4 I Discrete distributions 6 1 Classic discrete distribution 7 2 Not so-common discrete distribution 27 II Continuous distributions 34 3 Finite support distribution 35 4 The Gaussian family 47 5 Exponential distribution and its extensions 56 6 Chi-squared's ditribution and related extensions 75 7 Student and related distributions 84 8 Pareto family 88 9 Logistic distribution and related extensions 108 10 Extrem Value Theory distributions 111 3 4 CONTENTS III Multivariate and generalized distributions 116 11 Generalization of common distributions 117 12 Multivariate distributions 133 13 Misc 135 Conclusion 137 Bibliography 137 A Mathematical tools 141 Introduction This guide is intended to provide a quite exhaustive (at least as I can) view on probability distri- butions. It is constructed in chapters of distribution family with a section for each distribution. Each section focuses on the tryptic: definition - estimation - application. Ultimate bibles for probability distributions are Wimmer & Altmann (1999) which lists 750 univariate discrete distributions and Johnson et al. (1994) which details continuous distributions. In the appendix, we recall the basics of probability distributions as well as \common" mathe- matical functions, cf. section A.2. And for all distribution, we use the following notations • X a random variable following a given distribution, • x a realization of this random variable, • f the density function (if it exists), • F the (cumulative) distribution function, • P (X = k) the mass probability function in k, • M the moment generating function (if it exists), • G the probability generating function (if it exists), • φ the characteristic function (if it exists), Finally all graphics are done the open source statistical software R and its numerous packages available on the Comprehensive R Archive Network (CRAN∗).
    [Show full text]
  • Tables for Exam STAM the Reading Material for Exam STAM Includes a Variety of Textbooks
    Tables for Exam STAM The reading material for Exam STAM includes a variety of textbooks. Each text has a set of probability distributions that are used in its readings. For those distributions used in more than one text, the choices of parameterization may not be the same in all of the books. This may be of educational value while you study, but could add a layer of uncertainty in the examination. For this latter reason, we have adopted one set of parameterizations to be used in examinations. This set will be based on Appendices A & B of Loss Models: From Data to Decisions by Klugman, Panjer and Willmot. A slightly revised version of these appendices is included in this note. A copy of this note will also be distributed to each candidate at the examination. Each text also has its own system of dedicated notation and terminology. Sometimes these may conflict. If alternative meanings could apply in an examination question, the symbols will be defined. For Exam STAM, in addition to the abridged table from Loss Models, sets of values from the standard normal and chi-square distributions will be available for use in examinations. These are also included in this note. When using the normal distribution, choose the nearest z-value to find the probability, or if the probability is given, choose the nearest z-value. No interpolation should be used. Example: If the given z-value is 0.759, and you need to find Pr(Z < 0.759) from the normal distribution table, then choose the probability for z-value = 0.76: Pr(Z < 0.76) = 0.7764.
    [Show full text]
  • A New Generalized Burr Family of Distributions Based on Quantile Function -.:: Natural Sciences Publishing
    J. Stat. Appl. Pro. 6, No. 3, 499-504 (2017) 499 Journal of Statistics Applications & Probability An International Journal http://dx.doi.org/10.18576/jsap/060306 A New Generalized Burr Family of Distributions Based on Quantile Function Arslan Nasir1, Mohammad Aljarrah2, Farrukh Jamal3,∗ and M. H. Tahir4 1 Goverment Degree College Lodhran, Punjab, Pakistan. 2 Department of Mathematics, Tafila Technical University, At Tafila 66110, Jordan. 3 Goverment S.A. Post Graduate College Dera Nawab, Punjab, Pakistan. 4 Department of Statistics, The Islamia University Bahawalpur, Bahawalpur-63100, Pakistan. Received: 29 May 2017, Revised: 20 Aug. 2017, Accepted: 25 Aug. 2017 Published online: 1 Nov. 2017 Abstract: We propose a new family from Burr XII distribution, called T-Burr family of distributions based on the T-R{Y} framework. For this family, we consider the quantile functions of three well-known distributions, namely, Lomax, logistic and Weibull, and further developed three sub-families T-Burr{Lomax}, T-Burr{Log-logistic} and T-Burr{Weibull}. Some mathematical properties such as quantile function, mode, Shannon entropy, moments, and mean deviations, of T-R{Y} family are obtained. One special model, namely, Weibull-Burr{Lomax} from T-Burr{Lomax} family is considered and its properties are obtained. This model is flexible and can produce the shapes of the density such as left-skewed, right-skewed, symmetrical, J, and reversed-J, and can have constant, increasing and decreasing hazard rate shapes. The usefulness of this model is demonstrated through applications to censored and complete data sets. Keywords: Burr XII distribution, generalization, quantile function, T-X family, T-R{Y} family.
    [Show full text]
  • Use of the Three-Parameter Burr XII Distribution for Modelling Ambient Daily Maximum Nitrogen Dioxide Concentrations in the Gaborone Fire Brigade
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by American Scientific Research Journal for Engineering, Technology, and Sciences... American Scientific Research Journal for Engineering, Technology, and Sciences (ASRJETS) ISSN (Print) 2313-4410, ISSN (Online) 2313-4402 © Global Society of Scientific Research and Researchers http://asrjetsjournal.org/ Use of the Three-parameter Burr XII Distribution for Modelling Ambient Daily Maximum Nitrogen Dioxide Concentrations in the Gaborone Fire Brigade Wilson Moseki Thupeng* Department of Statistics, University of Botswana, Private Bag 00705, Gaborone, Botswana Email: [email protected] Abstract Air pollution constitutes one of the major problems in urban areas where many sources of airborne pollutants are concentrated. Identifying an appropriate probability model to describe the stochastic behaviour of extreme ambient air pollution level for a specific site or multiple sites forms an integral part of environmental management and pollution control. This paper proposes the use of the three-parameter Burr Type XII distribution to model maximum levels of nitrogen dioxide at a specific site. The study focuses on the daily maximum ambient nitrogen dioxide concentrations for Gaborone Fire Brigade in the winter season, May-July, 2014. The fit of the three-parameter Burr Type XII is compared to that of the three-parameter Dagum and three- parameter log-logistic by using the Kolmogorov-Smirnov and Anderson-Darling criterion for model selection. It is found that the three-parameter Burr Type XII gives the best fit, followed by the log-logistic and the Dagum I, in that order. The results justify the use of the three-parameter Burr Type XII to model ambient air pollution extreme values.
    [Show full text]
  • University of Tartu DISTRIBUTIONS in FINANCIAL MATHEMATICS (MTMS.02.023)
    University of Tartu Faculty of Mathematics and Computer Science Institute of Mathematical Statistics DISTRIBUTIONS IN FINANCIAL MATHEMATICS (MTMS.02.023) Lecture notes Ants Kaasik, Meelis K¨a¨arik MTMS.02.023. Distributions in Financial Mathematics Contents 1 Introduction1 1.1 The data. Prediction problem..................1 1.2 The models............................2 1.3 Summary of the introduction..................2 2 Exploration and visualization of the data2 2.1 Motivating example.......................2 3 Heavy-tailed probability distributions7 4 Detecting the heaviness of the tail 10 4.1 Visual tests............................ 10 4.2 Maximum-sum ratio test..................... 13 4.3 Records test............................ 14 4.4 Alternative definition of heavy tails............... 15 5 Creating new probability distributions. Mixture distribu- tions 17 5.1 Combining distributions. Discrete mixtures.......... 17 5.2 Continuous mixtures....................... 19 5.3 Completely new parts...................... 20 5.4 Parameters of a distribution................... 21 6 Empirical distribution. Kernel density estimation 22 7 Subexponential distributions 27 7.1 Preliminaries. Definition..................... 27 7.2 Properties of subexponential class............... 28 8 Well-known distributions in financial and insurance mathe- matics 30 8.1 Exponential distribution..................... 30 8.2 Pareto distribution........................ 31 8.3 Weibull Distribution....................... 33 i MTMS.02.023. Distributions in Financial Mathematics
    [Show full text]
  • Parametric Representation of the Top of Income Distributions: Options, Historical Evidence and Model Selection
    PARAMETRIC REPRESENTATION OF THE TOP OF INCOME DISTRIBUTIONS: OPTIONS, HISTORICAL EVIDENCE AND MODEL SELECTION Vladimir Hlasny Working Paper 90 March, 2020 The CEQ Working Paper Series The CEQ Institute at Tulane University works to reduce inequality and poverty through rigorous tax and benefit incidence analysis and active engagement with the policy community. The studies published in the CEQ Working Paper series are pre-publication versions of peer-reviewed or scholarly articles, book chapters, and reports produced by the Institute. The papers mainly include empirical studies based on the CEQ methodology and theoretical analysis of the impact of fiscal policy on poverty and inequality. The content of the papers published in this series is entirely the responsibility of the author or authors. Although all the results of empirical studies are reviewed according to the protocol of quality control established by the CEQ Institute, the papers are not subject to a formal arbitration process. Moreover, national and international agencies often update their data series, the information included here may be subject to change. For updates, the reader is referred to the CEQ Standard Indicators available online in the CEQ Institute’s website www.commitmentoequity.org/datacenter. The CEQ Working Paper series is possible thanks to the generous support of the Bill & Melinda Gates Foundation. For more information, visit www.commitmentoequity.org. The CEQ logo is a stylized graphical representation of a Lorenz curve for a fairly unequal distribution of income (the bottom part of the C, below the diagonal) and a concentration curve for a very progressive transfer (the top part of the C).
    [Show full text]
  • On the Robustness of Epsilon Skew Extension for Burr III Distribution on Real Line
    On The Robustness of Epsilon Skew Extension for Burr III Distribution on Real Line ∗ Authors: Mehmet Niyazi C¸ankaya – Department of Statistics, Faculty of Arts and Science U¸sak University, U¸sak ([email protected]) Abdullah Yalc¸ınkaya – Department of Statistics, Faculty of Science Ankara University, ([email protected]) Omer¨ Altındagˇ – Department of Statistics, Faculty of Arts and Science Bilecik S¸eyh Edebali University, ([email protected]) Olcay Arslan – Department of Statistics, Faculty of Science Ankara University , ([email protected]) Abstract: • The Burr III distribution is used in a wide variety of fields of lifetime data analysis, reliability theory, and financial literature, etc. It is defined on the positive axis and has two shape parameters, say c and k. These shape parameters make the distribution quite flexible. They also control the tail behavior of the distribution. In this study, we extent the Burr III distribution to the real axis and also add a skewness parameter, say ε, with epsilon-skew extension approach. When the parameters c and k have a relation such that ck ≈ 1 or ck < 1, it is skewed unimodal. Otherwise, it is skewed bimodal with the same level of peaks on the negative and positive sides of real line. Thus, ESBIII distribution can capture fitting the various data sets even when the number of parameters are three. Location and scale form of this distribution are also given. arXiv:1705.09988v1 [math.ST] 28 May 2017 Some distributional properties of the new distribution are investigated. The maximum likelihood (ML) estimation method for the parameters of ESBIII is considered.
    [Show full text]