У меня есть сомнение: я корректирую кратность с помощью BH (у меня довольно много сравнений и я немного более мягок), но когда я корректирую значения, результаты мне кажутся странными. Код для тестирования (и получение значений в виде вектора, чтобы я мог затем использовать его для корректировки)
prandom <- data.frame(p.value= sapply(random[,5:77], function(i)
t.test(i ~ random$Group)$p.value))
print(prandom)
which(prandom$p.value < 0.05)
x <- (prandom$p.value)
corrected <-p.adjust(x, method = "BH")
результат x равен:
[1] 0.007319452 0.044712307 0.059936099 0.081492910 0.085423931 0.108844961 0.134041535 0.216181963 0.225669850
[10] 0.229112552 0.264940706 0.266180569 0.316797423 0.321418692 0.338581157 0.358322392 0.361364009 0.381420023
Результат исправления равен :
[1] 0.5343200 0.9737560 0.9737560 0.9737560 0.9737560 0.9737560 0.9737560 0.9737560 0.9737560 0.9737560 0.9737560
[12] 0.9737560 0.9737560 0.9737560 0.9737560 0.9737560 0.9737560 0.9737560 0.9737560 0.9737560 0.9737560 0.9737560
У меня такое чувство, что я делаю что-то не так Большое спасибо