Worked examples of OLS point estimates

[21/76%]
Show
LaTeX

The variables xx and yy satisfy the linear equation y=ax+by=ax+b, whose coefficients need to be determined. From a series of equally precise measurements, the values xix_{i}, yiy_{i} (i=1,2,…,n)(i=1,2,\ldots ,n) are obtained for the quantities xx and yy. Using the method of least squares, find the most reliable values of the coefficients aa and bb.

?

In the plane, nn points Mi(xi,yi)M_{i}(x_{i},y_{i}) (i=1,2,…,n)(i=1,2,\ldots ,n) are given. In what position must the line xcos⁡α+ysin⁡α−p=0x\cos \alpha +y\sin \alpha -p=0 be so that the sum of squared deviations of the given points from this line is a minimum?

?

From the two observations X1=2,X2=4X_{1} = 2, X_{2} = 4 with h1=1,h2=2h_{1} = 1, h_{2} = 2, construct a realization of the OLS estimate of the unknown parameters in the simple linear regression model. Compute the coefficient of determination.

?

Using the data of Example 10.6, construct the dependence X=θ1+θ2hX = \theta_{1}+\theta_{2} h, assuming that:

?
(1)

h=[1,…,12]⊤h = [1, \ldots , 12]^{\top };

(2)

h=[1954,…,1967]⊤h = [1954, \ldots , 1967]^{\top }.

In both cases, estimate the error of the obtained realization of the estimate and compute the coefficient of determination.

Estimate the unknown parameters α,β\alpha , \beta in the model of a person's consumption of some product v=α+βln⁡wv = \alpha +\beta \ln w, where vv is the expenditure on the consumer product, ww is the weekly income (Table 10.4).

Table 10.4

kk1234567
wkw_{k}0,81,21,52,22,33,63,1
vkv_{k}1,72,73,65,76,78,112,0

Estimate the elasticity coefficient σ=dvdw⋅wv\sigma = \frac{d v}{d w} \cdot \frac{w}{v}. Assuming that the observation errors have the distribution N(0;σ2)\mathcal{N}\left(0 ; \sigma^{2}\right), use Fisher's test to test the hypothesis H0:β=1H_{0}: \beta = 1 at the significance level 0,05.

?

Table 14.4 gives the indices of output YkY_{k}, capital inputs KkK_{k} and labor inputs LkL_{k} for n=10n = 10 years (USA, 1913-1922).

Table 14.4

kkYkY_{k}KkK_{k}LkL_{k}kkYkY_{k}KkK_{k}LkL_{k}
11842361546223366200
21692441497218387193
31892661548231407193
42252981829179417147
522733519610240431161

The relationship between these variables is assumed to have the form

ln⁡(YkLk)=θ1+θ2ln⁡(KkLk)+εk,k=1,…,10 \ln \left(\frac{Y_{k}}{L_{k}}\right) = \theta _{1}+\theta _{2} \ln \left(\frac{K_{k}}{L_{k}}\right)+\varepsilon _{k}, \quad k = 1, \ldots , 10

Find estimates of the parameters θ1\theta_{1} and θ2\theta_{2} by OLS and LAD. Find a robust estimate of the parameter σ2=Var⁡[εk]\sigma^{2} = \operatorname {Var}\left[\varepsilon_{k}\right]. Analyze the data for outliers and reject them.

?
?
(a)

In Example 4.7, determine α^,β^\hat{\alpha }, \hat{\beta }, and hence ξ^i\hat{\xi }_{i} by minimizing Σ(Xi−α−βti)2\Sigma \left(X_{i}-\alpha -\beta t_{i}\right)^{2}.

(b)

Verify the expressions (4.12) for α\alpha and β\beta, and the corresponding expressions for α^\hat{\alpha } and β^\hat{\beta }.

In quadrilateral ABCDA B C D, independent measurements of equal precision of the angles ABD,DBC,ABC,BCD,CDBA B D, D B C, A B C, B C D, C D B, BDA,CDA,DABB D A, C D A, D A B (in degrees) gave the results 50.7850.78, 30.25,78.29,99.57,50.42,40.59,88.87,89.8630.25,78.29,99.57,50.42,40.59,88.87,89.86 respectively. Assuming that measurement errors are normally distributed according to the law N(0,σ2)\mathscr {N}\left(0, \sigma^{2}\right), find the optimal estimators of the angles β1=ABD,β2=DBC,β3=CDB,β4=BDA\beta_{1}=A B D, \beta_{2}=D B C, \beta_{3}=C D B, \beta_{4}=B D A and of the unknown variance σ2\sigma^{2}.

?

Let

Xi=β1+iβ2+ε0+…+εi, X_{i}=\beta _{1}+i \beta _{2}+\varepsilon _{0}+\ldots +\varepsilon _{i},

i=0,1,…,ni=0,1, \ldots , n, where β1,β2\beta_{1}, \beta_{2} are unknown parameters, and ε0,…,εn\varepsilon_{0}, \ldots , \varepsilon_{n} are independent random variables distributed according to the law N(0,σ2)\mathscr {N}\left(0, \sigma^{2}\right). Reduce the problem to a linear model and find the least-squares estimators for β1\beta_{1} and β2\beta_{2}, as well as an unbiased estimator for σ2\sigma^{2}.

