The Gauss-Markov theorem states that with full-column-rank design, mean-zero errors and covariance , the ordinary least squares estimator is the best linear unbiased estimator: its covariance is no larger than any other such estimator's in the Loewner order. Normality is unnecessary for that theorem.
In the given normal linear model, maximizing the likelihood function is equivalent to minimizing . The normal equations and full rank give
For a direct variance comparison, write . Unbiasedness for every requires . Set and , so . The cross covariance terms vanish because . Therefore
This proves the requested comparison for every contrast without relying only on the theorem's name.
For this model, the Gauss-Markov theorem states that is the unique minimum-variance linear unbiased estimator of : if
is unbiased for every , then .
Indeed,
so unbiasedness is equivalent to . Write
The constraint becomes . Since the errors are independent with common variance ,
Equality holds exactly when every , which gives and .
Write , where
Any other linear unbiased estimator is with , so and . The cross covariance vanishes because
Therefore
By the Gauss-Markov theorem,