In statistics, a linear probability model (LPM) is a special case of a binary regression model. Here the dependent variable for each observation takes values which are either 0 or 1. The probability of observing a 0 or 1 in any one case is treated as depending on one or more explanatory variables. For the "linear probability model", this relationship is a particularly simple one, and allows the model to be fitted by linear regression.
The model assumes that, for a binary outcome (Bernoulli trial),
Y
X
\Pr(Y=1|X=x)=x'\beta.
For this model,
E[Y|X]=0 ⋅ \Pr(Y=0|X)+1 ⋅ \Pr(Y=1|X)=\Pr(Y=1|X)=x'\beta,
and hence the vector of parameters β can be estimated using least squares. This method of fitting would be inefficient, and can be improved by adopting an iterative scheme based on weighted least squares,[1] in which the model from the previous iteration is used to supply estimates of the conditional variances,
\operatorname{Var}(Y|X=x)
A drawback of this model is that, unless restrictions are placed on
\beta
[0,1]
More formally, the LPM can arise from a latent-variable formulation (usually to be found in the econometrics literature[2]), as follows: assume the following regression model with a latent (unobservable) dependent variable:
y*=b0+x'b+\varepsilon, \varepsilon\midx\simU(-a,a).
The critical assumption here is that the error term of this regression is a symmetric around zero uniform random variable, and hence, of mean zero. The cumulative distribution function of
\varepsilon
F\varepsilon|x(\varepsilon\midx)=
\varepsilon+a | |
2a |
.
Define the indicator variable
y=1
y*>0
{\rmPr}(y=1\midx)={\rmPr}(y*>0\midx)={\rmPr}(b0+x'b+\varepsilon>0\midx)
={\rmPr}(\varepsilon>-b0-x'b\midx)=1-{\rmPr}(\varepsilon\leq-b0-x'b\midx)
=1-F\varepsilon|x(-b0-x'b\midx)=1-
-b0-x'b+a | |
2a |
=
b0+a | + | |
2a |
x'b | |
2a |
.
But this is the Linear Probability Model,
P(y=1\midx)=\beta0+x'\beta
with the mapping
\beta0=
b0+a | |
2a |
, \beta=
b | |
2a |
.
This method is a general device to obtain a conditional probability model of a binary variable: if we assume that the distribution of the error term is Logistic, we obtain the logit model, while if we assume that it is the Normal, we obtain the probit model and, if we assume that it is the logarithm of a Weibull distribution, the complementary log-log model.