如何在R中处理crossprod中的缺失值？

Question

我想计算两个向量的总和乘积，但其中之一有NA。

crossprod(c(4.5,4,4,3,NA), c(1,5,4,1,1)) 

在计算，我想对待NA作为一个零。但在原始数据框架中，我仍然希望将其保留为NA以进行进一步的数据处理。任何暗示将不胜感激。

Answer 1

如果你的十字架产品比之更复杂，你可以定义一个新的功能，它取代NA的参数：

crossprod.replacena <- function(x, y, val=0) { 
            crossprod(replace(x, is.na(x), val), 
              replace(y, is.na(y), val) 
           ) 
         } 

或者你可以拨打你的论点replace没有定义的功能，如replace不修改其参数。

Answer 2

可以使用sum功能，并告诉它忽略遗漏值：

sum(c(4.5,4,4,3,NA)*c(1,5,4,1,1), na.rm = TRUE)