The Non-Central Nakagami Distribution

Total Page:16

File Type:pdf, Size:1020Kb

The Non-Central Nakagami Distribution The non-central Nakagami distribution Søren Hauberg Section for Cognitive Systems Technical University of Denmark [email protected] Abstract The Nakagami distribution describe the square root of a random variable drawn from a Gamma distribution. Equivalently, the Gamma distribution can be seen as a one-dimensional Wishart distribution. In this note, we consider the distribution of the square root of a random variable drawn from a non-central one-dimensional Wishart distribution. We present the probability density function of this distribution along with closed-form expressions for its moments. This document is just a short technical note. Comments can be directed at [email protected] 1 The Nakagami distribution 2 Let zd ∼ N (0; σ ); d = 1;:::;D be iid samples from a zero-mean normal distribution. Then D X 2 + x = zd 2 R (1.1) d=1 follows a Gamma distribution. Equivalently, and more suitable for our purposes, x also follows a one-dimensional Wishart distribution [1], 2 x ∼ W1(D; σ ) (1.2) r D−2 2 x −1 x p(x) = W (x j D; σ ) = Γ(D=2) exp − : (1.3) 1 2D σ2D 2σ2 This distribution can be seen as describing the squaredp norm of a zero-mean normally distributed vector. The norm of this vector is then given by y = x, which follows a Nakagami distribution [2], 2 y ∼ Nakagami D=2; Dσ (1.4) 2 2 2 D−1 y D 2 p p(y) = Nakagami y j = ; Dσ = y exp − 2 : (1.5) Γ(D=2)( 2σ)D 2σ This expression is easily derived by the change of value theorem; see next section. The expectation and variance of a Nakagami distributed variable is D+1 p Γ 2 E [y] = D 2σ (1.6) Γ 2 2! Γ((D + 1)=2) var [y] = σ2 D − 2 : (1.7) Γ(D=2) We will derive a more general version of these results in the next section. Preprint. Work in progress. 2 The non-central Nakagami distribution 2 Now, let zd ∼ N (µd; σ ); d = 1;:::;D be independent samples from normal distributions with variance σ2. Then D X 2 + x = zd 2 R (2.1) d=1 follows a one-dimensional non-central Wishart distribution [1], 2 x ∼ W1(D; σ ; Ω) (2.2) PD µ2 Ω = d=1 d (2.3) σ2 2 p(x) = W1(x j D; σ ; Ω) r D−2 (2.4) x −1 −2 x Ω = Γ(D=2) F D=2; 1=4Ωσ x exp − exp − : 2D σ2D 0 1 2σ2 2 Here 0F1 is a generalized hypergeometric function; see Sec. 7.3 of Muirhead’s book [1]. This distribution can be seen as describing the squared norm ofp a normal distributed vector with isotropic variance. The norm of this vector is then given by y = x, which we say follows a non-central Nakagami distribution. This distribution can be derived by change-of-variables to be 2 2 p (y) = W1(y j D; σ ; Ω) · 2y (2.5) yD−1 D y2 y2 Ω = (D−2) 0F1 ; Ω exp − exp − : (2.6) 2 =2 σD Γ(D=2) 2 4σ2 2σ2 2 p To compute the moments of x we recall the following result from Muirhead [1, Theorem 10.3.7]. q×q Let X 2 R and X ∼ Wq(D; Σ; Ω) then D k k qk Γq( =2 + k) (det(X)) = (det(Σ)) 2 F (−k; D=2; −1=2Ω); (2.7) E D 1 1 Γq( =2) p where 1F1 is a generalized hypergeometric function. Since y = x is a positive scalar, its determinant is merely y, and we get hp 2ki Γ(D=2 + k) 2k k D 1 E x = σ 2 1F1(−k; =2; − =2Ω): (2.8) Γ(D=2) p From this we see that the mean and the variance of y = x is p Γ((D+1)=2) 1=2 1 D 1 E x = σ2 1F1(− =2; =2; − =2Ω) (2.9) Γ(D=2) Γ((D+2)=2) 2 D 1 E [x] = σ 2 1F1(−1; =2; − =2Ω) (2.10) Γ(D=2) p p 2 var x = E x − E x (2.11) σ22 D + 2 D 1 = Γ 1F1(−1; =2; − =2Ω) Γ(D=2) 2 (2.12) 2 ! Γ((D+1)=2) 1 D 1 2 − 1F1(− =2; =2; − =2Ω) : Γ(D=2) Acknowledgments. The author was supported by a research grant (15334) from VILLUM FONDEN. This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement no 757360). References [1] R. J. Muirhead. Aspects of Multivariate Statistical Theory. John Wiley & Sons, 2005. [2] M. Nakagami. The m-distribution, a general formula of intensity of rapid fading, iwc hoffman, ed, 1960. 2.
Recommended publications
  • A New Family of Odd Generalized Nakagami (Nak-G) Distributions
    TURKISH JOURNAL OF SCIENCE http:/dergipark.gov.tr/tjos VOLUME 5, ISSUE 2, 85-101 ISSN: 2587–0971 A New Family of Odd Generalized Nakagami (Nak-G) Distributions Ibrahim Abdullahia, Obalowu Jobb aYobe State University, Department of Mathematics and Statistics bUniversity of Ilorin, Department of Statistics Abstract. In this article, we proposed a new family of generalized Nak-G distributions and study some of its statistical properties, such as moments, moment generating function, quantile function, and prob- ability Weighted Moments. The Renyi entropy, expression of distribution order statistic and parameters of the model are estimated by means of maximum likelihood technique. We prove, by providing three applications to real-life data, that Nakagami Exponential (Nak-E) distribution could give a better fit when compared to its competitors. 1. Introduction There has been recent developments focus on generalized classes of continuous distributions by adding at least one shape parameters to the baseline distribution, studying the properties of these distributions and using these distributions to model data in many applied areas which include engineering, biological studies, environmental sciences and economics. Numerous methods for generating new families of distributions have been proposed [8] many researchers. The beta-generalized family of distribution was developed , Kumaraswamy generated family of distributions [5], Beta-Nakagami distribution [19], Weibull generalized family of distributions [4], Additive weibull generated distributions [12], Kummer beta generalized family of distributions [17], the Exponentiated-G family [6], the Gamma-G (type I) [21], the Gamma-G family (type II) [18], the McDonald-G [1], the Log-Gamma-G [3], A new beta generated Kumaraswamy Marshall-Olkin- G family of distributions with applications [11], Beta Marshall-Olkin-G family [2] and Logistic-G family [20].
    [Show full text]
  • Dictionary-Based Stochastic Expectation– Maximization for SAR
    188 IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING, VOL. 44, NO. 1, JANUARY 2006 Dictionary-Based Stochastic Expectation– Maximization for SAR Amplitude Probability Density Function Estimation Gabriele Moser, Member, IEEE, Josiane Zerubia, Fellow, IEEE, and Sebastiano B. Serpico, Senior Member, IEEE Abstract—In remotely sensed data analysis, a crucial problem problem as a parameter estimation problem. Several strategies is represented by the need to develop accurate models for the have been proposed in the literature to deal with parameter statistics of the pixel intensities. This paper deals with the problem estimation, e.g., the maximum-likelihood methodology [18] of probability density function (pdf) estimation in the context of and the “method of moments” [41], [46], [53]. On the contrary, synthetic aperture radar (SAR) amplitude data analysis. Several theoretical and heuristic models for the pdfs of SAR data have nonparametric pdf estimation approaches do not assume any been proposed in the literature, which have been proved to be specific analytical model for the unknown pdf, thus providing effective for different land-cover typologies, thus making the a higher flexibility, although usually involving internal archi- choice of a single optimal parametric pdf a hard task, especially tecture parameters to be set by the user [18]. In particular, when dealing with heterogeneous SAR data. In this paper, an several nonparametric kernel-based estimation and regression innovative estimation algorithm is described, which faces such a problem by adopting a finite mixture model for the amplitude architectures have been described in the literature, that have pdf, with mixture components belonging to a given dictionary of proved to be effective estimation tools, such as standard Parzen SAR-specific pdfs.
    [Show full text]
  • Estimation and Testing Procedures of the Reliability Functions of Nakagami Distribution
    Austrian Journal of Statistics January 2019, Volume 48, 15{34. AJShttp://www.ajs.or.at/ doi:10.17713/ajs.v48i3.827 Estimation and Testing Procedures of the Reliability Functions of Nakagami Distribution Ajit Chaturvedi Bhagwati Devi Rahul Gupta Dept. of Statistics, Dept. of Statistics, Dept. of Statistics, University of Delhi University of Jammu University of Jammu Abstract A very important distribution called Nakagami distribution is taken into considera- tion. Reliability measures R(t) = P r(X > t) and P = P r(X > Y ) are considered. Point as well as interval procedures are obtained for estimation of parameters. Uniformly Mini- mum Variance Unbiased Estimators (U.M.V.U.Es) and Maximum Likelihood Estimators (M.L.Es) are developed for the said parameters. A new technique of obtaining these esti- mators is introduced. Moment estimators for the parameters of the distribution have been found. Asymptotic confidence intervals of the parameter based on M.L.E and log(M.L.E) are also constructed. Then, testing procedures for various hypotheses are developed. At the end, Monte Carlo simulation is performed for comparing the results obtained. A real data analysis is performed to describe the procedure clearly. Keywords: Nakagami distribution, point estimation, testing procedures, confidence interval, Markov chain Monte Carlo (MCMC) procedure.. 1. Introduction and preliminaries R(t), the reliability function is the probability that a system performs its intended function without any failure at time t under the prescribed conditions. So, if we suppose that life- time of an item or any system is denoted by the random variate X then the reliability is R(t) = P r(X > t).
    [Show full text]
  • Newdistns: an R Package for New Families of Distributions
    JSS Journal of Statistical Software March 2016, Volume 69, Issue 10. doi: 10.18637/jss.v069.i10 Newdistns: An R Package for New Families of Distributions Saralees Nadarajah Ricardo Rocha University of Manchester Universidade Federal de São Carlos Abstract The contributed R package Newdistns written by the authors is introduced. This pack- age computes the probability density function, cumulative distribution function, quantile function, random numbers and some measures of inference for nineteen families of distri- butions. Each family is flexible enough to encompass a large number of structures. The use of the package is illustrated using a real data set. Also robustness of random number generation is checked by simulation. Keywords: cumulative distribution function, probability density function, quantile function, random numbers. 1. Introduction Let G be any valid cumulative distribution function defined on the real line. The last decade or so has seen many approaches proposed for generating new distributions based on G. All of these approaches can be put in the form F (x) = B (G(x)) , (1) where B : [0, 1] → [0, 1] and F is a valid cumulative distribution function. So, for every G one can use (1) to generate a new distribution. The first approach of the kind of (1) proposed in recent years was that due to Marshall and Olkin(1997). In Marshall and Olkin(1997), B was taken to be B(p) = βp/ {1 − (1 − β)p} for β > 0. The distributions generated in this way using (1) will be referred to as Marshall Olkin G distributions. Since Marshall and Olkin(1997), many other approaches have been proposed.
    [Show full text]
  • University of Nevada, Reno Generalized Univariate
    University of Nevada, Reno Generalized Univariate Distributions and a New Asymmetric Laplace Model A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Mathematics By Palash Sharma Dr. Tomasz Kozubowski/Thesis Advisor August, 2017 c 2017 Palash Sharma ALL RIGHTS RESERVED THE GRADUATE SCHOOL We recommend that the thesis prepared under our supervision by Palash Sharma entitled Generalized Univariate Distributions and a New Asymmetric Laplace Model be accepted in partial fulfillment of the requirements for the degree of MASTER OF SCIENCE Tomasz J. Kozubowski, Ph.D., Advisor Anna Panorska, Ph.D., Committee Member Minggen Lu, Ph.D., Graduate School Representative David Zeh, Ph.D., Dean, Graduate School August, 2017 i ABSTRACT Generalized Univariate Distributions and a New Asymmetric Laplace Model By Palash Sharma This work provides a survey of general class of distributions generated from a mixture of beta random variables. We provide an extensive review of the literature, concern- ing generating new distributions via the inverse CDF transformation. In particular, we account for beta generated and Kumaraswamy generated families of distributions. We provide a brief summary of each of their families of distributions. We also propose a new asymmetric mixture distribution, which is an alternative to beta generated dis- tributions. We provide basic properties of this new class of distributions generated from the Laplace model. We also address the issue of parameter estimation of this new skew generalized Laplace model. ii ACKNOWLEDGMENTS At first, I would like to thank my honorable thesis advisor, Professor Tomasz J. Kozubowski, who showed me a great interest in the field of statistics and probability theory.
    [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]
  • Probability Models and Statistical Tests for Extreme Precipitation Based on Generalized Negative Binomial Distributions
    mathematics Article Probability Models and Statistical Tests for Extreme Precipitation Based on Generalized Negative Binomial Distributions Victor Korolev 1,2,3,4 and Andrey Gorshenin 1,2,3,* 1 Moscow Center for Fundamental and Applied Mathematics, Lomonosov Moscow State University, 119991 Moscow, Russia; [email protected] 2 Faculty of Computational Mathematics and Cybernetics, Lomonosov Moscow State University, 119991 Moscow, Russia 3 Federal Research Center “Computer Science and Control” of the Russian Academy of Sciences, 119333 Moscow, Russia 4 Department of Mathematics, School of Science, Hangzhou Dianzi University, Hangzhou 310018, China * Correspondence: [email protected] Received: 4 April 2020; Accepted: 14 April 2020; Published: 16 April 2020 Abstract: Mathematical models are proposed for statistical regularities of maximum daily precipitation within a wet period and total precipitation volume per wet period. The proposed models are based on the generalized negative binomial (GNB) distribution of the duration of a wet period. The GNB distribution is a mixed Poisson distribution, the mixing distribution being generalized gamma (GG). The GNB distribution demonstrates excellent fit with real data of durations of wet periods measured in days. By means of limit theorems for statistics constructed from samples with random sizes having the GNB distribution, asymptotic approximations are proposed for the distributions of maximum daily precipitation volume within a wet period and total precipitation volume for a wet period. It is shown that the exponent power parameter in the mixing GG distribution matches slow global climate trends. The bounds for the accuracy of the proposed approximations are presented. Several tests for daily precipitation, total precipitation volume and precipitation intensities to be abnormally extremal are proposed and compared to the traditional PoT-method.
    [Show full text]
  • Concentration Inequalities from Likelihood Ratio Method
    Concentration Inequalities from Likelihood Ratio Method ∗ Xinjia Chen September 2014 Abstract We explore the applications of our previously established likelihood-ratio method for deriving con- centration inequalities for a wide variety of univariate and multivariate distributions. New concentration inequalities for various distributions are developed without the idea of minimizing moment generating functions. Contents 1 Introduction 3 2 Likelihood Ratio Method 4 2.1 GeneralPrinciple.................................... ...... 4 2.2 Construction of Parameterized Distributions . ............. 5 2.2.1 WeightFunction .................................... .. 5 2.2.2 ParameterRestriction .............................. ..... 6 3 Concentration Inequalities for Univariate Distributions 7 3.1 BetaDistribution.................................... ...... 7 3.2 Beta Negative Binomial Distribution . ........ 7 3.3 Beta-Prime Distribution . ....... 8 arXiv:1409.6276v1 [math.ST] 1 Sep 2014 3.4 BorelDistribution ................................... ...... 8 3.5 ConsulDistribution .................................. ...... 8 3.6 GeetaDistribution ................................... ...... 9 3.7 GumbelDistribution.................................. ...... 9 3.8 InverseGammaDistribution. ........ 9 3.9 Inverse Gaussian Distribution . ......... 10 3.10 Lagrangian Logarithmic Distribution . .......... 10 3.11 Lagrangian Negative Binomial Distribution . .......... 10 3.12 Laplace Distribution . ....... 11 ∗The author is afflicted with the Department of Electrical
    [Show full text]
  • MPS: an R Package for Modelling New Families of Distributions
    MPS: An R package for modelling new families of distributions Mahdi Teimouri Department of Statistics Gonbad Kavous University Gonbad Kavous, IRAN Abstract: We introduce an R package, called MPS, for computing the probability density func- tion, computing the cumulative distribution function, computing the quantile function, simulating random variables, and estimating the parameters of 24 new shifted families of distributions. By considering an extra shift (location) parameter for each family more flexibility yields. Under some situations, since the maximum likelihood estimators may fail to exist, we adopt the well-known maximum product spacings approach to estimate the parameters of shifted 24 new families of distributions. The performance of the MPS package for computing the cdf, pdf, and simulating random samples will be checked by examples. The performance of the maximum product spacings approach is demonstrated by executing MPS package for three sets of real data. As it will be shown, for the first set, the maximum likelihood estimators break down but MPS package find them. For the second set, adding the location parameter leads to acceptance the model while absence of the location parameter makes the model quite inappropriate. For the third set, presence of the location parameter yields a better fit. Keywords: Cumulative distribution function; Maximum likelihood estimation; Method of maxi- mum product spacings; Probability density function; Quantile function; R package; Simulation; 1 Introduction Over the last two decades, generalization of the statistical distributions has attracted much at- tention in the literature. Most of these extensions have been spawned by applications found in arXiv:1809.02959v1 [stat.CO] 9 Sep 2018 analyzing lifetime data.
    [Show full text]
  • The Estimation of the M Parameter of the Nakagami Distribution
    WSEAS TRANSACTIONS on BIOLOGY and BIOMEDICINE Li-Fei Huang, Jen-Jen Lin The estimation of the m parameter of the Nakagami distribution Li-Fei Huang Jen-Jen Lin Ming Chuan University Ming Chuan University Dept. of Applied Stat. and Info. Science Dept. of Applied Stat. and Info. Science 5 Teh-Ming Rd., Gwei-Shan Taoyuan City 5 Teh-Ming Rd., Gwei-Shan Taoyuan City Taiwan Taiwan [email protected] [email protected] Abstract: This paper introduces the Nakagami-m distribution which is usually used to simulate the ultrasound image. The gamma distribution is used to derive the moment estimator and the maximum likelihood estimator because the Nakagami distribution has no moment generation function and too complicate likelihood function. The moment estimator of m is normally distributed with a smaller bias and a larger standard deviation. The second order maximum likelihood estimator of m is also normally distributed with a larger bias and a smaller standard deviation. The confidence interval for the ratio of medians from two independent distributions of Nakagami-m estimators is constructed. The moment estimator provides a quick understanding about the m parameter, while the second order maximum likelihood estimator provides a full understanding about the m parameter. Key–Words: The Nakagami distribution, The moment estimator, The maximum likelihood estimator, The distribu- tion of the Nakagami-m, The ratio of medians. 1 Introduction agami distribution is as follows. 2 m m 2m−1 − m n2 f(n) = ( ) n e Ω The Nakagami distribution is usually used to simu- Γ(m) Ω late the ultrasound image.
    [Show full text]
  • Assessing Probabilistic Modelling for Wind Speed from Numerical Weather Prediction Model and Observation in the Arctic
    www.nature.com/scientificreports OPEN Assessing probabilistic modelling for wind speed from numerical weather prediction model and observation in the Arctic Hao Chen1*, Yngve Birkelund1, Stian Normann Anfnsen2, Reidar Staupe‑Delgado1 & Fuqing Yuan1 Mapping Arctic renewable energy resources, particularly wind, is important to ensure the transition into renewable energy in this environmentally vulnerable region. The statistical characterisation of wind is critical for efectively assessing energy potential and planning wind park sites and is, therefore, an important input for wind power policymaking. In this article, diferent probability density functions are used to model wind speed for fve wind parks in the Norwegian Arctic region. A comparison between wind speed data from numerical weather prediction models and measurements is made, and a probability analysis for the wind speed interval corresponding to the rated power, which is largely absent in the existing literature, is presented. The results of the present study suggest that no single probability function outperforms across all scenarios. However, some diferences emerged from the models when applied to diferent wind parks. The Nakagami and Generalised extreme value distributions were chosen for the numerical weather predicted prediction and the observed wind speed modelling, respectively, due to their superiority and stability compared with other methods. This paper, therefore, provides a novel direction for understanding the numerical weather prediction wind model and shows that its speed statistical features are better captured than those of real wind. With the growing reliance on renewable energy resources in many regions of the world, studying the predictabil- ity of renewable energy is becoming progressively important1. As one of the cleanest renewable energy sources, wind energy has attracted growing attention worldwide 2.
    [Show full text]
  • The Generalized Odd Nakagami-G Family of Distributions: Properties and Applications
    NATURENGS, MTU Journal of Engineering and Natural Sciences 1:2 (2020) 1-16 NATURENGS MTU Journal of Engineering and Natural Sciences https://dergipark.org.tr/tr/pub/naturengs DOI: 10.46572/nat.2020.12 The Generalized Odd Nakagami-G Family of Distributions: Properties and Applications Ibrahim ABDULLAHI1*, Job OBALOWU2 1Department of Mathematics and Statistics, Faculty of Science, Yobe State University, 620, Yobe, Nigeria. 2Department of Statistics, Faculty of Physical Science, University of Ilorin, 240003, Kwara, Nigeria. (Received: 15.08.2020; Accepted: 01.10.2020) ABSTRACT: In this study, the Generalized Odd Nakagami-G distribution has been empirically investigated using Frechet distribution as the baseline distribution. Some mathematical statistics properties viz: quantile function, moments, probability-weighted moment, entropies, and order statistics were derived among others. The method of maximum likelihood estimation is suitable for the derivation of estimators of the distribution parameters. And the Generalized odd Nakagami Frechet (GONak-Fr) distribution gives the best fit vis-a-vis its competitors via application to real life data sets. Keywords: GONak-G, GONak-Fr, Probability-weighted moment, Entropies, Order statistics. 1. INTRODUCTION There have been recent developments focus on generalized odd classes of continuous distributions by adding at least one shape parameters to the baseline distribution, has led mathematical statistician researcher to make of the new flexible distributions and studying the properties of these distributions
    [Show full text]