inference Archives • Statisticelle

Parametric statistics assume that the unknown CDF $F$ belongs to a family of CDFs characterized by a parameter (vector) $\theta$ . As the form of $F$ is assumed, the target of estimation is its parameters $\theta$ . Thus, all uncertainty about $F$ is comprised of uncertainty about its parameters. Parameters are estimated by $\hat{\theta}$ , and estimates are be substituted into the assumed distribution to conduct inference for the quantities of interest. If the assumed distribution $F$ is incorrect, inference may also be inaccurate, or trends in the data may be missed.

To demonstrate the parametric approach, consider $n = 100$ independent and identically distributed random variables $X_1, …, X_n$ generated from an exponential distribution with rate $\lambda = 2$ . Investigators wish to estimate the 75 $^{th}$ percentile and erroneously assume that their data is normally distributed. Thus, $F$ is assumed to be the Normal CDF but $\mu$ and $\sigma^2$ are unknown. The parameters $\mu$ and $\sigma$ are estimated in their typical way by $\bar{x}$ and $\sigma^2$ , respectively. Since the normal distribution belongs to the location-scale family, an estimate of the $p^{th}$ percentile is provided by,

$x_p = \bar{x} + s\Phi^{-1}(p)$

where $\Phi^{-1}$ is the standard normal quantile function, also known as the probit.

set.seed(12345)
library(tidyverse, quietly = T)

# Generate data from Exp(2)
x <- rexp(n = 100, rate = 2)

# True value of 75th percentile with rate = 2
true <- qexp(p = 0.75, rate = 2) 
true

## [1] 0.6931472

# Estimate mu and sigma
xbar <- mean(x)
s    <- sd(x)

# Estimate 75th percentile assuming mu = xbar and sigma = s
param_est <- xbar + s * qnorm(p = 0.75)
param_est

## [1] 0.8792925

The true value of the 75 $^{th}$ percentile of $\text{Exp}(2)$ is 0.69 while the parametric estimate is 0.88.

Nonparametric statistics make fewer distributions about the unknown distribution $F$ , requiring only mild assumptions such as continuity or the existence of specific moments. Instead of estimating parameters of $F$ , $F$ itself is the target of estimation. $F$ is commonly estimated by the empirical cumulative distribution function (ECDF) $\hat{F}$ ,

$\hat{F}(x) = \frac{1}{n} \sum_{i=1}^{n} \mathbb{I}(X_i \leq x).$

Any statistic that can be expressed as a function of the CDF, known as a statistical functional and denoted $\theta = T(F)$ , can be estimated by substituting $\hat{F}$ for $F$ . That is, plug-in estimators can be obtained as $\hat{\theta} = T(\hat{F})$ .

Continue reading Parametric vs. Nonparametric Approach to Estimations