我有一個矩陣(5,000 x 5,000),這是我的因變量隨着時間的推移,我有幾個自變量矩陣隨着時間的推移,在相同的格式。這兩個矩陣不時包含NA,所以這些必須通過na.exclude來照顧。 我做了一些樣本數據來說明我的問題: y <- matrix(rnorm(25000),5000,5000)
x1 <- matrix(rnorm(25000),5000,5000)
x2 <- matri
library(lmPerm)
x <- lmp(formula = a ~ b * c + d + e, data = df, perm = "Prob")
summary(x) # truncated output, I can see `NA` rows here!
#Coefficients: (1 not defined because of singularities)
#