Functions are provided for the density function,
	     distribution function, quantiles and random number
	     generation for the skew hyperbolic
	     t-distribution. There are also functions that fit
	     the distribution to data. There are functions for the
	     mean, variance, skewness, kurtosis and mode of a given
	     distribution and to calculate moments of any order about
	     any centre. To assess goodness of fit, there are
	     functions to generate a Q-Q plot, a P-P plot and a tail plot.

It fits scale mixture of skew-normal linear mixed models using either an expectation–maximization (EM) type algorithm or its accelerated version (Damped Anderson Acceleration with Epsilon Monotonicity, DAAREM), including some possibilities for modeling the within-subject dependence <doi:10.18637/jss.v115.i07>.

Provide data generation and estimation tools for the multivariate scale mixtures of normal 
             presented in Lange and Sinsheimer (1993) <doi:10.2307/1390698>, the multivariate scale 
             mixtures of skew-normal presented in Zeller, Lachos and Vilca (2011) 
             <doi:10.1080/02664760903406504>, the multivariate skew scale mixtures of normal 
             presented in Louredo, Zeller and Ferreira (2021) <doi:10.1007/s13571-021-00257-y>
             and the multivariate scale mixtures of skew-normal-Cauchy presented in Kahrari et al. (2020)
	     <doi:10.1080/03610918.2020.1804582>.

Estimate necessary sample sizes for comparing the location
    of data from two groups or categories when the distribution of the 
    data is skewed. The package 
    offers a non-parametric method for a Wilcoxon Mann-Whitney test of 
    location shift as well as methods for several generalized linear 
    models, for instance, Gamma regression.

Density, distribution function, quantile function and
        random generation for the skewed t distribution of Fernandez
        and Steel.

Provide estimation and data generation tools for the skew-unit family
		discussed based on Mukhopadhyay and Brani (1995) <doi:10.2307/2348710>. 
		The family contains extensions for popular distributions such as the ArcSin 
		discussed in Arnold and Groeneveld  (1980) <doi:10.1080/01621459.1980.10477449>, 
		triangular, U-quadratic and Johnson-SB proposed in Cortina-Borja (2006) <doi:10.1111/j.1467-985X.2006.00446_12.x> 
		distributions, among others.