Neyman construction

Neyman construction, named after Jerzy Neyman, is a frequentist method to construct an interval at a confidence level ${\displaystyle C,\,}$ such that if we repeat the experiment many times the interval will contain the true value of some parameter a fraction ${\displaystyle C\,}$ of the time.

Assume ${\displaystyle X_{1},X_{2},...X_{n}}$ are random variables with joint pdf ${\displaystyle f(x_{1},x_{2},...x_{n}|\theta _{1},\theta _{2},...,\theta _{k})}$, which depends on k unknown parameters. For convenience, let ${\displaystyle \Theta }$ be the sample space defined by the n random variables and subsequently define a sample point in the sample space as ${\displaystyle X=(X_{1},X_{2},...X_{n})}$
Neyman originally proposed defining two functions ${\displaystyle L(x)}$ and ${\displaystyle U(x)}$ such that for any sample point,${\displaystyle X}$,

• ${\displaystyle L(X)\leq U(X)}$ ${\displaystyle \forall X\in \Theta }$
• L and U are single valued and defined.

Given an observation, ${\displaystyle X^{'}}$, the probability that ${\displaystyle \theta _{1}}$ lies between ${\displaystyle L(X^{'})}$ and ${\displaystyle U(X^{'})}$ is defined as ${\displaystyle P(L(X^{'})\leq \theta _{1}\leq U(X^{'})|X^{'})}$ with probability of ${\displaystyle 0}$ or ${\displaystyle 1}$. These calculated probabilities fail to draw meaningful inference about ${\displaystyle \theta _{1}}$ since the probability is simply zero or unity. Furthermore, under the frequentist construct the model parameters are unknown constants and not permitted to be random variables.^[1] For example if ${\displaystyle \theta _{1}=5}$, then ${\displaystyle P(2\leq 5\leq 10)=1}$. Likewise, if ${\displaystyle \theta _{1}=11}$, then ${\displaystyle P(2\leq 11\leq 10)=0}$

As Neyman describes in his 1937 paper, suppose that we consider all points in the sample space, that is, ${\displaystyle \forall X\in \Theta }$, which are a system of random variables defined by the joint pdf described above. Since ${\displaystyle L}$ and ${\displaystyle U}$ are functions of ${\displaystyle X}$ they too are random variables and one can examine the meaning of the following probability statement:

Under the frequentist construct the model parameters are unknown constants and not permitted to be random variables. Considering all the sample points in the sample space as random variables defined by the joint pdf above, that is all ${\displaystyle X\in \Theta }$ it can be shown that ${\displaystyle L}$ and ${\displaystyle U}$ are functions of random variables and hence random variables. Therefore one can look at the probability of ${\displaystyle L(X)}$ and ${\displaystyle U(X)}$ for some ${\displaystyle X\in \Theta }$. If ${\displaystyle \theta _{1}^{'}}$ is the true value of ${\displaystyle \theta _{1}}$, we can define ${\displaystyle L}$ and ${\displaystyle U}$ such that the probability ${\displaystyle L(X)\leq \theta _{1}^{'}}$ and ${\displaystyle \theta _{1}^{'}\leq U(X)}$ is equal to pre-specified confidence level${\displaystyle ,C}$.

That is, ${\displaystyle P(L(X)\leq \theta _{1}^{'}\leq U(X)|\theta _{1}^{'})=C}$ where ${\displaystyle 0\leq C\leq 1}$ and ${\displaystyle L(X)}$ and ${\displaystyle U(X)}$ are the upper and lower confidence limits for ${\displaystyle \theta _{1}}$^[1]

The coverage probability, ${\displaystyle C}$, for Neyman construction is the frequency of experiments in which the confidence interval contains the actual value of interest. Generally, the coverage probability is set to a ${\displaystyle 95\%}$ confidence. For Neyman construction, the coverage probability is set to some value ${\displaystyle C}$ where ${\displaystyle 0<C<1}$.

A Neyman construction can be carried out by performing multiple experiments that construct data sets corresponding to a given value of the parameter. The experiments are fitted with conventional methods, and the space of fitted parameter values constitutes the band which the confidence interval can be selected from.

Suppose ${\displaystyle X\sim N(\theta ,\sigma ^{2})}$, where ${\displaystyle \theta }$ and ${\displaystyle \sigma ^{2}}$ are unknown constants where we wish to estimate ${\displaystyle \theta }$. We can define (2) single value functions, ${\displaystyle L}$ and ${\displaystyle U}$, defined by the process above such that given a pre-specified confidence level, ${\displaystyle C}$, and random sample ${\displaystyle X^{*}=(x_{1},x_{2},...x_{n})}$

${\displaystyle L(X^{*})={\bar {x}}-t{\frac {s}{\sqrt {n}}}}$

${\displaystyle U(X^{*})={\bar {x}}+t{\frac {s}{\sqrt {n}}}}$

where ${\displaystyle s/{\sqrt {n}}}$ is the standard error, and the sample mean and standard deviation are:

${\displaystyle {\bar {x}}={\frac {1}{n}}\sum _{i=1}^{n}x_{i}={\frac {1}{n}}(x_{1},x_{2},...x_{n})}$

${\displaystyle s={\sqrt {{\frac {1}{n-1}}\sum _{i=1}^{n}(x_{i}-{\bar {x}})^{2}}}}$

The factor ${\displaystyle t}$ follows a t distribution with (n-1) degrees of freedom, ${\displaystyle t}$~t_{${\displaystyle ({1-C}/2,n-1)}$} ^[2]

${\displaystyle X_{1},X_{2},...,X_{n}}$ are iid random variables, and let ${\displaystyle T=(X_{1},X_{2},...,X_{n})}$. Suppose ${\displaystyle T\sim N(\mu ,\sigma ^{2})}$. Now to construct a confidence interval with ${\displaystyle C}$ level of confidence. We know ${\displaystyle {\bar {x}}}$ is sufficient for ${\displaystyle \mu }$. So,

${\displaystyle p(-Z_{\frac {\alpha }{2}}\leq {\frac {{\bar {x}}-\mu }{\sigma ^{2}}}\leq Z_{\frac {\alpha }{2}})=C}$

${\displaystyle p(-Z_{\frac {\alpha }{2}}\sigma ^{2}\leq {\bar {x}}-\mu \leq Z_{\frac {\alpha }{2}}\sigma ^{2})=C}$

${\displaystyle p({\bar {x}}-Z_{\frac {\alpha }{2}}\sigma ^{2}\leq \mu \leq {\bar {x}}+Z_{\frac {\alpha }{2}}\sigma ^{2})=C}$

This produces a ${\displaystyle 100(C)\%}$ confidence interval for ${\displaystyle \mu }$ where,

${\displaystyle L(T)={\bar {x}}-Z_{\frac {\alpha }{2}}\sigma ^{2}}$

${\displaystyle U(T)={\bar {x}}+Z_{\frac {\alpha }{2}}\sigma ^{2}}$.

^[3]

• Probability interpretations