Add function to compute and save ind-level joint likelihoods
Modify plotting range in plot_joint_ind