我试图进行简单的矩阵乘法，之前从未给过我这个问题。
我有一个X矩阵（nxfe）和一个观察矩阵y（nx1）。出于本示例的目的，我们只说
X <- matrix(c(1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1), nrow=4)
y <- matrix(c(1, 2, 3, 4), nrow=4)
这种乘法的简单代码是
Xt <- t(X)
Xy <- Xt%*%y
然而，当我这样做时，我得到一个NA＆＃39s（nx1）的矩阵。
这从未发生过。我已经检查过并且在我的所有矩阵中都有适当的值。任何帮助，将不胜感激。