DOCSLIB.ORG

On the Scale Parameter of Exponential Distribution

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
On the Scale Parameter of Exponential Distribution Review of the Air Force Academy No.2 (34)/2017 ON THE SCALE PARAMETER OF EXPONENTIAL DISTRIBUTION Anca Ileana LUPAŞ Military Technical Academy, Bucharest, Romania ([email protected]) DOI: 10.19062/1842-9238.2017.15.2.16 Abstract: Exponential distribution is one of the widely used continuous distributions in various fields for statistical applications. In this paper we study the exact and asymptotical distribution of the scale parameter for this distribution. We will also define the confidence intervals for the studied parameter as well as the fixed length confidence intervals. 1. INTRODUCTION Exponential distribution is used in various statistical applications. Therefore, we often encounter exponential distribution in applications such as: life tables, reliability studies, extreme values analysis and others. In the following paper, we focus our attention on the exact and asymptotical repartition of the exponential distribution scale parameter estimator. 2. SCALE PARAMETER ESTIMATOR OF THE EXPONENTIAL DISTRIBUTION We will consider the random variable X with the following cumulative distribution function: x F(x ; ) 1 e ( x 0 , 0) (1) where is an unknown scale parameter Using the relationships between MXXX( ) ; 22( ) ; ( ) , we obtain ()X a theoretical variation coefficient 1. This is a useful indicator, especially if MX() you have observational data which seems to be exponential and with variation coefficient of the selection closed to 1. If we consider x12, x ,... xn as a part of a population that follows an exponential distribution, then by using the maximum likelihood estimation method we obtain the following estimate n ˆ 1 xi (2) n i1 119 On the Scale Parameter of Exponential Distribution Since M ˆ , it follows that ˆ is an unbiased estimator for . Similarly, because 2 2 D2 ()ˆ and lim 0 , we obtain that ˆ is an absolutely correct estimator. n n n The efficiency of the estimator can be calculated using the following formula: 1 2 nlnfx ( , ) f(,) x dx (3) e ()ˆ 0 n D()ˆ 1 x where f(;) x e ( x 0 , 0) is the probability density function of the 2 n exponential distribution. After the calculation we obtain nlnfx ( ; ) f(;) x dx 2 0 ˆ and thus en ( ) 1 which implies that is also an efficient estimator. Taking into account the reproducibility property of the exponential distribution we can calculate the exact distribution of the random variable ˆ . A sum XXX12 ... n , consisting of n randomly selected variables, all exponentially distribuited, has a distribution density function equal to n1 x 1 x fn (). x e (4) ()n This leads us to the cumulative distribution function nx u 1 n1 HxP()()(...ˆ xPX Xnx ) u eduFnx () n1 n n (5) ()n 0 from which we obtain that the distribution density function of the estimator ˆ is equal to n1 nx n nx hnn()() x f nx n e (6) ()n From the above we obtain that M()()ˆ xh x dx and n 0 2 D22()()()ˆ x h x which is consistent with previously known results. n 0 n By using the central limit theorem we can state that the random variable X follows a 2 2 normal distribution with the parameters and , which means that XN ; . n n 120 Review of the Air Force Academy No.2 (34)/2017 Based on this we can calculate the asymptotical cumulative distribution function for the estimator ˆ , ˆ Gn ( x ) lim P ( x ) lim P ( X x ) . (7) nn This shows that the estimator has the same asymptotical distribution as and a distribution density function equal to: 2 ()xn 1 n 2 g() x e (8) n 2 Using the asymptotical distribution of the estimator ˆ we can determine the confidence intervals for the scale parameter . To do this we need to consider the reduced ()ˆ n normal random variable Z and the significance level . By definition we will get PZZ( ) 1 1 (9) 2 where Z is the 1 quartile of the standard normal distribution. Equation (9) can be 1 2 2 rewritten such as ()ˆ n PZZ( ) 1 . (10) 11 22 Z 1 We may assume that 102 because we deal with asymptotical distribution and thus n n is sufficiently high. After carrying out the calculation we obtain for the parameter the following confidence interval which depends on the maximum likelihood estimation ˆ ˆˆ . ZZ 11 (11) 1122 nn 3. SET LENGTH CONFIDENCE INTERVALS FOR THE SCALE PARAMETER Using a similar method to that which Stein proposed for the double selection method we find a set length confidence interval for the parameter . Let p be a system of independent selections, each with a volume n, taken from a population which has the same exponential distribution as the random variableX X 121 On the Scale Parameter of Exponential Distribution xx11,..., 1n ... (12) xxp1,..., pn These selections allow us to obtain independent and identical assigned estimators 2 ˆ i N; i 1,..., p . (13) n 1 p 1 p Let ˆˆ and s22()ˆˆ. Because ˆ is the selection mean for the i ˆ i p i1 p 1 i1 2 random variables ˆ ; i 1,...,p we know that ˆ N(;) . i np Next let us consider a second system with m independent selections, each of volume n, taken from a population which has the same exponential distribution as the random variable X xxp1;1,..., p 1; n ... (14) xx,..., p m;1 p m ; n These new selections allow us to calculate m independent and identically assigned estimators 2 ˆ N( ; ) j1,..., m . (15) pj n Let ˆ be the calculated average of both selection systems: pm ˆˆ1 i (16) pm i1 2 ˆ ˆˆ It is know that N ; . Note that for m=0 we obtain . Due to the n() m p 2 fact that ˆ N( ; ) and s2 is the variance of random variable ˆ considered in n equation 2, we deduce the following equation 2 2 s ˆ p1 (17) 2 (p 1) n 122 Review of the Air Force Academy No.2 (34)/2017 Also the random variable matches with the reduced random variable through the following equation: n() m p (ˆ )N (0,1) (18) The random variables from equations (17) and (18) are independent and thus we can consider the following ratio n() m p ()ˆ ()mp ()ˆ t (19) sn s p1 ˆ ˆ where t p1 is a Student random variable with p-1 degrees of freedom. For the significance level we have the following equation: mp ˆ Pt( ) 1 (20) s ;1p ˆ 2 from which we deduce the following confidence interval for s t s t ˆˆ;pp 1 ; 1 ˆˆ22 (21) m p m p The length of this interval can be easily calculated and is equal to 2st ˆ ˆ ;1p l 2 . (22) mp The length of this interval must not be greater than the considered length . We will start with m=0 and compare the length l with . If l then the interval ˆˆ11 ; is a confidence interval for of size , that has a confidence 22 coefficient greater or at least equal to 1 , build only with the help of the first selection system. If l , we need to carry out the second selection system, where m is equal with the 2st ˆ ;1p 11 2 ˆˆ smallest integer for which . In this case ; will be the mp 22 confidence interval for of size , build on both selection systems and having a confidence coefficient greater then 1-. 123 On the Scale Parameter of Exponential Distribution REFERENCES [1] Gheorghe, F.C. , Obreja , G. Teste statistice pentru selecţii de volum mic ( n 50) Ed. Univ. Piteşti 1999. [2] Margolin,H.B., Maurer,W. Tests of the Kolmogorov Smirnov type for exponential data with unknow scale and related problems. Biometrica (1976), 63,1 pg.149-160 124 .
Load more

