Trying to get marginal residuals but not working yet.
Added prediction function to predict fluoroescence at stations that weren't
fitted in the model to test extrapolation.