finished debugging MCMC pair, need to speed up a little if possible
pairwise warping finished and debugged, need plotting functions