LL_model {seedreg}

Analysis: Logistic regression

Description

Logistic models with three (LL.3) or four (LL.4) continuous data parameters. This model was extracted from the drc package and adapted for temperature analysis in seed germination.

Usage

LL_model(
  trat,
  resp,
  npar = "LL.3",
  error = "SE",
  ylab = "Germination (%)",
  xlab = expression("Temperature ("^"o" * "C)"),
  theme = theme_classic(),
  legend.position = "top",
  cardinal = 0,
  r2 = "all",
  width.bar = NA,
  scale = "none",
  textsize = 12,
  pointsize = 4.5,
  linesize = 0.8,
  pointshape = 21,
  font.family = "sans"
)

Arguments

trat

Numerical or complex vector with treatments

resp

Numerical vector containing the response of the experiment.

npar

Number of model parameters

error

Error bar (It can be SE - default, SD or FALSE)

ylab

Variable response name (Accepts the expression() function)

xlab

Treatments name (Accepts the expression() function)

theme

ggplot2 theme (default is theme_bw())

legend.position

Legend position (default is c(0.3,0.8))

cardinal

Defines the value of y considered extreme (default considers 0 germination)

r2

Coefficient of determination of the mean or all values (default is all)

width.bar

Bar width

scale

Sets x scale (default is none, can be "log")

textsize

Font size

pointsize

shape size

linesize

line size

pointshape

format point (default is 21)

font.family

Font family (default is sans)

Details

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).

Value

Coefficients

Coefficients and their p values

Optimum temperature

Optimum temperature (equivalent to the maximum point)

Optimum temperature response

Response at the optimal temperature (equivalent to the maximum point)

Minimal temperature

Temperature that has the lowest response

Minimal temperature response

Lowest predicted response

Predicted maximum basal value

Lower basal limit temperature based on the value set by the user (default is 0)

Predicted minimum basal value

Upper basal limit temperature based on the value set by the user (default is 0)

AIC

Akaike information criterion

BIC

Bayesian Inference Criterion

r-squared

Determination coefficient

RMSE

Root mean square error

grafico

Graph in ggplot2 with equation

Note

if the maximum predicted value is equal to the maximum x, the curve does not have a maximum point within the studied range. If the minimum value is less than the lowest point studied, disregard the value.

Author(s)

Model imported from the drc package (Ritz et al., 2016)

Gabriel Danilo Shimizu

Leandro Simoes Azeredo Goncalves

References

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.

Examples

library(seedreg)
data("aristolochia")
attach(aristolochia)

#================================
# Germination
#================================
LL_model(trat,germ)

#================================
# Germination speed
#================================
LL_model(trat, vel, ylab=expression("v"~(dias^-1)))