我有一个协变量矩阵(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
,但没有运气。检查函数代码或小插图没有帮助,搜索 StackOverflow 也没有帮助。level
homals()
homals
通过构造,我的矩阵没有缺失值或无限值,所以我真的无法解释错误和警告。
有什么建议吗?谢谢!