<div dir="ltr">Thanks for the response Steve. Thanks for catching that—it isn't exactly 1000. I ran the code you gave me and now this pops up. Thank you so much for the help with this.<div><br></div><div><div style="font-size:12.8000001907349px">pop_epoch():64 epochs selected</div><div style="font-size:12.8000001907349px">Epoching...</div><div style="font-size:12.8000001907349px">pop_epoch():64 epochs generated</div><div style="font-size:12.8000001907349px">pop_epoch(): checking epochs for data discontinuity</div><div style="font-size:12.8000001907349px">pop_rmbase(): Removing baseline...</div><div style="font-size:12.8000001907349px">Warning: Out of range or non-integer values truncated during conversion to character.</div></div><div style="font-size:12.8000001907349px"><br></div><div style="font-size:12.8000001907349px">Then I get:</div><div style="font-size:12.8000001907349px"><br></div><div style="font-size:12.8000001907349px"><u>Error using horzcat</u></div><div style="font-size:12.8000001907349px">CAT arguments dimensions are not consistent.</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Sep 5, 2015 at 4:03 PM, Stephen Politzer-Ahles <span dir="ltr"><<a href="mailto:stephen.politzer-ahles@ling-phil.ox.ac.uk" target="_blank">stephen.politzer-ahles@ling-phil.ox.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div><div><div>Hello Tara,<br><br></div>It's hard to tell without seeing the data, but this sort of thing might happen when the timepoints in your data are not exactly e.g. 20 and 100 ms (which might be the case if your sampling rate was not 1000 Hz). Is the issue fixed if you replace the line in question with the following?<br><br>STLResponse = squeeze(max(EEG.data(1,find(EEG.times>=20 & EEG.times<=100),:)));<br><br></div>Best,<br></div>Steve<br></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Sat, Sep 5, 2015 at 8:17 PM, Tara Miskovich <span dir="ltr"><<a href="mailto:miskovi2@uwm.edu" target="_blank">miskovi2@uwm.edu</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr">Hello all,<div><br></div><div>I am having difficulty with a script I obtained for processing startle data in EEGlab:<br clear="all"><div><br></div><div><div>%% High pass filter</div><div> [HP28_1,HP28_2] = butter(2,28/(EEG.srate/2),'high'); % generate filter coefficient for 4th order 28Hz High pass Butterworth filter</div><div> EEG.data = filtfilt(HP28_1,HP28_2,EEG.data); %apply that filter to ORB channel; filter order is divided by 2 because we use the filtfilt function to avoid phase shifts in the data</div><div><br></div><div>%% Rectify </div><div> EEG.data = abs(EEG.data); %the 'abs' function rectifies the data</div><div> </div><div>%% Low pass filter </div><div> [LP30_1,LP30_2] = butter(2,30/(EEG.srate/2),'low'); % generate filter coefficient for 4th order 30Hz Low pass Butterworth filter</div><div> EEG.data = filtfilt(LP30_1,LP30_2,EEG.data); %%apply that filter to ORB channel; filter order is divided by 2 because we use the filtfilt function to avoid phase shifts in the data</div><div><br></div><div>%% Epoch and baseline correct</div><div> EventCodes = num2cell(unique([EEG.event.type])); %Extract event codes for epoching </div><div> EEG = pop_epoch( EEG, EventCodes, [-0.05 0.25], 'epochinfo', 'yes'); %Epoch file for events that indicate startle probes</div><div> EEG = pop_rmbase( EEG, [-50 0]); %Baseline correct</div><div><br></div><div>%% Mark and reject trials with excessive artificat </div><div> Artifact = squeeze(max(EEG.data(1,1:find(EEG.times==0),:)) - min( EEG.data(1,1:find(EEG.times==0),:)) > 40); %Identify trials with greater than 40 microvolt deflections in baseline</div><div> EEGNA = pop_select(EEG,'notrial', find(Artifact)); %Reject these trials </div><div> </div><div>%% Score startle peaks </div><div> STLResponse = squeeze(max(EEG.data(1,find(EEG.times==20):find(EEG.times==100),:))); %Identify peak response of startle trials in scoring window</div><div> STLResponse = [[EEG.event.type]' STLResponse]; %Store peak values by event type </div></div><div><br></div><div>I get the following output:</div><div><br></div><div>pop_epoch():64 epochs selected</div><div>Epoching...</div><div>pop_epoch():64 epochs generated</div><div>pop_epoch(): checking epochs for data discontinuity</div><div>pop_rmbase(): Removing baseline...</div><div>Warning: Concatenation involves an empty array with an incorrect number of rows.</div><div>This may not be allowed in a future release.</div><div><br></div><div>I am not sure what is the issue here, but after "STLResponse = squeeze(max(EEG.data(1,find(EEG.times==20):find(EEG.times==100),:))); %Identify peak response of startle trials in scoring window" I end up with an empty 0 by 64 matrix. </div><div><br></div><div>I am just a little new to this and can't seem to figure out what is going on here.</div><div><br></div><div>Thank you!</div><span><font color="#888888"><div>Tara</div><div><br></div>
</font></span></div></div>
<br></div></div>_______________________________________________<br>
Eeglablist page: <a href="http://sccn.ucsd.edu/eeglab/eeglabmail.html" rel="noreferrer" target="_blank">http://sccn.ucsd.edu/eeglab/eeglabmail.html</a><br>
To unsubscribe, send an empty email to <a href="mailto:eeglablist-unsubscribe@sccn.ucsd.edu" target="_blank">eeglablist-unsubscribe@sccn.ucsd.edu</a><br>
For digest mode, send an email with the subject "set digest mime" to <a href="mailto:eeglablist-request@sccn.ucsd.edu" target="_blank">eeglablist-request@sccn.ucsd.edu</a><br></blockquote></div><br></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><font style="font-size:small;font-family:arial,helvetica,sans-serif">Thank you,<br><br></font><font size="2"><span style="font-family:arial,helvetica,sans-serif;color:rgb(33,33,33)">Tara A. Miskovich, M.S.</span><br style="font-family:arial,helvetica,sans-serif;color:rgb(33,33,33)"><span style="font-family:arial,helvetica,sans-serif;color:rgb(33,33,33)">Affective Neuroscience Laboratory</span><br style="font-family:arial,helvetica,sans-serif;color:rgb(33,33,33)"><span style="font-family:arial,helvetica,sans-serif;color:rgb(33,33,33)">Department of Psychology</span><br style="font-family:arial,helvetica,sans-serif;color:rgb(33,33,33)"><span style="font-family:arial,helvetica,sans-serif;color:rgb(33,33,33)">University of Wisconsin-Milwaukee</span><br style="font-family:arial,helvetica,sans-serif;color:rgb(33,33,33)"><span style="font-family:arial,helvetica,sans-serif;color:rgb(33,33,33)">334 Garland Hall</span><br style="font-family:arial,helvetica,sans-serif;color:rgb(33,33,33)"><span style="font-family:arial,helvetica,sans-serif;color:rgb(33,33,33)"><a href="mailto:miskovi2@uwm.edu" target="_blank">miskovi2@uwm.edu</a></span></font><br></div></div></div></div>
</div>