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 社区中对于分配首选= 或<- 的意见不同）。