Changes in version 2.1.3 (2026-02-24) o In msn.mle and its auxiliary functions, exact symmetry of the estimated Omega matrix is being forced, to avoid a numerical issue arising in particular instances. o Working of the show method for selm and mselm classes has been reinstated. This had effectively been inhibited in a previous version of sn by unfortunate re-definitions of the methods. o In function st.prelimFit, the call to quantreg::rq.wfit now has method="fn" to reduce the chance of warning messages “Solution may be nonunique”. o Fixed a minor bug in convolutionSUNdistr. o Various scattered improvements in the documentation. Changes in version 2.1.2 (2026-01-25) o Fixed a bug which affected usage of ... in pst and qst, and improved documentation of these functions. o Improved documentation in "SUNdistr-base". Changes in version 2.1.1 (2023-04-04) o Coding of pst is improved for extreme x arguments by employing asymptotic expressions of tail probabilites. As an implication, this change improves the working of qst. Also, new arguments lower.tail and log.p are added. o Checking on the dimensionality (d,m) of calls to {d,p}sun, sunMean, sunVcov, sunMardia takes place before a call to a function of mnormt is issued. If the maximal dimensionality is exceeded, either a stop is generated, or NAs are returned, depending on the value of the silent argument. o Fixed bugs affecting calls to psn, pst when length(x) was shorter than the length of the parameters (among xi, omega, alpha). Changes in version 2.1.0 (2022-08-11) o New function fitdistr.grouped, with pertaining class and methods. o More informative messages are issued when trace=TRUE, at several places. o An experimental way for initializing the parameters of sn.mple fitting is introduced; it only takes effect when it improves over the existing scheme. o Some adjustments in the documenation on request of the CRAN group. o Fixed a bug affecting st.prelimFit when the argument w was used. o The NEWS file is now in Rd format. Changes in version 2.0.2 (2022-03-08) o Calling rmsn, rmst, rmsc with parameters for distributions with dimension 1 now works. o The same feature applies to mst.prelimFit with one column matrix y. o In psn, improved handling of non-finite x values (NA, NaN, Inf), which led to crashes; also, streamlined coding of qsn. o Minor fixes in the documentation of zeta and qsn; improved wording elsewhere. Changes in version 2.0.1 (2021-11-26) o Change of the tuning arguments when numDeriv::hessian is called, to improve computation of the information matrix following a ST model fitting. o Fixed bug which in some cases prevented the use of the 'start' argument. o Re-organization of some internal functions computing SUN summary quantities. o Improved documentation for selm, selm.fit and SUNdistr-base. Changes in version 2.0.0 (2021-03-28) o Support for the SUN family is introduced, as for probability distribution operations. Two modes of working are envisaged: (i) using classic-style functions for probability distributions, plus some functions of similar style; (ii) using S4 objects of the new class SUNdistr. For mode (i), the functions {d,p,r}sun, sun{Mean,Vcov,Mardia} are provided. For mode (ii), there are {make, marginal, conditional, affineTrans, join, convolution, summary}SUNdistr, and related S4 methods. o Additional facilities include convertSN2SUNdistr, convertCSN2SUNpar and two matrix operations (tr and blockDiag). o Fixed a bug in plot.SECdistr affecting the plot of a subset of the variables. Changes in version 1.6-2 (2020-05-26) o Fixed wrong computation of standard errors when a multivariate ST model was fitted with the constraint alpha=0; similar fix of function confint. o Change of the algorithm used in rsn: the additive representation is now used both fo r 0 and for non-0 values of tau. o Use of the vignette builder R.rsp. Changes in version 1.6-1 (2020-04-01) o Some minor changes in the selm function documentation. Changes in version 1.6-0 (2020-03-28) o New intialization technique for numerical MLE search when function selm is called with family="ST". This is performed by the new functions st.prelimFit and mst.prelimFit, with the aid of galton_moors2alpha_nu. o In st.prelimFit, the package quantreg is employed for preliminary linear predictor estimation. o Related new functions of more general interest are: fournum, pprodn2, pprodt2, qprodt2. o Additional facts: an improved version of function profile.selm. o Improved coding of some internal functions. Changes in version 1.5-5 (2020-01-30) o In sn.infoMv, new argument at.MLE is introduced. o Arguments of selm are updated to match changes in R. o Fix a bug in dsn when called with a non-scalar argument alpha. o Changed internal checks on 'try' output, to accomplish R changes. Changes in version 1.5-4 (2019-05-14) o On request from the CRAN group, a modification is inserted to get around numerical problems arising in connection with use of OpenBLAS (version 0.3.5). o For plotting of a multivariate SECdistr when 'range' is not supplied and 'data' is non-NULL, calculation of the plotting range has been modified. o Proper handling of the 'name' argument of marginalSECdistr when the result is a univariate distribution. Changes in version 1.5-3 (2018-11-08) o In pst, improved implementation of Method 2; slight modification of the automatic selection method when 'method=0'. o Internal function qst_bounds introduced for better initial bracketing of ST quantiles; improved qst coding. o In profile.selm the selected parameter area does no longer need to include the MLE/MPLE point. o In sn.infoMv removed bugs in the stage of parameter parsing, which prevented computing the expected information matrix. o In internal functions st.pdev.gh and mst.pdev.grad, improved computation of 'nu' component of logLik gradient, yielding faster fitting of ST models. o Improved handling of arguments of sn.infoMv. o Improved documentation of modeSECdistr and coding of modeSECdistrMv. Changes in version 1.5-2 (2018-04-24) o Improved checking of input arguments to lower level fitting procedures (those below selm) and improved handling for those of pst. o Fixed improper handling when not positive-definite information in st.infoUv occurs, pointed out by the CRAN group. Changes in version 1.5-1 (2017-11-22) o More extensive documentation: - addition of 'overview' entry in standard documentation; - in directory 'doc', inclusion of a PDF file providing a tutorial introduction to the package; - additions and improvements at various places in Rd files. o Fix a minor bug in sn.infoUv causing crash when the DP information matrix is not invertible. o Improved numerical inversion of st.cp2dp in extreme situations. Changes in version 1.5-0 (2017-02-10) o Tools for user-defined symmetry-modulated (AKA skew-symmetric) distributions are introduced: {d,r}[m]SymmModulated and its bivariate density plotting. o Fixed a bug in dsn affecting the cases (a) x=Inf, alpha=0, and (b) omega<=0. Changes in version 1.4-0 (2016-06-30) o Introduce methods confint and predict for selm-class objects. o Fix bug in rmst causing some dependence among subsequent samples. o Fix bug of modeSECdistrMv affecting certain ST cases. o Now plot.SECdistrBv allows to overlap plots. o Improved naming of output. o profile.selm can now be called with vector(s) param.values of length 1. Changes in version 1.3-0 (2015-11-11) o Method profile.selm is introduced. o The object returned by plot.SECdistrMv now includes the coordinates of the contour curves. o Fixed a bug affecting rmsn when called using dp=, also dp[[1]] is named beta instead of xi. Changes in version 1.2-5 o Not released Changes in version 1.2-4 (2015-08-25) o Output of plot.SECdistr is better structured and documented. o In pmst, handling of case nu=Inf required a fix. o Corrected a bug of internal function msn.dp2dp when called with aux=TRUE and d=1; this affected rmsn and rmst if d=1. Changes in version 1.2-3 (2015-07-16) o Fixed a bug in evaluation of the feasible CP parameter space of univariate ST. o Fixed a bug which crashed pmst when called with fractional degrees of freedom. o Functions dmsn, pmsn and dmst now expand a single value supplied as 'xi' into a vector or matrix of suitable dimension. Changes in version 1.2-2 (2015-06-05) o Fixed a bug in extractSECdistr from mselm-class objects. o Fixed a bug that prevented calling low level fitting functions with non-null 'penalty' argument. o Improved documentation of selm.fit and related functions. Changes in version 1.2-1 (2015-04-29) o Optimization parameters are now passed from selm to sn.mple and st.mple as indicated in the documentation. o Plotting of selm-class and mselm-class objects avoids clash of par('cex') parameters. o Computation of sn.infoMv now takes into account whether method="MPLE" was used at the estimation stage. Changes in version 1.2-0 (2015-03-25) o Created new functions extractSECdistr and modeSECdistr. o New methods mean and sd for class SECdistrUv, and new methods mean and vcov for class SECdistrMv. o Computation of qst switches to qsn if nu>1e4, instead of nu=Inf as before. o Fixed a bug in st.pdev.hessian (correction in args sequence). o Improved detection of singular distributions in selm output. o Improved handling of component names of SECdistr. Changes in version 1.1-2 (2014-12-01) o Fixed a bug affecting plotting of mselm objects under certain circumstances. o Fixed a bug affecting function selm when the weights argument contained 0's. o Improved coding in some functions. o More functions are exported and their documentation is added. Changes in version 1.1-1 (2014-10-30) o Function qsn has an additional argument 'solver'. o Functions pmsn and pmst can now be called with argument 'xi' of matrix type. o More functions are now exported in NAMESPACE. o Fixed a bug about selm.control argument of selm.fit. o Improved documentation of various functions. Changes in version 1.1-0 (2014-08-04) o Main few feature is the possibility to set the constraint alpha=0 in function selm and in lower level fitting functions. o Other additions or changes are: introduction of OP parameterization; fix a bug in qst; more efficient coding of dmsn and dmst; pmsn can now be called with argument 'x' of matrix type; in pst and qst, new argument method allows to select the algorithm employed. o More detailed documentation of pst and other functions and methods. Changes in version 1.0-0 (2014-01-07) o This is a major upgrade of the package, with much of the code completely new or largely re-written, leading to changes in the syntax and the user interface. o The key new functions are selm and makeSECdistr, with various related functions and methods. o S4 methods are adopted. o Many existing functions are updated, a few are dropped. o See help(SN) for more information. o (Development of "version 1" was started in June 2007.) Changes in version 0.4-18 (2013-05-01) o Various minor adjustments, many of them to fulfill CRAN programming standards. Changes in version 0.4-0 (2006-04-11) o many routines allow use of composite parameter 'dp'; o multivariate normal and t probabilities are now computed by 'mnormt'; o use of NAMESPACE introduced; o some more routines introduced, eg. st.cumulants.inversion; o various fixes/improvements in documentation. Changes in version 0.3-5 (2006-01-16) o Added some new functions (these include msn.affine, sn.mmle, sn.Einfo, sn.mle.grouped), fix various errors, and other improvements (eg. improved pst). Changes in version 0.3-0 o The main change is the addition of routines for (multivariate) skew-t distribution; also some other routines, e.g. mle for grouped data Changes in version 0.2-2 o Fix error in sn.dev.gh, improved qsn Changes in version 0.2-0 o This is the first public release and distribution, originally from in October 1998, shortly afterwards from CRAN.