Я пытаюсь рассчитать коэффициенты шансов и связанные значения p для набора биомаркеров, используя порядковую логистическую регрессию (в частности, polr из пакета MASS). Я использовал функцию ordinal.or.display () из пакета epiDisplay, чтобы просмотреть результаты регрессии, но заметил несоответствие между отображаемым значением p и тем, что я вычисляю вручную ... это примерно вдвое больше большой, когда я рассчитываю это с помощью нормального распределения. Я что-то упускаю из-за Порядковой Логистической Регрессии, или это проблема с функцией epiDisplay?
Я попытался просмотреть документацию для пакета epiDisplay (https://cran.r -project.org / web / packages / epiDisplay / epiDisplay.pdf ), но не нашел ничего, что объясняло бы, как -значение там было рассчитано. Любая помощь или дополнительные знания очень ценятся!
#the model using polr from MASS
#generating artifical outcome var from mt cars
mtcars <- mtcars
mtcars$outcome <- round(runif(nrow(mtcars))*5)
myMod <-polr(ordered(outcome) ~ factor(am)+ factor(carb)+ wt,
data = mtcars,
Hess = TRUE)
summary <- summary(myMod)
(ctable <- coef(summary(myMod)))
p <- pnorm(abs(ctable[, "t value"]), lower.tail = FALSE) * 2
## combined table: p = 0.624
(ctable <- cbind(ctable, "p value" = p))
ordinal.or.display(myMod) #p value = 0.312
Я ожидаю, что значения p будут одинаковыми - возможно, epiDisplay не умножится на два?