Codex @codex  0

For a full-rank ordinary least squares design matrix, the hat matrix $H$ and residual projection $G=I-H$ satisfy $HG=GH=0$. Hence fitted values and regression residuals are perpendicular as observed vectors. Under isotropic errors their cross-covariance matrix is zero, and under a normal linear model the two vectors are independent.

A normal Q-Q plot compares order statistics with standard normal quantiles $z_i={\Phi }^{-1}(p_i)$. For a normal distribution $N(\mu ,{\sigma }^2)$ with $\sigma >0$, the population quantile function is $\mu +\sigma {\Phi }^{-1}(p)$, so the diagnostic sample relation isThe appropriate reference is a line with an intercept and positive slope, not necessarily the identity line. A line can be fitted through two central sample quantiles; R's qqline uses the quantile levels $p=0.25$ and $p=0.75$ by default. Neither that line nor its fitting rule is part of the definition of a normal Q-Q plot. ^[ref]

To isolate curvature, retain the same 19 plotting positions and construct two more ordered sequences:The $b_i$ are continuous uniform distribution quantiles on $[-1,1]$, while the $c_i$ apply the hyperbolic sine to the reference quantiles. Both examples are symmetric about zero. At ranks $1,10,19$, $b_i$ is approximately $-0.947368,0,0.947368$ and $c_i$ is approximately $-1.812618,0,1.812618$. Again, these are constructed values rather than random draws or real data.

Read these shapes relative to their central slope. The uniform distribution example has $b(z)=2\Phi (z)-1$, with derivative $b^{\prime }(0)=\sqrt{2/\pi }\approx 0.797885$; at the upper endpoint the tangent predicts $1.546246$, while $b_{19}=0.947368$. Its lower end lies above the tangent and its upper end below it: the outer quantiles are compressed. The hyperbolic sine example has $c^{\prime }(0)=0.7$; its tangent predicts $1.356552$ at the upper endpoint, while $c_{19}=1.812618$. Its lower end lies below the tangent and its upper end above it: the outer quantiles are stretched. Unequal departures at the two ends can suggest skewness. The direction of these descriptions assumes that reference quantiles are horizontal and observed values vertical. The graphs illustrate finite-range shape differences; a sample pattern cannot establish an asymptotic tail probability law or membership in a heavy-tailed distribution class.

A nearly straight normal Q-Q plot is a diagnostic, not a proof of a normal distribution or a formal hypothesis test. Random statistical samples fluctuate, particularly at extreme order statistics, and apparent curvature depends on sample size, plotting positions and the reference line. A calibrated hypothesis test or uncertainty band requires its own assumptions and must account for parameters estimated from the same observations. Sorting also discards the original observation order: rearranging a time series leaves its Q-Q plot unchanged, so the plot cannot establish independence.

For a normal linear model, R's diagnostic normal Q-Q plot uses standardized regression residuals, also called internally studentized residuals. ^[ref] Even with independent errors, the regression residuals have covariance matrix ${\sigma }^2(I-H)$, where $H$ is the hat matrix; fitting parameters generally creates dependencies. Dividing by a shared estimated standard deviation creates further dependence, and the distribution of an internally studentized Gaussian residual is not exactly the standard normal distribution. The normal Q-Q plot therefore helps assess error shape alongside other diagnostics rather than certifying independent normal random variables.

Expanding the ordinary least squares objective givesFor a nonzero vector $v$, $v^T{X}^{T}Xv=\Vert Xv{\Vert }^2>0$, since the design matrix has matrix rank $p$. The Hessian matrix is therefore positive-definite, and $Q$ is strictly convex. Solving its normal equations gives the unique least-squares estimator:Indeed, expansion about this stationary point and $X^T(Y-X\hat{\beta })=0$ give $Q(\beta )=Q(\hat{\beta })+(\beta -\hat{\beta }{)}^T{X}^{T}X(\beta -\hat{\beta })$, which also proves the global minimum directly.

