[Eeglablist] Applying ICA weights from 1Hz filtered data to 0.1Hz filtered data

[chao wang]

Hi, I am new in ICA and have a question related to this issue. Each column
of ICA weight matrix is corresponding to a specific ICA component, but the
sequences of components resulting from two different ICA processes (e.g.
ica of 1 Hz and 0.01 Hz filtered datasets) might be very different. So, I
just wonder how it assures the same correspondence across the 2 datasets
when simply copying ICA weights from one to another.

Thanks,
Chao


On Sat, Aug 22, 2015 at 7:43 AM, <eeglablist-request at sccn.ucsd.edu> wrote:

>
> ---------- Forwarded message ----------
> From: Makoto Miyakoshi <mmiyakoshi at ucsd.edu>
> To: Neil Bailey <neil.bailey at monash.edu>
> Cc: EEGLAB List <eeglablist at sccn.ucsd.edu>
> Date: Fri, 21 Aug 2015 16:43:24 -0700
> Subject: Re: [Eeglablist] Applying ICA weights from 1Hz filtered data to
> 0.1Hz filtered data
> Dear Neil,
>
> I saw your code. I don't know what's wrong with it.
> Maybe after copying the ICA matrices, you want to try
>
> EEG = eeg_checkset(EEG, 'ica')'
>
> to let EEGLAB re-calculate ICA activations etc. I'm not sure though if it
> will solve the problem.
>
> Makoto
>
> On Wed, Aug 19, 2015 at 8:29 PM, Neil Bailey <neil.bailey at monash.edu>
> wrote:
>
>> Hi all
>>
>> I'm having trouble applying my ICA weights from 1Hz filtered data to my
>> 0.1Hz filtered data.
>>
>> I've been trying to use a script adapted from a previous post by Stefan
>> Debener to take the ICA weights out of the 1Hz data and put them into the
>> 0.1Hz data. The script is effective for transfering the ICA weights (they
>> show up in the EEG.ica... sections of the 0.1Hz dataset), but when I plot
>> the 0.1Hz dataset the components that I removed are still in the trace (ie.
>> eye blinks remain in the trace). An example section of the script I'm using
>> is listed below.
>>
>> I'm applying the ICAs to epoched data rather than continuous data, but I
>> can't see that that would be the problem. Am I missing a required step?
>> Like perhaps removing the artifact components?
>>
>> Load 1Hz filtered data set after ICA has been conducted and bad
>> components have been removed, then:
>>
>> >>* TMP.icawinv = EEG.icawinv;
>> *>>* TMP.icasphere = EEG.icasphere;
>> *>>* TMP.icaweights = EEG.icaweights;
>> *>>* TMP.icachansind = EEG.icachansind;
>> *>>>>* % apply to epoched dataset
>> *>>* clear EEG;
>> *>>* EEG = pop_loadset('filename', [SUBJ{s}, '_bcgrem.set'], 'filepath', PATHIN);
>> *>>* EEG.icawinv = TMP.icawinv;
>> *>>* EEG.icasphere = TMP.icasphere;
>> *>>* EEG.icaweights = TMP.icaweights;
>> *>>* EEG.icachansind = TMP.icachansind;
>> *>>* clear TMP;
>> *>>* EEG = pop_saveset(EEG, 'filename',[SUBJ{s}, '_epoched.set'],
>> *>>* 'filepath',PATHOUT);
>> *>>* [ALLEEG, EEG, CURRENTSET] = eeg_store(ALLEEG, EEG, 0);
>> *>>>>* end
>> *>>* eeglab redraw;*
>>
>> Kind regards,
>>
>> --
>> Neil Bailey (PhD)
>> Monash Alfred Psychiatry Research Centre
>> Neil.Bailey at Monash.edu
>> (03) 9076 5032
>>
>> _______________________________________________
>> Eeglablist page: http://sccn.ucsd.edu/eeglab/eeglabmail.html
>> To unsubscribe, send an empty email to
>> eeglablist-unsubscribe at sccn.ucsd.edu
>> For digest mode, send an email with the subject "set digest mime" to
>> eeglablist-request at sccn.ucsd.edu
>>
>
>
>
> --
> Makoto Miyakoshi
> Swartz Center for Computational Neuroscience
> Institute for Neural Computation, University of California San Diego
>
> _______________________________________________
> eeglablist mailing list eeglablist at sccn.ucsd.edu
> Eeglablist page: http://www.sccn.ucsd.edu/eeglab/eeglabmail.html
> To unsubscribe, send an empty email to eeglablist-unsub at sccn.ucsd.edu
> To switch to non-digest mode, send an empty email to
> eeglablist-nodigest at sccn.ucsd.edu
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://sccn.ucsd.edu/pipermail/eeglablist/attachments/20150822/44bb81eb/attachment.html>