Я должен выполнить линейную регрессию (лм) для двух модифицированных типов p-значения: q-значения и p-значения, скорректированного по Бенджамини-Хохбергу. Результаты дают два астрономических выброса, однако после их удаления всегда присутствуют новые выбросы. Вот полный код:
p = 5000
m = 1000
pval = c(rbeta(m,1,100), runif(p-m,0,1))
BHpval <- p.adjust(pval,method="BH")
qval_ <- qvalue(pval)
print(qval_$pi0)
fit2 <- lm(qval_$qvalues ~ BHpval)
plot(fit2)
Их невозможно удалить, поэтому я подозреваю, что в моем понимании есть пробел в линейной регрессии или в этих двух модификациях p-значения. Может кто-нибудь объяснить, что я делаю не так?