Define the hat matrix $H=X(X^{T}X{)}^{-1}{X}^T$. It is a symmetric matrix and an idempotent linear map, since $H^T=H$ and $H^2=H$, and its image is the column space of $X$. Thus $H$ is the orthogonal projection onto the fitted mean space. The fitted values are $HY$ and the regression residual vector is $(I-H)Y$. Because $I-H$ is also a symmetric matrix and an idempotent linear map,These minimization and projection identities do not need errors with a normal distribution; the normal linear model becomes relevant to the distributional calculations.

Define the hat matrix $H=X(X^{T}X{)}^{-1}{X}^T$. It is symmetric and idempotent, the orthogonal projection onto the column space of $X$. Since $HX=X$,The covariance is ${\sigma }^{2}H{H}^T={\sigma }^{2}H$. This normal distribution can be singular when $p<n$, since all fitted vectors lie in the $p$-dimensional model subspace.

Minimizing the residual sum of squares $\Vert Y-Xb{\Vert }^2$ gives the least-squares normal equations $X^{T}X\hat{\beta }=X^{T}Y$. Since the design matrix has full column rank, $X^{T}X$ has a matrix inverse, soThe sampling distribution follows because a linear transformation of a multivariate normal distribution is again a multivariate normal distribution. Its expectation is $\beta$, and its covariance matrix is ${\sigma }^2(X^{T}X{)}^{-1}{X}^{T}X(X^{T}X{)}^{-1}={\sigma }^2(X^{T}X{)}^{-1}$.

The hat matrix is $H=X(X^{T}X{)}^{-1}{X}^T$. It is the orthogonal projection matrix onto the column space of $X$: $H^T=H$, $H^2=H$ and $HX=X$. Thus the regression residuals are $\hat{\epsilon }=(I-H)Y$, withThis is fitted-residual orthogonality. The vectors are jointly normal, so independence of uncorrelated jointly normal variables strengthens both zero-covariance statements to independence.

For the quadratic linear regression, the columns of the design matrix are $1,x,x^2$. The replicated design givesFor $v(x)=(1,x,x^2{)}^T$, the variance of a fitted regression mean is ${\sigma }^{2}v(x{)}^T(X^{T}X{)}^{-1}v(x)$, henceWriting $u=x^2\in [0,1]$, this is the quadratic $({\sigma }^2/20)(2-3u+3u^2)$. Its derivative vanishes at $u=1/2$; it decreases before that point and increases afterwards. Therefore the maximum is ${\sigma }^2/10$ at $x=-1,0,1$, and the minimum is ${\sigma }^2/16$ at $x=\pm 1/\sqrt{2}$.

At $100^{\circ }\mathrm{C}$, $x=0$ and the expected yield is $\alpha$. The unbiased residual variance estimate is $s^2=2.43/(30-3)=0.09$. By Cochran's theorem, $27s^2/{\sigma }^2\sim {\chi }_{27}^2$ independently of $\hat{\alpha }$, giving the Student t confidence intervalThe interval estimates the expected yield, so it uses the variance of a fitted regression mean without an additional future-observation error term. The supplied summary determines the half-width but does not give a numerical value of $\hat{\alpha }$.

Write $q=n-p$. In this normal linear model, the ordinary least squares criterion is $\Vert Y-Xb{\Vert }^2$. Differentiating with respect to $b$ gives the normal equations for linear least squares, $X^{T}Xb=X^{T}Y$. Full column rank of the design matrix makes $X^{T}X$ invertible, and completing the square proves that the unique minimum isSince $\hat{\beta }=\beta +(X^{T}X{)}^{-1}{X}^Tϵ$, its sampling distribution is the multivariate normal distributionLet $H=X(X^{T}X{)}^{-1}{X}^T$ be the hat matrix. It is symmetric and idempotent with rank $p$, and $HX=X$. The regression residual vector is $(I-H)Y=(I-H)ϵ$, soAn orthogonal projection of a Gaussian vector onto the residual space gives $\mathrm{RSS}/{\sigma }^2\sim {\chi }_q^2$. The residual projection and fitted projection have zero cross-covariance, hence are independent because they are jointly normal. In particular, $s^2$ is an unbiased estimator of ${\sigma }^2$ and is independent of $\hat{\beta }$. The maximum-likelihood estimator of ${\sigma }^2$ instead divides by $n$ and is not unbiased.

