r - 来自 R 博客网页的 EM 算法中的 R 代码错误

Question

我参考了这个页面的 EM 代码： http ://www.r-bloggers.com/imputing-missing-data-with-expectation-maximization/

最初，当我运行代码时raw < - replicate(10, rpois(50,100))，我收到了这个错误：

Error in raw < -replicate(10, rpois(50, 100)) : 
  comparison (3) is possible only for atomic and list types

然后我将该随机数据更改为raw<-faithful$waiting. 但随后，我收到另一个错误：

Error in solve.default(sig[!pick.miss, !pick.miss]) : 'a' is 0-diml

这里有人熟悉这段代码以及如何解决这些错误吗？

Answer 1

复制代码时要小心，尤其是从 R-Blogger 等低质量网站复制代码时：R 中的赋值是由 完成的<-，而不是由完成的< -。两个字符之间的空格是错误的，你需要删除它。

通常，您也可以替换<-为=以避免此问题（R 社区内部对于分配是否=或<-首选分配的意见不同）。