Displaying 1 result from an estimated 1 matches for "richinit".
Did you mean:
mic_init
2006 May 24
1
problem-nlme
...onBertInit)
class(vonBert)
#Richards
Rich <- function(x, Linf, K, t0, m)
Linf*(1-exp(-K*(x-t0)))^(1/(1-m))
size ~ Rich(age, Linf, K, t0, m)
Rich <-deriv(~ Linf*(1-exp(-K*(x-t0)))^(1/(1-m)),
c("Linf","K","t0","m"),function(x,Linf,K,t0,m){})
RichInit <- function(mCall, LHS, data)
{
xy <- sortedXyData(mCall[["x"]], LHS, data)
Linf <- 900
if (Linf != max(xy[,"y"])) Linf <- -Linf
K <- 0.3
t0<-0
m <- 0.3
value <- c(Linf, K, t0, m)
names(value) <- mCall[c("Linf",...