2015-02-11 36 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