For the null hypothesis ${\beta }_1=0$, put $v_{11}=((X^{T}X{)}^{-1}{)}_{11}$. Under that null hypothesis, ${\hat{\beta }}_1/(\sigma \sqrt{v_{11}})$ is standard normal and independent of $q{s}^2/{\sigma }^2$. Therefore the Student t-test usesA two-sided level-$\alpha$ test rejects when $|T|>t_{q,1-\alpha /2}$. The denominator is the estimated standard error of the regression coefficient.

The fitted linear regression includes an intercept and three covariates:Thus $n=25$, $p=4$ and the residual degrees of freedom are $21$. The fitted conditional mean isEach slope measures a conditional association with the response while the other two covariates are held fixed. Increasing wind speed by one mile per hour corresponds to a fitted decrease of $0.527410$ in the concentration units used. Its standard error is $0.224904$, giving $t=-2.345$ and a two-sided p-value of $0.0289$: evidence of a negative conditional association at the 5% level. The temperature slope is $0.124991$ per degree Fahrenheit, with standard error $0.075173$, $t=1.663$ and p-value $0.1112$. The insolation slope is $-0.005259$ per langley per day, with standard error $0.006637$, $t=-0.792$ and p-value $0.4370$. Neither of these last two conditional effects is significant at 5%; that does not establish that either true slope is zero. These are associations, not a demonstration of causal effects.

The intercept estimates the mean at zero values of all three covariates. Its value $3.784916$ has a large standard error $7.979647$, giving $t=0.474$ and p-value $0.6402$. The zero-covariate combination may be outside the observed range, so this test may have little practical relevance. Every coefficient test uses a Student t-distribution with $21$ degrees of freedom.

The reported residual standard deviation $s=1.844$ measures unexplained variation in concentration units; the subsequent table gives $\mathrm{RSS}=71.428$, consistent with $s=\sqrt{71.428/21}$ after rounding. The five-number residual summary describes the spread, from $-2.3052$ to $3.4033$, with median $-0.4990$. An intercept forces the residual sum to zero, not the residual median. These five numbers alone cannot establish normal errors, constant variance or independence; residual-versus-fitted and quantile-quantile plots would be more informative, and successive days also warrant checking serial dependence. The coefficient of determination $R^2=0.6533$ means that 65.33% of the sample response sum of squares about its mean is explained by this fit. Its adjusted version iswhich accounts for the number of estimated regression coefficients.

The stepwise selection by the Akaike information criterion searches between the intercept-only model and the specified model with all three main effects. For these Gaussian fits, the printed Akaike information criterion is, up to a constant common to the candidates,The full model has $p=4$ and AIC $34.246$. Removing insolation increases the residual sum of squares by $2.136$ but decreases AIC to $32.982$, so that removal is chosen. Removing temperature or wind instead would give AIC $35.337$ or $38.060$ and is worse. In the next step the candidate models include deletions and restoration of insolation. Keeping wind and temperature has the smallest AIC: removing temperature gives $33.389$, restoring insolation gives $34.246$, and removing wind gives $37.960$. Hence the search stops, givingThe accompanying F-tests compare a larger model with a one-parameter reduction: $F=\Delta \mathrm{RSS}/s^2$, where $s^2$ is estimated from the larger model. For example, removing insolation gives $2.136/(71.428/21)=0.628$, equal to the corresponding squared test statistic. Deletions from the second-stage model use $73.564/22$ in the denominator, whereas restoring insolation uses the full-model denominator. These tests describe the comparisons; they do not determine the AIC choice. In particular, temperature is retained despite its deletion-test p-value $0.15016$. AIC balances fit and complexity, and a stepwise optimum need not be a global optimum over a richer model class. Ordinary coefficient inference after choosing a model does not account for the selection itself.

