Logistic regression is a very popular analysis in agrarian sciences, such as in fruit growth curves, seed germination, etc...The logistic function performs the analysis using 3 or 4 parameters of the logistic model, being imported from the LL function .3 or LL.4 of the drc package (Ritz & Ritz, 2016).
logistic(
trat,
resp,
npar = "LL.3",
error = "SE",
ylab = "Dependent",
xlab = expression("Independent"),
theme = theme_classic(),
legend.position = "top",
r2 = "all",
width.bar = NA,
scale = "none",
textsize = 12,
font.family = "sans"
)
The function allows the automatic graph and equation construction of the logistic model, provides important statistics, such as the Akaike (AIC) and Bayesian (BIC) inference criteria, coefficient of determination (r2), square root of the mean error ( RMSE).
Numerical or complex vector with treatments
Numerical vector containing the response of the experiment.
Number of model parameters
Error bar (It can be SE - default, SD or FALSE)
Variable response name (Accepts the expression() function)
Treatments name (Accepts the expression() function)
ggplot2 theme (default is theme_bw())
Legend position (default is c(0.3,0.8))
Coefficient of determination of the mean or all values (default is all)
Bar width
Sets x scale (default is none, can be "log")
Font size
Font family (default is sans)
Model imported from the drc package (Ritz et al., 2016)
Gabriel Danilo Shimizu
Leandro Simoes Azeredo Goncalves
The three-parameter log-logistic function with lower limit 0 is $$f(x) = 0 + \frac{d}{1+\exp(b(\log(x)-\log(e)))}$$ The four-parameter log-logistic function is given by the expression $$f(x) = c + \frac{d-c}{1+\exp(b(\log(x)-\log(e)))}$$ The function is symmetric about the inflection point (e).
Seber, G. A. F. and Wild, C. J (1989) Nonlinear Regression, New York: Wiley and Sons (p. 330).
Ritz, C.; Strebig, J.C.; Ritz, M.C. Package ‘drc’. Creative Commons: Mountain View, CA, USA, 2016.
data("emerg")
with(emerg, logistic(time, resp,xlab="Time (days)",ylab="Emergence (%)"))
with(emerg, logistic(time, resp,npar="LL.4",xlab="Time (days)",ylab="Emergence (%)"))
Run the code above in your browser using DataLab