Implementation of Taylor Regression Estimator (TRE), 
   Tulip Extreme Finding Estimator (TEFE), Bell Extreme Finding Estimator (BEFE),
   Integration Extreme Finding Estimator (IEFE) and 
   Integration Root Finding Estimator (IRFE) for roots, extrema and inflections of a curve .     
   Christopoulos, DT (2019) <doi:10.13140/RG.2.2.17158.32324> .
   Christopoulos, DT (2016) <doi:10.2139/ssrn.3043076> .
   Christopoulos, DT (2016) <https://demovtu.veltech.edu.in/wp-content/uploads/2016/04/Paper-04-2016.pdf> .
   Christopoulos, DT (2014) <doi:10.48550/arXiv.1206.5478> .

Routines to find the root of nonlinear functions, and to perform steady-state and equilibrium analysis of ordinary differential equations (ODE). 
  Includes routines that: (1) generate gradient and jacobian matrices (full and banded),
  (2) find roots of non-linear equations by the 'Newton-Raphson' method, 
  (3) estimate steady-state conditions of a system of (differential) equations in full, banded or sparse form, using the 'Newton-Raphson' method, or by dynamically running,
  (4) solve the steady-state conditions for uni-and multicomponent 1-D, 2-D, and 3-D partial differential equations, that have been converted to ordinary differential equations
    by numerical differencing (using the method-of-lines approach).
  Includes fortran code.

Provides distribution functions and log-scale tail probabilities for Roy's largest root
  in single and double Wishart (Jacobi ensemble) settings. This package is derived from the
  'rootWishart' package by Maxime Turgeon and extends it with numerically robust log-CDF/log-survival
  evaluation, tail-aware adaptive precision, and high-dimensional validation utilities, based on
  Chiani (2014) <DOI:10.1016/j.jmva.2014.04.002> and Chiani (2016) <DOI:10.1016/j.jmva.2015.10.007>.

Selects one model with variable selection FDR controlled at a
    specified level. A q-value for each potential variable is also returned. The
    input, variable selection counts over many bootstraps for several levels of
    penalization, is modeled as coming from a beta-binomial mixture
    distribution.