eTSMSN {TSMSN} | R Documentation |
Estimate the parameters of Truncated Scale Mixtures of Skew-Normal Distributions
Description
This function obtains the maximum likelihood estimators of the TSMSN (Skew-Normal, Skew-t, Skew-Slash and Skew-Contaminated Normal) distribution parameters by direct maximization.
Usage
eTSMSN(x, mu = 0.01, sigma2 = 1.01, lambda = 0, nu = 5, a = -Inf,
b = Inf, dist = "SN", shape = FALSE)
Arguments
x |
Dataset. |
mu |
Initial location parameter (Optional). |
sigma2 |
Initial scale parameter (Optional). |
lambda |
Initial skewness parameter (Optional). |
nu |
Fixed shape parameter. Must be NULL in case of Skew-Normal distribution. Must be a bidimensional vector in case of skew-contaminated normal distribution (SCN) and contaminated normal distribution (CN). |
a |
Lower bound. |
b |
Upper bound. |
dist |
Distribution to be used: "SN" for Skew-Normal model, "ST" for Skew-t model, "SSL" for Skew-slash model and "SCN" for Skew-contaminated Normal model. |
shape |
For ST, SSL and SCN distribution. Consider the parameter nu as fixed and known. If TRUE nu must be provided. |
Details
For the SMN family, consider lambda = 0. For the Skew-contaminated Normal and Contaminated Normal distribution, each component of the bidimensional vector "nu" must be on (0,1). For the estimation in the cases of distributions ST, SSL and SCN nu is considered fixed, but may be known or unknown. The shape parameter is the one that regulates if nu is known or unknown.
Value
Returns a vector with the maximum likelihood estimators of the distribution parameters.
References
Lachos, V. H.; Garay, A. M.; Cabral, C. R. "Moments of truncated scale mixtures of skew-normal distributions." Brazilian Journal of Probability and Statistics (In press).
Basso, Rodrigo M., et al. "Robust mixture modeling based on scale mixtures of skew-normal distributions." Computational Statistics & Data Analysis 54.12 (2010): 2926-2941.
Examples
x <- rTSMSN(n = 100, mu = 0, sigma2 = 1, lambda = 0, nu = NULL, a = -2, b = 2, dist = "SN")
eTSMSN(x, a = -2, b = 2, dist = "SN")