Блог

У меня есть простая логистическая регрессия, такая как

  logit(p[i]) lt;- b0   b1*x[i]   

но данные для p[i] даются как «да», «нет». Как кодировать («да»/»нет») в (1,0) ?

например, smt, как это

  logit(encode(p[i])) lt;- b0   b1*x[i]   

Попробуйте старомодный способ! Я предполагаю, что p и x-это столбцы данных, которые я объединяю в фрейм данных с именем df.

 library(tidyverse) library(stats) # for glm df lt;- as.dataframe(cbind(p,x)) df %gt;% mutate(p=ifelse(p=="yes",1,0) logit.model lt;- glm(p~x, data=df, family=binomial(link="logit")) summary(logit.model)