?

Using the linear regression method, construct an approximation of the function f(x)f(x) by a cubic polynomial based on the following data:

f(xi)f\left(x_{i}\right)3.93.95.05.05.75.76.56.57.17.17.67.67.87.88.18.18.48.4
xix_{i}4.04.05.25.26.16.17.07.07.97.98.68.68.98.99.59.59.99.9
?

Let

Y1=θ+ε1Y2=2θ−ϕ+ε2Y3=θ+2ϕ+ε3 \begin{aligned} & Y_{1} = \theta +\varepsilon _{1} \\ & Y_{2} = 2 \theta -\phi +\varepsilon _{2} \\ & Y_{3} = \theta +2 \phi +\varepsilon _{3} \end{aligned}

where E[εi]=0(i=1,2,3)\mathbb {E}\left[\varepsilon_{i}\right] = 0(i = 1,2,3). Find the least squares estimates of θ\theta and ϕ\phi.

?

The tension TT observed in a nonextensible string required to maintain a body of unknown weight ww in equilibrium on a smooth inclined plane of angle θ(0<θ<π/2)\theta (0 < \theta < \pi / 2) is a random variable with mean E[T]=wsin⁡θ\mathbb {E}\left[T\right] = w \sin \theta. If for θ=θi(i=1,2,…,n)\theta = \theta_{i}(i = 1,2, \ldots , n) the corresponding values of TT are TiT_{i} (i=1,2,…,ni = 1,2, \ldots , n), find the least squares estimate of ww.

?

Let Yi=β0+β1xi+εi(i=1,2,…,n)Y_{i} = \beta_{0}+\beta_{1} x_{i}+\varepsilon_{i}(i = 1,2, \ldots , n), where E[ε]=0\mathbb {E}\left[\varepsilon \right] = 0 and Var⁡[ε]=σ2In\operatorname {Var}\left[\varepsilon \right] = \sigma^{2} \mathbf{I}_{n}. Find the least squares estimates of β0\beta_{0} and β1\beta_{1}. Prove that they are uncorrelated if and only if xˉ=0\bar{x} = 0.

?

In order to estimate two parameters θ\theta and ϕ\phi it is possible to make observations of three types: (a) the first type have expectation θ\theta, (b) the second type have expectation θ+ϕ\theta +\phi, and (c) the third type have expectation θ−2ϕ\theta -2 \phi. All observations are subject to uncorrelated errors of mean zero and constant variance. If mm observations of type (a), mm observations of (b), and nn observations of type (c) are made, find the least squares estimates θ^\widehat{\theta } and ϕ^\widehat{\phi }. Prove that these estimates are uncorrelated if m=2nm = 2 n.

?

Let

Yi=β0+β1(xi1−xˉ1)+β2(xi2−xˉ2)+εi(i=1,2,…,n), Y_{i} = \beta _{0}+\beta _{1}\left(x_{i 1}-\bar{x}_{1}\right)+\beta _{2}\left(x_{i 2}-\bar{x}_{2}\right)+\varepsilon _{i} \quad (i = 1,2, \ldots , n),

where xˉj=∑i=1nxij/n,E[ε]=0\bar{x}_{j} = \sum_{i = 1}^{n} x_{i j} / n, \mathbb {E}\left[\varepsilon \right] = \mathbf{0}, and Var⁡[ε]=σ2In\operatorname {Var}\left[\varepsilon \right] = \sigma^{2} \mathbf{I}_{n}. If β^1\widehat{\beta }_{1} is the least squares estimate of β1\beta_{1}, show that

Var⁡[β^1]=σ2∑i(xi1−xˉ1)2(1−r2) \operatorname {Var}\left[\widehat{\beta }_{1}\right] = \frac{\sigma ^{2}}{\sum _{i}\left(x_{i 1}-\bar{x}_{1}\right)^{2}\left(1-r^{2}\right)}

where rr is the correlation coefficient of the nn pairs (xi1,xi2)\left(x_{i 1}, x_{i 2}\right).

?

Obtain the least squares estimates of the μi\mu_{i} using β^=(XTX)−1XTY\widehat{\mathbf{\beta }} = \left(\mathbf{X}^{T} \mathbf{X}\right)^{-1} \mathbf{X}^{T} \mathbf{Y} directly, where X\mathbf{X} is given by (8.2).

?

Let Yi=aiβ1+biβ2+εi(i=1,2,…,n)Y_{i} = a_{i} \beta_{1}+b_{i} \beta_{2}+\varepsilon_{i}(i = 1,2, \ldots , n), where the ai,bia_{i}, b_{i} are known and the εi\varepsilon_{i} are independently and identically distributed as N(0,σ2)N\left(0, \sigma^{2}\right). Find a necessary and sufficient condition for the least squares estimates of β1\beta_{1} and β2\beta_{2} to be independent.

?

Let

