= Solution
Let
$$
H=X(X^TX)^{-1}X^T
$$
be the <hat matrix>. Then
$$
\widehat\beta=(X^TX)^{-1}X^T(HY)
$$
depends only on the orthogonal projection $HY$, while
$$
n\widehat\sigma^2=\|(I-H)Y\|^2
$$
depends only on the orthogonal residual projection $(I-H)Y$.
The random vectors $HY$ and $(I-H)Y$ are jointly <multivariate normal distribution>[multivariate normal], and their cross-covariance is
$$
\sigma^2H(I-H)=0.
$$
Uncorrelated jointly Gaussian vectors are <independent random variables>[independent]. Therefore
$$
\boxed{\widehat\beta\ \text{and}\ \widehat\sigma^2
\text{ are independent}}.
$$
Back to article page