Normal linear model

$Y=X\beta +\epsilon$, $\epsilon \sim N(0,{\sigma }^{2}I)$, and full-rank $X$ gives $\hat{\beta }=(X^{T}X{)}^{-1}{X}^{T}Y$.