|
刚有学员问用循环做了很多回归,每次拷贝回归模型的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应该是迎刃而解!
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-27 06:24
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社