赖江山的博客分享 http://blog.sciencenet.cn/u/laijiangshan 生态、统计与R语言

博文

如何提取多元回归lm里面的全模型的p值

已有 571 次阅读 2021-1-25 23:02 |系统分类:科研笔记

刚有学员问用循环做了很多回归,每次拷贝回归模型的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应该是迎刃而解!




http://blog.sciencenet.cn/blog-267448-1268968.html

上一篇:为什么SEM中通径系数会大于1?
下一篇:通过置换检验(也就是 Monte Carlo模拟检验)获得相关系数的p值

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2021-4-15 02:29

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部