我有一个协变量矩阵(n = 759,p = 175),有大量缺失的数据,我想看看样本中缺失模式是否有一些相似之处。
由于我不习惯使用二进制数据,因此我阅读了一些方法并遇到了非常有趣的homals
R包。
从is.na(data_covariates)
我删除了非信息性特征(没有样本是一个完整的案例,所以我留下所有样本)并留下了759 x 164矩阵data_covariates_missing
,我作为参数传递到homals()。
(编辑:我还删除了两个重复的行)
命令homals(data=data_covariates_missing)
导致错误
Error in qr.default((1/sqrt(w)) * x) :
NA/NaN/Inf in foreign function call (arg 1)
和许多警告,都有相同的信息:
In is.na(x) : is.na() applied to non-(list or vector) of type 'NULL'
我在ndim, rank
参数中玩level
和homals()
,但没有运气。检查homals
功能代码或插图没有帮助,也没有搜索StackOverflow。
通过构造我的矩阵没有遗漏或无限值,所以我真的无法解释错误和警告。
有什么建议吗?谢谢!