Estimate and confidence/credible intervals for an unknown
    regressor x0 given an observed y0.

Perform two linear combination methods for biomarkers: (1) Empirical performance optimization for specificity (or sensitivity) at a controlled sensitivity (or specificity) level of Huang and Sanda (2022) <doi:10.1214/22-aos2210>, and (2) weighted maximum score estimator with empirical minimization of averaged false positive rate and false negative rate. Both adopt the algorithms of Huang and Sanda (2022) <doi:10.1214/22-aos2210>. 'MOSEK' solver is used and needs to be installed; an academic license for 'MOSEK' is free.

Sample truncated multivariate Normal distribution following Gessner, A., Kanjilal, O., & Hennig, P. (2019). Integrals over Gaussians under Linear Domain Constraints. 108. <arxiv:1910.09328>.

L-systems or Lindenmayer systems are parallel rewriting systems which can
    be used to simulate biological forms and certain kinds of fractals.
    Briefly, in an L-system a series of symbols in a string are replaced
    iteratively according to rules to give a more complex string. Eventually,
    the symbols are translated into turtle graphics for plotting. Wikipedia has
    a very good introduction: en.wikipedia.org/wiki/L-system This package
    provides basic functions for exploring L-systems.

Provides a set of streamlined functions that allow
    easy generation of linear regression diagnostic plots necessarily 
    for checking linear model assumptions.
    This package is meant for easy scheming of linear 
    regression diagnostics, while preserving merits of 
    "The Grammar of Graphics" as implemented in 'ggplot2'.
    See the 'ggplot2' website for more information regarding the
    specific capability of graphics.

Computes the probability density function, the cumulative distribution function, the hazard rate function, the quantile function and random generation for Lindley Power Series distributions, see Nadarajah and Si (2018) <doi:10.1007/s13171-018-0150-x>.