add_index function

Adds an indirect effect index to a fitted mediation model