Come già accennato, il teorema di Gauss-Markov dice che, considerando un modello lineare in notazione matriciale y=X?+? dove E[?]=0 e E[??']=(?^2) I, essendo ?°=[(XX')^-1]X'y uno stimatore qualsiasi ottenuto come combinazione lineare degli "y": b=Ly è tale per cui var (b)-var (?)=E (b-?)(b-?)'-E (?°-?)(?°-?)'. Quindi, il teorema può essere dimostrato prendendo in considerazione che uno stimatore lineare generico b=Ly in cui si decomponga la matrice "L" come L=D+[(XX')^-1]X'. Dunque, è necessario scrivere: E[b]=E{Dy+[(XX')^-1]X'y}=?. Tuttavia, questa dimostrazione non è sempre possibile, ma lo diventa solo se DX=0 e anche E[L?]=E[D?]=0. Per ottenere la matrice varianze - covarianze, bisogna scrivere: E (b-?)(b-?)'={D+[(XX')^-1]X'}??'{D+[(XX')^-1]X'}'=(?^2)(DD'[(XX')^-1]), perché la correttezza di "b" impone che DX=0.