The Poisson Burr X Inverse Rayleigh Distribution and Its Applications

Total Page:16

File Type:pdf, Size:1020Kb

The Poisson Burr X Inverse Rayleigh Distribution and Its Applications Journal of Data Science,18(1). P. 56 – 77,2020 DOI:10.6339/JDS.202001_18(1).0003 The Poisson Burr X Inverse Rayleigh Distribution And Its Applications Rania H. M. Abdelkhalek* Department of Statistics, Mathematics and Insurance, Benha University, Egypt. ABSTRACT A new flexible extension of the inverse Rayleigh model is proposed and studied. Some of its fundamental statistical properties are derived. We assessed the performance of the maximum likelihood method via a simulation study. The importance of the new model is shown via three applications to real data sets. The new model is much better than other important competitive models. Keywords: Inverse Rayleigh; Poisson; Simulation; Modeling. * [email protected] Rania H. M. Abdelkhalek 57 1. Introduction and physical motivation The well-known inverse Rayleigh (IR) model is considered as a distribution for a life time random variable (r.v.). The IR distribution has many applications in the area of reliability studies. Voda (1972) proved that the distribution of lifetimes of several types of experimental (Exp) units can be approximated by the IR distribution and studied some properties of the maximum likelihood estimation (MLE) of the its parameter. Mukerjee and Saran (1984) studied the failure rate of an IR distribution. Aslam and Jun (2009) introduced a group acceptance sampling plans for truncated life tests based on the IR. Soliman et al. (2010) studied the Bayesian and non- Bayesian estimation of the parameter of the IR model Dey (2012) mentioned that the IR model has also been used as a failure time distribution although the variance and higher order moments does not exist for it. The probability density function (PDF) and cumulative distribution function (CDF) of the IR distribution are given by (for 푥 ≥ 0) (훼) 2 −3 2 −2 (훼) 2 −2 푔퐼푅 (푥) = 2훼 푥 exp[−훼 푥 ] and 퐺퐼푅 (푥) = exp[−훼 푥 ], respectively, where 훼 > 0 is a scale parameter. It is worth mentioning that the IR model is a special case for the Inverse Weibull model. In the statistical literature, there are many useful extensions of the IR model, for example: beta IR by Barreto-Souza et al. (2011), Marshall-Olkin IR by Krishna et al. (2013), Transmuted IR by Mahmoud and Mandouh (2013), Weibull IR by Afify et al. (2016), odd Lindley IR by Korkmaz et al. (2017), Burr X IR by Yousof et al. (2017a), transmuted Topp- Leone IR by Yousof et al. (2017b), Topp-Leone Generated IR by Yousof et al. (2018b) and odd log-logistic IR Yousof et al. (2018b), among others. In this paper we propose and study a new extension of the IR distribution using the zero truncated Poisson (ZTP) distribution. Suppose that a system has 푁 subsystems functioning independently at a given time where 푁 has ZTP distribution with parameter 휆. It is the conditional probability distribution of a Poisson-distributed r.v., given that the value of the r.v. is not zero. The probability mass function (PMF) of 푁 is given by 휆 푛 푃푀퐹푍푇푃(푁 = 푛) = [exp(−휆) 휆 ]/{푛! [휏휆]}|푛=1,2,… Note that for ZTP r.v., the expected value 퐄(푁|휆) and variance 푉푎푟(푁|휆) are, respectively, given by 퐄(푁|휆) = 휆/휏(휆), where 휏(휆) = −exp(−휆) + 1 and 휆2 휆 + 휆2 푉푎푟(푁|휆) = − 2 + . 휏( ) [휏(휆)] 휆 Suppose that the failure time of each subsystem has the Burr X IR (BX-IR(휃, 훼) for short) defined by the CDF and PDF given by 58 The Poisson Burr X Inverse Rayleigh Distribution And Its Applications 2 휃 exp⁡[−훼2푥−2] 퐻휃,훼 (푥) = [1 − exp⁡(− { } )] 퐵푋−퐼푅 1 − exp⁡[−훼2푥−2] and 2 휃−1 exp[−훼2푥−2] ℎ휃,훼 (푥) = 4휃훼2푥−(2+1){1 − exp[−훼2푥−2]}−3 × [1 − exp (− { } )] 퐵푋−퐼푅 1 − exp[−훼2푥−2] 2 exp[−훼2푥−2] ⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡× exp⁡[−2(훼푥−1)2]exp⁡(− { } ) 1 − exp[−훼2푥−2] respectively, where 훼 > 0 is a scale parameter and 휃 > 0 is the shape parameter. Let 푌 denote the failure time of the ith subsystem and let 푋 = min{푌1, 푌2, ⋯ , 푌푁}. Then the conditional CDF of 푋 given 푁 is 휃,훼 푁 퐹(푥|푁) = 1 − Pr(푋 > 푥|푁) = 1 − [1 − 퐻퐵푋−퐼푅(푥)] Therefore, the unconditional CDF of the Poisson BX-IR (PBX-IR) density function, as described in Ristić and Nadarajah (2014), can be expressed as 휃,훼 휆,휃,훼 1 − exp⁡(−휆퐻퐵푋−퐼푅(푥)) 퐹푃퐵푋−퐼푅(푥) = 휏휆 then we have 2 휃 exp[−훼2푥−2] 1 − exp⁡(−휆 [1 − exp (− { } )] ) ( 2 −2) 휆,휃,훼 1 − exp −훼 푥 퐹푃퐵푋−퐼푅(푥) = 휏휆 the corresponding PDF is 2 (휆,휃,훼) 4휃휆2훼 −(2+1) 2 −2 −3 푓푃퐵푟푋퐼푊(푥) = 푥 {1 − exp[−훼 푥 ]} 휏(휆) 2 exp[−훼2푥−2] ⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡× exp (−2(훼푥−1)2 − { } ) 1 − exp[−훼2푥−2] 2 휃−1 exp[−훼2푥−2] ⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡× [1 − exp (− { } )] 1 − exp[−훼2푥−2] 2 휃 exp[−훼2푥−2] ⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡× exp⁡{−휆 [1 − exp (− { } )] } 1 − exp[−훼2푥−2] The transformation of Ristić and Nadarajah (2014) used by many authors such as Ramos et al. (2015), Aryal and Yousof (2017) and Korkmaz et al. (2018), among others. Rania H. M. Abdelkhalek 59 (휆,휃,훼) The hazard rate function (HRF) can be easily calculated via 푓푃퐵푟푋퐼푊(푥)/[1 − (휆,휃,훼) 퐹푃퐵푋−퐼푅(푥)]. The PBX-IR density can be right-skewed, unimodal and symmetric (see Figure 1) whereas the PBX-IR HRF can be unimodal then bathtub, increasing and bathtub (see Figure 2). Figure 1: Plots of the PBX-IR PDF for selected parameter values. 60 The Poisson Burr X Inverse Rayleigh Distribution And Its Applications Figure 2: Plots of the PBX-IR HRF for selected parameter values. This article is organized as follows: In Section 2, we derive some mathematical properties of the new model. Maximum likelihood method for the model parameters is addressed in Section 3. Sections 4 presents the simulation studies. In Section 5, The potentiality of the proposed model is illustrated by means of three real data sets. Section 6 provides some concluding remarks. 2. Mathematical properties 2.1 Useful expansions Upon the power series ∞ 휁푞 exp(휁) = ∑ ( ) , 푞! 푞=0 Rania H. M. Abdelkhalek 61 the PDF in (6) can be written as ∞ 1+ℎ 2 ℎ (휆,휃,훼) 2휃휆 2훼 (−1) ⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡푓푃퐵푟푋퐼푊(푥) = ∑ ℎ! [휏(휆)] 푞=0 exp[−훼2푥−2] ⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡× exp (푥−(2+1) ) {1 − exp[−훼2푥−2]}3 2 exp[−훼2푥−2] ⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡× exp[− ( ) ] 1 − exp[−훼2푥−2] 2 휃(ℎ+1)−1 exp[−훼2푥−2] ⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡× [1 − exp (− { } )] 1 − exp[−훼2푥−2] If |휁| < 1 and 휏 > 0 is a real non-integer, the following power series holds ∞ (−1)Γ(1 + 휏) (1 − 휁)휏 = ∑ 휁 Γ(1 + 휏 − ) =0 Applying (8) to (7) we have 2 −(2+1) 2 −2 (휆,휃,훼) 2휃2훼 푥 exp[−2훼 푥 ] ⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡푓푃퐵푟−퐼푅(푥) = [휏(휆)] ∞ 휆1+ℎ(−1)ℎ+Γ(휃(⁡ℎ + 1)) ⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡× ∑ i! Γ(휃(ℎ + 1) − ) ℎ,=0 2 exp[−훼2푥−2] 푒푥푝[−( + 1) ( ) ] 1 − exp[−훼2푥−2] ⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡× [1 − exp[−훼2푥−2]]3 Via applying the power series to the term 2 exp[−훼2푥−2] exp [−( + 1) ( ) ], 1 − exp[−훼2푥−2] equation (9) becomes ∞ 4훼2휃휆1+ℎ(−1)ℎ++휔( + 1)휔푥−3 ⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡푓(휆,휃,훼) (푥) = ∑ 푃퐵푋−퐼푅 ! 휔! [휏 ] ℎ,,휔=0 (휆) Γ(θ(h + 1)) exp[−훼2푥−2] {exp[−훼2푥−2]}2휔+1 ⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡⁡× Γ(θ(h + 1) − i) {1 − exp[−훼2푥−2]}2휔+3 Then consider the series expansion 62 The Poisson Burr X Inverse Rayleigh Distribution And Its Applications ∞ Γ(푐 + 휅) (1 − 휁)−푐퐼 = ∑ 휁휅 |휁|<1,푐>0 휅! Γ(푐) =0 Applying the expansion in (11) to (10) for the term {1 − exp[−훼2푥−2]}2휔+3, equation (10) becomes 1+ℎ ℎ++휔 휔 (휆,휃,훼) ∞ 2휃휆 (−1) ⁡(+1) 푓푃퐵푋−퐼푅(푥) = ⁡ℎ,,휔,휅=0 ⁡ !⁡휔!휅![휏(휆)][2(1+휔)+휅] Γ(휃(ℎ+1))Γ(3+2휔+휅) × [2(1 + 휔) + 휅] Γ(휃(ℎ+1)−)Γ(2휔+3) × 2훼2푥−(2+1)exp{−[2(1 + 휔) + 휅](훼푥−1)2}. This can be written as ∞ (휆,휃,훼) 푓푃퐵푋−퐼푅(푥) = ∑ 휐휔,휅⁡ℎ[2(1+휔)+휅](푥; 훼), 휔,휅=0 where ∞ 2휃휆1+ℎ(−1)휔Γ(3 + 2휔 + 휅) (−1)ℎ+⁡Γ(휃(ℎ + 1))( + 1)휔 휐 = ⁡ ∑ ,⁡ 휔,휅 휔! 휅! [휏 ]Γ(2휔 + 3)[2(1 + 휔) + 휅] ! ⁡Γ(휃(ℎ + 1) − ) (휆) ℎ,=0 1 and ℎ[2(1+휔)+휅](푥; 훼) is the IR density with scale parameter 훼[2(1 + 휔) + 휅]2. Similarly, the CDF of the PBX-IR model can also be expressed as ∞ (휆,휃,훼) 퐹푃퐵푋−퐼푅(푥) = ∑ 휐휔,휅퐻[2(1+휔)+휅](푥; 훼), 휔,휅=0 1 where 퐻[2(1+휔)+휅](푥; 훼) is the the IR density with scale parameter 훼[2(1 + 휔) + 휅]2. 2.2 Quantile and random number generation The quantile function (QF) of 푋, where 푋 ∼PBX-IR (휆, 휃, 훼), is obtained by inverting (5) as 1 − 1 2 1 2 휃 −ln{1 − 푢[휏(휆)]} 푄(푢) = 훼 −ln 1 + {−ln [1 − ( ) ]} , 휆 { [( )]} so, simulating the PBX-IR r.v. is straightforward (see Section 4). If 푈 is a uniform variate on the unit interval (0,1), then the r.v. 푋 = 푄(푈) follows (6). 2.3 Moments 푡ℎ ′ The 푟 ordinary moment of 푋, say 휇푟, follows from (12) as Rania H. M. Abdelkhalek 63 ∞ 푟 푟 휇′ | = 퐄(푋푟) = ∑ 휐 훼푟[2(1 + 휔) + 휅]2⁡Γ (1 − ), 푟 (푟<2) 휔,휅 2 휔,휅=0 where ∞ 훕−1 훕 Γ(1 + 훕)|(훕∈ℝ+) =푤=0 (훕 − 푤) = ∫ ⁡푥 ⁡exp(−푡)푑푥 = 훕!, 0 upon setting 푟 = 1 in (15) gives the mean of 푋 as ∞ 1 1 퐄(푋) = ∑ 휐 훼[2(1 + 휔) + 휅]2⁡Γ (1 − ).
Recommended publications
  • New Proposed Length-Biased Weighted Exponential and Rayleigh Distribution with Application
    CORE Metadata, citation and similar papers at core.ac.uk Provided by International Institute for Science, Technology and Education (IISTE): E-Journals Mathematical Theory and Modeling www.iiste.org ISSN 2224-5804 (Paper) ISSN 2225-0522 (Online) Vol.4, No.7, 2014 New Proposed Length-Biased Weighted Exponential and Rayleigh Distribution with Application Kareema Abed Al-kadim1 Nabeel Ali Hussein 2 1,2.College of Education of Pure Sciences, University of Babylon/Hilla 1.E-mail of the Corresponding author: [email protected] 2.E-mail: [email protected] Abstract The concept of length-biased distribution can be employed in development of proper models for lifetime data. Length-biased distribution is a special case of the more general form known as weighted distribution. In this paper we introduce a new class of length-biased of weighted exponential and Rayleigh distributions(LBW1E1D), (LBWRD).This paper surveys some of the possible uses of Length - biased distribution We study the some of its statistical properties with application of these new distribution. Keywords: length- biased weighted Rayleigh distribution, length- biased weighted exponential distribution, maximum likelihood estimation. 1. Introduction Fisher (1934) [1] introduced the concept of weighted distributions Rao (1965) [2]developed this concept in general terms by in connection with modeling statistical data where the usual practi ce of using standard distributions were not found to be appropriate, in which various situations that can be modeled by weighted distributions, where non-observe ability of some events or damage caused to the original observation resulting in a reduced value, or in a sampling procedure which gives unequal chances to the units in the original That means the recorded observations cannot be considered as a original random sample.
    [Show full text]
  • Study of Generalized Lomax Distribution and Change Point Problem
    STUDY OF GENERALIZED LOMAX DISTRIBUTION AND CHANGE POINT PROBLEM Amani Alghamdi A Dissertation Submitted to the Graduate College of Bowling Green State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY August 2018 Committee: Arjun K. Gupta, Committee Co-Chair Wei Ning, Committee Co-Chair Jane Chang, Graduate Faculty Representative John Chen Copyright c 2018 Amani Alghamdi All rights reserved iii ABSTRACT Arjun K. Gupta and Wei Ning, Committee Co-Chair Generalizations of univariate distributions are often of interest to serve for real life phenomena. These generalized distributions are very useful in many fields such as medicine, physics, engineer- ing and biology. Lomax distribution (Pareto-II) is one of the well known univariate distributions that is considered as an alternative to the exponential, gamma, and Weibull distributions for heavy tailed data. However, this distribution does not grant great flexibility in modeling data. In this dissertation, we introduce a generalization of the Lomax distribution called Rayleigh Lo- max (RL) distribution using the form obtained by El-Bassiouny et al. (2015). This distribution provides great fit in modeling wide range of real data sets. It is a very flexible distribution that is related to some of the useful univariate distributions such as exponential, Weibull and Rayleigh dis- tributions. Moreover, this new distribution can also be transformed to a lifetime distribution which is applicable in many situations. For example, we obtain the inverse estimation and confidence intervals in the case of progressively Type-II right censored situation. We also apply Schwartz information approach (SIC) and modified information approach (MIC) to detect the changes in parameters of the RL distribution.
    [Show full text]
  • A New Generalization of the Generalized Inverse Rayleigh Distribution with Applications
    S S symmetry Article A New Generalization of the Generalized Inverse Rayleigh Distribution with Applications Rana Ali Bakoban * and Ashwaq Mohammad Al-Shehri Department of Statistics, College of Science, University of Jeddah, Jeddah 22254, Saudi Arabia; [email protected] * Correspondence: [email protected] Abstract: In this article, a new four-parameter lifetime model called the beta generalized inverse Rayleigh distribution (BGIRD) is defined and studied. Mixture representation of this model is derived. Curve’s behavior of probability density function, reliability function, and hazard function are studied. Next, we derived the quantile function, median, mode, moments, harmonic mean, skewness, and kurtosis. In addition, the order statistics and the mean deviations about the mean and median are found. Other important properties including entropy (Rényi and Shannon), which is a measure of the uncertainty for this distribution, are also investigated. Maximum likelihood estimation is adopted to the model. A simulation study is conducted to estimate the parameters. Four real-life data sets from difference fields were applied on this model. In addition, a comparison between the new model and some competitive models is done via information criteria. Our model shows the best fitting for the real data. Keywords: beta generalized inverse Rayleigh distribution; statistical properties; mean deviations; quantile; Rényi entropy; Shannon entropy; incomplete beta function; Montecarlo Simulation Citation: Bakoban, R.A.; Al-Shehri, MSC: Primary 62E10; Secondary 60E05 A.M. A New Generalization of the Generalized Inverse Rayleigh Distribution with Applications. Symmetry 2021, 13, 711. https:// 1. Introduction doi.org/10.3390/sym13040711 Modeling and analysis of lifetime phenomena are important aspects of statistical Academic Editor: Sergei D.
    [Show full text]
  • Location-Scale Distributions
    Location–Scale Distributions Linear Estimation and Probability Plotting Using MATLAB Horst Rinne Copyright: Prof. em. Dr. Horst Rinne Department of Economics and Management Science Justus–Liebig–University, Giessen, Germany Contents Preface VII List of Figures IX List of Tables XII 1 The family of location–scale distributions 1 1.1 Properties of location–scale distributions . 1 1.2 Genuine location–scale distributions — A short listing . 5 1.3 Distributions transformable to location–scale type . 11 2 Order statistics 18 2.1 Distributional concepts . 18 2.2 Moments of order statistics . 21 2.2.1 Definitions and basic formulas . 21 2.2.2 Identities, recurrence relations and approximations . 26 2.3 Functions of order statistics . 32 3 Statistical graphics 36 3.1 Some historical remarks . 36 3.2 The role of graphical methods in statistics . 38 3.2.1 Graphical versus numerical techniques . 38 3.2.2 Manipulation with graphs and graphical perception . 39 3.2.3 Graphical displays in statistics . 41 3.3 Distribution assessment by graphs . 43 3.3.1 PP–plots and QQ–plots . 43 3.3.2 Probability paper and plotting positions . 47 3.3.3 Hazard plot . 54 3.3.4 TTT–plot . 56 4 Linear estimation — Theory and methods 59 4.1 Types of sampling data . 59 IV Contents 4.2 Estimators based on moments of order statistics . 63 4.2.1 GLS estimators . 64 4.2.1.1 GLS for a general location–scale distribution . 65 4.2.1.2 GLS for a symmetric location–scale distribution . 71 4.2.1.3 GLS and censored samples .
    [Show full text]
  • Slashed Rayleigh Distribution
    Revista Colombiana de Estadística January 2015, Volume 38, Issue 1, pp. 31 to 44 DOI: http://dx.doi.org/10.15446/rce.v38n1.48800 Slashed Rayleigh Distribution Distribución Slashed Rayleigh Yuri A. Iriarte1;a, Héctor W. Gómez2;b, Héctor Varela2;c, Heleno Bolfarine3;d 1Instituto Tecnológico, Universidad de Atacama, Copiapó, Chile 2Departamento de Matemáticas, Facultad de Ciencias Básicas, Universidad de Antofagasta, Antofagasta, Chile 3Departamento de Estatística, Instituto de Matemática y Estatística, Universidad de Sao Paulo, Sao Paulo, Brasil Abstract In this article we study a subfamily of the slashed-Weibull family. This subfamily can be seen as an extension of the Rayleigh distribution with more flexibility in terms of the kurtosis of distribution. This special feature makes the extension suitable for fitting atypical observations. It arises as the ratio of two independent random variables, the one in the numerator being a Rayleigh distribution and a power of the uniform distribution in the denominator. We study some probability properties, discuss maximum likelihood estimation and present real data applications indicating that the slashed-Rayleigh distribution can improve the ordinary Rayleigh distribution in fitting real data. Key words: Kurtosis, Rayleigh Distribution, Slashed-elliptical Distribu- tions, Slashed-Rayleigh Distribution, Slashed-Weibull Distribution, Weibull Distribution. Resumen En este artículo estudiamos una subfamilia de la familia slashed-Weibull. Esta subfamilia puede ser vista como una extensión de la distribución Ray- leigh con más flexibilidad en cuanto a la kurtosis de la distribución. Esta particularidad hace que la extensión sea adecuada para ajustar observa- ciones atípicas. Esto surge como la razón de dos variables aleatorias in- dependientes, una en el numerador siendo una distribución Rayleigh y una aLecturer.
    [Show full text]
  • Hand-Book on STATISTICAL DISTRIBUTIONS for Experimentalists
    Internal Report SUF–PFY/96–01 Stockholm, 11 December 1996 1st revision, 31 October 1998 last modification 10 September 2007 Hand-book on STATISTICAL DISTRIBUTIONS for experimentalists by Christian Walck Particle Physics Group Fysikum University of Stockholm (e-mail: [email protected]) Contents 1 Introduction 1 1.1 Random Number Generation .............................. 1 2 Probability Density Functions 3 2.1 Introduction ........................................ 3 2.2 Moments ......................................... 3 2.2.1 Errors of Moments ................................ 4 2.3 Characteristic Function ................................. 4 2.4 Probability Generating Function ............................ 5 2.5 Cumulants ......................................... 6 2.6 Random Number Generation .............................. 7 2.6.1 Cumulative Technique .............................. 7 2.6.2 Accept-Reject technique ............................. 7 2.6.3 Composition Techniques ............................. 8 2.7 Multivariate Distributions ................................ 9 2.7.1 Multivariate Moments .............................. 9 2.7.2 Errors of Bivariate Moments .......................... 9 2.7.3 Joint Characteristic Function .......................... 10 2.7.4 Random Number Generation .......................... 11 3 Bernoulli Distribution 12 3.1 Introduction ........................................ 12 3.2 Relation to Other Distributions ............................. 12 4 Beta distribution 13 4.1 Introduction .......................................
    [Show full text]
  • Journal of Science Log-Gamma
    Research Article GU J Sci 31(3): 967-983 (2018) Gazi University Journal of Science http://dergipark.gov.tr/gujs Log-Gamma - Rayleigh Distribution: Properties and Applications Samir ASHOUR 1, Darwish DARWISH 1, Ahmad SAAD 1,* 1 Institute of Statistical Studies & Research, Cairo University, 5 Ahmed Zowil St., Orman, Giza, Egypt Article Info Abstract In this paper, for the first time we defined and studied a new two parameter lifetime model by Received: 14/08/2017 using the T-X method, called the Log-Gamma Rayleigh distribution. This distribution can be Accepted: 20/04/2018 considered as a new generalization of the Gamma distribution and the Rayleigh distribution. We obtain some of its mathematical properties. Some structural properties of the new distribution are studied. Maximum likelihood estimation method is used for estimating the model Keywords parameters. An application to real data set is given to show the flexibility and potentiality of the T-X method new model. Gamma distribution Rayleigh distribution Quantile function Hazard function Moments Maximum likelihood estimation 1. INTRODUCTION Eugene et al. [6] for the first time introduced the beta-generated family of distributions. They noted that the probability density function pdf of the beta random variable and the cumulative distribution function CDF of any distribution are between 0 and 1. The beta-generated random variable X is defined with the following CDF and pdf Fx() G()() x b t dt (1) 0 and 1 g() x f () x F11 ()[1 x F ()] x B(,) where b(t) is the pdf of the beta random variable with parameters and , F(x) and f(x) are the CDF and the pdf of any random variable.
    [Show full text]
  • Slashed Exponentiated Rayleigh Distribution
    Revista Colombiana de Estadística July 2015, Volume 38, Issue 2, pp. 453 to 466 DOI: http://dx.doi.org/10.15446/rce.v38n2.51673 Slashed Exponentiated Rayleigh Distribution Distribución recortada Rayleigh exponenciada Hugo S. Salinas1;a, Yuri A. Iriarte2;b, Heleno Bolfarine3;c 1Departamento de Matemática, Facultad de Ingeniería, Universidad de Atacama, Copiapó, Chile 2Instituto Tecnológico, Universidad de Atacama, Copiapó, Chile 3Departamento de Estatística, IME, Universidad de Sao Paulo, Sao Paulo, Brasil Abstract In this paper we introduce a new distribution for modeling positive data with high kurtosis. This distribution can be seen as an extension of the exponentiated Rayleigh distribution. This extension builds on the quotient of two independent random variables, one exponentiated Rayleigh in the numerator and Beta(q; 1) in the denominator with q > 0. It is called the slashed exponentiated Rayleigh random variable. There is evidence that the distribution of this new variable can be more flexible in terms of modeling the kurtosis regarding the exponentiated Rayleigh distribution. The properties of this distribution are studied and the parameter estimates are calculated using the maximum likelihood method. An application with real data re- veals good performance of this new distribution. Key words: Exponentiated Rayleigh Distribution, Kurtosis, Maximum Likelihood, Rayleigh Distribution, Slash Distribution. Resumen En este trabajo presentamos una nueva distribución para modelizar datos positivos con alta curtosis. Esta distribución puede ser vista como una exten- sión de la distribución Rayleigh exponenciada. Esta extensión se construye en base al cuociente de dos variables aleatorias independientes, una Raileigh exponenciada en el numerador y una Beta(q; 1) en el denominador con q > 0.
    [Show full text]
  • Field Guide to Continuous Probability Distributions
    Field Guide to Continuous Probability Distributions Gavin E. Crooks v 1.0.0 2019 G. E. Crooks – Field Guide to Probability Distributions v 1.0.0 Copyright © 2010-2019 Gavin E. Crooks ISBN: 978-1-7339381-0-5 http://threeplusone.com/fieldguide Berkeley Institute for Theoretical Sciences (BITS) typeset on 2019-04-10 with XeTeX version 0.99999 fonts: Trump Mediaeval (text), Euler (math) 271828182845904 2 G. E. Crooks – Field Guide to Probability Distributions Preface: The search for GUD A common problem is that of describing the probability distribution of a single, continuous variable. A few distributions, such as the normal and exponential, were discovered in the 1800’s or earlier. But about a century ago the great statistician, Karl Pearson, realized that the known probabil- ity distributions were not sufficient to handle all of the phenomena then under investigation, and set out to create new distributions with useful properties. During the 20th century this process continued with abandon and a vast menagerie of distinct mathematical forms were discovered and invented, investigated, analyzed, rediscovered and renamed, all for the purpose of de- scribing the probability of some interesting variable. There are hundreds of named distributions and synonyms in current usage. The apparent diver- sity is unending and disorienting. Fortunately, the situation is less confused than it might at first appear. Most common, continuous, univariate, unimodal distributions can be orga- nized into a small number of distinct families, which are all special cases of a single Grand Unified Distribution. This compendium details these hun- dred or so simple distributions, their properties and their interrelations.
    [Show full text]
  • On the Arcsecant Hyperbolic Normal Distribution. Properties, Quantile Regression Modeling and Applications
    S S symmetry Article On the Arcsecant Hyperbolic Normal Distribution. Properties, Quantile Regression Modeling and Applications Mustafa Ç. Korkmaz 1 , Christophe Chesneau 2,* and Zehra Sedef Korkmaz 3 1 Department of Measurement and Evaluation, Artvin Çoruh University, City Campus, Artvin 08000, Turkey; [email protected] 2 Laboratoire de Mathématiques Nicolas Oresme, University of Caen-Normandie, 14032 Caen, France 3 Department of Curriculum and Instruction Program, Artvin Çoruh University, City Campus, Artvin 08000, Turkey; [email protected] * Correspondence: [email protected] Abstract: This work proposes a new distribution defined on the unit interval. It is obtained by a novel transformation of a normal random variable involving the hyperbolic secant function and its inverse. The use of such a function in distribution theory has not received much attention in the literature, and may be of interest for theoretical and practical purposes. Basic statistical properties of the newly defined distribution are derived, including moments, skewness, kurtosis and order statistics. For the related model, the parametric estimation is examined through different methods. We assess the performance of the obtained estimates by two complementary simulation studies. Also, the quantile regression model based on the proposed distribution is introduced. Applications to three real datasets show that the proposed models are quite competitive in comparison to well-established models. Keywords: bounded distribution; unit hyperbolic normal distribution; hyperbolic secant function; normal distribution; point estimates; quantile regression; better life index; dyslexia; IQ; reading accu- racy modeling Citation: Korkmaz, M.Ç.; Chesneau, C.; Korkmaz, Z.S. On the Arcsecant Hyperbolic Normal Distribution. 1. Introduction Properties, Quantile Regression Over the past twenty years, many statisticians and researchers have focused on propos- Modeling and Applications.
    [Show full text]
  • Weighted Inverse Rayleigh Distribution
    International Journal of Statistics and Systems ISSN 0973-2675 Volume 12, Number 1 (2017), pp. 119-137 © Research India Publications http://www.ripublication.com Weighted Inverse Rayleigh Distribution Kawsar Fatima and S.P Ahmad Department of Statistics, University of Kashmir, Srinagar, India. Corresponding author Abstract In this paper, we have introduced weighted inverse Rayleigh (WIR) distribution and investigated its different statistical properties. Expressions for the Mode and entropy have also been derived. In addition, it also contains some special cases that are well known. Moreover, we apply the maximum likelihood method to estimate the parameter , and applications to two real data sets show the superiority of this new distribution by comparing the fitness with its special cases. Keywords: Inverse Rayleigh distribution, weighted distribution, Reliability Analysis, Entropy, Maximum likelihood estimation, Real life data sets. 1. INTRODUCTION The IR distribution was proposed by Voda (1972). He studies some properties of the MLE of the scale parameter of inverse Rayleigh distribution which is also being used in lifetime experiments. If X has IR distribution, its probability density function (pdf) takes the following form: 2 2 g() x ex ; x 0, 0 )1.1( x3 The corresponding cumulative distribution function is 2 G( x ) ex ; x 0, 0 )2.1( Where x 0 the scale parameter 0 The kth moment of the IR distribution is given as the following: 120 Kawsar Fatima and S.P Ahmad 2 k k 2 x Eg () x x e dx. 0 x3 1 1 2 Making the substitution y , dx dy , so that x 2/1 ,we obtain x 2 x3 y y E() xk y(1k / 2) 1 e dy.
    [Show full text]
  • Lecture 17: Bivariate Normal Distribution F (X, Y) Is Given By
    Section 5.3 Bivariate Unit Normal Bivariate Unit Normal If X and Y have independent unit normal distributions then their joint Lecture 17: Bivariate Normal distribution f (x; y) is given by 2 − 1 (x2+y 2) f (x; y) = φ(x)φ(y) = c e 2 Sta230 / Mth230 Colin Rundel April 2, 2014 Sta230 / Mth230 (Colin Rundel) Lecture 17 April 2, 2014 1 / 20 Section 5.3 Bivariate Unit Normal Section 5.3 Bivariate Unit Normal Bivariate Unit Normal, cont. Bivariate Unit Normal - Normalizing Constant We can rewrite the joint distribution in terms of the distance r from the First, rewriting the joint distribution in this way also makes it possible to origin evaluate the constant of integration, c. From our previous discussions of the normal distribution we know that p 2 2 r = x + y c = p1 , but we never saw how to derive it explicitly. 2π 2 − 1 (x2+y 2) 2 − 1 r 2 f (x; y) = c e 2 = c e 2 This tells us something useful about this special case of the bivariate normal distributions: it is rotationally symmetric about the origin. r r + eps This particular fact is incredibly powerful and helps us solve a variety of Y problems. Sta230 / Mth230 (Colin Rundel) Lecture 17 April 2, 2014 2 / 20 Sta230 / Mth230 (Colin Rundel)X Lecture 17 April 2, 2014 3 / 20 Section 5.3 Bivariate Unit Normal Section 5.3 Bivariate Unit Normal Bivariate Unit Normal - Normalizing Constant Rayleigh Distribution The distance from the origin of a point (x; y) derived from X ; Y ∼ N (0; 1) is called the Rayleigh distribution.
    [Show full text]