交叉产品矩阵意外地充满了NA

Question

所以我从CSV中获取了一些信息,将其存储为矩阵,并试图在结果上计算以下操作,但它给了我一个2x2的NA数组.没有在这里看到问题.

data <- read.csv('qog.csv', sep=';')
X <- matrix( log( data$wdi_gnipc ) )
X <- cbind(X, data$ciri_empinx_new)
t(X) %*% X

当我看到X和t(X)时,它们看起来就像我期望的那样,所以我将2xn矩阵与nx2矩阵(n是像193这样的大数字)进行矩阵乘法,因此矩阵乘积应该很好 - 定义并给出有意义的2x2答案.

什么想法可能会出错？

注意:我试试的时候

a < - rbind(c(1,2),c(3,4))t(a)%*%a

它给出了期望的结果.不确定它与我正在对数据做什么之间的重要区别.

Answer 1

让我们回答一下吧.对于要填充的叉积,每列内部NA必须至少有一NA列X.您可以NA通过运行找到每列的s数:

colSums(is.na(X))

并且很有可能你会拥有它

all(colSums(is.na(X)) > 0)
# [1] TRUE