刚有学员问用循环做了很多回归，每次拷贝回归模型的R2和p值很麻烦，能否循环自动做成表格。我再次查看一下summary.lm函数，http://127.0.0.1:19929/library/stats/html/summary.lm.html

发现虽然直接显示是有全模型的p值，但输出的对象中无法直接提取模型p值（可能可以，但是我是不知道怎么提取，我说的不是系数p值，是全模型F检验的p值），而是给一个元素叫fstatistic，包括F值和两个自由度，其实有个这些值，就可以自己算p值。代码如下

fvalue <- summary(lm(y~x1+x2+x3+...xi,data))$fstatistic

pvalue <- 1-pf(fvalue[1],fvalue[2],fvalue[3])

有了求p值的方法，那么自己写循环提取p值和R2应该是迎刃而解！