stat_fit_augment function

Augment data with fitted values and statistics