R/funcion_vb.R

In PabloMBooster/CAMPtools:

# funcion para calcular los parametros de crecimiento

ab_fun <- function(data, Peso_medio, t){

  model   = nls(Peso_medio ~ Winf*(1-exp(-k*(t)))^3, data = datos, control = list(maxiter = 1000),
                start = c(Winf = 227, k = 0.02))
  Winf <- as.numeric(coefficients(model)[1])
  k    <- as.numeric(coefficients(model)[2])
  vb_list = list()
  vb_list$k = k
  vb_list$Winf = Winf

  return(vb_list)
}