yth_glm function

Fits Hamilton's alternative model