Correct parameters to gamma distribution
Sample from prediction distribution instead of using mean