Я определил, что наиболее применимым статистическим анализом после разборки после одностороннего ANOVA моих данных является тест Даннетта.Я выполнял это с помощью R в прошлом, однако теперь я ограничен Python из-за пакетов, которые я использовал для автоматизации моего рабочего процесса (автоматический анализ больших объемов данных).
Я нашел несколько пакетов (например,sci-kit, tukeyHSD), которые предоставляют несколько различных специальных тестов, однако ни один из них не включает тесты Даннетта.Например, я могу легко выполнить одностороннюю ANOVA в scipy:
import scipy.stats as stats
# Made up data
a = [10, 12, 10, 14, 18] # Control
b = [15, 14, 18, 10, 38]
c = [20, 22, 23, 34, 20]
d = [50, 48, 42, 51, 51]
stats.f_oneway(a, b, c, d)
> F_onewayResult(statistic=26.92639734366354, pvalue=1.7207487532445122e-06)
Однако я хочу провести многократный сравнительный анализ после этого по сравнению с одной контрольной группой с нормально распределенными данными (n ~ 1000).Я знаю о Rpy2, однако я хочу выполнить это на машине без использования докера.Любые предложения по способным пакетам?
(Я также биолог с довольно базовыми знаниями в области сценариев, поэтому есть большая вероятность, что мне здесь не хватает чего-то простого)