I want to fit a linear model by R with family=binomial(link="identity")
, however, binomial family do not have identity link. What should I do?
Solved – R binomial family with identity link
binomial distributiongeneralized linear modellink-functionr
Best Answer
See Wikipedia on the linear probability model, & CV posts here & here for the statistical background. Though not "wrong", you'd want a good reason for using an identity link to model a Bernoulli probability.
According to the
family
manualBut
So
family=binomial(link="identity")
works butfamily=binomial(link=identity)
doesn't. (If you find differently it might be to do with the R version.) To allow for over-dispersion, then usefamily=quasi(link="identity", variance = "mu(1-mu)")
.