Watch this Author's files
Nicely done. It is useful to have a time rounding fn to determine if times are equal within a tolerance, especially if they may have lost significant digits.
Nice work. Much faster than factoran. I got another speed boost by calling dsyevr from Lapack instead of svd.
Nicely done; simple and efficient code.
This is an interesting form of robust regression. I checked the likelihood function, and it looks correct.