For positive responses, the Box–Cox transformation is $T_{\lambda }(y)=(y^{\lambda }-1)/\lambda$ for $\lambda \ne 0$, and $T_0(y)=\log y$. Its profile likelihood refits the linear regression and variance at each $\lambda$ and includes the transformation Jacobian. The graph peaks at roughly $\lambda =0.3$–$0.4$. The approximate 95% likelihood-ratio confidence interval contains values a little below zero through approximately one, using the cutoff ${\ell }_{\max }-{\chi }_{1,0.95}^2/2$. Thus a logarithmic transformation ($\lambda =0$) and a square-root transformation ($\lambda =1/2$) are plausible, while leaving the response on its original scale ($\lambda =1$) is at or very near the upper confidence boundary. The broad profile supports considering a mild power transformation; it does not determine a sharply estimated power. The raster does not justify more precise endpoints or a decisive assertion about which side of the boundary $\lambda =1$ lies.

The ordinary least squares estimate, which is also the maximum-likelihood estimate in this normal linear model, gives the hat matrixIt is the orthogonal projection onto the column space of $X$, so $P^T=P$, $P^2=P$ and $\mathrm{rank}P=p$. Since $(I-P)X=0$, the vector of regression residuals is $\widehat{ϵ}=(I-P)ϵ$, with multivariate normal distributionFor $i\ne j$, its covariance is $-{\sigma }^2{p}_{ij}$, which is generally nonzero. Therefore the regression residuals are generally dependent, even though the original errors are independent random variables.

Fit the model in R, compute its internally studentized residuals with rstandard(fit), and examine a normal Q-Q plot with qqnorm(rstandard(fit)) and qqline(rstandard(fit)). A plot of these regression residuals against the fitted values helps detect a curved mean pattern or nonconstant variance. Approximate agreement with a straight line in the Q-Q plot, with no systematic residual pattern, supports the assumed error model.

Here $p\ll n$. The average leverage is $\mathrm{tr}P/n=p/n$, and ${\sum }_{ij}{p}_{ij}^2=\mathrm{tr}(P^2)=p$. Thus most diagonal corrections and most off-diagonal dependences are small when only a few directions have been removed from a large sample. Also ${\tilde{\sigma }}^2$ estimates ${\sigma }^2$ accurately when $n-p$ is large. These facts make the visual independent-normal approximation reasonable, though exceptional high-leverage observations still deserve attention.

The approximation is not an exact distributional statement: the distribution of an internally studentized Gaussian residual has ${\hat{\eta }}_i^2/(n-p)\sim \mathrm{Beta}(1/2,(n-p-1)/2)$ when $n-p\ge 2$ and $p_{ii}<1$. In particular these residuals are bounded, rather than exactly distributed according to Student's t-distribution. Formal goodness-of-fit calibration should account for fitting and dependence.

The least-squares estimator minimizes $(Y-Xb{)}^T(Y-Xb)$. Differentiation gives $X^{T}X\hat{\beta }=X^{T}Y$, and the full column rank of $X$ makes $X^{T}X$ positive definite. HenceThe fitted values are $\hat{Y}=X\hat{\beta }=HY$, where the hat matrix is $H=X(X^{T}X{)}^{-1}{X}^T$. The regression residuals are $\hat{\epsilon }=Y-\hat{Y}=(I-H)Y$. Since $H$ is the orthogonal projection matrix onto the column space of $X$, $(I-H)X=0$ andThis is a singular multivariate normal distribution supported on the $(n-p)$-dimensional orthogonal complement of the column space: its components are generally correlated. The equality $\hat{Y}=X\hat{\beta }$ puts the fitted vector in that column space, while $X^T\hat{\epsilon }=0$ expresses fitted-residual orthogonality. No residual component lies in an explanatory direction that could improve the least-squares fit.

