POP_JOINTPROB Reject artifacts in an EEG dataset using joint probability of the recorded electrode or component activities observed at each time point. e.g., Observing large absoluate values at most electrodes or components is improbable and may well mark the presence of artifact.


Usage: >> pop_jointprob( INEEG, typerej) % pop-up interative window mode
>> [OUTEEG, locthresh, globthresh, nrej] = = pop_jointprob( INEEG, typerej, elec_comp, locthresh, globthresh, superpose, reject, vistype);

Graphic interface:
"Electrode"   
[edit box] electrode|component number(s) to take into consideration for rejection. Sets the 'elec_comp' parameter in the command line call (see below).
"Single-channel limit(s)"   
[edit box] activity probability limit(s) (in std. dev.) Sets the 'locthresh' command line parameter. If more than one, defined individual electrode|channel limits. If fewer values than the number of electrodes | components specified above, the last input value is used for all remaining electrodes|components.
"All-channel limit(s)"   
[edit box] activity probability limit(s) (in std. dev.) for all channels (grouped). Sets the 'globthresh' command line parameter.
"Display with previously marked rejections?"   
[edit box] either YES or NO. Sets the command line option 'superpose'.
"Reject marked trial(s)?"   
[edit box] either YES or NO. Sets the command line option 'reject'.
"visualization mode"   
[edit box] either REJECTRIALS or EEGPLOT. Sets the command line option 'vistype'.

Inputs:
INEEG   
input dataset
typerej   
[1|0] data to reject on (0 = component activations;
1   
electrode data). {Default: 1 = electrode data}.
elec_comp   
[n1 n2 ...] electrode|component number(s) to take into consideration for rejection
locthresh   
activity probability limit(s) (in std. dev.) See "Single- channel limit(s)" above.
globthresh   
global limit(s) (all activities grouped) (in std. dev.)
superpose   
[0|1] 0 = Do not superpose rejection marks on previously marks stored in the dataset: 1 = Show both current and previous marks using different colors. {Default: 0}.

reject - 0 = do not reject marked trials (but store the marks:
1   
reject marked trials {Default: 1}.
vistype   
visualization mode: 0 = rejstatepoch(); 1 = eegplot() {Default: 0}.

Outputs:
OUTEEG   
output dataset with updated joint probability array
locthresh   
electrodes probability of activity thresholds in terms of standard-dev.
globthresh   
global threshold (where all electrode activity are regrouped).
nrej   
number of rejected sweeps

Author: Arnaud Delorme, CNL / Salk Institute, 2001

See also: jointprob(), rejstatepoch(), eegplot(), eeglab(), pop_rejepoch()

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

Back to functions