Note that this reference documentation is identical to the help that is displayed in Matlab when you type “help ft_sourcestatistics”.

  FT_SOURCESTATISTICS computes the probability for a given null-hypothesis using
  a parametric statistical test or using a non-parametric randomization test.
  
  Use as
    [stat] = ft_sourcestatistics(cfg, source1, source2, ...)
  where the input data is the result from FT_SOURCEANALYSIS, FT_SOURCEDESCRIPTIVES
  or FT_SOURCEGRANDAVERAGE.  The source structures should be spatially alligned
  to each other and should have the same positions for the source grid.
 
  The configuration should contain the following option for data selection
    cfg.parameter  = string, describing the functional data to be processed, e.g. 'pow', 'nai' or 'coh'
 
  Furthermore, the configuration should contain:
    cfg.method       = different methods for calculating the probability of the null-hypothesis,
                     'montecarlo'    uses a non-parametric randomization test to get a Monte-Carlo estimate of the probability,
                     'analytic'      uses a parametric test that results in analytic probability,
                     'glm'           uses a general linear model approach,
                     'stats'         uses a parametric test from the Matlab statistics toolbox,
                     'parametric'    uses the Matlab statistics toolbox (very similar to 'stats'),
                     'randomization' uses randomization of the data prior to source reconstruction,
                     'randcluster'   uses randomization of the data prior to source reconstruction 
                                     in combination with spatial clusters.
 
  You can restrict the statistical analysis to regions of interest (ROIs)
  or to the average value inside ROIs using the following options:
    cfg.atlas        = filename of the atlas
    cfg.roi          = string or cell of strings, region(s) of interest from anatomical atlas
    cfg.avgoverroi   = 'yes' or 'no' (default = 'no')
    cfg.hemisphere   = 'left', 'right', 'both', 'combined', specifying this is
                       required when averaging over regions
    cfg.inputcoord   = 'mni' or 'tal', the coordinate system in which your source 
                       reconstruction is expressed
 
  The other cfg options depend on the method that you select. You
  should read the help of the respective subfunction STATISTICS_XXX
  for the corresponding configuration options and for a detailed
  explanation of each method.
 
 
  See also FT_SOURCEANALYSIS, FT_SOURCEDESCRIPTIVES, FT_SOURCEGRANDAVERAGE