handle missing values
the downdate may be inaccurate if the cross product of design matrix is close 
to singular. Potentially check this doing estimation and re-do the QR 
decomposition if this is the case
