REJKURT Calculation of kutosis of a 1D, 2D or 3D array and rejection of outliers values of the input data array using the discrete kutosis of the values in that dimension.

Usage: >> [kurtosis rej] = rejkurt( signal, threshold, kurtosis, normalize);

Inputs:
signal   
one dimensional column vector of data values, two dimensional column vector of values of size sweeps x frames or three dimensional array of size component x sweeps x frames. If three dimensional, all components are treated independantly.
threshold   
Absolute threshold. If normalization is used then the threshold is expressed in standard deviation of the mean. 0 means no threshold.
kurtosis   
pre-computed kurtosis (only perform thresholding). Default is the empty array [].
normalize   
0 = do not not normalize kurtosis. 1 = normalize kurtosis. Default is 0.

Outputs:
kurtosis   
normalized joint probability of the single trials (same size as signal without the last dimension)
rej   
rejected matrix (0 and 1, size: 1 x sweeps)

Remarks: The exact values of kurtosis depend on the size of a time
step and thus can not be considered as absolute.
This function uses the kurtosis function from the statistival
matlab toolbox. If the statistical toolbox is not installed,
it uses the 'kurt' function of the ICA/EEG toolbox.

See also: realproba

See the matlab file rejkurt.m (may require other functions)

Back to functions