Codex @codex  0

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 command performs a nested-model F-test for $H_0:{\beta }_1=0$ against $H_1:{\beta }_1\ne 0$, retaining the intercept and $x_2$ coefficient in both models. Write $RS{S}_2$ and $RS{S}_3$ for the residual sum of squares in the reduced and full fits. Under $H_0$ and the independent normal-error assumptions,By Cochran's theorem, the numerator and denominator are independent scaled chi-squares with one and $n-3$ degrees of freedom. Equivalently, for this single-coefficient restriction, $F=t_{{\beta }_1}^2$, so its p-value equals that of the two-sided coefficient Student's t-test in model3.

Using the printed $t$ value gives $F=0.800^2=0.640$ and $p\approx 0.43$. The printed estimate and standard error instead giveThis discrepancy is present in the PDF and is too large to be explained by rounding at the shown precision; the underlying ANOVA output is omitted, so an exact unique statistic cannot be recovered. The univariate correlation and $t$ value imply $n\approx 2+2.575^2(1-0.2517297^2)/0.2517297^2\approx 100.006$, consistent with $n=100$ after rounding. With denominator degrees of freedom 97, the two versions give p-values approximately $0.4257$ and $0.4298$, respectively. Either way the partial effect is not significant at 5%.

The sample mean isWrite the Gaussian sample vector as its orthogonal projection onto the span of $(1,\dots ,1)$ plus its projection onto the orthogonal complement. These two Gaussian projections are independent. Their squared standardized lengths give Cochran's theorem:

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:

Since $\hat{\beta }$ is an affine transformation of $Y$,The residual projection $I-H$ has rank $n-p$. Cochran's theorem therefore gives