在計算點積時跳過NA's

Question

我正在通過將它們的投影減去前1-2個主成分來調整數據矩陣中的測量值。問題是，如果數據矩陣中甚至有一個單獨的NA（幾千次測量幾乎不可避免），內積操作x%*%y（我也試過sum(x*y)，對於矢量x,y）返回NA。有沒有一種簡單的方法（即避免條件語句和循環）在非NA值上計算內積，以便操作實際返回某些內容？順便說一句，我想避免用0來代替NA's，因爲那樣我就不得不在每個階段重新歸一化向量。

Answer 1

你可以試試這個命令：

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