p_mediation function

p-value from three-variable mediation analysis simulation