risk_mod function

Fit an Integer Risk Score Model