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

[Neil Bailey]

Dear Makoto

Thanks for the advice. After advice from someone on the list I've tried
loading the ICA weights onto the 0.1Hz data and then removing the
artifactual components from there. It seems to work well.

Kind regards,

Neil

On 22 August 2015 at 09:43, Makoto Miyakoshi <mmiyakoshi at ucsd.edu> wrote:

> 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
>



-- 
Neil Bailey (PhD)
Monash Alfred Psychiatry Research Centre
Neil.Bailey at Monash.edu
(03) 9076 5032
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://sccn.ucsd.edu/pipermail/eeglablist/attachments/20150824/51175fd3/attachment.html>