em算法的r语言实现