JOINTPROB Rejection of odd columns of a data array using joint probability of the values in that column (and using the probability distribution of all columns).

Usage: >> [rej jp] = jointprob( signal );
>> [rej jp] = jointprob( signal, threshold, jp, normalize, discret);

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.
jp   
pre-computed joint probability (only perform thresholding). Default is the empty array [].
normalize   
0 = do not not normalize entropy. 1 = normalize entropy. Default is 0.
discret   
discretization variable for calculation of the discrete probability density. Default is 1000 points.

Outputs:
jp   
normalized joint probability of the single trials (size component x sweeps)
rej   
rejected matrix (0 and 1, size comp x sweeps)

Remark: The exact values of joint-probability depend on the size of a time
step and thus can not be considered as absolute.

See also: realproba

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

Back to functions