Я сгруппировал данные, в которых я выполняю тест хи-квадрат, и хотел бы вернуть сводную таблицу, содержащую несколько значений из объекта htest
.Например ( из предыдущего вопроса ),
library(dplyr)
set.seed(1)
foo <- data.frame(
partido=sample(c("PRI", "PAN"), 100, 0.6),
genero=sample(c("H", "M"), 100, 0.7),
GM=sample(c("Bajo", "Muy bajo"), 100, 0.8)
)
foo %>%
group_by(GM) %>%
summarise(p.value=chisq.test(partido, genero)$p.value))
возвращает значение p, но вместо этого я хотел бы получить несколько значений (скажем, p.value
и statistic
) из htest
объект, который будет возвращен в виде различных столбцов в сводной таблице.
Я пытался
foo %>%
group_by(GM) %>%
summarise(htest=chisq.test(partido, genero)) %>%
mutate(p.value=htest$p.value, statistic=htest$statistic)
, но это выдает ошибку
Ошибка в summarise_impl(.data, точки):
Столбец htest
должен иметь длину 1 (итоговое значение), а не 9
Как этого добиться с помощью инструментов Tidyverse?