R/get_skewness.R

In padpadpadpad/rTPC: Fitting and Analysing Thermal Performance Curves

#' Estimates skewness of a thermal performance curve
#'
#' @param model nls model object that contains a model of a thermal performance curve
#' @details Skewness is calculated from the values of activation energy (e) and deactivation energy (eh) as: skewness = e - eh. A negative skewness indicates the TPC is left skewed, the drop after the optimum is steeper than the rise up to the optimum. A positive skewness means that the TPC is right skewed and a value of 0 would mean the curve is symmetrical around the optimum.
#' @return Numeric estimate of skewness
#' @export get_skewness

get_skewness <- function(model){

  e <- rTPC::get_e(model)
  eh <- rTPC::get_eh(model)

  return(e - eh)
}