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)