loess_model {seedreg} | R Documentation |
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"
)
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) |
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
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)
Gabriel Danilo Shimizu
Leandro Simoes Azeredo Goncalves
See Also
Examples
library(seedreg)
data("aristolochia")
attach(aristolochia)
#================================
# Germination
#================================
loess_model(trat,germ)
#================================
# Germination speed
#================================
loess_model(trat, vel, ylab=expression("v"~(dias^-1)))