Persamaan model regresi yang sering kita jumpai adalah sebagai berikut:
y=Xβ+ε (1)
dimana y adalah vektor T x 1, X adalah matrix T x p, β adalah vektor p x 1, dan ε adalah vektor T x 1. Sistem OLS estimator-nya didapat dengan minimisasi error sum of squares, ESS=ε’ε sebagai fungsi β yang tidak diketahui:
ESS = (y-Xβ)^’ (y-Xβ)=y^’ y-2y^’ Xβ+β’X’Xβ (2)
Kalau dibuat kodenya, kira-kira seperti ini: