r_score_IM function

Simulation from the interaction model