Встреча с "Ошибка в eval(predvars, data, env): объект" Возраст "не найден" при использовании функции прогнозирования в логистической регрессии "

Question

#Logistic Regression
glm.fit <- glm(recent_cannabis_use~.,data = drug_use_train, family = binomial)
summary(glm.fit)


predict(glm.fit, with(drug_use_train, data.frame(Gender = "Male")), type = "response")

Попытка найти прогнозируемую вероятность недавнего_каньяка для мужчины.

Answer 1

Вы должны использовать predict(glm.fit, newdata = data.frame(Gender = "Male")), С помощью with в этом случае не гарантируется, так как вы не обращаетесь ни к одной из переменных в drug_use_train,

Обратите внимание, что это предполагает, что ваша формула, после расширения, recent_cannabis_use ~ Gender, Если у вас есть другая переменная, и вы хотите изучить только эффект Gender, вам нужно будет установить (предварительно рассчитать или составить) все другие переменные на некоторое фиксированное значение (помните, как интерпретируются коэффициенты - меняйте y с одной единицей изменения x, если все остальное остается прежним). Смотрите, например, этот пост.