Methods for calculating (usually) accurate
	numerical first and second order derivatives. Accurate calculations 
	are done using 'Richardson''s' extrapolation or, when applicable, a
	complex step derivative is available. A simple difference 
	method is also provided. Simple difference is (usually) less accurate
	but is much quicker than 'Richardson''s' extrapolation and provides a 
	useful cross-check. 
	Methods are provided for real scalar and vector valued functions.

High-dimensional datasets that do not exhibit a clear intrinsic clustered structure pose a challenge to conventional clustering algorithms. For this reason, we developed an unsupervised framework that helps scientists to better subgroup their datasets based on visual cues, please see Gao S, Mutter S, Casey A, Makinen V-P (2019) Numero: a statistical framework to define multivariable subgroups in complex population-based datasets, Int J Epidemiology, 48:369-37, <doi:10.1093/ije/dyy113>. The framework includes the necessary functions to construct a self-organizing map of the data, to evaluate the statistical significance of the observed data patterns, and to visualize the results.

Format numbers and plots for publication; includes the removal of leading zeros,
           standardization of number of digits, addition of affixes, and a p-value formatter. These
           tools combine the functionality of several 'base' functions such as 'paste()',
           'format()', and 'sprintf()' into specific use case functions that are named in a way
           that is consistent with usage, making their names easy to remember and easy to deploy.

To add the table of numbers at risk below the Kaplan-Meier plot.

Optimizing regular numeric problems in optically stimulated luminescence dating, 
        such as: equivalent dose calculation, dose rate determination, growth curve fitting, decay 
        curve decomposition, statistical age model optimization, and statistical plot visualization.

Tools for detecting numeric persistence ("spells") and rigidity patterns
    in time-ordered numeric data. The package identifies periods of stability,
    computes spell-based rigidity metrics, and provides plain-language
    interpretations suitable for policy and applied analysis.