R/fitplotsize.R

In biotools: Tools for Biometry and Applied Statistics in Agricultural Science

Defines functions fitplotsize

Documented in fitplotsize

fitplotsize <-
function(plotsize, CV)
{
   if (length(plotsize) != length(CV))
      stop ("incompatible dimensions!")

   ini <- lm(log(CV) ~ log(plotsize))$coef

   fit <- nls(CV ~ a * plotsize^(-b),
      start = list(a = ini[1], b = ini[2]))
   fit
}

Try the biotools package in your browser