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 at all frequencies
  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"
    3drefit -sublabel 0 FFTr{$r}f5 -subrepkey 0 FFTr{$r}f5 {$ec}r{$r}vrdtblFFTf5+orig


  foreach r (1 2 3 4 5 6 7 8)
  # Retinotopic stimulus has periodicity at frequency 5, calculate phase at only that frequency
  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
  # Calculate magnitude squared at all frequencies
   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