Algorithms for iteratively refining the parameter estimates and residuals from the fitting of a regression model using QR decomposition are described. It is shown that, if algorithms for performing the QR decomposition are free from calculation of square root, the related iterative refinement algorithms can also be square root free. Verification of the algorithms is carried out and comments are made on the accuracy of parameter estimates.