Как построить нечетное соотношение из множественной логистической регрессии с CI в r - PullRequest
Новая
       14

Как построить нечетное соотношение из множественной логистической регрессии с CI в r

0 голосов
/ 04 июля 2019

У меня есть модель с двоичным результатом 0,1 и четырьмя независимыми переменными, и все они являются категориальными переменными.Я хочу построить странное соотношение моего GLM.

My model is glm(outcome ~ v1 + v2 + v3, data = data, family = binomial)
All the variables v1, v2 and v3 have multiple categories.

Я также обнаружил, что некоторые люди предлагают использовать sjp.glm для построения отношения шансов для логистической регрессии.Но я не могу скачать этот пакет в версии 3.6.0.

Я нашел ниже код с этого сайта http://www.jscarlton.net/post/2015-10-24VisualizingLogistic/ Исход сюжета очень похож на то, что я хочу получить.Однако я не хочу делать это вручную.Есть ли способ построить график отношения шансов модели glm напрямую?

library(ggplot2)

# Create labels
boxLabels = c("Past encounter", "Mutualism WVO:\nNon-hunting wildlife experiences", "Domination WVO:\nEthics of hunting", "Domination WVO:\nWildlife management", "Age", "Education", "Negative attitudes toward box turtles", "Gender")

# Enter summary data. boxOdds are the odds ratios (calculated elsewhere), boxCILow is the lower bound of the CI, boxCIHigh is the upper bound.

df <- data.frame(
  yAxis = length(boxLabels):1,
  boxOdds = c(2.23189,1.315737,1.22866,.8197413,.9802449,.9786673,.6559005,.5929812),
  boxCILow = c(.7543566,1.016,.9674772,.6463458,.9643047,.864922,.4965308,.3572142),
  boxCIHigh = c(6.603418,1.703902,1.560353,1.039654,.9964486,1.107371,.8664225,.9843584)
)

# Plot
p <- ggplot(df, aes(x = boxOdds, y = yAxis))
p + geom_vline(aes(xintercept = 1), size = .25, linetype = "dashed") +
  geom_errorbarh(aes(xmax = boxCIHigh, xmin = boxCILow), size = .5, height = .2, color = "gray50") +
  geom_point(size = 3.5, color = "orange") +
  theme_bw() +
  theme(panel.grid.minor = element_blank()) +
  scale_y_continuous(breaks = yAxis, labels = boxLabels) +
  scale_x_continuous(breaks = seq(0,7,1) ) +
  coord_trans(x = "log10") +
  ylab("") +
  xlab("Odds ratio (log scale)") +
  annotate(geom = "text", y =1.1, x = 3.5, label ="Model p < 0.001\nPseudo R^2 = 0.10", size = 3.5, hjust = 0) + ggtitle("Intention to remove box turtles from the road")

пакет 'sjp.glm' недоступен (для версии R 3.6.0)

...