一共16个因素x1~x16，分别对y进行一元回归，r1=lm(y~x1);r2=lm(y~x2);……r16=lm(y~x16);
然后summary(r1);summary(r2);……summary(r16);
r程序运行结果是16个格式如下的结果
> summary(r13);
Call:
lm(formula = y ~ x13)
Residuals:
    Min      1Q  Median      3Q     Max
-42.416 -10.935  -4.618  12.309  82.417
Coefficients:
            Estimate Std. Error t value Pr(>|t|)  
(Intercept)   25.163     14.010   1.796   0.0862 .
x13           -4.780      3.202  -1.493   0.1497  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 28.65 on 22 degrees of freedom
Multiple R-squared: 0.09197,    Adjusted R-squared: 0.0507
F-statistic: 2.228 on 1 and 22 DF,  p-value: 0.1497

目标就是上面蓝色的数字P值0.1497
请问有没有一段程序能够直接得到16个P值，而不用一个个复制粘贴到excel中得到？
谢谢！