Я пытаюсь провести анализ логистической регрессии в R
. Я посещал курсы по этому материалу с использованием STATA. Мне очень трудно копировать функциональность в R
. Это зрелый в этой области? Там, кажется, мало документации или руководства доступны. Кажется, что для получения отношения шансов требуется установка epicalc
и / или epitools
и / или другие, ни один из которых я не могу приступить к работе, устарел или не имеет документации. Я привык glm
делать логистическую регрессию. Любые предложения будут приветствоваться.
Я бы лучше сделал это реальным вопросом. Как я могу провести логистическую регрессию и получить соотношение шансов R
?
Вот что я сделал для одномерного анализа:
x = glm(Outcome ~ Age, family=binomial(link="logit"))
И для многомерного:
y = glm(Outcome ~ Age + B + C, family=binomial(link="logit"))
Я посмотрел на x
, y
, summary(x)
и summary(y)
.
Имеет ли x$coefficients
какое-либо значение?
источник
cbind( exp(coef(x)), exp(summary(x)$coefficients[,1] - 1.96*summary(x)$coefficients[,2]), exp(summary(x)$coefficients[,1] + 1.96*summary(x)$coefficients[,2]) )
. Существует также дельта-метод: ats.ucla.edu/stat/r/faq/deltamethod.htmВы правы в том, что выходные данные R обычно содержат только существенную информацию, и больше нужно рассчитывать отдельно.
coefficients()
дает оценочные параметры регрессии . Хотя проще интерпретировать e x p ( b j ) (кроме перехвата).Чтобы получить отношение шансов, нам нужна таблица классификации исходного дихотомического DV и прогнозируемая классификация в соответствии с некоторым порогом вероятности, который должен быть выбран первым. Вы также можете увидеть функцию
ClassLog()
в пакетеQuantPsyc
(как упоминалось в смежном вопросе ).Для отношения шансов вы можете использовать пакет
vcd
или выполнить расчет вручную.источник
logit
иlogistic
получать коэффициенты шансов и доверительные интервалы легко. Я несколько разочарован тем, что это кажется настолько сложным и нестандартнымR
. Могу ли я использоватьexp(cbind(coef(x), confint(x)))
ответ Фабиана ниже, чтобы получить ОД и КИ? Мне не понятно, что дает ваш ответ?cTab
На странице статистики UCLA есть хороший обзор выполнения логистической регрессии в R. Она включает в себя краткий раздел по вычислению коэффициентов шансов.
источник
Пакет epiDisplay делает это очень легко.
источник
outreg
илиxtable
?