函数kappa()
df<-data.frame()df_cor=cor(df)kappa(df_cor, exact=T)
当 κ<100 , 说明共线性程度小;
当 100<κ<1000 , 有较强的多重共线性;
当 κ>1000 ,存在严重的多重共线性。
函数qr()
x<-matrix()qr(x)$rank
qr(X)$rank 计算X矩阵的秩,如果不是满秩的,说明其中有 xi 可以用其他 xj 的线性组合表示;此时可以进行逐步回归,用step()命令。
fm<-lm()step(fm)