Как извлечь `Pr (> F)` из `leveneTest ()`?

Question

> leveneTest(Sepal.Length~Species,iris,na.action='na.omit')
Levene's Test for Homogeneity of Variance (center = median: "na.omit")
       Df F value   Pr(>F)   
group   2  6.3527 0.002259 **
      147                    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

После leveneTest() я проверил структуру результата следующим образом:

> leveneTest(Sepal.Length~Species,iris,na.action='na.omit') %>% str
Classes ‘anova’ and 'data.frame':   2 obs. of  3 variables:
 $ Df     : int  2 147
 $ F value: num  6.35 NA
 $ Pr(>F) : num  0.00226 NA
 - attr(*, "heading")= chr "Levene's Test for Homogeneity of Variance (center = median: \"na.omit\")"

Затем я извлек Pr(>F), но получил ошибку, как показано ниже:

> leveneTest(Sepal.Length~Species,iris,na.action='na.omit')$Pr(>F)
Error: unexpected '>' in "leveneTest(Sepal.Length~Species,iris,na.action='na.omit')$Pr(>"

В чем проблема? Как извлечь Pr(>F)?

Answer 1

Обычный x$Pr(>F) не будет работать, потому что > - это специальный символ, поэтому вам нужно x[["Pr(>F)"]] или x$`Pr(>F)`

Обратите внимание, что у вектора, который вы получаете, есть значение, за которым вы хотите следовать на NA (то есть c(0.002258528,NA)