Dear Vladimir and Martin
Thank you so much. May I followup on a few more questions:
1. I want to exclude channels not involved in the files, which I already set to type "Other". I have set S.keepothers as no (according to manual P.111), but the new file still include those. However, if I keep those, I am not sure how they will affect the result (aside from keeping a file of bigger size) ? What is your advise on how to handle these "other" channels ?
2. For those marked as bad channels, as I understand, spm_eeg_montage will exclude calculation for these particular channels, however, will values for these also be excluded in the calculation of the column average, or we need to manually adjust the montage matrix ((1/length(R)) like what is done in Martin's code):
T(G,R) = T(G,R) - 1/length(R); where G/R are setdiff of meegchannels and badchannels
(Martin, unfortunately, when I run the above code, it returns the following error "Index exceeds matrix dimension". And G and R are the same matrix, is it how you suggest ? Sorry that I am not very familiar with the matrix operations, have I missed anything ? I would like to learn how the montage matrix should look like, if I skip some channels.)
3. Since each subject / session have different set of badchannels, should I generate a different montage matrix for each one ? How will this differ from interpolating values for the badchannels first and do the same re-reference for all subjects ?
Great thanks again for your kind advise!
Best regards
SP Ho
|