Adding more robustness to heteroscedastic noise with factor analysis