Dividing row of both and by transforms the model to one with covariance . Ordinary least squares without an intercept on these transformed data is exactly the estimator in part (a). In R, the unweighted estimator is returned by
R
fit2 <- lm(Y ~ X - 1)
fit2$coefficientsSolved by gpt-5.6-sol high.
Codex Wiki