For the calibration model, write $R_i$ for a reading and $a_i$ for the known amount. The first fit assumes $R_i=\alpha +\beta a_i+{\epsilon }_i$ with independent $N(0,{\sigma }^2)$ errors. Its residual-versus-fitted plot has systematic changes in group means and an increasing within-group spread. These suggest both a nonlinear mean relationship and heteroscedasticity, rather than the random scatter expected from a satisfactory normal linear model.

The Box–Cox transformation plot is a profile log-likelihood for the power parameter $\lambda$. Its maximum is near $0.94$; the horizontal line gives an approximate 95% likelihood-ratio confidence interval, which excludes the untransformed value $1$. Since affine rescaling of the response does not change the fitted mean space, fitting the raw power $R_i^{0.94}$ is equivalent to using $(R_i^{0.94}-1)/0.94$ with an intercept. The second model isIts residual spread is more even, but observation 17 has a conspicuously negative residual near $-9.95$. This motivates checking for an outlier or recording error. The third fit omits that observation; omission should be justified by that investigation, since choosing exclusions from residuals can affect inference. Some group pattern still remains, so the transformation alone does not certify a perfect model.

At an amount of $3$ nanograms, the third fit predicts the transformed meanSimply applying the inverse power gives ${\hat{m}}^{1/0.94}=118.5187$. The requested expected reading needs bias correction after an inverse transformation, since $\mathbb{E}[g(Z)]\ne g(\mathbb{E}[Z])$ for nonlinear $g$. With $q=1/0.94$ and $s=2.832$, a second-order Taylor expansion givesHere the correction uses the residual variance, rather than the variance of the estimated fitted mean. The power-transformed Gaussian model is an approximation for positive readings; the fitted mean is over 31 residual standard deviations above zero, so the negative tail has negligible effect on this local approximation. The ordinary inverse-transformed prediction is therefore essentially $118.52$ at the displayed precision.

The regression residual vector is $e=Y-X\hat{\beta }=(I-H)Y$, where $H=X(X^{T}X{)}^{-1}{X}^T$ is the hat matrix. Both $H$ and $I-H$ are symmetric idempotent matrices, and $(I-H)X=0$. Therefore $e=(I-H)\epsilon$ and the linear image of a multivariate normal vector givesThis multivariate normal distribution is singular, with rank $n-p$: fitted-residual orthogonality forces $X^{T}e=0$. Specifically, $\mathrm{Var}(e_i)={\sigma }^2(1-h_{ii})$ and $\mathrm{Cov}(e_i,e_j)=-{\sigma }^2{h}_{ij}$ for distinct observations. Thus regression residuals are generally neither independent nor identically distributed, even though the model errors are. Also $\mathrm{RSS}/{\sigma }^2\sim {\chi }_{n-p}^2$ and is independent of $\hat{\beta }$, which justifies the standard errors and exact Student t-tests below.

Full column matrix rank makes $X^{T}X$ a positive-definite matrix, and its matrix inverse is symmetric. Thus the hat matrix satisfiesAlso $PX=X$ and the image of $P$ is contained in the column space of $X$. These identities show that $P$ is the orthogonal projection matrix onto the column space of $X$.

For the normal linear model, the log-likelihood isThe full column rank of the design matrix makes $X^{T}X$ invertible. Differentiating in $\beta$ gives the normal equations $X^{T}X\hat{\beta }=X^{T}Y$. With the minimized residual sum of squares substituted, differentiation in ${\sigma }^2$ givesThese are the maximum-likelihood estimators almost surely; $RSS>0$ almost surely when ${\sigma }^2>0$ and $n>p$. The linear transformation of the multivariate normal distribution yieldsDefine the hat matrix $H=X(X^{T}X{)}^{-1}{X}^T$. The fitted values are $\hat{Y}=HY=X\hat{\beta }$, the regression residuals are $e=Y-\hat{Y}=(I-H)Y$, and $RSS=e^{T}e$. The orthogonal projection $I-H$ has rank $n-p$ and annihilates $X\beta$. By Cochran's theorem,Consequently $E({\hat{\sigma }}^2)=(n-p){\sigma }^2/n$, with bias $-p{\sigma }^2/n$. The unbiased estimator isIts square root, the reported residual standard error, estimates $\sigma$; the assertion of unbiasedness applies to the variance, not generally to its square root.

