Add EM.q_init() and EM.t_q_estimate() methods. This works but uses a ton of memory.