Y1=θ1+θ2+ε1Y2=θ1−2θ2+ε2 \begin{aligned} & Y_{1} = \theta _{1}+\theta _{2}+\varepsilon _{1} \\ & Y_{2} = \theta _{1}-2 \theta _{2}+\varepsilon _{2} \end{aligned}

and

Y3=2θ1−θ2+ε3 Y_{3} = 2 \theta _{1}-\theta _{2}+\varepsilon _{3}

where E[εi]=0(i=1,2,3)\mathbb {E}\left[\varepsilon_{i}\right] = 0(i = 1,2,3). Find the least squares estimates of ∂1\partial_{1} and θ2\theta_{2}. If the equations above are augmented to

Y1=θ1+θ2+θ3+ε1,Y2=θ1−2θ2+θ3+ε2,Y3=2θ1−θ2+θ3+ε3 \begin{aligned} & Y_{1} = \theta _{1}+\theta _{2}+\theta _{3}+\varepsilon _{1}, \\ & Y_{2} = \theta _{1}-2 \theta _{2}+\theta _{3}+\varepsilon _{2}, \\ & Y_{3} = 2 \theta _{1}-\theta _{2}+\theta _{3}+\varepsilon _{3} \end{aligned}

find the least squares estimate of θ3\theta_{3}.

?

Suppose that E[Yt]=β0+β1cos⁡(2πk1t/n)+β2sin⁡(2πk2t/n)\mathbb {E}\left[Y_{t}\right] = \beta_{0}+\beta_{1} \cos \left(2 \pi k_{1} t / n\right)+\beta_{2} \sin \left(2 \pi k_{2} t / n\right), where t=1,2,…,nt = 1,2, \ldots , n, and k1k_{1} and k2k_{2} are positive integers. Find the least squares estimates of β0,β1\beta_{0}, \beta_{1}, and β2\beta_{2}.

?

Suppose that

Xij=αi+θtij+εij,i=1,…,a, j=1,…,b, X_{ij}=\alpha _i+\theta t_{ij}+\varepsilon _{ij},\quad i=1,\ldots ,a,\ j=1,\ldots ,b,

where αi\alpha_i and θ\theta are unknown parameters, tijt_{ij} are known constants, and εij\varepsilon_{ij} are independent and identically distributed random variables with mean 00. Find explicit forms for the least squares estimators (LSE's) of θ\theta, αi\alpha_i, i=1,…,ai=1,\ldots ,a.

?
Remark.
?

Corresponding exercise in Shao (2003): 3.58.

Consider the linear model X=Zβ+εX=Z\beta +\varepsilon, where ZZ is a known n×pn\times p matrix, β\beta is a pp-vector of unknown parameters, and ε\varepsilon is a random nn-vector whose components are independent and identically distributed with mean 00 and variance σ2\sigma^2. Let Z=(Z1,Z2)Z=(Z_1,Z_2) and β=(β1,β2)\beta =(\beta_1,\beta_2), where ZjZ_j is n×pjn\times p_j and βj\beta_j is a pjp_j-vector, j=1,2j=1,2. Assume that (Z1τZ1)−1(Z_1^\tau Z_1)^{-1} and [Z2τZ2−Z2τZ1(Z1τZ1)−1Z1τZ2]−1[Z_2^\tau Z_2-Z_2^\tau Z_1(Z_1^\tau Z_1)^{-1}Z_1^\tau Z_2]^{-1} exist.

?
(i)

Derive the LSE of β\beta in terms of Z1Z_1, Z2Z_2, and XX.

(ii)

Let β^=(β^1,β^2)\widehat{\beta }=(\widehat{\beta }_1,\widehat{\beta }_2) be the LSE in (i). Calculate the covariance between β^1\widehat{\beta }_1 and β^2\widehat{\beta }_2.

(iii)

Suppose that it is known that β2=0\beta_2=0. Let β~1\widetilde\beta_1 be the LSE of β1\beta_1 under the reduced model X=Z1β1+εX=Z_1\beta_1+\varepsilon. Show that, for any l∈Rp1l\in \mathbb {R}^{p_1}, lτβ~1l^\tau \widetilde\beta_1 is better than lτβ^1l^\tau \widehat{\beta }_1 in terms of their variances.

Remark.
?

Corresponding exercise in Shao (2003): 3.70.

References

[Dem97] Demidovich B.P. Сборник задач и упражнений по математическому анализу. 13-е изд., Издательство Московского университета / Издательство "ЧеРо", 1997.
[ГПП12] Горяинова Е.Р., Панков А.Р., Платонов Е.Н. Прикладные методы анализа статистических данных. Издательский дом Высшей школы экономики, 2012.
[LeCa98] Lehmann E.L., Casella G. Theory of Point Estimation. 2nd ed., Springer, 1998.
[ЖРШ16] Жуковский М.Е., Родионов И.В., Шабанов Д.А. Введение в математическую статистику. МФТИ, 2016.
[SeLe03] Seber G.A., Lee A.J. Linear Regression Analysis. 2nd ed., John Wiley & Sons, 2003.
[Sha05] Shao J. Mathematical Statistics: Exercises and Solutions. Springer, 2005.