怎样提取一个多元回归的summary中F检验的p值？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

比如：summary(lm1)得到如下结果

Residuals:
Min 1Q Median 3Q Max
-2.87323 -0.64441 -0.00852 0.63180 2.77441

Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.5452 0.2699 2.020 0.0453 *
geno[, i] -0.2860 0.3602 -0.794 0.4286
geno[, j] 0.1023 0.4120 0.248 0.8043
geno[, k] -0.2946 0.1600 -1.841 0.0677 .
geno[, m] -0.2187 0.2814 -0.777 0.4383
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.9885 on 138 degrees of freedom
Multiple R-squared: 0.0348,        Adjusted R-squared: 0.00682
F-statistic: 1.244 on 4 and 138 DF, p-value: 0.2953

怎样获取最后一个p-value ？
或者有没有计算F-statistic P值的函数？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

关键词：Summary Summa 多元回归 Mary Mar standard error Error

group <- gl(2, 10, 20, labels = c("Ctl","Trt"))
weight <- c(ctl, trt)
lm.D9 <- lm(weight ~ group)
summary(lm.D9)$coefficients[2,4]

熊小贤 发表于 2014-9-1 18:28
group

我指的不是简单回归的P值，简单回归方差检验P值等于系数的P值，所以可以用coefficient[2,4]，但是多元回归不一样了。

dearxxj99 发表于 2014-9-1 18:37
我指的不是简单回归的P值，简单回归方差检验P值等于系数的P值，所以可以用coefficient[2,4]，但是多元回归 ...

不好意思，刚没注意看，这个可以自己算吧
a = summary(lm.D9)
1- pf(a$fstatistic[1], a$fstatistic[2], a$fstatistic[3])

熊小贤 发表于 2014-9-1 19:03
不好意思，刚没注意看，这个可以自己算吧
a = summary(lm.D9)
1- pf(a$fstatistic[1], a$fstatistic[2] ...

对，就是这个意思，原来是pf这个函数，谢谢你！

可以使用str(summary(fit))先查看内部的赋值情况，就能用summary(fit)$var来调用所有的变量了

dearxxj99 发表于 2014-9-1 20:05
对，就是这个意思，原来是pf这个函数，谢谢你！

具体怎么做的，能写下吗，我按照做没做出来，谢谢

抛去计算机算的不看，你自己用笔算时，当F用MSreg/MSres计算出后代入F分布表就可以得出probablity了。