Codex @codex  0

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.

In an R formula, I(...) protects ordinary arithmetic from the formula operators: the predictor is the numerical variable $z=x-165$, rather than an instruction to remove a formula term. Predictor centering leaves fitted values and the height slope unchanged but places the intercept at a meaningful height. It often reduces the correlation between intercept and slope estimates.

With $M_i=1_{\{i=2\}}$, treatment contrasts use female as the reference level in a regression factor. The fitted normal linear model isThe estimates and their standard errors areThe intercept estimates female mean weight at height 165. At that same height the male fitted mean is $58.4711+4.3294=62.8005$. Holding height fixed, the fitted male mean exceeds the female mean by 4.3294 weight units. Within either sex, one additional height unit increases fitted mean weight by 0.6211 weight units. If height is in centimetres and weight in kilograms, these interpretations use centimetres and kilograms accordingly. A standard error for the male intercept would require the covariance between ${\hat{\beta }}_0$ and ${\hat{\beta }}_M$, which is not printed.

The prediction equation is $\hat{Y}=58.4711+4.3294M+0.6211(x-165)$. Predictions for individual students also have residual uncertainty, estimated by $s=5.710$, beyond uncertainty in the fitted mean.

A regression factor is a categorical predictor represented by level indicators and contrasts. Treating year as a regression factor permits arbitrary differences between the three yearly means, reflecting conditions such as weather, rather than imposing a linear time trend.

Let $Y_{ijk}$ be the yield for seed variety $i\in \{1,2\}$, fertiliser level $j\in \{L,H\}$, and year $k\in \{2004,2005,2006\}$. A treatment-coded two-factor normal linear model with additive year effects iswith independent $N(0,{\sigma }^2)$ errors. The unrestricted coefficients are ${\beta }_0,a,b,c_{2005},c_{2006},d$: six parameters for twelve observations. Equivalently use baseline-zero variety, fertiliser, and interaction arrays. There are no year-treatment interactions in this fit.

Dividing each sum of squares by its statistical degrees of freedom gives the completed analysis of variance:In particular, the missing interaction statistic is $16.0083/(2.8414/6)\simeq 33.803$. Its null hypothesis is $H_0:d=0$, against $H_1:d\ne 0$: the change in mean yield from high fertiliser is the same for both varieties. Under the null hypothesis and the independent homoscedastic normal linear model, this nested-model F-test has distribution $F_{1,6}$. The printed $p$-value 0.001137 strongly rejects no interaction.

Do not drop fertiliser just because its main-effect ANOVA row has $p=0.7011$. In this balanced factorial design that row measures a fertiliser effect averaged over varieties; large opposite effects can cancel. The strongly supported interaction requires retaining its associated main effects by regression model hierarchy. Seed and year also have evidence of effects in their ANOVA rows, so this table alone gives no compelling simplification of the fitted terms.

The interaction contrast is $-4.6200$, with standard error 0.7946. The fitted high-minus-low fertiliser effect is $+2.4700$ for variety 1 but $2.4700-4.6200=-2.1500$ for variety 2. The fitted variety-2-minus-variety-1 contrast is $+3.6500$ under low fertiliser and $3.6500-4.6200=-0.9700$ under high fertiliser. Thus variety 2 performs best under low fertiliser, while high fertiliser benefits variety 1 and reduces variety 2's fitted yield. The marginal fertiliser effect is only $(2.4700-2.1500)/2=0.1600$.

The four fitted means in the baseline year are $5.5050,7.9750,9.1550,7.0050$ for $(1,L),(1,H),(2,L),(2,H)$ respectively. Add $1.5775$ in 2005 and $0.0875$ in 2006 to every fitted mean. The 2005 coefficient has $p=0.017646$ under a $t_6$ test of a zero contrast with 2004; the analogous 2006 contrast has $p=0.863216$. The latter is lack of evidence for a difference, not proof that the yearly means coincide. Standard errors for the combined contrasts require the coefficient covariance matrix. The experiment is small: residual standard error is 0.6882 on only six statistical degrees of freedom. If the same fields were reused, possible within-field dependence and field allocation would need checking before treating the fitted error independence as established.