robbieyeah 发表于 2016-4-19 22:06 
给个例子
> cor(dt)
Y X1 X2 X3 X4 X5 X6 X7 X8 X9
Y 1.00000 0.91321 0.39624 0.62225 0.65082 0.32599 -0.29950 -0.23087 0.26687 -0.08998
X1 0.91321 1.00000 0.53889 0.59236 0.73131 0.38184 -0.15522 -0.35249 0.18238 -0.04305
X2 0.39624 0.53889 1.00000 0.32984 0.56939 0.20025 -0.09374 -0.43817 0.02971 0.12010
X3 0.62225 0.59236 0.32984 1.00000 0.52144 -0.03112 -0.47090 -0.08988 0.27215 0.06434
X4 0.65082 0.73131 0.56939 0.52144 1.00000 0.07544 -0.28853 -0.03899 0.31786 -0.26465
X5 0.32599 0.38184 0.20025 -0.03112 0.07544 1.00000 0.35186 -0.53474 -0.24747 0.05003
X6 -0.29950 -0.15522 -0.09374 -0.47090 -0.28853 0.35186 1.00000 -0.12701 -0.54097 -0.06658
X7 -0.23087 -0.35249 -0.43817 -0.08988 -0.03899 -0.53474 -0.12701 1.00000 0.37833 -0.62926
X8 0.26687 0.18238 0.02971 0.27215 0.31786 -0.24747 -0.54097 0.37833 1.00000 -0.34441
X9 -0.08998 -0.04305 0.12010 0.06434 -0.26465 0.05003 -0.06658 -0.62926 -0.34441 1.00000
> lmxy<-lm(Y~X1+X2+X3+X4+X5+X6+X7+X8+X9,data=dt)
> lmxy
Call:
lm(formula = Y ~ X1 + X2 + X3 + X4 + X5 + X6 + X7 + X8 + X9,
data = dt)
Coefficients:
(Intercept) X1 X2 X3 X4 X5 X6 X7
59.5867 0.5528 -0.1461 2.8870 -0.1616 0.0438 -0.5116 0.0593
X8 X9
-0.0537 -0.1095
> summary(lmxy)
Call:
lm(formula = Y ~ X1 + X2 + X3 + X4 + X5 + X6 + X7 + X8 + X9,
data = dt)
Residuals:
Min 1Q Median 3Q Max
-18.450 -2.642 -0.033 2.370 25.743
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 59.58669 58.97967 1.01 0.31266
X1 0.55281 0.01373 40.25 < 2e-16 ***
X2 -0.14609 0.03173 -4.60 4.8e-06 ***
X3 2.88702 0.82991 3.48 0.00053 ***
X4 -0.16160 0.03472 -4.66 3.8e-06 ***
X5 0.04384 0.00816 5.37 1.0e-07 ***
X6 -0.51165 0.04537 -11.28 < 2e-16 ***
X7 0.05933 0.01870 3.17 0.00157 **
X8 -0.05370 0.05757 -0.93 0.35127
X9 -0.10946 0.04652 -2.35 0.01886 *
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 4.67 on 811 degrees of freedom
Multiple R-squared: 0.884, Adjusted R-squared: 0.883
F-statistic: 689 on 9 and 811 DF, p-value: <2e-16
> vif(lmxy,digit=3)#vif(方差膨胀因子)的值大于10则存在共线性
X1 X2 X3 X4 X5 X6 X7 X8 X9
3.96 2.03 2.12 3.45 2.11 2.08 3.98 1.89 2.67