ENH: Fix numerical issues in links and families.

Closes #2024.
Should this use `t = np.exp(-z) and `return 1 / (1 + t)**2` ?