全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
22444 7
2014-09-01

比如: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值的函数?


二维码

扫码加我 拉你入群

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

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

全部回复
2014-9-1 18:28:00
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:37:40
熊小贤 发表于 2014-9-1 18:28
group
我指的不是简单回归的P值,简单回归方差检验P值等于系数的P值,所以可以用coefficient[2,4],但是多元回归不一样了。
二维码

扫码加我 拉你入群

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

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

2014-9-1 19:03:06
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 20:05:54
熊小贤 发表于 2014-9-1 19:03
不好意思,刚没注意看,这个可以自己算吧
a = summary(lm.D9)
1- pf(a$fstatistic[1], a$fstatistic[2] ...
对,就是这个意思,原来是pf这个函数,谢谢你!
二维码

扫码加我 拉你入群

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

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

2015-9-25 10:05:32
可以使用str(summary(fit))先查看内部的赋值情况,就能用summary(fit)$var来调用所有的变量了
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群