У меня вопрос по поводу нескольких t-тестов для одной групповой переменной.
У меня есть фрейм данных генной экспрессии, просто для примера предположим, что вот так:
df <- data.frame(a=runif(100), b=runif(100), c=runif(100)+0.5, d=runif(100)+0.5)
d <- melt(df)
Тогда у меня есть фрейм данных с переменной группировки. Как это:
dg <- data.frame(c(rep("A",50), rep("B",50)))
names(dg)<-"Group"
Я нашел несколько похожих вопросов, и мне удалось выполнить несколько t-тестов с lapply для каждого столбца первого фрейма данных и переменной группировки.
Однако я думаю, что необходимо отрегулировать значение p для множественных сравнений (у меня всего 15 столбцов). Вот моя проблема, как я могу это сделать? Например с Бонферрони? Могу ли я реализовать это в этом коде?
Кроме того, я думаю, что было бы хорошо проверить предположения для t-тестов также с этими несколькими столбцами данных (нормальное распределение / равные отклонения).
Есть здесь кто-нибудь, кто знает, как все это комбинировать? Идеальным будет сценарий, в котором статистика также преобразуется во фрейм данных или что-то вроде этого ...
Спасибо за любую помощь !!