Я хотел бы выполнить несколько bartlett.test между значениями в двух столбцах после идентификатора в другом столбце.
Данные
structure(list(moda = structure(c(3L, 3L, 3L, 3L, 3L, 3L, 3L,
3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L,
1L, 1L, 1L, 1L, 1L), .Label = c("BS1", "HW1", "PG"), class = "factor"),
t0 = c(1.5, 1.5, 2, 1, 1.5, 1.2, 1, 2.4, 1.3, 1.4, 1.7, 2,
1.8, 2.3, 2.5, 2.5, 1.5, 1.5, 2, 2.1, 1.8, 1.3, 2, 1.5, 2,
3.5, 1.5, 1.7), t14 = c(1.5, 1.5, 2, 1, 1.5, 1.2, 1, 2.4,
1.3, 1.4, 5, 7, 2.5, 6.5, 5.4, 5, 6, 5.7, 7, 3.5, 4.7, 4.5,
5.9, 5, 6, 7, 6, 5.5)), class = "data.frame", row.names = c(NA,
-28L))
Код
library(dplyr)
group <- function(data1, var1, var2){
var1 <- enquo(var1); var2 <- enquo(data1$t0), var3 <- enquo(data1$t14);
data1 %>%
group_by(!!data1$moda) %>%
summarise(n = n(),
bartlett.test(var2, var3))
}
group(df = data1, var1 = data1$moda, var2 = data1$t0, var3=data1$t14)
Если можно поместить результат в фрейм данных и сохранить идентификатор, он будет идеальным.РЕДАКТИРОВАТЬ: может быть, лучше использовать сплит раньше.И bartlett.test нужен один столбец со значением и один столбец с идентификатором, поэтому я не знаю, как это сделать в моем случае.