HEAVYmodel function

HEAVY model estimation