Запуск множественной логистической регрессии с одной и той же зависимой переменной в Rmarkdown - PullRequest
0 голосов
/ 13 мая 2019

Я создаю двоичные зависимые переменные, используя пороговое правило

# here a simple dataset such as 

x <- rnorm(n = 50, mean = 1, sd = 3)
names1 <- c(rep("gold", 10), rep("silver", 20), rep("bronze", 20))
df <- data.frame(names1, x)

thresholds <- c("L1", "L3")

logit_conditions <- c(function(x) ifelse(x >= 1, 1, 0), 
                      function(x) ifelse(x >= 3, 1, 0))

dat <- df[0]

for (i in 1:length(logit_conditions)) dat[i] <- logit_conditions[[i]](df$x)

colnames(dat) <- c("L1", "L3")
df <- data.frame(df, dat)

Затем я запускаю простую логистическую регрессию, которая отлично работает в консоли R, а также в Rmarkdown.

fit <- glm(L1 ~ x, data = df, family = "binomial")

Когда я пытаюсь создать много регрессий с одной и той же формулой и разными DV через lpply, в духе аналогичных сообщений SO ( Повторить регрессию с изменяющейся зависимой переменной , R Несколько циклов регрессии и коэффициенты извлечения ) приведенный ниже код запускается только из консоли R, в то время как зависимая переменная не найдена в Rmarkdown.Что я делаю неправильно?

model1 <- lapply(thresholds, function(dv) {
  fit <- glm(get(dv) ~ x, data = df, family = "binomial")
             })
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...