<div dir="ltr"><div>I'm trying to import and analyze some Neuromag 306 data, and I keep
running into the following error message whenever I try to detect blink
artifacts with the threshold method:<br><br><b>??? Error using ==> filter_with_correction at 44<br>
Calculated filter coefficients have poles on or outside the unit circle
and will not be stable. Try a higher cutoff frequency or a different
type/order of filter.<br><br>Error in ==> filter_with_correction at 44<br>
error('Calculated filter coefficients have poles on or outside the
unit circle and will not be stable. Try a higher cutoff frequency or a
different type/order of filter.');<br><br>Error in ==> ft_preproc_bandpassfilter at 148<br>
filt = filter_with_correction(B,A,dat,dir);<br><br>Error in ==> preproc at 299<br>if
strcmp(cfg.bpfilter, 'yes'), dat = ft_preproc_bandpassfilter(dat,
fsample, cfg.bpfreq, cfg.bpfiltord, cfg.bpfilttype, cfg.bpfiltdir, cfg.bpinstabilityfix); end<br><br>Error in ==> ft_artifact_threshold at 148<br> dat = preproc(dat, channel, offset2time(cfg.trl(trlop,3), hdr.Fs, size(dat,2)), artfctdef);<br>
</b><br></div>Based
on previous advice in these threads, I have tried changing the cutoff
frequency, order, direction, and type of filtering of my data using cfg
fields like bpfreq, bpfilttype, bpfiltdir, and bpinstabilityfix. Long
story short, I've tried every logically possible combination of filter
settings (including no filtering at all), and I get the same error
regardless. Does anyone have some insight into this? <br clear="all">
<br>Take care,<br>Tristan Davenport<br>UC San Diego</div>