Ordinary differential equations (Equation 1–5) were solved using the stiff system solver ode15s in Matlab (MathWorks, Natick, MA). The standard errors in the model parameters were determined using Markov chain Monte Carlo (Gilks et al., 1996). The Markov chain Monte Carlo procedure was also designed to ensure convergence to the global minimum and ensure that the model parameters were optimal.