Я пытаюсь запустить двусторонний t-тест Уэлча на следующем примере фрейма данных:
structure(list(CellType = structure(c(1L, 2L, 3L, 3L, 1L, 1L, 2L, 3L, 3L, 1L), .Label = c("CellA", "CellB", "CellC"), class = "factor"), Gene1 = structure(c(2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L), .Label = c("Mut", "Wt"), class = "factor")), class = "data.frame", row.names = c(NA,-10L))
Я думаю, что код может выглядеть примерно так, поэтому я могу запустить этот t-тест между speci c Типы ячеек и Gene1 (имеющие только двоичные значения типа Mutant vs Wild):
t.test(Test$Gene1[Test$CellTypeT=="CellA"],
Test$Gene1[Test$CellType=="CellB"])$p.value
Однако этот код не работает как есть и может потребовать преобразования текстовых значений (например, «Mut»). "и" Wt ") в числовые значения c, которые могут быть распознаны функцией t-теста. Любые предложения о том, как я могу исправить эту ситуацию? Спасибо!