In the Installation Notes (Sep 2018) we tried to document all that you have to expect when FSDA is installed manually by unpacking the compressed tar file FSDA.tar.gz, or automatically with our setup program for Windows platforms.


FSDA TOOLBOX DOWNLOAD (Release 2018b available from September 28th, 2018)

DOWNLOAD FSDA. A setup executable for MS Windows platforms will install the toolbox and update the search path of your local MATLAB installation.

DOWNLOAD A COMPRESSED TAR FILE of the toolbox, suitable for Unix platforms installation. In this case, you have to add manually FSDA folder and sub-folders to the MATLAB path or use our routine addFSDA2path.

Download a working paper describing the main characteristics of the FSDA toolbox or see Riani, Perrotta, Torti (2012).

Highlights of the last releases

FSDA follows typical MATLAB timetable in the sense that there are two releases per year. The first typically in May/June and the second around October/November

We decided to label FSDA release with the most recent MATLAB version to which it is associated.

September 2018 Release 2018b.

(1) New function qqplotFS that enables to create a qqplot of residuals with confidence bands

(2) New function  mtR which generates the same random numbers produced by R software with Mersenne Twister mt19937ar

(3) New functions associated with Rocke biweght estimator. See for example RKrho, RKpsi, RKpsider, RKwei, RKbdp, RKeff.

(4) Routines FSR, FSRmdr, FSRbsb extended to time series (see new functions FSRts, FSRtsmdr, FSRtsbsb and regressts)

(5) New function verlessthanFS.  It is a faster version of MATLAB function verlessthan.

(6) New datasets added to the collection.

(7) New routine publishBibliography to create in a automatic way the bibliography from the citations present inside the .m files.

June 2018 Release 2018a (Download .tar .exe)

(1) New function tclusteda that helps choosing the best tclust model. It computes tclustfor different values of the trimming factor and produces plots that allow to find the optimal level of trimming. This function uses the parallel processing toolbox, if available. 

(2) Extension of the score test. New function ScoreYJpn  that computes the score test for Yeo Johnson transformation separately for positive and negative observations. FSRfan now accepts the new option family "YJpn" and it is possible to monitor the score test for both positive and negative observations (output arguments out.Scorep and out.Scoren).

(4) New functions for time series analysis. simulateTS simulates a time series with trend (up to third order), seasonality (constant or of varying amplitude) with a different number of  harmonics and a level shift. forecastTS produces forecasts with confidence bands for a time series estimated with function LTSts.   

(4) CorAna has an improved display of results. New function CorAnaplot draws a rich Correspondence Analysis graph with different types of confidence ellipses for selected rows and columns.

(5) New function verlessthanExt.  It is a faster version of MATLAB function verlessthan.

(6) Documentation of yXplot considerably improved. New options added (xlimx, ylimy, namey, nameX).

(7) MixSimreg extended to account for multiple parameter distribution (betadistrib option) 

(8) histFS has a new optional argument (weights) for plotting a weighted histogram.

(9) options labenv has been added to mmdrsplot.

(10) option axesellipse added to ellipse

(11) New output argument OldAndNewIndexes used in function UnitsSameCluster, to track the indexes permutations used to rearch a desired cluster labelling.


November 2017 Release 2017b.

Major statistical release. Highlights:

FSDA has introduced two new categories of tools, one for (robust) time series analysis; another for analyzing categorical data and contingency tables. More precisely:

(1) Function LTSts  extends LTS estimator to time series. A related new graphical plot associated to a time series, wedgeplot, provides information on the presence of outliers and level shifts.

(2) CorAna performs correspondence analysis; SparseTableTest computes independence test for large and sparse contingency tables; CressieRead computes the power divergence family of tests, to check the discrepancy/distance between observed and expected frequencies in a contingency table; rcontFS generates a random two-way table with given marginal totals; barnardtest computes the Barnard test,  corrNominal measures strength of association between two unordered (nominal) categorical variables. Similarly for ordinal data with corrOrdinal. crosstab2datamatrix recreates the original data matrix X from contingency table N. This group of functions is complemented by file examples_categorical.m as in style of FSDA.

The two categories of functions will be progressively enriched.

Other new functions which are included are boxtest (test of equality of covariance matrices used for example in tkmeans), GYfilt (Gervini and Yohai, univariate outlier identifier), mmdrsplot (interactive plot of the trajectories of minimum Mahalanobis distances from different starting points), overlapmap to plot the ordered pairwise overlap values between components, dempk to perform a merging of components found by tkmeans, ncpci to compute a non centrality parameter confidence interval. 

Finally, spmplot has been enriched to superimpose ellipses, density and contour functions to data and extract single panels from the scatter matrix.

