假设我有一个数据矩阵 d
pc = prcomp(d)
# pc1 and pc2 are the principal components
pc1 = pc$rotation[,1]
pc2 = pc$rotation[,2]
那么这应该适合线性回归模型吧?
r = lm(y ~ pc1+pc2)
但后来我得到这个错误:
Errormodel.frame.default(formula = y ~ pc1+pc2, drop.unused.levels = TRUE) :
unequal dimensions('pc1')
我想那里有一个包可以自动执行此操作,但这也应该有效吗?