Как перейти от OLS к логической регрессии?

Question

Я пытаюсь вычислить свою регрессию логита из OLS. Итак, на первом шаге я использовал функцию lm (), затем вычислил yhat и вероятность. Затем появляется сообщение об ошибке, что создаются NaN. Что я сделал не так?

start<-lm(invest~profit+gender,data=data)
summary(start)

data$yhat<-intercept+b1*data$profit+b2*data$gender

data$w<-sqrt(1/(data$yhat*(1-data$yhat)))

logit_2<-lm(w~profit+gender,data=data)
summary(logit_2)

Спасибо за вашу помощь!

Answer 1

Я чувствую, что так и должно быть, а не sqrt. Попробуй.

data$w<-ln(1/(data$yhat*(1-data$yhat)))

Примечание:

p = a0 + a1X1 + a2X2 + … + akXk    (linear)
ln[p/(1-p)] = b0 + b1X1 + b2X2 + … + bkXk       (logistic)