Извлечь 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-value). Я пробовал это:

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

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

Очевидно, без везения.

Первоначальная помощь, которую я получил от профессора Google, привела меня к тому, что я описал выше.

Искренне признателен за любую помощь. Спасибо!

PS. Обратите внимание, что данные приведены только для примера. Я не тестировал нормальное распределение и не хочу делать это с этими данными. Мои исходные рабочие данные не распространяются нормально, но я не буду использовать их или их части в качестве примера из-за факторов конфиденциальности. Данные примера ведут себя так же, как мои исходные рабочие данные.


person Olli J    schedule 11.02.2015    source источник
comment
names(dfKW) и str(dfKW) полезны, чтобы увидеть, какие объекты можно извлечь   -  person user20650    schedule 11.02.2015


Ответы (1)


Это даст вам значение p как numeric:

> dfKW$p.value
[1] 0.1929473

Посмотрите на раздел «Значение» в ?kruskal.test.

person Stephan Kolassa    schedule 11.02.2015
comment
Спасибо! У меня есть печально известная склонность пролистывать руководства. Теперь я хотел бы дать это еще раз. - person Olli J; 11.02.2015