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

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

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]?
这是我试了两个小时试出来的。

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：Summary Summa Mary 数据提取 Mar library

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中第一个元素的第五个元素的第一个元素。（额，好拗口）

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

1. library(multcomp)
   
2. attach(cholesterol)
   
3. fit<-aov(response~trt)
   
4. a<-summary(fit)
   
5. edit(a);
   
6. #关掉弹出来的窗口
   
7. a[[1]];
   
8. #关掉弹出来的窗口
   
9. a[[1]][1,5];
   
10. #然后我想你应该已经懂了

复制代码

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

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

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

谢谢你了

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

不客气~