This data was acquired in 2002 or 2003.

Dimensions are 128x128x5x272 complex.
Use magnitude of first image as anatomical.
Remove first 3 of 272 from fMRI analysis.


Data Slice 5 Used in: 
-ROWE DB, Logan BR. A complex way to compute fMRI activation. Neuroimage 23(3):1078-1092 (2004).
-ROWE DB. Modeling both the magnitude and phase of complex-valued fMRI data. Neuroimage 25(4):1310-1324 (2005).
-ROWE DB, Meller CP, Hoffmann RG. Characterizing phase-only fMRI data with an angular regression model. J. Neurosci. Methods 161(2):331-341 (2007).
-Adrian DW, Maitra R, ROWE DB: Ricean versus Gaussian modelling in magnitude fMRI Analysis – Added Complexity with Few Practical Benefits. Stat, 2(1):303–316 (2013). 

Data description from RL04.
A bilateral sequential finger-tapping experiment was performed in a block design with 16-s off followed by eight epochs of 16-s on and 16-s off. Scanning was performed using a 1.5 T GE Signa in which five axial slices of size 96x96 wereacquired. In image reconstruction, the acquired data was zero filled to 128x128. After Fourier reconstruction, each voxel has dimensions in mm of 1.5625x1.5625x5, with TE = 47 ms. Observations were taken every TR = 1000 ms, so that there are 272 in each voxel. Data from a single axial slice through the motor cortex were selected for analysis. Preprocessing using an ideal 0/1 frequency filter (Gonzales and Woods, 1992; Press et al., 1992) was performed to remove respiration and low frequency physiological noise in addition to the removal of the first three points to omit machine warm-up effects.

Reference Waveform and Design Matrix:
q=2; r=1;
temp=[-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1]; % makes ref function
ref=[temp,temp,temp,temp,temp,temp,temp,temp]';
ref=[ref;-ones(16,1)]
X=[ones(n,1),(1:n)'-mean(1:n),ref];
trshift=5;
tempX=[-ones(trshift,q+1);X];
tempX(n+1:n+trshift,:)=[];
X(1:n,q+1)=tempX(1:n,q+1);