I'm getting inconsistent output about the link function from the
summary() command when fitting an inverse Gaussian GLM:
> summary(glm(ig.formula, family=inverse.gaussian(link = "log"),
+ data=mydata, start=start.vals))$call
glm(formula = ig.formula, family = inverse.gaussian(link = "log"),
data = mydata, start = start.vals)
> summary(glm(ig.formula, family=inverse.gaussian(link = "log"),
+ data=mydata, start=start.vals))$family
Family: inverse.gaussian
Link function: 1/mu^2
Has anyone else run into this problem?
I'm running v1.9.1 in Windows.
E.J.