R/fitcond3.R

In jtominaga/dirci: Direct Ci Calculations

#' Fits cuticular and intercellular conductance to data
#' \code{fitcond3} Fits cuticular and intercellular conductances
#'
#' @param data A dataframe.
#' @param rtc..rtc Difference in total and stomatal resistances to CO2 diffusion
#' @param gsw  Stomatal conductance to water
#' @param gcw  Cuticular conductance to water
#' @param gias Intercellular conductance to water?
#' @param a Initial value for gcw in curve fitting
#' @param b Initial value for gias in curve fitting
#'
#' @return fitcond3 fits cuticular and intercellular conductances
#' @importFrom stats nls
#' @importFrom stats coef
#' @export
#'
#'
#'
fitcond3 <- function(data, rtc..rtc, gsw, gcw, gias, a, b){
  a <- ifelse(missing(a) == TRUE, 2, a)
  b <- ifelse(missing(b) == TRUE, 0.5, b)
  gcw <- ifelse(missing(gcw) == FALSE,
                ifelse(missing(gias) == FALSE,
                m      <- nls(rtc..rtc ~ 1.6 * (1 / (gsw - gcw) - 1 / gsw)
                              + 1 / gias, data = data),
                m      <- nls(rtc..rtc ~ 1.6 * (1 / (gsw - gcw) - 1 / gsw)
                              + 1 / gias, data = data,
                              start = list(gias = b))
                ),
                ifelse(missing(gias) == FALSE,
                       m      <- nls(rtc..rtc ~ 1.6 * (1 / (gsw - gcw) - 1 / gsw)
                                     + 1 / gias, data = data,
                                     start = list(gcw = a)),
                       m      <- nls(rtc..rtc ~ 1.6 * (1 / (gsw - gcw) - 1 / gsw)
                                     + 1 / gias, data = data,
                                     start = list(gcw = a, gias = b))))
  cond   <- coef(m)
  output <- as.list(c(m, cond, summary(m)))
  return(output)
}