May 2017 Release 2017a.

Major statistical release. Highlights:
    Function tclustreg now includes trimmed Cluster Weighted Restricted Models.
    New function tclustIC for the automatic selection of the best number of groups.
    New function tlclustICsol to extract a set of relevant solutions (and associated tclustICplot).
    New function UnitsSameCluster to to control the labels of the clusters which contain predefined units.
    New function to compare two partitions (Fowlkes and Mallows index)
    Updated routine simdatasetreg to generate new outlier patterns


    New routines for density estimation and thinning, for univariate and bivariate data (used in tclustreg).
    bwe, rthin, wthin, WNChygepdf.

    New functions wraptextFS, removeextraspacesLF

    New functions mveeda, MMmulteda
    Smulteda, Sregeda

     New function ScoreYJ which implements the score test for the Yeo and Johnson transformation. This new transformation has also been embedded inside function FSRfan.

    updated functions randsampleFS and subsets.
    New routines for thinning.

  • NEW mlx files
    .mlx files introduced for examples_MixSim

    New function to create the car-bike plot to find the most relevant solutions (carbikeplot).
    Functions resfwdplot, malfwdplot generalized in order to take as input the output of procedures which monitor robust estimators.

  • The FSDA help folder now contains XML files associated to the functions documentation. This is in view of generating/updating automatically or using a GUI the functions documentation, in html as well as in the function head


October 2016 Release 2016b.

Major statistical release. Highlights:
  • New functions for  bivariate density estimation and random thinning (kdebiv.m, rthin.m) used to extend tclustreg.m features.
  • The FSDA help folder now contains XML files associated to the functions documentation. This is in view of generating/updating automatically or using a GUI the functions documentation, in html as well as in the function head


June 2016 Release 2016a.

Major statistical release. Highlights:
  • New features added to the tclust function, including determinant restriction and new adjusted BIC criterion for the estimation of the number of groups.
  • Added functions for reweighting FSR and FSRB (FSRr and FSRBr).
  • Functions FSR, FSRB and FSRH redesigned; a routing implementing the core of the Forward Search algorithm (FSRcore) introduced to avoid code redundancies.
  • New function, winsor, to winsor data.
  • New function FSMbsb, which will replace FSMbbm.
  • New function randindexFS, to evaluate the quality of different clusterings.
  • New routines poolClose and poolPrepare introduced to conveniently open and close a pool of parallel workers.
  • Several new robust functions to generate, for example, the Tukey Biweigh rho function (HUrho), the tuning constant associated to a certain efficiency (HUeff), the psi functions (HUpsi), its derivative (HUpsider), etc. For a full list, see functions under utilities_stats folder.

Major highlights in the documentations and examples:

  • New function, publishFS, introduced to generate documentation pages directly from the .m files.
  • New function, makecontentsfileFS, introduced to create a the list of files present in a FSDA folder and/or subfolders. It extends MATLAB function makecontentsfile.
  • .mlx files introduced for examples_multivariate and examples_regression

September 2015 Release 2015b.

  • Function simdataset.m modified to allow the user to simulate outliers from different distributions and contamination schemes and/or contaminate existing datasets.
  • New Bayesian regression analysis routines: FSRB.m, FSRBeda.m, FSRBmdr.m, regressB.m.
  • In FSReda.m: monitoring of confidence intervals of beta and sigma2.
  • In FSRBeda.m: monitoring of HPD (highest posterior density regions) of beta and sigma2.
  • New functions for inverse gamma computation: inversegampdf.m, inversegamcdf.m, inversegaminv.m.
  • Added functions to monitor units forming subset in heterosckedastic and Bayesian regression: FSRHbsb.m, FSRBbsb.m.
  • Added new datasets for Bayesian examples.
  • Added option for the robust transformation in the Yeo-Johnson family.
  • addFSDA2path.m: modified for compatibility with unix platforms and to address changes in the folder organization of FSDA functions.
  • Added routines to compute and visualize robust bivariate boxplot (function boxplotb.m).
  • New routine for automatic outlier detection in heteroskedastic regression (FSRH.m).

Please, do not hesitate to contact us for any bug you might find and for any suggestion you might have!!! Thanks to all those who have already contacted us and have helped us to correct several bugs and improve the performance of the code.

Marco Riani, Andrea Cerioli and Aldo Corbellini (University of Parma)

Domenico Perrotta and Francesca Torti (EC, Joint Research Centre)

Patrizia Calcaterra, Andrea Cerasa, Emmanuele Sordini, Daniele Palermo (EC, Joint Research Centre)

FSDA Toolbox Download