loess_model: Analysis: loess regression

Description

Fit a polynomial surface determined by one or more numerical predictors, using local fitting.

Usage

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

Value

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)

grafico

Graph in ggplot2 with equation

Arguments

trat: Numerical or complex vector with treatments
resp: Numerical vector containing the response of the experiment.
ylab: Variable response name (Accepts the expression() function)
xlab: treatments name (Accepts the expression() function)
theme: ggplot2 theme (default is theme_bw())
error: Error bar (It can be SE - default, SD or FALSE)
cardinal: defines the value of y considered extreme (default considers 0 germination)
width.bar: bar width
legend.position: legend position (default is c(0.3,0.8))
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)

Author

Gabriel Danilo Shimizu

Leandro Simoes Azeredo Goncalves

Examples

Run this code

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

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

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