4.II.19H

Statistics
Part IB, 2008

(i) Consider the linear model

Yi=α+βxi+εiY_{i}=\alpha+\beta x_{i}+\varepsilon_{i}

where observations Yi,i=1,,nY_{i}, i=1, \ldots, n, depend on known explanatory variables xix_{i}, i=1,,ni=1, \ldots, n, and independent N(0,σ2)N\left(0, \sigma^{2}\right) random variables εi,i=1,,n\varepsilon_{i}, i=1, \ldots, n.

Derive the maximum-likelihood estimators of α,β\alpha, \beta and σ2\sigma^{2}.

Stating clearly any results you require about the distribution of the maximum-likelihood estimators of α,β\alpha, \beta and σ2\sigma^{2}, explain how to construct a test of the hypothesis that α=0\alpha=0 against an unrestricted alternative.

(ii) A simple ballistic theory predicts that the range of a gun fired at angle of elevation θ\theta should be given by the formula

Y=V2gsin2θY=\frac{V^{2}}{g} \sin 2 \theta

where VV is the muzzle velocity, and gg is the gravitational acceleration. Shells are fired at 9 different elevations, and the ranges observed are as follows:

θ (degrees) 51525354555657585sin2θ0.17360.50.76600.939710.93970.76600.50.1736Y( m)4322118981748520664212961949115572100273458\begin{array}{cccccccccc}\theta \text { (degrees) } & 5 & 15 & 25 & 35 & 45 & 55 & 65 & 75 & 85 \\ \sin 2 \theta & 0.1736 & 0.5 & 0.7660 & 0.9397 & 1 & 0.9397 & 0.7660 & 0.5 & 0.1736 \\ Y(\mathrm{~m}) & 4322 & 11898 & 17485 & 20664 & 21296 & 19491 & 15572 & 10027 & 3458\end{array}

The model

Yi=α+βsin2θi+εiY_{i}=\alpha+\beta \sin 2 \theta_{i}+\varepsilon_{i}

is proposed. Using the theory of part (i) above, find expressions for the maximumlikelihood estimators of α\alpha and β\beta.

The tt-test of the null hypothesis that α=0\alpha=0 against an unrestricted alternative does not reject the null hypothesis. Would you be willing to accept the model ()(*) ? Briefly explain your answer.

[You may need the following summary statistics of the data. If xi=sin2θix_{i}=\sin 2 \theta_{i}, then xˉn1xi=0.63986,Yˉ=13802,Sxx(xixˉ)2=0.81517,Sxy=Yi(xixˉ)=\bar{x} \equiv n^{-1} \sum x_{i}=0.63986, \bar{Y}=13802, S_{x x} \equiv \sum\left(x_{i}-\bar{x}\right)^{2}=0.81517, S_{x y}=\sum Y_{i}\left(x_{i}-\bar{x}\right)= 17186. ]