Wikipedia / Parametric statistics
Wikipedia / Parametric statistics
Parametric statistics is a branch of statistics that is concerned with the analysis of and inference from data assuming that the underlying distribution, from which the observed data was drawn, can be described by a finite set of (unknown) parameters. In contrast, nonparametric statistics does not assume explicit (finite-parametric) mathematical forms for distributions when modeling data. However, it may make some assumptions about that distribution, such as continuity or symmetry, or even an explicit mathematical shape but have a model for a distributional parameter that is not itself finite-parametric.
Most well-known statistical methods are parametric. Regarding nonparametric (and semiparametric) models, Sir David Cox has said, "These typically involve fewer assumptions of structure and distributional form but usually contain strong assumptions about independencies".
Overview
The main goal of statistical inference is to provide methods to systematically analyse data and infer properties of the probability distribution from which the data was drawn. The fundamental assumption underlying parametric statistics is that the data distribution belongs to a more general family of distributions that can be parameterized by a finite number of parameters, which are unknown. Such a family of distributions is called a parametric model. Typical questions in this setting are:
Parameter estimation: Which choice of parameters best explains the observed data or leads to best predictions?
Interval estimation with high probability?
Hypothesis testing: Is it plausible that the "best" parameter has a certain value or lies within a certain range?
Model (mis-)specification: What is an appropriate family of distributions that can explain the observed data?
Parameter estimation
, called the likelihood (or a function thereof) based on the observed data.
Frequentist approaches
from the observed data. The most common methods for parameter estimation are the following.
Maximum Likelihood estimation (MLE): The model parameters are chosen such that the probability (or probability density) of the given observation is maximal. Under certain assumptions on the likelihood, the MLE is (strongly) consistent, asymptotically efficient and asymptotically normal.
Method of moments of the moments of the distribution are continuous, then the moment estimator is strongly consistent (due to the strong law of large numbers are even differentiable, then the moment estimator is also asymptotically normal (due to the central limit theorem and the delta method).
Least square estimation (LSE): This method applies to a regression is normally distributed.
Bayesian approaches
In a Bayesian reflecting the statistician's prior belief. Given the observed data, the parameter distribution is updated via Bayes' rule. Bayesian estimators therefore rather give a best estimate given the statisticians beliefs.
Posterior mean estimation: The estimator takes the average of the posterior distribution.
Posterior median estimation: The estimator takes the median of the posterior distribution.
Maximum à-posteriori estimation (MAP): The estimator takes the parameter that maximizes the posterior distribution. If the prior is the uniform distribution, then MAP estimation is equivalent to maximum likelihood estimation.
Optimality of estimators
Uniformly minimum-variance unbiased estimators
Uniformly minimum-variance unbiased estimators (UMVUE), sometimes called best unbiased estimators as well, are estimators that have minimum variance among all unbiased estimators. Due to the bias-variance decomposition, they are optimal in the sense that they minimise the mean squared error among all unbiased estimators.
with finite variance exists, then it must be unique due to the Rao–Blackwell theorem. If there exists a complete sufficient statistic in the chosen statistical model, which is often easy to determine using the Neyman–Fisher theorem, then it holds the Lehmann–Scheffé theorem, then it is also the unique UMVUE.
Best linear unbiased estimator
. In a sense, BLUEs are UMVU estimators with the additional restriction of being linear.
In the regression setting, the Gauss–Markov theorem states that the least squares estimator is the unique BLUE. This theorem also shows (considering the special case of constant regression) that the empirical mean is a BLUE for estimating the mean.
Information inequality and exponential families
In regular is bounded from below by the universal bound is the Fisher-information matrix of the statistical model. The right hand-side is the so-called Cramér–Rao bound.
Estimators for which the variance equals the Cramér–Rao bound are called efficient and are almost UMVU, since they minimize the variance over all regular is efficient if and only if (i) the statistical model is an exponential family is the natural sufficient statistic.
Parametric models
The choice of the model in regression/classification problems, lies at the core of parametric procedures. Here is a list of common models used in practice.
Density estimation
Exponential families (e.g. normal distribution, exponential distribution, log-normal distribution, Gamma distribution, Chi-squared distribution, Erlang distribution, Beta distribution, Gumbel distribution, Pareto distribution, (Negative-)Binomial distribution, Poisson distribution, geometric distribution)
Regression
linear model (special cases thereof are ANOVA and ANCOVA)
generalized linear model (GLM)
Classification
Example
The normal family of distributions all have the same general shape and are parameterized by mean and standard deviation. That means that if the mean and standard deviation are known and if the distribution is normal, the probability of any future observation lying in a given range is known.
Suppose that we have a sample of 99 test scores with a mean of 100 and a standard deviation of 1. If we assume all 99 test scores are random observations from a normal distribution, then we predict there is a 1% chance that the 100th test score will be higher than 102.33 (that is, the mean plus 2.33 standard deviations), assuming that the 100th test score comes from the same distribution as the others. Parametric statistical methods are used to compute the 2.33 value above, given 99 independent observations from the same normal distribution.
A non-parametric estimate of the same thing is the maximum of the first 99 scores. We don't need to assume anything about the distribution of test scores to reason that before we gave the test it was equally likely that the highest score would be any of the first 100. Thus there is a 1% chance that the 100th score is higher than any of the 99 that preceded it.
History
Parametric statistics was mentioned by R. A. Fisher in his work Statistical Methods for Research Workers in 1925, which created the foundation for modern statistics.
See also
Notes
References