关于summary(aov(...))中的数据提取问题

lww1993

10041

收藏 2013-10-24

library(multcomp)
attach(cholesterol)
fit<-aov(response~trt)
a<-summary(fit)
detach(cholesterol)
a
Df Sum Sq Mean Sq F value Pr(>F)
trt       4 1351.4 337.8 32.43 9.82e-13 ***
Residuals 45  468.8 10.4
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1''1
我们要提取Pr(>F)下的9.82e-13,
为什么提取方式是a[[1]]["Pr(>F)"][[1]][1]?
这是我试了两个小时试出来的。

扫码加我拉你入群

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

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

全部回复

namgalsip

2013-10-24 13:08:48

a为list格式，a[[1]]表示引用该列表的第一个元素，
定义 b <- a[[1]]["Pr(>F)"]，它表示引用a第一个元素的"Pr(>F)"列，
b还是list格式，所以使用b[[1]]引用b中第一个元素，b[[1]]就是一个向量了，
b[[1]][1]则为选其中的第一个元素，即trt的 Pr(>F) 值。

应该和a[[1]][[5]][1]效果相同，a中第一个元素的第五个元素的第一个元素。（额，好拗口）

扫码加我拉你入群

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

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

求证1加1

2013-10-24 14:11:37

这其实就是一个嵌套list，lz用一下edit(a)看一下他的结构就明白了，第一步用a[[1]]的原因其实也就是把外面套的这个list给去掉，虽然根据输出来看貌似没啥区别，但你这时在用edit()时候你就知道区别了。附上代码，lz可以照着试试

复制代码

扫码加我拉你入群

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

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

lww1993

2013-10-24 20:58:50

求证1加1 发表于 2013-10-24 14:11
这其实就是一个嵌套list，lz用一下edit(a)看一下他的结构就明白了，第一步用a[[1]]的原因其实也就是把外面套 ...

我明白了，谢谢了，edit真好用

扫码加我拉你入群

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

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

lww1993

2013-10-24 20:59:55

namgalsip 发表于 2013-10-24 13:08
a为list格式，a[[1]]表示引用该列表的第一个元素，
定义 b F) 值。

谢谢你了

扫码加我拉你入群

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

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

求证1加1

2013-10-24 21:01:33

lww1993 发表于 2013-10-24 20:58
我明白了，谢谢了，edit真好用

不客气~

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群