Estimation of RegArima models

We consider the linear model $y = X \beta + \varepsilon$ where $\varepsilon~N(0, \Omega)$

QR residuals