In regard to the code for the plots, the function is:
mv_mean_contour(self, mu1_l, mu1_u, mu2_l, mu2_u, step1, step2,
levs=[.2, .1, .05, .01, .001], plot_dta=False)
The user enters the mean ranges of both variables and the interval in which to conduct a joint hypothesis test of the means. Before, these had to be entered by the user because there were problems with hypothesis tests with very large likelihood ratios. however, I think I fixed this problem and can now make all of the input parameters optional. It will take some fiddling first.