Я хотел бы выполнить bartlett.test между двумя столбцами значений, следуя специальному коду внутри другого столбца, и я хочу применить эту функцию ко всем моим группам, запустить ее в Rmarkdown и получить все мои результаты в данные .Рамка. Итак, вот краткий пример моих данных (в моем первом столбце больше групп).
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))
result_list <- by(data1, data1$moda, function(sub_data) {
res <- t.test(sub_data$t0, sub_data$t14, mu=0, alt="two.sided",
paired = TRUE, conf.level = 0.95)
output <- c(t_stat=unname(res$estimate), p_value=res$p.value)
})
result_matrix <- do.call(rbind, result_list)
result_matrix
Кто-то дал способ выполнить t.test, но я бы знал, что выполнить bartlett.test. И, может быть, wilcox.test я понимаю, как это работает для подмножества.
t.test между двумя столбцами после определенного кода в третьем столбце