Determination of the Weibull Parameters from the Mean Value and the Coefficient of Variation of the Measured Strength for Brittle Ceramics

Total Page:16

File Type:pdf, Size:1020Kb

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]. In general, a two-parameter form of the value, E( ) , and the variance, D( ) , of the two- Weibull function is adopted to give the cumulative parameter Weibull distribution can be expressed as [19]: failure probability, P, of a component or specimen at a 1 E()0 1 (2) given applied stress, , as [2]: m m 22 21 P 1exp (1) D()0 1 1 (3) mm 0 Thus, the coefficient of variation for this distribution, where m is called the Weibull modulus and 0 is the ()cvw, is scale parameter with the same dimension as . 1/2 Different methods have been used to estimate the 212 11 D() mm Weibull parameters, m and 0 , the most popular being ()cvw (4) the least-square (LS) method and the maximum- E() 1 1 m Equation (4) shows that ()c is dependent only on * Corresponding author. vw E-mail: [email protected] the Weibull modulus, m, and decreases as m increases www.springer.com/journal/40145 150 J Adv Ceram 2017, 6(2): 149–156 when m > 1. Based on this relation, Gong et al. [20,21] the strength value, i , with the prescribed mtr : proposed an empirical equation for determining the 1/mtr 1 Weibull modulus, mes , from the coefficient of variation, i ln (6) 1RND i ()cvf, of measured fracture strength: Thus, a sample containing N strength values, , m (5) 1 es ()cvf 2, …, N, is obtained. where and are both constants dependent on the (ii) The yielded strength values are ordered from the sample size, and were determined and formulated by lowest to the highest and the ith result in the set of N analyzing the data produced by Monte Carlo simulation strength data, , is assigned a cumulative fracture in the original works of Gong et al. [20,21]. i probability, Pi , which is calculated with [22–24]: It should be noted that Gong et al.’s method is only i 0.5 P (7) for the determination of Weibull modulus, and the other i N Weibull parameter, the scale parameter 0 , should be The strengths and the probabilities are then analyzed, determined still with the conventional methods, i.e., the by using a simple, least-square regression method, LS method or the ML method. Undoubtedly, this according to the alternative form of Eq. (1): drawback would impede the application of this simple 1 method. ln lnmmes lnes ln( 0 ) es (8) 1 Pi The aim of this paper is to contribute supplementary to give the estimated values of the Weibull parameters, information towards a more accurate determination of m and () . the Weibull modulus based on Eq. (5) and then establish es 0es an approach for the direct determination of the scale (iii) The mean values, f , and the standard deviation, parameter 0 . Considering its versatility and simplicity, Sf , of the N strength data are calculated as the LS method was employed in the present work to 1 N yield the referenced Weibull parameters for the analyses f i (9) N i1 and comparison. N 221 Sff()i (10) N i1 2 Monte Carlo simulation procedure Then coefficient of variation, ()cvf, is given by Sf ()cvf (11) The basic idea of the Monte Carlo simulation can be f described as: suppose we have a material whose To study the effect of sample size, N is increased fracture strength variation follows a two-parameter progressively from 10 to 200. Similarly, mtr is distribution of known parameters, mtr and () 0tr (for increased progressively from 5 to 25 to study the effect the sake of simplicity, 0 is set to be unity throughout of Weibull modulus. For each set of given sample size, the whole study); thus we can choose N random N, and the prescribed Weibull modulus, mtr , the specimens of this material to yield a sample, “measure” above-mentioned procedure is repeated for 100 000 the fracture strength of each specimen by assigning a times in order to ensure statistical convergence of the random number in the internal of 0–1 as its results. corresponding fracture probability, P, and then estimate the Weibull parameters. In continuation of this idea, a computer program is 3 Results and discussion written to establish the relationships between the 3. 1 Refinement of Gong et al.’s method estimated Weibull parameters, mes and () 0es, the mean values, f , and the coefficient of variation, ()cvf, In the original works of Gong et al. [20,21], by of the “measured” fracture strengths. repeating the above-mentioned Monte Carlo simulation (i) Firstly, a sample of random number, RND (i = 1, i procedure 4000 times for a given N and mtr = 2, 3, 4, …, 2, …, N) in the interval of 0–1 is produced to calculate 25, 26, 100 000 (= 4000 25) samples of strength data www.springer.com/journal/40145 J Adv Ceram 2017, 6(2): 149–156 151 were produced. Then, a least-square regression analysis was applied with all the resultant 100 000 estimated parameters, mes and (),cvf to yield the parameters and in Eq. (5). Such a treatment, however, seems to be questioned. Figure 1 shows, as an example, some data produced with N = 30 and mtr = 2.5, 10, 20. It can be seen from the insert, which shows the enlargement of the small mtr region, that the data points produced with mtr= 2.5 show large scatter, making it unreliable to fit these data to a straight line. The dashed line in Fig. 1 is obtained by least-square regression analysis of the data produced with mtr 20. Clearly, most of the data points corresponding to mtr 2.5 locate above the dashed line. If all the data shown in Fig. 1 were fitted according to Eq. (5), however, most of the data points Fig. 2 Best-fit value of as a function of the prescribed mtr . corresponding to mtr 2.5 locate below the best-fit line (solid line in Fig. 1). of statistics, fitting the data produced with different mtr To get a clearer understanding on the experimental to a single straight line would result in some phenomena shown in Fig. 1, 100 000 samples of uncertainties in the best-fit values of and . strength data were produced with given N and mtr . As mentioned above, in the original works of Gong et Then the resultant data were analyzed according to Eq. al. [20,21], the regression analyses for determining (5) to yield and . Figure 2 shows the best-fit value of and were conducted within mtr ranging from 2 to 26 the parameter as a function of the prescribed mtr for and only 4000 data produced with each prescribed mtr some selected N. As can be seen, strongly depends on were used. Based on the above discussion, such a the prescribed mtr value, especially when the sample simulation procedure seems to be insufficient to ensure size N is smaller. On the other hand, the best-fit value of the statistical convergence and warrant a further the parameter is found to be very small, less than improvement. 5 510 in each case. These results imply that, each set of Therefore, three key issues were considered here to the mes – ()cvf data points produced with a given mtr improve Gong et al.’s simulation procedure: would fall along a straight line starting from the origin (i) For each given set of N and mtr , 100 000 samples and with a slope dependent of mtr . From the viewpoint were analyzed to ensure the statistical convergence. (ii) Noting that the Weibull modulus has typical values ranging 5 to 20 for brittle ceramics, the smallest and the largest mtr were set to be 5 and 25, respectively. (iii) For a given N, the mean values of mes and ()cvf corresponding to each prescribed mtr , were employed for the regression analyses for determining the parameters and . The relationships between the mean values of mes and ()cvf for different sample size, N, are now shown in Fig.
Recommended publications
  • 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 .
    [Show full text]
  • 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]
  • 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]