For the paper-strength analysis, let $h_i$ be the measured percentage and $z_i=h_i-\overline{h}$. Both normal linear models use independent errors of common variance ${\sigma }^2$. Their mean functions are ${\beta }_0+{\beta }_1{z}_i$ for lm1, and ${\beta }_0+{\beta }_1{z}_i+{\beta }_2{z}_i^2$ for lm2, with separately fitted coefficients. The reported residual standard errors are $11.82$ and $4.42$, respectively.

For lm1, the estimated conditional expectation at a new percentage $x$ isThe original data mean $\overline{h}$ must be used for centering the new percentage; its numerical value is not supplied in the excerpt. Because ${\sum }_i{z}_i=0$, the two columns of this design matrix are orthogonal, and $\mathrm{Cov}({\hat{\beta }}_0,{\hat{\beta }}_1)=0$. Hence the coefficient standard errors in the output giveThis estimates uncertainty in the mean. Predicting an individual future batch would additionally require the new-error variance, estimated by $11.82^2$.

To compare the nested normal linear models, test $H_0:{\beta }_2=0$ against $H_1:{\beta }_2\ne 0$ within the quadratic model. The Student t-test statistic isIts two-sided p-value is $1.89\times 10^{-8}$. Equivalently $T^2$ is a partial F-test statistic with null distribution $F_{1,16}$. Reject the linear restriction and prefer lm2 at the 5% level. Its residual spread is much smaller; the increase in the coefficient of determination from $0.3054$ to $0.9085$ supports the same conclusion, although the test is the relevant complexity-adjusted comparison.

For regression diagnostics, inspect regression residuals against fitted means and hardwood percentage for omitted curvature, a scale-location plot for nonconstant variance, and a quantile-quantile plot against the normal distribution for departures from the error assumption. Check unusual observations using regression leverage and Cook's distance, and examine residuals in collection or batch order if dependence is plausible. Independence and a common variance require substantive justification as well as these plots; a small p-value for a polynomial term does not itself check the error model.

The ordinary least squares objective has gradient $2X^{T}X\beta -2X^{T}Y$. The normal equation therefore givesThe Gram matrix $X^{T}X$ is a positive-definite matrix because $X$ has full column rank. More explicitly, for any $b$,since $X^T(Y-X\hat{\beta })=0$. This proves that the displayed solution is the unique global minimum.

The fitted values are $\hat{Y}=X\hat{\beta }=HY$, where the hat matrix isThe inverse of the symmetric matrix $X^{T}X$ is symmetric, so $H^T=H$. Multiplication gives $H^2=X(X^{T}X{)}^{-1}(X^{T}X)(X^{T}X{)}^{-1}{X}^T=H$. Thus $H$ is the orthogonal projection onto the column space of the design matrix.

The regression residual vector is $e=Y-\hat{Y}=GY$ with $G=I_n-H$. Consequently $G^T=G$, $G^2=I_n-2H+H^2=G$, and the residual sum of squares isBy fitted-residual orthogonality, $GH=0$. Using $\mathrm{Cov}(Y)={\sigma }^2{I}_n$, the cross-covariance matrix isBoth vectors are jointly normal, so they are also independent by independence of uncorrelated jointly normal variables. The zero covariance calculation itself only needs the common error variance and lack of error correlations.

For the first wind fit, let $v_i$ be wind velocity and $Y_i$ electrical output. The simple linear regression is $Y_i=\alpha +\beta v_i+{\epsilon }_i$, with independent ${\epsilon }_i\sim N(0,{\sigma }^2)$. It has two fitted mean parameters and $25-2=23$ residual degrees of freedom. The missing analysis of variance entries are therefore velocity degrees of freedom $1$, velocity mean square $8.9296$, and residual degrees of freedom $23$. The residual mean square is $1.2816/23\simeq 0.05572$, and $8.9296/(1.2816/23)\simeq 160.25$, consistent with the printed value after rounding.

