在计算点积时跳过NA

Question

我正在通过将它们的投影减去前1-2个主成分来调整数据矩阵中的测量值。问题是，如果数据矩阵中只有一个NA（几千次测量几乎不可避免），内部产品操作x%*%y（我也试过sum(x*y)，对于矢量{ {1}}）返回x,y。是否有一种简单的方法（即避免条件语句和循环）计算非NA值的内积，以便操作实际返回一些东西？

顺便说一句，我想避免只用NA替换NA，因为那时我必须在每个阶段重新规范化向量。

Answer 1

您可以尝试以下命令：

sum(x*y, na.rm = TRUE)