2015-02-11 39 views
2

比方说,我有一个数据帧提取p值

> col1<-c(1,5,2,6,8,1,3,8,9,1,8) 
> col2<-c(1,2,1,1,2,2,1,2,2,1,1) 
> df<-data.frame(col1,col2) 
> df 

    col1 col2 
1  1 1 
2  5 2 
3  2 1 
4  6 1 
5  8 2 
6  1 2 
7  3 1 
8  8 2 
9  9 2 
10 1 1 
11 8 1 

我已经跑了秩和检验测试的数据我在df

> dfKW<-kruskal.test(col1~col2, data=df) 
> dfKW 

Kruskal-Wallis rank sum test 

data: col1 by col2 
Kruskal-Wallis chi-squared = 1.695, df = 1, p-value = 0.1929 

我想什么要做的是将p值提取到一个向量中(只有没有标签'p值'的值)。我试过这个:

> dfKWx<-sapply(dfKW, '[', 'p.value') 
> dfKWx 

statistic.NA parameter.NA  p.value  method data.name 
      NA   NA   NA   NA   NA 

显然没有运气。

我从Google教授那里获得的初步帮助让我明白了上面所述的观点。

衷心感谢所有帮助。谢谢!

PS。请注意,这些数据仅用于举例说明。我还没有测试过正常分布,我也不想用这些数据来做。我的原始工作数据没有正常分发,但由于保密性因素,我不会将其用作示例 - 或其中的一部分。示例数据的行为与我原始的工作数据类似。

+0

'name(dfKW)'和'str(dfKW)'用于查看可以提取哪些对象 – user20650 2015-02-11 12:42:57

回答

2

这会给你的p值一个numeric

> dfKW$p.value 
[1] 0.1929473 

查找范围?kruskal.test “值” 部分。

+0

谢谢!我有一个臭名昭着的倾向,通过手册浏览。现在我希望我能再给一次这个时间。 – 2015-02-11 12:39:36