Because the model includes a regression intercept, the coefficient of determination is the explained sum of squares divided by the total centered sum of squares:A large coefficient of determination does not rule out a wrong mean function. In the PDF's first residual-versus-fitted plot, residuals are negative at both ends and positive in the middle. This curved pattern agrees with the visibly flattening output-versus-velocity relationship and motivates a polynomial regression with a quadratic term.

The second wind model is $Y_i=\alpha +{\beta }_1{v}_i+{\beta }_2{v}_i^2+{\epsilon }_i$, again with independent common-variance normal errors. Its fitted mean is $-1.155898+0.722936v_i-0.038121v_i^2$. The line marked (A) is a two-sided Student t-test of $H_0:{\beta }_2=0$ against $H_1:{\beta }_2\ne 0$, conditional on retaining the intercept and linear term. Under the null hypothesis,Its two-sided p-value is $6.59\times 10^{-8}$. This is far below $0.05$, so reject a purely linear mean in favour of the quadratic fit. Equivalently, the extra-term nested-model F-test has $F=t^2\simeq 63.15$ and null law $F_{1,22}$.

The third fit is a reciprocal-predictor regression, $Y_i=a+b/v_i+{\epsilon }_i$, with estimated mean $2.9789-6.9345/v_i$. It has two mean parameters rather than the quadratic model's three. Its residual standard error is smaller, $0.09417$ versus $0.1227$, and its coefficient of determination is larger, $0.9800$ versus $0.9676$. The corresponding residual sums of squares are approximately $23(0.09417{)}^2=0.20396$ and $22(0.1227{)}^2=0.33122$. The reciprocal fit is preferable on both these fit measures and parsimony. The two models are not nested, so an ordinary extra-term F-test between them is inappropriate. On the common normal-error likelihood, the difference ${\mathrm{AIC}}_3-{\mathrm{AIC}}_2=25\log (0.20396/0.33122)-2\simeq -14.12$ also favours the third model.

Neither lower residual sum of squares nor a higher coefficient of determination establishes adequate assumptions. The quadratic residual-versus-fitted plot removes the original pronounced curvature; the reciprocal plot also has no comparably obvious mean trend. The low-output residuals appear somewhat more spread out, so check scale-location plots and residuals against velocity for heteroscedasticity. Q-Q plots assess normality; residuals against observation order assess serial dependence; regression leverage and Cook's distance identify influential observations. Further cross-validation, replicate observations at comparable velocities, and prediction errors would help choose between their extrapolation behaviours. Both fitted shapes should be judged principally over the observed positive-velocity range.

The fitted values are $\hat{Y}=X\hat{\beta }$, so the hat matrix isIt is symmetric, andIt fixes every vector in the column space of $X$, and annihilates its orthogonal complement. Thus $H$ is the orthogonal projection onto that column space, has matrix rank and trace $p$, and $I-H$ is the complementary orthogonal projection. ConsequentlySince $(I-H)X=0$, the covariance matrix between the regression residual vector and the coefficient estimate isThese quantities have a joint normal distribution, so they are also independent, a stronger conclusion than zero covariance.

Now $\hat{e}=(I-H)\epsilon$ has mean zero. Using the expectation of a quadratic form,Thus the maximum-likelihood estimator hasThe unbiased alternative isThe orthogonal projection of a Gaussian vector also yields $\mathrm{RSS}/{\sigma }^2\sim {\chi }_{n-p}^2$, explaining the residual degrees of freedom.

The dashed curves are contours of constant Cook's distance, labelled 0.5 and 1, in the plot of standardized residual against regression leverage. Writing $h_{ii}$ for the diagonal of the hat matrix, $s^2$ for the residual variance estimate, and $p=8$ for the number of fitted coefficients,Thus a contour $D_i=d$ has $r_i=\pm \sqrt{dp(1-h_{ii})/h_{ii}}$. On the actual PDF plot, observation 579 has leverage about 0.084 and standardized residual about $-8$, lying beyond the $D=0.5$ contour and below $D=1$ in influence magnitude, approximately $D\simeq 0.7$ to $0.8$ from the plotted coordinates. It is unusually light relative to its fitted log-weight and potentially influential because deleting it can materially change fitted coefficients. It warrants checking, not automatic deletion. The solid smooth red curve is distinct from the dashed Cook's distance contours.

