Function factor_analysis 

pub fn factor_analysis(
    trains: &[&[i32]],
    n_factors: usize,
    bin_size: usize,
    n_iter: usize,
) -> (Vec<f64>, Vec<f64>)

Factor analysis via EM. Rubin & Thayer 1982.