Как извлечь эту пару из pairwise.t.test с определенным значением?

Question

Я думаю, у меня есть несколько простой вопрос.

q1 <- with(airquality, pairwise.t.test(Ozone, Month))

и я хочу найти минимум этого, я знаю, что могу использовать

min(q1$p.value, na.rm = TRUE)

и скажите, что это 0.001. У меня вопрос, как мне найти пару, которая дала мне это конкретное значение из моего q1? Извините, если это немного расплывчато. Любая помощь будет оценена. Спасибо!

Answer 1

Будет проще, если вы предоставите пример, который мы можем использовать, поэтому вот пример из документации pairwise.t.test.

Я выполняю t-тест, а затем печатаю матрицу p-значения. Я нахожу индексы элемента, равные минимальному p-значению, и использую эти индексы, чтобы вытащить имена строк и столбцов.

# Air quality example
attach(airquality)
Month <- factor(Month, labels = month.abb[5:9])
res <- pairwise.t.test(Ozone, Month)

# Print resulting p-value matrix
print(res$p.value)
#>              May        Jun         Jul         Aug
#> Jun 1.0000000000         NA          NA          NA
#> Jul 0.0002638036 0.05112741          NA          NA
#> Aug 0.0001949061 0.04987333 1.000000000          NA
#> Sep 1.0000000000 1.00000000 0.004878798 0.003878108

# Find indices of minimum value
ind <- which(res$p.value == min(res$p.value, na.rm = TRUE), arr.ind = TRUE)

# Pull row and column names
rownames(res$p.value)[ind[1]]
#> [1] "Aug"
colnames(res$p.value)[ind[2]]
#> [1] "May"

^{Создано в 2019-03-15 пакетом Представления (v0.2.1)}