Может ли кто-нибудь помочь мне понять, как glm обрабатывает бесконечность?
(log) отношение шансов колеблется от -Inf до + Inf.
log(((0)/(4))/((4-0)/(4)))
#-Inf
log(((2)/(4))/((4-2)/(4)))
#0
log(((4)/(4))/((4-4)/(4)))
#Inf
Когда вы подходите к модели, как показано ниже в glm перехватчик показывает некоторые значения (например, 23,56606889) вместо Inf
d<-data.frame(ar=c(1,1,1,1))
m<-glm(ar~1,d,family=binomial)
coef(m)[1]
#(Intercept)
#23.56606889
Это означает, что glm принимает приближенное значение: например,
log(((3.999999999)/(4))/((4-3.999999999)/(4)))
#22.10956012
log(((3.9999999999)/(4))/((4-3.9999999999)/(4)))
#24.41214521
Я хотел бы знать как glm принимает приближенные значения.