Pre-release for review

Review

Change to using mcmc matrix for diagnostic plots instead of fit object