How do I add channels to exclude while keeping the channels already excluded?

Hi all, I'm having trouble with this piece of code. I'm trying to exclude a whole bunch of channels from a data set, but whenever I try and exclude extras, it ends up undoing the exclusionary criteria. Here's the code:

%% Channel Exclusion

        figure;bar(ChanImpedances)
        thres = 4*1e6; % set the threshold for the bad channel
        hold on,plot(ones(1,length(ChanImpedances))*thres,'r--','LineWidth',2)


        bad_channel = [];
        bad_channel_o = [];
        bad_channel_o = find(ChanImpedances > thres)
        bad_channel = find(ChanImpedances > thres)

        bad_channel_o = [bad_channel_o;32;128];
        bad_channel = [bad_channel;32;128];

So it's already excluding a bunch of channels based upon impedance, but I also want to exclude others based upon visual criteria. Whenever I add:

        bad_channel = [46];

It undoes the bad channel selection it just ran. How do I tell it to also exclude channel 46 along with the ones it's already excluded for impedance?

Thank you!

Answer 1

[Star Strider]

If you want to add 46 to the bad channels, probably the easiest way is to ‘grow’ the vector:

bad_channel = [bad_channel; 46];

(This is not great programming practice. In some instances, it works best.)

[Alexandra Rowan O'Donnell]

Thank you! I actually ended up adding it to this line and it worked.

        bad_channel_o = [bad_channel_o;32;46;128];
        bad_channel = [bad_channel;32;46;128];

It wasn't running properly because bad_channel is recalled further down for interpolation and that's where it screws up.

[Star Strider]

As always, my pleasure!

If the order is important, efficiency is secondary, and you want them sorted ascending (lowest first), you can add them whenever you want, then use:

bad_channel = sort(bad_channel);