Я создаю двоичные зависимые переменные, используя пороговое правило
# 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")
})