R/mean.R

In csgillespie/voltagefit: Voltage Curve Fitting

get_curve = function(pars, dev_curve, x = seq(-4, 6, length.out=1000)) {
  pars = as.numeric(pars)
  y = dev_curve(x, pars)
  data.frame(VG = x, ID = exp(y), stringsAsFactors = FALSE)
}


#' @export
mean.fit_manova = function(x, direction="Forward", ...){
  pars = get_params(x, combine=TRUE)
  #pars = pars[pars$direction == direction, , drop=FALSE]
  
  par_loc = get_par_loc(pars)
  dd = NULL; i = 1
  for(i in seq_len(nrow(pars))) {
    dd_tmp = get_curve(pars[i, par_loc], dev_curve = attr(x, "dev_curve"))
    dd_tmp$type = pars$type[i]
    dd_tmp$direction = pars$direction[i]
    dd = rbind(dd, dd_tmp)
  }
  dd = trans_device(dd)
  dd
}