Hi Ziad,
Thanks for the response. I've succeeded in getting past the first error, but have now come across another (although they vary based on whether I set the Opt.Quiet setting to 1 or 0; see below).
Thanks for your help.
Brandon
>> RetroTS(Opt)
Opt =
PhysFS: 40
Nslices: 30
VolTR: 2
Respfile: '/Volumes/HD_2/RSVP/fMRI/Physio/P3_Run1_Resp.txt'
Cardfile: '/Volumes/HD_2/RSVP/fMRI/Physio/P3_Run1_Card.txt'
Quiet: 0
err: 1
zerophaseoffset: 0
SliceMajor: 1
RVTshifts: [0 5 10 15 20]
ResampFS: 40
RespCutoffFreq: 3
CardCutoffFreq: 3
ResamKernel: 'linear'
FIROrder: 40
Demo: 0
Prefix: 'oba'
Resp_out: 1
Card_out: 1
RVT_out: 1
SepDups: 0
SliceOffset: [30x1 double]
SliceOrder: 'alt+z'
Slice timing:
ans =
0
1.0000
0.0667
1.0667
0.1333
1.1333
0.2000
1.2000
0.2667
1.2667
0.3333
1.3333
0.4000
1.4000
0.4667
1.4667
0.5333
1.5333
0.6000
1.6000
0.6667
1.6667
0.7333
1.7333
0.8000
1.8000
0.8667
1.8667
0.9333
1.9333
File /Volumes/HD_2/RSVP/fMRI/Physio/P3_Run1_Resp.txt exists and will be read.
Reading file
--> Load signal
--> Smooth signal
--> Calculate analytic signal Z
--> Find zero crossing of imag(Z)
--> Improved peak location
--> Removed duplicates
--> Calculated the period (from beat to beat)
??? Index exceeds matrix dimensions.
Error in ==> plotsign2 at 78
if (CallPath(1)=='.') CallPath = pwd; end
Error in ==> PeakFinder at 281
if (~Opt.Quiet), plotsign2(1); end
Error in ==> RetroTS at 274
[R,e]= PeakFinder(Opt.Respfile,OptR);
>> Opt.Quiet = 1
Opt =
PhysFS: 40
Nslices: 30
VolTR: 2
Respfile: '/Volumes/HD_2/RSVP/fMRI/Physio/P3_Run1_Resp.txt'
Cardfile: '/Volumes/HD_2/RSVP/fMRI/Physio/P3_Run1_Card.txt'
Quiet: 1
>> RetroTS(Opt)
Opt =
PhysFS: 40
Nslices: 30
VolTR: 2
Respfile: '/Volumes/HD_2/RSVP/fMRI/Physio/P3_Run1_Resp.txt'
Cardfile: '/Volumes/HD_2/RSVP/fMRI/Physio/P3_Run1_Card.txt'
Quiet: 1
err: 1
zerophaseoffset: 0
SliceMajor: 1
RVTshifts: [0 5 10 15 20]
ResampFS: 40
RespCutoffFreq: 3
CardCutoffFreq: 3
ResamKernel: 'linear'
FIROrder: 40
Demo: 0
Prefix: 'oba'
Resp_out: 1
Card_out: 1
RVT_out: 1
SepDups: 0
SliceOffset: [30x1 double]
SliceOrder: 'alt+z'
File /Volumes/HD_2/RSVP/fMRI/Physio/P3_Run1_Resp.txt exists and will be read.
Reading file
File /Volumes/HD_2/RSVP/fMRI/Physio/P3_Run1_Card.txt exists and will be read.
Reading file
Estimating phase for R
Estimating phase for E
Computing RVT from peaks
Showing RVT Peaks for R
??? Index exceeds matrix dimensions.
Error in ==> plotsign2 at 78
if (CallPath(1)=='.') CallPath = pwd; end
Error in ==> Show_RVT_Peak at 58
plotsign2(fg);
Error in ==> RetroTS at 305
Show_RVT_Peak(R,1);