Why is the canonical link function used so frequently with GLMs? What makes it "natural"?
Is there any reason to think that, $Q(\theta _i)$ (where $Q$ is the canonical link function, and $\theta _i$ is the parameter of interest) is better described by a linear combination of predictor variables than some other function of $\theta_i$.
That is, is there any reason to believe that:
$$Q(\theta _i) = \sum_j \beta _j x_{ij}$$
is superior to:
$$f(\theta _i) = \sum_j \beta _j x_{ij}$$ where $Q \neq f$
If not, is there anything else that makes $Q$ better than $f$? The text book I am using just mentions what a canonical link function is, and makes use of it (pretty much exclusively), but does not explain what distinguishes it from any arbitrary link function.