R/GetBinNum.R

In functionaldata/tPACE: Functional Data Analysis and Empirical Dynamics

GetBinNum = function(n, m, dataType, verbose ){
  
  # Get the number of bins 
  # n : number of curves
  # m : median or max value of number of time-points
  # dataType : indicator about structure of the data 
  #     (dense (2), or  dataType data with missing values (1) or sparse (0))
  # verbose : outpit diagnostics/progress
  
  numBin = NULL;
  if (m <= 20){
    if (dataType =='Sparse'){
      str = 'Median of ni';
    } else {
      str = 'Maximum of ni';
    }
    if (verbose){
      message(paste0(str, 'is no more than 20! No binning is performed!\n')) 
    }
    return(NULL)
  }

  
  if (m >400){
    numBin = 400;
  }
  
  if (n > 5000){
    mstar = max(20,(((5000-n)*19)/2250)+400);
    if (mstar < m){
      numBin = ceiling(mstar);
    } else {
      if (verbose){
        message('No binning is needed!\n');    
      }
      return(NULL)
    }
  }
  
  if( verbose && is.null(numBin) ) {   
    message('No binning is needed!\n');
  }

  return(numBin)
}