处理单变量拟合时，mclust-package中出现错误

Question

在使用Mclust进行单变量拟合时，我得到了以下错误：

Error in mstepE(data = as.matrix(data)[initialization$subset, ], z = z,  : 
  row dimension of z should equal data length

我正在使用中提到的代码：https://cran.r-project.org/web/packages/mclust/vignettes/mclust.html#initialisation

这是我遇到错误的代码部分：

df1 <- dataSample
BIC <- NULL
for(j in 1:20){
  rBIC <- mclustBIC(df1, verbose = T,
                    initialization = list(hcPairs = randomPairs(df1)))
  BIC <- mclustBICupdate(BIC, rBIC)
}
summary(BIC)

下面的链接包含要传递给变量'df1‘(文件name:dataSample.csv) https://drive.google.com/open?id=0Bzau9RsRnQreYk9XOWVBSm91b2o4NTQ4RlA2UFdWbDBVOVpR的数据

Answer 1

这是我从作者之一( Luca Scrucca教授)那里得到的'mclust‘库的解决方案：

“在提供hcPairs时，由于自动子集的使用会发生冲突，所以出现了一个错误。我已经在mclust的当前开发版本中修复了这个问题。由于提交给CRAN不会很快发生，您可以使用以下代码来避免当前版本的mclust出现错误：

rBIC <- mclustBIC(df1, verbose = T,
                  initialization = list(hcPairs = randomPairs(df1),
                                        subset = 1:NROW(df1)))

当bug修复程序发布时，可以省略subset参数，因为它是多余的。“

现在，代码运行得很好。