Confidence intervals with linear_contrast()

This release includes a new function linear_contrast(), which calculates robust confidence intervals and p-values for linear contrasts of regression coefficients from a fitted model. The function works with constrain_pairwise() and other constrain_*() helper functions.