161
edits
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" | |||
end | end | ||
edits