Beauchamp:Retinotopy: Difference between revisions

From OpenWetWare
Jump to navigationJump to search
Line 50: Line 50:
     3dDFT -prefix {$ec}r{$r}vrdtblFFT -detrend junkr{$r}+orig
     3dDFT -prefix {$ec}r{$r}vrdtblFFT -detrend junkr{$r}+orig
   end
   end
The retinotopic stimulus has a periodicity at frequency 5, so calculate phase at only that frequency.  Then, calculate amplitude squared (power) at all frequencies.
The retinotopic stimulus has a periodicity at frequency 5, so calculate phase at only that frequency for all runs.  Then, calculate amplitude squared (power) at all frequencies for all runs.
   foreach r (1 2 3 4 5 6 7 8)
   foreach r (1 2 3 4 5 6 7 8)
     3dcalc -datum short -prefix {$ec}r{$r}vrdtblFFTf5 -cx2r PHASE -a {$ec}r{$r}vrdtblFFT+orig'[5]'  -expr "a"
     3dcalc -datum short -prefix {$ec}r{$r}vrdtblFFTf5 -cx2r PHASE -a {$ec}r{$r}vrdtblFFT+orig'[5]'  -expr "a"
     3drefit -sublabel 0 FFTr{$r}f5 -subrepkey 0 FFTr{$r}f5 {$ec}r{$r}vrdtblFFTf5+orig
     3drefit -sublabel 0 FFTr{$r}f5 -subrepkey 0 FFTr{$r}f5 {$ec}r{$r}vrdtblFFTf5+orig
 
    3dcalc -datum short -prefix {$ec}r{$r}vrdtblFFTMAG -cx2r ABS -a {$ec}r{$r}vrdtblFFT+orig  -expr "a*a"
  3dcalc -datum short -prefix {$ec}r{$r}vrdtblFFTMAG -cx2r ABS -a {$ec}r{$r}vrdtblFFT+orig  -expr "a*a"
   end
   end


Navigation menu