[pymvpa] question about detrend and zscore
John Clithero
john.clithero at gmail.com
Mon Oct 5 15:05:34 UTC 2009
Hi all,
I have a naive question about some of the preprocessing steps in PyMVPA.
I am loading in my data and detrending as follows (similar to some
examples listed):
##Load Data##
dataset = NiftiDataset((wb_file),
labels=attr.labels,
chunks=attr.chunks,
mask=os.path.join(roidir,'wb.nii.gz'))
##Detrend Data##
detrend(dataset, perchunk=True, model='linear')
zscore(dataset, targetdtype='float32')
I have two types of trials (A and B from the labels).
If I plot the average voxel value of A trials versus B trials, I get a
perfectly negatively correlated line.
In other words, if mean(sample voxel on A trials) = .5, then
mean(sample voxel on B trials) = -.5. This is true for all voxels.
I have looked over miscfx.py, but I thought would send an email to see
if (1) this is what "should happen" and (2) if so, what the idea is
for making such a split before running a classifier.
Best,
John
More information about the Pkg-ExpPsy-PyMVPA
mailing list