Code to (i) get the MLE of the cupid effect, (ii) impute and rerandomize the data and (iii) plot the pvalues