R/growthfd.evaluate.R

In ondrej-klima/growthfd: Fitting FPCA-based growth curve model

#' Generate a Discrete Growth Curve
#'
#' This function evaluates a curve function for given ages. Depending 
#' on a degree of derivation, the function produces stature, velocity 
#' or acceleration curve.
#' 
#' @param x Ages to be evaluated
#' @param par Parameters of the model
#' @param model FPCA growth model
#' @param deriv Path to the input file 
#' @return Y-values of the evaluated curve
#' @export
growthfd.evaluate <- function(x, par, model, deriv=0) {
  nscores <- sum(model$scores.elements);
  scores <- par[1:sum(nscores)];
  
  #message('eval begin')
  
  suppressWarnings(growthfd <- growthfd.std(scores, model))
  #message(sprintf('%f ', x))
  suppressWarnings(r <- fda::eval.fd(x, growthfd, deriv))
  
  #message('eval end')
  
  
  return(r);
}