joint_trait_prob_gradient function

Computes partial response curves of joint probabilities