Commit e35f64af authored by patrick.lambert's avatar patrick.lambert
Browse files

to allow multiple calculations

parent 23842291
......@@ -61,7 +61,7 @@ vonBertalanffyWithNextIncrement = function(L, L0, Linf, K, timeStepDuration, sig
}
else {
mu = log((Linf - L) * (1 - exp(-K * timeStepDuration))) - (sigma * sigma) / 2
increment = exp(rnorm(1, mu, sigma))
increment = exp(rnorm(length(mu), mu, sigma))
}
L = L + increment * tempEffect
return(L)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment