Enhancing shiftplot - can now use colour inside aes() and can pass a theme into the function. Aslo adding margarita.simulate.simple to simulate similar data to the original observed data using all components of the fitted margarita model