estimate_hat_vals function

Estimate hat (fitted) values