Of course, there exists statistics library in Julia which do the job. But which library is the â€śbestâ€ť with respect to stability, CPU time, etc ? Especially for these functions :

RGIVN fits a multivariate linear regression model via fast Givens transformations.

RLEQU fits a multivariate linear regression model with linear equality restrictions imposed on the regression parameter.

RSTAT computes statistics related to a regression fit given the coefficient estimates.

(WRRL is a basic output for matrix, forget this function)