Receives two vectors, computes appropriate function for group comparison (i.e., t-test, Mann-Whitney; equality of variances), and reports the findings (mean/median, standard deviation, test statistic, p-value, effect size) in APA format (Fay, M.P., & Proschan, M.A. (2010)<DOI: 10.1214/09-SS051>).

Methods for dividing data into groups. 
    Create balanced partitions and cross-validation folds. 
    Perform time series windowing and general grouping and splitting of data. 
    Balance existing groups with up- and downsampling or collapse them to fewer groups.

These 'Rcpp'-based functions compute the efficient score statistics for grouped time-to-event data (Prentice and Gloeckler, 1978), with the optional inclusion of baseline covariates. Functions for estimating the parameter of interest and nuisance parameters, including baseline hazards, using maximum likelihood are also provided. A parallel set of functions allow for the incorporation of family structure of related individuals (e.g., trios). Note that the current implementation of the frailty model (Ripatti and Palmgren, 2000) is sensitive to departures from model assumptions, and should be considered experimental. For these data, the exact proportional-hazards-model-based likelihood is computed by evaluating multiple variable integration. The integration is accomplished using the 'Cuba' library (Hahn, 2005), and the source files are included in this package. The maximization process is carried out using Brent's algorithm, with the C++ code file from John Burkardt and John Denker (Brent, 2002).

The 'Grouphmap' was implemented in R, an open-source programming environment, and was released under the provided website. The difference analysis is based on the 'limma' package, which can cover gene and protein expression profiles (Reference: Matthew E Ritchie , Belinda Phipson , Di Wu , Yifang Hu , Charity W Law , Wei Shi , Gordon K Smyth (2015) <doi:10.1093/nar/gkv007>). The GO enrichment analysis is based on the 'clusterProfiler' package and supports three common species: human, mouse, and yeast (Reference: Guangchuang Yu, Li-Gen Wang, Yanyan Han, Qing-Yu He (2012) <doi:10.1089/omi.2011.0118>). The results of batch difference analysis and enrichment analysis are output in separate folders for easy viewing and further visualization of the results during the process. The results returned a heatmap in R and exported to 3 folders named DEG, go, and merge.