Add argument to mjoint() that allows for a seed to be specified. Need to be careful that it is controlled properly in stepEM() and C++ functions.