added code for using refund package functions

Added code from another old simulation that fits mediation with
different functions from the refund package