Recommended publications
  • A Family of Skew-Normal Distributions for Modeling Proportions and Rates with Zeros/Ones Excess
    S S symmetry Article A Family of Skew-Normal Distributions for Modeling Proportions and Rates with Zeros/Ones Excess Guillermo Martínez-Flórez 1, Víctor Leiva 2,* , Emilio Gómez-Déniz 3 and Carolina Marchant 4 1 Departamento de Matemáticas y Estadística, Facultad de Ciencias Básicas, Universidad de Córdoba, Montería 14014, Colombia; [email protected] 2 Escuela de Ingeniería Industrial, Pontificia Universidad Católica de Valparaíso, 2362807 Valparaíso, Chile 3 Facultad de Economía, Empresa y Turismo, Universidad de Las Palmas de Gran Canaria and TIDES Institute, 35001 Canarias, Spain; [email protected] 4 Facultad de Ciencias Básicas, Universidad Católica del Maule, 3466706 Talca, Chile; [email protected] * Correspondence: [email protected] or [email protected] Received: 30 June 2020; Accepted: 19 August 2020; Published: 1 September 2020 Abstract: In this paper, we consider skew-normal distributions for constructing new a distribution which allows us to model proportions and rates with zero/one inflation as an alternative to the inflated beta distributions. The new distribution is a mixture between a Bernoulli distribution for explaining the zero/one excess and a censored skew-normal distribution for the continuous variable. The maximum likelihood method is used for parameter estimation. Observed and expected Fisher information matrices are derived to conduct likelihood-based inference in this new type skew-normal distribution. Given the flexibility of the new distributions, we are able to show, in real data scenarios, the good performance of our proposal. Keywords: beta distribution; centered skew-normal distribution; maximum-likelihood methods; Monte Carlo simulations; proportions; R software; rates; zero/one inflated data 1.
    [Show full text]
  • Generalized Inferences for the Common Scale Parameter of Several Pareto Populations∗
    InternationalInternational JournalJournal of of Statistics Statistics and and Management Management System System Vol.Vol. 4 No. 12 (January-June,(July-December, 2019) 2019) International Journal of Statistics and Management System, 2010, Vol. 5, No. 1–2, pp. 118–126. c 2010 Serials Publications Generalized inferences for the common scale parameter of several Pareto populations∗ Sumith Gunasekera†and Malwane M. A. Ananda‡ Received: 13th August 2018 Revised: 24th December 2018 Accepted: 10th March 2019 Abstract A problem of interest in this article is statistical inferences concerning the com- mon scale parameter of several Pareto distributions. Using the generalized p-value approach, exact confidence intervals and the exact tests for testing the common scale parameter are given. Examples are given in order to illustrate our procedures. A limited simulation study is given to demonstrate the performance of the proposed procedures. 1 Introduction In this paper, we consider k (k ≥ 2) independent Pareto distributions with an unknown common scale parameter θ (sometimes referred to as the “location pa- rameter” and also as the “truncation parameter”) and unknown possibly unequal shape parameters αi’s (i = 1, 2, ..., k). Using the generalized variable approach (Tsui and Weer- ahandi [8]), we construct an exact test for testing θ. Furthermore, using the generalized confidence interval (Weerahandi [11]), we construct an exact confidence interval for θ as well. A limited simulation study was carried out to compare the performance of these gen- eralized procedures with the approximate procedures based on the large sample method as well as with the other test procedures based on the combination of p-values.
    [Show full text]
  • Procedures for Estimation of Weibull Parameters James W
    United States Department of Agriculture Procedures for Estimation of Weibull Parameters James W. Evans David E. Kretschmann David W. Green Forest Forest Products General Technical Report February Service Laboratory FPL–GTR–264 2019 Abstract Contents The primary purpose of this publication is to provide an 1 Introduction .................................................................. 1 overview of the information in the statistical literature on 2 Background .................................................................. 1 the different methods developed for fitting a Weibull distribution to an uncensored set of data and on any 3 Estimation Procedures .................................................. 1 comparisons between methods that have been studied in the 4 Historical Comparisons of Individual statistics literature. This should help the person using a Estimator Types ........................................................ 8 Weibull distribution to represent a data set realize some advantages and disadvantages of some basic methods. It 5 Other Methods of Estimating Parameters of should also help both in evaluating other studies using the Weibull Distribution .......................................... 11 different methods of Weibull parameter estimation and in 6 Discussion .................................................................. 12 discussions on American Society for Testing and Materials Standard D5457, which appears to allow a choice for the 7 Conclusion ................................................................
    [Show full text]
  • Volatility Modeling Using the Student's T Distribution
    Volatility Modeling Using the Student’s t Distribution Maria S. Heracleous Dissertation submitted to the Faculty of the Virginia Polytechnic Institute and State University in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Economics Aris Spanos, Chair Richard Ashley Raman Kumar Anya McGuirk Dennis Yang August 29, 2003 Blacksburg, Virginia Keywords: Student’s t Distribution, Multivariate GARCH, VAR, Exchange Rates Copyright 2003, Maria S. Heracleous Volatility Modeling Using the Student’s t Distribution Maria S. Heracleous (ABSTRACT) Over the last twenty years or so the Dynamic Volatility literature has produced a wealth of uni- variateandmultivariateGARCHtypemodels.Whiletheunivariatemodelshavebeenrelatively successful in empirical studies, they suffer from a number of weaknesses, such as unverifiable param- eter restrictions, existence of moment conditions and the retention of Normality. These problems are naturally more acute in the multivariate GARCH type models, which in addition have the problem of overparameterization. This dissertation uses the Student’s t distribution and follows the Probabilistic Reduction (PR) methodology to modify and extend the univariate and multivariate volatility models viewed as alternative to the GARCH models. Its most important advantage is that it gives rise to internally consistent statistical models that do not require ad hoc parameter restrictions unlike the GARCH formulations. Chapters 1 and 2 provide an overview of my dissertation and recent developments in the volatil- ity literature. In Chapter 3 we provide an empirical illustration of the PR approach for modeling univariate volatility. Estimation results suggest that the Student’s t AR model is a parsimonious and statistically adequate representation of exchange rate returns and Dow Jones returns data.
    [Show full text]
  • Determination of the Weibull Parameters from the Mean Value and the Coefficient of Variation of the Measured Strength for Brittle Ceramics
    Journal of Advanced Ceramics 2017, 6(2): 149–156 ISSN 2226-4108 DOI: 10.1007/s40145-017-0227-3 CN 10-1154/TQ Research Article Determination of the Weibull parameters from the mean value and the coefficient of variation of the measured strength for brittle ceramics Bin DENGa, Danyu JIANGb,* aDepartment of the Prosthodontics, Chinese PLA General Hospital, Beijing 100853, China bAnalysis and Testing Center for Inorganic Materials, State Key Laboratory of High Performance Ceramics and Superfine Microstructure, Shanghai Institute of Ceramics, Shanghai 200050, China Received: January 11, 2017; Accepted: April 7, 2017 © The Author(s) 2017. This article is published with open access at Springerlink.com Abstract: Accurate estimation of Weibull parameters is an important issue for the characterization of the strength variability of brittle ceramics with Weibull statistics. In this paper, a simple method was established for the determination of the Weibull parameters, Weibull modulus m and scale parameter 0 , based on Monte Carlo simulation. It was shown that an unbiased estimation for Weibull modulus can be yielded directly from the coefficient of variation of the considered strength sample. Furthermore, using the yielded Weibull estimator and the mean value of the strength in the considered sample, the scale parameter 0 can also be estimated accurately. Keywords: Weibull distribution; Weibull parameters; strength variability; unbiased estimation; coefficient of variation 1 Introduction likelihood (ML) method. Many studies [9–18] based on Monte Carlo simulation have shown that each of these Weibull statistics [1] has been widely employed to methods has its benefits and drawbacks. model the variability in the fracture strength of brittle According to the theory of statistics, the expected ceramics [2–8].
    [Show full text]
  • Power Birnbaum-Saunders Student T Distribution
    Revista Integración Escuela de Matemáticas Universidad Industrial de Santander H Vol. 35, No. 1, 2017, pág. 51–70 Power Birnbaum-Saunders Student t distribution Germán Moreno-Arenas a ∗, Guillermo Martínez-Flórez b Heleno Bolfarine c a Universidad Industrial de Santander, Escuela de Matemáticas, Bucaramanga, Colombia. b Universidad de Córdoba, Departamento de Matemáticas y Estadística, Montería, Colombia. c Universidade de São Paulo, Departamento de Estatística, São Paulo, Brazil. Abstract. The fatigue life distribution proposed by Birnbaum and Saunders has been used quite effectively to model times to failure for materials sub- ject to fatigue. In this article, we introduce an extension of the classical Birnbaum-Saunders distribution substituting the normal distribution by the power Student t distribution. The new distribution is more flexible than the classical Birnbaum-Saunders distribution in terms of asymmetry and kurto- sis. We discuss maximum likelihood estimation of the model parameters and associated regression model. Two real data set are analysed and the results reveal that the proposed model better some other models proposed in the literature. Keywords: Birnbaum-Saunders distribution, alpha-power distribution, power Student t distribution. MSC2010: 62-07, 62F10, 62J02, 62N86. Distribución Birnbaum-Saunders Potencia t de Student Resumen. La distribución de probabilidad propuesta por Birnbaum y Saun- ders se ha usado con bastante eficacia para modelar tiempos de falla de ma- teriales sujetos a la fátiga. En este artículo definimos una extensión de la distribución Birnbaum-Saunders clásica sustituyendo la distribución normal por la distribución potencia t de Student. La nueva distribución es más flexi- ble que la distribución Birnbaum-Saunders clásica en términos de asimetría y curtosis.
    [Show full text]
  • Bayesian Modelling of Skewness and Kurtosis with Two-Piece Scale and Shape Transformations
    Bayesian modelling of skewness and kurtosis with two-piece scale and shape transformations F.J. Rubio and M.F.J. Steel∗ Department of Statistics, University of Warwick, Coventry, CV4 7AL, U.K. Abstract We introduce the family of univariate double two–piece distributions, obtained by using a density–based transformation of unimodal symmetric continuous distributions with a shape parameter. The resulting distributions contain five interpretable parameters that control the mode, as well as the scale and shape in each direction. Four-parameter subfamilies of this class of distributions that capture different types of asymmetry are presented. We propose interpretable scale and location-invariant benchmark priors and derive conditions for the pro- priety of the corresponding posterior distribution. The prior structures used allow for mean- ingful comparisons through Bayes factors within flexible families of distributions. These distributions are applied to data from finance, internet traffic data, and medicine, comparing them with appropriate competitors. 1 Introduction In the theory of statistical distributions, skewness and kurtosis are features of interest since they provide information about the shape of a distribution. Definitions and quantitative measures of these features have been widely discussed in the statistical literature (see e.g. van Zwet, 1964; Groeneveld and Meeden, 1984; Critchley and Jones, 2008). Distributions containing parameters that control skewness and/or kurtosis are attractive since they can lead to robust models. This sort of flexible distributions are typically obtained by adding ∗keywords: model comparison; posterior existence; prior elicitation; scale mixtures of normals; unimodal con- tinuous distributions 1 parameters to a known symmetric distribution through a parametric transformation.
    [Show full text]
  • 3.4 Exponential Families
    3.4 Exponential Families A family of pdfs or pmfs is called an exponential family if it can be expressed as ¡ Xk ¢ f(x|θ) = h(x)c(θ) exp wi(θ)ti(x) . (1) i=1 Here h(x) ≥ 0 and t1(x), . , tk(x) are real-valued functions of the observation x (they cannot depend on θ), and c(θ) ≥ 0 and w1(θ), . , wk(θ) are real-valued functions of the possibly vector-valued parameter θ (they cannot depend on x). Many common families introduced in the previous section are exponential families. They include the continuous families—normal, gamma, and beta, and the discrete families—binomial, Poisson, and negative binomial. Example 3.4.1 (Binomial exponential family) Let n be a positive integer and consider the binomial(n, p) family with 0 < p < 1. Then the pmf for this family, for x = 0, 1, . , n and 0 < p < 1, is µ ¶ n f(x|p) = px(1 − p)n−x x µ ¶ n ¡ p ¢ = (1 − p)n x x 1 − p µ ¶ n ¡ ¡ p ¢ ¢ = (1 − p)n exp log x . x 1 − p Define 8 >¡ ¢ < n x = 0, 1, . , n h(x) = x > :0 otherwise, p c(p) = (1 − p)n, 0 < p < 1, w (p) = log( ), 0 < p < 1, 1 1 − p and t1(x) = x. Then we have f(x|p) = h(x)c(p) exp{w1(p)t1(x)}. 1 Example 3.4.4 (Normal exponential family) Let f(x|µ, σ2) be the N(µ, σ2) family of pdfs, where θ = (µ, σ2), −∞ < µ < ∞, σ > 0.
    [Show full text]
  • A Primer on the Exponential Family of Distributions
    A Primer on the Exponential Family of Distributions David R. Clark, FCAS, MAAA, and Charles A. Thayer 117 A PRIMER ON THE EXPONENTIAL FAMILY OF DISTRIBUTIONS David R. Clark and Charles A. Thayer 2004 Call Paper Program on Generalized Linear Models Abstract Generahzed Linear Model (GLM) theory represents a significant advance beyond linear regression theor,], specifically in expanding the choice of probability distributions from the Normal to the Natural Exponential Famdy. This Primer is intended for GLM users seeking a hand)' reference on the model's d]smbutional assumptions. The Exponential Faintly of D,smbutions is introduced, with an emphasis on variance structures that may be suitable for aggregate loss models m property casualty insurance. 118 A PRIMER ON THE EXPONENTIAL FAMILY OF DISTRIBUTIONS INTRODUCTION Generalized Linear Model (GLM) theory is a signtficant advance beyond linear regression theory. A major part of this advance comes from allowmg a broader famdy of distributions to be used for the error term, rather than just the Normal (Gausstan) distributton as required m hnear regression. More specifically, GLM allows the user to select a distribution from the Exponentzal Family, which gives much greater flexibility in specifying the vanance structure of the variable being forecast (the "'response variable"). For insurance apphcations, this is a big step towards more realistic modeling of loss distributions, while preserving the advantages of regresston theory such as the ability to calculate standard errors for estimated parameters. The Exponentml family also includes several d~screte distributions that are attractive candtdates for modehng clatm counts and other events, but such models will not be considered here The purpose of this Primer is to give the practicmg actuary a basic introduction to the Exponential Family of distributions, so that GLM models can be designed to best approximate the behavior of the insurance phenomenon.
    [Show full text]
  • Optimal Segmentation Scale Parameter, Feature Subset and Classification Algorithm for Geographic Object-Based Crop Recognition Using Multisource Satellite Imagery
    remote sensing Article Optimal Segmentation Scale Parameter, Feature Subset and Classification Algorithm for Geographic Object-Based Crop Recognition Using Multisource Satellite Imagery Lingbo Yang 1,2 , Lamin R. Mansaray 1,2,3 , Jingfeng Huang 1,2,* and Limin Wang 4,5 1 Institute of Applied Remote Sensing and Information Technology, Zhejiang University, Hangzhou 310058, China; [email protected] (L.Y.); [email protected] (L.R.M.) 2 Key Laboratory of Agricultural Remote Sensing and Information Systems, Zhejiang University, Hangzhou 310058, China 3 Laboratory of Agro-meteorology and Geo-informatics, Magbosi Land, Water and Environment Research Centre (MLWERC), Sierra Leone Agricultural Research Institute (SLARI), Tower Hill, Freetown PMB 1313, Sierra Leone 4 Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing 100081, China; [email protected] 5 Key Laboratory of Agri-informatics, Ministry of Agriculture, Beijing 100081, China * Correspondence: [email protected]; Tel.: +86-571-8898-2830 Received: 4 February 2019; Accepted: 27 February 2019; Published: 3 March 2019 Abstract: Geographic object-based image analysis (GEOBIA) has been widely used in the remote sensing of agricultural crops. However, issues related to image segmentation, data redundancy and performance of different classification algorithms with GEOBIA have not been properly addressed in previous studies, thereby compromising the accuracy of subsequent thematic products. It is in this regard that the current study investigates the optimal scale parameter (SP) in multi-resolution segmentation, feature subset, and classification algorithm for use in GEOBIA based on multisource satellite imagery. For this purpose, a novel supervised optimal SP selection method was proposed based on information gain ratio, and was then compared with a preexisting unsupervised optimal SP selection method.
    [Show full text]
  • POLS 7050 the Exponential Family
    POLS 7050 Spring 2009 April 23, 2009 Pulling it All Together: Generalized Linear Models Throughout the course, we've considered the various models we've discussed discretely. In fact, however, most of the models we've discussed can actually be considered under a unified framework, as examples of a class of models known as generalized linear models (usually abbreviated GLMs). To understand GLMs, we first have to get a handle on the exponential family of probability distributions. Once we have that, we can show how a bunch of the models we've been learning about are special cases of the more general class of GLMs. The Exponential Family We'll start with a random variable Z, which as realizations called z. We are typically interested in the conditional density (PDF) of Z, where the conditioning is on some parameter or parameters : f(zjθ) = Pr(Z = zj ) A PDF is a member of the exponential family of probability functions if it can be written in the form:1 f(zjθ) = r(z)s( ) exp[q(z)h( )] (1) where r(·) and q(·) are functions of z that do not depend on , and (conversely) s(·) and h(·) are functions of that do not depend on z. For reasons that will be apparent in a minute, it is also necessary that r(z) > 0 and s( ) > 0. With a little bit of algebra, we can rewrite (1) as: f(zj ) = exp[ln r(z) + ln s( ) + q(z)h( )]: (2) Gill dubs the first two components the \additive" part, and the second the \interactive" bit.
    [Show full text]
  • Statistics of the Spectral Kurtosis Estimator
    Statistics of the Spectral Kurtosis Estimator Gelu M. Nita1 and Dale E. Gary1 ABSTRACT Spectral Kurtosis (SK; defined by Nita et al. 2007) is a statistical approach for detecting and removing radio frequency interference (RFI) in radio astronomy data. In this paper, the statistical properties of the SK estimator are investigated and all moments of its probability density function are analytically determined. These moments provide a means to determine the tail probabilities of the esti- mator that are essential to defining the thresholds for RFI discrimination. It is shown that, for a number of accumulated spectra M ≥ 24, the first SK standard moments satisfy the conditions required by a Pearson Type IV (Pearson 1985) probability distribution function (PDF), which is shown to accurately reproduce the observed distributions. The cumulative function (CF) of the Pearson Type IV, in both analytical and numerical form, is then found suitable for accurate estimation of the tail probabilities of the SK estimator. This same framework is also shown to be applicable to the related Time Domain Kurtosis (TDK) es- timator (Ruf, Gross, & Misra 2006), whose PDF corresponds to Pearson Type IV when the number of time-domain samples is M ≥ 46. The PDF and CF are determined for this case also. Subject headings: SK- Spectral Kurtosis, TDK- Time Domain Kurtosis, RFI- Radio Frequency Interference 1. Introduction Given the expansion of radio astronomy instrumentation to ever broader bandwidths, and the simultaneous increase in usage of the radio spectrum for wireless communication, radio frequency interference (RFI) has become a limiting factor in the design of a new generation of radio telescopes.
    [Show full text]