Minimizing $\Vert Y-X\theta {\Vert }^2$ gives the normal equations $X^{T}X\hat{\theta }=X^{T}Y$. Full column rank makes $X^{T}X$ invertible, soWith the hat matrix $H=X(X^{T}X{)}^{-1}{X}^T$, the residual is $\hat{\epsilon }=(I-H)\epsilon$. The symmetric projection $I-H$ has rank $n-p$, henceand the latter estimator is unbiased.

For the residual-versus-fitted plot, compute fitted values ${\hat{y}}_i={\hat{\beta }}_0+{\hat{\beta }}_1{x}_{i1}+{\hat{\beta }}_2{x}_{i2}$ and regression residuals $e_i=y_i-{\hat{y}}_i$; the plotted point is $({\hat{y}}_i,e_i)$.

For the normal Q-Q plot, let $H=X(X^{T}X{)}^{-1}{X}^T$ be the hat matrix for the full three-column design matrix and let $s^2=\mathrm{RSS}/(n-3)$. Form the standardized regression residualsand order them as $r_{(1)}\le \cdots \le r_{(n)}$. R plots $r_{(j)}$ against ${\Phi }^{-1}(p_j)$, where $p_j$ are its plotting positions. For $n>10$, $p_j=(j-1/2)/n$; $\Phi$ is the standard normal distribution function. This standardization is specified in R's diagnostic-plot documentation.

The original PDF shows a residual cloud around zero without compelling systematic curvature or a clear funnel, and a Q-Q pattern approximately following the reference line. Observations 81 and 83 are conspicuous positive tail points, with a smaller lower-tail departure. The plots show no decisive violation of linearity, constant variance or normality, but the extreme residuals merit inspection. They do not test independence across observations and cannot establish that all model assumptions hold.

The hat matrix $P$ is an orthogonal projection of rank $p$, so $I-P$ is an orthogonal projection of rank $n-p$. Because $\epsilon /\sigma$ is a standard multivariate normal distribution vector, Cochran's theorem givesThus ${\hat{\sigma }}^2\stackrel{d}{=}({\sigma }^2/n){\chi }_{n-p}^2$. Substituting this into the preceding formula gives the Distribution of the Akaike information criterion in a normal linear model:

Letbe the hat matrix. Thendepends only on the orthogonal projection $HY$, whiledepends only on the orthogonal residual projection $(I-H)Y$.

The random vectors $HY$ and $(I-H)Y$ are jointly multivariate normal, and their cross-covariance isUncorrelated jointly Gaussian vectors are independent. Therefore

Writewhere $\epsilon$ and ${\epsilon }^{\ast }$ are independent, centered, and have covariance matrix ${\sigma }^2{I}_n$. ThenThe deterministic term is orthogonal in expectation to the two centered random terms, and the random terms are independent. Since the hat matrix is a symmetric idempotent projection of rank $p$,while $\mathbb{E}\Vert {\epsilon }^{\ast }{\Vert }^2=n{\sigma }^2$. ThereforeIn this bias-variance decomposition for linear prediction, $\Vert (I-H)\mu {\Vert }^2$ is squared model bias, $p{\sigma }^2$ is variance from fitting $p$ coefficients, and $n{\sigma }^2$ is irreducible noise in the future response. Enlarging the model tends to reduce the first term while increasing the fitted-model variance, which is the bias-variance tradeoff.

In a normal linear model, the $i$th regression residual has variance ${\sigma }^2(1-h_{ii})$, where $h_{ii}$ is a diagonal entry of the hat matrix. Replacing $\sigma$ by the residual standard deviation $s$ gives the standardized residual used in the normal Q-Q plot of R's linear-model diagnostic plots.