interaction_p_value function

Calculate interaction p-value