def test_load(tmpdir):
sim = Simulator()
sigma = 1
y = [0, 1]
n_reps = 3
output_dir = str(tmpdir)
dat = sim.create_data(y, sigma, reps=n_reps, output_dir=output_dir)
# if MNI_Template["resolution"] == '2mm':
# shape_3d = (91, 109, 91)
# shape_2d = (6, 238955)
# elif MNI_Template["resolution"] == '3mm':
# shape_3d = (60, 72, 60)
# shape_2d = (6, 71020)
y = pd.read_csv(os.path.join(str(tmpdir.join('y.csv'))), header=None, index_col=None)
# holdout = pd.read_csv(os.path.join(str(tmpdir.join('rep_id.csv'))), header=None, index_col=None)
# Test load list of 4D images
file_list = [str(tmpdir.join('data.nii.gz')), str(tmpdir.join('data.nii.gz'))]
dat = Brain_Data(file_list)
dat = Brain_Data([nb.load(x) for x in file_list])
# Test load list
dat = Brain_Data(data=str(tmpdir.join('data.nii.gz')), Y=y)
# Test Write
dat.write(os.path.join(str(tmpdir.join('test_write.nii'))))
assert Brain_Data(os.path.join(str(tmpdir.join('test_write.nii'))))
def test_roc(tmpdir):
sim = Simulator()
sigma = .1
y = [0, 1]
n_reps = 10
# output_dir = str(tmpdir)
dat = sim.create_data(y, sigma, reps=n_reps, output_dir=None)
# dat = Brain_Data(data=sim.data, Y=sim.y)
algorithm = 'svm'
# output_dir = str(tmpdir)
# cv = {'type': 'kfolds', 'n_folds': 5, 'subject_id': sim.rep_id}
extra = {'kernel': 'linear'}
output = dat.predict(algorithm=algorithm, plot=False, **extra)
# Single-Interval
roc = Roc(input_values=output['yfit_all'], binary_outcome=output['Y'] == 1)
roc.calculate()
roc.summary()
assert roc.accuracy == 1
# Forced Choice
binary_outcome = output['Y'] == 1
forced_choice = list(range(int(len(binary_outcome)/2))) + list(range(int(len(binary_outcome)/2)))
forced_choice = forced_choice.sort()
roc_fc = Roc(input_values=output['yfit_all'], binary_outcome=binary_outcome, forced_choice=forced_choice)
roc_fc.calculate()
assert roc_fc.accuracy == 1
assert roc_fc.accuracy == roc_fc.auc == roc_fc.sensitivity == roc_fc.specificity
def test_load(tmpdir):
sim = Simulator()
sigma = 1
y = [0, 1]
n_reps = 3
output_dir = str(tmpdir)
dat = sim.create_data(y, sigma, reps=n_reps, output_dir=output_dir)
# if MNI_Template["resolution"] == '2mm':
# shape_3d = (91, 109, 91)
# shape_2d = (6, 238955)
# elif MNI_Template["resolution"] == '3mm':
# shape_3d = (60, 72, 60)
# shape_2d = (6, 71020)
y = pd.read_csv(os.path.join(str(tmpdir.join("y.csv"))),
header=None,
index_col=None)
# holdout = pd.read_csv(os.path.join(str(tmpdir.join('rep_id.csv'))), header=None, index_col=None)
# Test load list of 4D images
file_list = [
str(tmpdir.join("data.nii.gz")),
str(tmpdir.join("data.nii.gz"))
]
dat = Brain_Data(file_list)
dat = Brain_Data([nb.load(x) for x in file_list])
# Test load string and path
dat = Brain_Data(data=str(tmpdir.join("data.nii.gz")), Y=y)
dat = Brain_Data(data=Path(tmpdir.join("data.nii.gz")), Y=y)
# Test Write
dat.write(os.path.join(str(tmpdir.join("test_write.nii"))))
assert Brain_Data(os.path.join(str(tmpdir.join("test_write.nii"))))
# Test i/o for hdf5
dat.write(os.path.join(str(tmpdir.join("test_write.h5"))))
b = Brain_Data(os.path.join(tmpdir.join("test_write.h5")))
for k in ["X", "Y", "mask", "nifti_masker", "file_name", "data"]:
if k == "data":
assert np.allclose(b.__dict__[k], dat.__dict__[k])
elif k in ["X", "Y"]:
assert all(b.__dict__[k].eq(dat.__dict__[k]).values)
elif k == "mask":
assert np.allclose(b.__dict__[k].affine, dat.__dict__[k].affine)
assert np.allclose(b.__dict__[k].get_fdata(),
dat.__dict__[k].get_fdata())
assert b.__dict__[k].get_filename(
) == dat.__dict__[k].get_filename()
elif k == "nifti_masker":
assert np.allclose(b.__dict__[k].affine_, dat.__dict__[k].affine_)
assert np.allclose(b.__dict__[k].mask_img.get_fdata(),
dat.__dict__[k].mask_img.get_fdata())
else:
assert b.__dict__[k] == dat.__dict__[k]
def test_simulator(tmpdir):
sim = Simulator()
r = 10
sigma = 1
y = [0, 1]
n_reps = 3
output_dir = str(tmpdir)
shape = (91, 109, 91)
dat = sim.create_data(y, sigma, reps=n_reps, output_dir=None)
assert len(dat) == n_reps * len(y)
assert len(dat.Y) == n_reps * len(y)
def sim_brain_data():
# MNI_Template["resolution"] = request.params
sim = Simulator()
r = 10
sigma = 1
y = [0, 1]
n_reps = 3
dat = sim.create_data(y, sigma, reps=n_reps)
dat.X = pd.DataFrame({'Intercept':np.ones(len(dat.Y)),
'X1':np.array(dat.Y).flatten()}, index=None)
return dat
def test_simulator(tmpdir):
sim = Simulator()
r = 10
sigma = 1
y = [0, 1]
n_reps = 3
output_dir = str(tmpdir)
sim.create_data(y, sigma, reps=n_reps, output_dir=output_dir)
flist = glob.glob(str(tmpdir.join("centered*nii.gz")))
shape = (91, 109, 91)
sim_img = nb.concat_images(flist)
assert len(sim.data) == n_reps * len(y)
assert sim_img.shape[0:3] == shape
def test_find_spikes():
sim = Simulator()
y = [0, 1]
n_reps = 50
s1 = create_sphere([0, 0, 0], radius=3)
d1 = sim.create_data(y, 1, reps=n_reps, output_dir=None).apply_mask(s1)
spikes = find_spikes(d1)
assert isinstance(spikes, pd.DataFrame)
assert spikes.shape[0] == len(d1)
spikes = find_spikes(d1.to_nifti())
assert isinstance(spikes, pd.DataFrame)
assert spikes.shape[0] == len(d1)
def test_simulator(tmpdir):
sim = Simulator()
r = 10
sigma = 1
y = [0, 1]
n_reps = 3
output_dir = str(tmpdir)
sim.create_data(y, sigma, reps=n_reps, output_dir=output_dir)
flist = glob.glob(str(tmpdir.join('centered*nii.gz')))
shape = (91, 109, 91)
sim_img = nb.concat_images(flist)
assert len(sim.data) == n_reps * len(y)
assert sim_img.shape[0:3] == shape
def test_groupby(tmpdir):
# Simulate Brain Data
sim = Simulator()
r = 10
sigma = 1
y = [0, 1]
n_reps = 3
output_dir = str(tmpdir)
sim.create_data(y, sigma, reps=n_reps, output_dir=output_dir)
s1 = create_sphere([12, 10, -8], radius=r)
s2 = create_sphere([22, -2, -22], radius=r)
mask = Brain_Data([s1, s2])
y = pd.read_csv(os.path.join(str(tmpdir.join('y.csv'))),
header=None,
index_col=None)
data = Brain_Data(glob.glob(str(tmpdir.join('data.nii.gz'))), Y=y)
data.X = pd.DataFrame(
{
'Intercept': np.ones(len(data.Y)),
'X1': np.array(data.Y).flatten()
},
index=None)
dat = Groupby(data, mask)
# Test length
assert len(dat) == len(mask)
# Test Index
assert isinstance(dat[1], Brain_Data)
# Test apply
mn = dat.apply('mean')
assert len(dat) == len(mn)
# assert mn[0].mean() > mn[1].mean() #JC edit: it seems this check relies on chance from simulated data
assert mn[1].shape() == np.sum(mask[1].data == 1)
reg = dat.apply('regress')
assert len(dat) == len(mn)
# r = dict([(x,reg[x]['beta'][1]) for x in reg.iterkeys()])
# Test combine
combine_mn = dat.combine(mn)
assert len(combine_mn.shape()) == 1
def test_data(tmpdir):
sim = Simulator()
r = 10
sigma = 1
y = [0, 1]
n_reps = 3
output_dir = str(tmpdir)
sim.create_data(y, sigma, reps=n_reps, output_dir=output_dir)
shape_3d = (91, 109, 91)
shape_2d = (6, 238955)
y=pd.read_csv(os.path.join(str(tmpdir.join('y.csv'))),header=None,index_col=None).T
flist = glob.glob(str(tmpdir.join('centered*.nii.gz')))
dat = Brain_Data(data=flist,Y=y)
# Test shape
assert dat.shape() == shape_2d
# Test Mean
assert dat.mean().shape()[0] == shape_2d[1]
# Test Std
assert dat.std().shape()[0] == shape_2d[1]
# Test to_nifti
d = dat.to_nifti
assert d().shape[0:3] == shape_3d
# # Test T-test
out = dat.ttest()
assert out['t'].shape()[0]==shape_2d[1]
# Test Regress
dat.X = pd.DataFrame({'Intercept':np.ones(len(dat.Y)),'X1':np.array(dat.Y).flatten()},index=None)
out = dat.regress()
assert out['beta'].shape() == (2,shape_2d[1])
# Test indexing
assert out['t'][1].shape()[0] == shape_2d[1]
# Test threshold
i=1
tt = threshold(out['t'][i], out['p'][i], threshold_dict={'fdr':.05})
assert tt.shape()[0] == shape_2d[1]
def test_predict_multi():
# Simulate data 100 images worth
sim = Simulator()
sigma = 1
y = [0, 1]
n_reps = 50
output_dir = '.'
dat = sim.create_data(y, sigma, reps=n_reps, output_dir=output_dir)
y = pd.read_csv('y.csv', header=None, index_col=None)
dat = Brain_Data('data.nii.gz', Y=y)
# Predict within given ROIs
# Generate some "rois" (in reality non-contiguous, but also not overlapping)
roi_1 = dat[0].copy()
roi_1.data = np.zeros_like(roi_1.data, dtype=bool)
roi_2 = roi_1.copy()
roi_3 = roi_1.copy()
idx = np.random.choice(range(roi_1.shape()[-1]), size=9999, replace=False)
roi_1.data[idx[:3333]] = 1
roi_2.data[idx[3333:6666]] = 1
roi_3.data[idx[6666:]] = 1
rois = roi_1.append(roi_2).append(roi_3)
# Load in all 50 rois so we can "insert" signal into the first one
# rois = expand_mask(Brain_Data(os.path.join(get_resource_path(), 'k50.nii.gz')))
# roi = rois[0]
from sklearn.datasets import make_classification
X, Y = make_classification(n_samples=100, n_features=rois[0].data.sum(), n_informative=500, n_redundant=5, n_classes=2)
dat.data[:, rois[0].data.astype(bool)] = X
dat.Y = pd.Series(Y)
out = dat.predict_multi(algorithm='svm', cv_dict={'type': 'kfolds', 'n_folds': 3}, method='rois', n_jobs=-1, rois=rois[:3], kernel='linear')
assert len(out) == 3
assert np.sum([elem['weight_map'].data.shape for elem in out]) == rois.data.sum()
# Searchlight
roi_mask = rois[:2].sum()
out = dat.predict_multi(algorithm='svm', cv_dict={'type': 'kfolds', 'n_folds': 3}, method='searchlight', radius=4, verbose=50, n_jobs=-1, process_mask=roi_mask)
assert len(np.nonzero(out.data)[0]) == len(np.nonzero(roi_mask.data)[0])
def test_brain_data(tmpdir):
# Add 3mm to list to test that resolution as well
for resolution in ['2mm']:
MNI_Template["resolution"] = resolution
sim = Simulator()
r = 10
sigma = 1
y = [0, 1]
n_reps = 3
output_dir = str(tmpdir)
dat = sim.create_data(y, sigma, reps=n_reps, output_dir=output_dir)
if MNI_Template["resolution"] == '2mm':
shape_3d = (91, 109, 91)
shape_2d = (6, 238955)
elif MNI_Template["resolution"] == '3mm':
shape_3d = (60, 72, 60)
shape_2d = (6, 71020)
y = pd.read_csv(os.path.join(str(tmpdir.join('y.csv'))),header=None, index_col=None)
holdout = pd.read_csv(os.path.join(str(tmpdir.join('rep_id.csv'))),header=None,index_col=None)
# Test load list of 4D images
file_list = [str(tmpdir.join('data.nii.gz')), str(tmpdir.join('data.nii.gz'))]
dat = Brain_Data(file_list)
dat = Brain_Data([nb.load(x) for x in file_list])
# Test load list
dat = Brain_Data(data=str(tmpdir.join('data.nii.gz')), Y=y)
# Test concatenate
out = Brain_Data([x for x in dat])
assert isinstance(out, Brain_Data)
assert len(out)==len(dat)
# Test to_nifti
d = dat.to_nifti()
assert d.shape[0:3] == shape_3d
# Test load nibabel
assert Brain_Data(d)
# Test shape
assert dat.shape() == shape_2d
# Test Mean
assert dat.mean().shape()[0] == shape_2d[1]
# Test Std
assert dat.std().shape()[0] == shape_2d[1]
# Test add
new = dat + dat
assert new.shape() == shape_2d
# Test subtract
new = dat - dat
assert new.shape() == shape_2d
# Test multiply
new = dat * dat
assert new.shape() == shape_2d
# Test Indexing
index = [0, 3, 1]
assert len(dat[index]) == len(index)
index = range(4)
assert len(dat[index]) == len(index)
index = dat.Y == 1
assert len(dat[index.values.flatten()]) == index.values.sum()
assert len(dat[index]) == index.values.sum()
assert len(dat[:3]) == 3
# Test Iterator
x = [x for x in dat]
assert len(x) == len(dat)
assert len(x[0].data.shape) == 1
# # Test T-test
out = dat.ttest()
assert out['t'].shape()[0] == shape_2d[1]
# # # Test T-test - permutation method
# out = dat.ttest(threshold_dict={'permutation':'tfce','n_permutations':50,'n_jobs':1})
# assert out['t'].shape()[0]==shape_2d[1]
# Test Regress
dat.X = pd.DataFrame({'Intercept':np.ones(len(dat.Y)),
'X1':np.array(dat.Y).flatten()}, index=None)
# Standard OLS
out = dat.regress()
assert type(out['beta'].data) == np.ndarray
assert type(out['t'].data) == np.ndarray
assert type(out['p'].data) == np.ndarray
assert type(out['residual'].data) == np.ndarray
assert type(out['df'].data) == np.ndarray
assert out['beta'].shape() == (2, shape_2d[1])
assert out['t'][1].shape()[0] == shape_2d[1]
# Robust OLS
out = dat.regress(mode='robust')
assert type(out['beta'].data) == np.ndarray
assert type(out['t'].data) == np.ndarray
assert type(out['p'].data) == np.ndarray
assert type(out['residual'].data) == np.ndarray
assert type(out['df'].data) == np.ndarray
assert out['beta'].shape() == (2, shape_2d[1])
assert out['t'][1].shape()[0] == shape_2d[1]
# Test threshold
i=1
tt = threshold(out['t'][i], out['p'][i], .05)
assert isinstance(tt, Brain_Data)
# Test write
dat.write(os.path.join(str(tmpdir.join('test_write.nii'))))
assert Brain_Data(os.path.join(str(tmpdir.join('test_write.nii'))))
# Test append
assert dat.append(dat).shape()[0] == shape_2d[0]*2
# Test distance
distance = dat.distance(method='euclidean')
assert isinstance(distance, Adjacency)
assert distance.square_shape()[0] == shape_2d[0]
# Test predict
stats = dat.predict(algorithm='svm',
cv_dict={'type': 'kfolds', 'n_folds': 2},
plot=False, **{'kernel':"linear"})
# Support Vector Regression, with 5 fold cross-validation with Platt Scaling
# This will output probabilities of each class
stats = dat.predict(algorithm='svm',
cv_dict=None, plot=False,
**{'kernel':'linear', 'probability':True})
assert isinstance(stats['weight_map'], Brain_Data)
# Logistic classificiation, with 2 fold cross-validation.
stats = dat.predict(algorithm='logistic',
cv_dict={'type': 'kfolds', 'n_folds': 2},
plot=False)
assert isinstance(stats['weight_map'], Brain_Data)
# Ridge classificiation,
stats = dat.predict(algorithm='ridgeClassifier', cv_dict=None, plot=False)
assert isinstance(stats['weight_map'], Brain_Data)
# Ridge
stats = dat.predict(algorithm='ridge',
cv_dict={'type': 'kfolds', 'n_folds': 2,
'subject_id':holdout}, plot=False, **{'alpha':.1})
# Lasso
stats = dat.predict(algorithm='lasso',
cv_dict={'type': 'kfolds', 'n_folds': 2,
'stratified':dat.Y}, plot=False, **{'alpha':.1})
# PCR
stats = dat.predict(algorithm='pcr', cv_dict=None, plot=False)
# Test Similarity
r = dat.similarity(stats['weight_map'])
assert len(r) == shape_2d[0]
r2 = dat.similarity(stats['weight_map'].to_nifti())
assert len(r2) == shape_2d[0]
r = dat.similarity(stats['weight_map'], method='dot_product')
assert len(r) == shape_2d[0]
r = dat.similarity(stats['weight_map'], method='cosine')
assert len(r) == shape_2d[0]
r = dat.similarity(dat, method='correlation')
assert r.shape == (dat.shape()[0],dat.shape()[0])
r = dat.similarity(dat, method='dot_product')
assert r.shape == (dat.shape()[0],dat.shape()[0])
r = dat.similarity(dat, method='cosine')
assert r.shape == (dat.shape()[0],dat.shape()[0])
# Test apply_mask - might move part of this to test mask suite
s1 = create_sphere([12, 10, -8], radius=10)
assert isinstance(s1, nb.Nifti1Image)
masked_dat = dat.apply_mask(s1)
assert masked_dat.shape()[1] == np.sum(s1.get_data() != 0)
# Test extract_roi
mask = create_sphere([12, 10, -8], radius=10)
assert len(dat.extract_roi(mask)) == shape_2d[0]
# Test r_to_z
z = dat.r_to_z()
assert z.shape() == dat.shape()
# Test copy
d_copy = dat.copy()
assert d_copy.shape() == dat.shape()
# Test detrend
detrend = dat.detrend()
assert detrend.shape() == dat.shape()
# Test standardize
s = dat.standardize()
assert s.shape() == dat.shape()
assert np.isclose(np.sum(s.mean().data), 0, atol=.1)
s = dat.standardize(method='zscore')
assert s.shape() == dat.shape()
assert np.isclose(np.sum(s.mean().data), 0, atol=.1)
# Test Sum
s = dat.sum()
assert s.shape() == dat[1].shape()
# Test Groupby
s1 = create_sphere([12, 10, -8], radius=10)
s2 = create_sphere([22, -2, -22], radius=10)
mask = Brain_Data([s1, s2])
d = dat.groupby(mask)
assert isinstance(d, Groupby)
# Test Aggregate
mn = dat.aggregate(mask, 'mean')
assert isinstance(mn, Brain_Data)
assert len(mn.shape()) == 1
# Test Threshold
s1 = create_sphere([12, 10, -8], radius=10)
s2 = create_sphere([22, -2, -22], radius=10)
mask = Brain_Data(s1)*5
mask = mask + Brain_Data(s2)
m1 = mask.threshold(upper=.5)
m2 = mask.threshold(upper=3)
m3 = mask.threshold(upper='98%')
m4 = Brain_Data(s1)*5 + Brain_Data(s2)*-.5
m4 = mask.threshold(upper=.5,lower=-.3)
assert np.sum(m1.data > 0) > np.sum(m2.data > 0)
assert np.sum(m1.data > 0) == np.sum(m3.data > 0)
assert np.sum(m4.data[(m4.data > -.3) & (m4.data <.5)]) == 0
assert np.sum(m4.data[(m4.data < -.3) | (m4.data >.5)]) > 0
# Test Regions
r = mask.regions(min_region_size=10)
m1 = Brain_Data(s1)
m2 = r.threshold(1, binarize=True)
# assert len(r)==2
assert len(np.unique(r.to_nifti().get_data())) == 2
diff = m2-m1
assert np.sum(diff.data) == 0
# Test Bootstrap
masked = dat.apply_mask(create_sphere(radius=10, coordinates=[0, 0, 0]))
n_samples = 3
b = masked.bootstrap('mean', n_samples=n_samples)
assert isinstance(b['Z'], Brain_Data)
b = masked.bootstrap('std', n_samples=n_samples)
assert isinstance(b['Z'], Brain_Data)
b = masked.bootstrap('predict', n_samples=n_samples, plot=False)
assert isinstance(b['Z'], Brain_Data)
b = masked.bootstrap('predict', n_samples=n_samples,
plot=False, cv_dict={'type':'kfolds','n_folds':3})
assert isinstance(b['Z'], Brain_Data)
b = masked.bootstrap('predict', n_samples=n_samples,
save_weights=True, plot=False)
assert len(b['samples'])==n_samples
# Test decompose
n_components = 3
stats = dat.decompose(algorithm='pca', axis='voxels',
n_components=n_components)
assert n_components == len(stats['components'])
assert stats['weights'].shape == (len(dat), n_components)
stats = dat.decompose(algorithm='ica', axis='voxels',
n_components=n_components)
assert n_components == len(stats['components'])
assert stats['weights'].shape == (len(dat), n_components)
dat.data = dat.data + 2
dat.data[dat.data<0] = 0
stats = dat.decompose(algorithm='nnmf', axis='voxels',
n_components=n_components)
assert n_components == len(stats['components'])
assert stats['weights'].shape == (len(dat), n_components)
stats = dat.decompose(algorithm='fa', axis='voxels',
n_components=n_components)
assert n_components == len(stats['components'])
assert stats['weights'].shape == (len(dat), n_components)
stats = dat.decompose(algorithm='pca', axis='images',
n_components=n_components)
assert n_components == len(stats['components'])
assert stats['weights'].shape == (len(dat), n_components)
stats = dat.decompose(algorithm='ica', axis='images',
n_components=n_components)
assert n_components == len(stats['components'])
assert stats['weights'].shape == (len(dat), n_components)
dat.data = dat.data + 2
dat.data[dat.data<0] = 0
stats = dat.decompose(algorithm='nnmf', axis='images',
n_components=n_components)
assert n_components == len(stats['components'])
assert stats['weights'].shape == (len(dat), n_components)
stats = dat.decompose(algorithm='fa', axis='images',
n_components=n_components)
assert n_components == len(stats['components'])
assert stats['weights'].shape == (len(dat), n_components)
# Test Hyperalignment Method
sim = Simulator()
y = [0, 1]
n_reps = 10
s1 = create_sphere([0, 0, 0], radius=3)
d1 = sim.create_data(y, 1, reps=n_reps, output_dir=None).apply_mask(s1)
d2 = sim.create_data(y, 2, reps=n_reps, output_dir=None).apply_mask(s1)
d3 = sim.create_data(y, 3, reps=n_reps, output_dir=None).apply_mask(s1)
# Test procrustes using align
data = [d1, d2, d3]
out = align(data, method='procrustes')
assert len(data) == len(out['transformed'])
assert len(data) == len(out['transformation_matrix'])
assert data[0].shape() == out['common_model'].shape()
transformed = np.dot(d1.data, out['transformation_matrix'][0])
centered = d1.data - np.mean(d1.data, 0)
transformed = (np.dot(centered/np.linalg.norm(centered), out['transformation_matrix'][0])*out['scale'][0])
np.testing.assert_almost_equal(0, np.sum(out['transformed'][0].data - transformed), decimal=5)
# Test deterministic brain_data
bout = d1.align(out['common_model'], method='deterministic_srm')
assert d1.shape() == bout['transformed'].shape()
assert d1.shape() == bout['common_model'].shape()
assert d1.shape()[1] == bout['transformation_matrix'].shape[0]
btransformed = np.dot(d1.data, bout['transformation_matrix'])
np.testing.assert_almost_equal(0, np.sum(bout['transformed'].data - btransformed))
# Test deterministic brain_data
bout = d1.align(out['common_model'], method='probabilistic_srm')
assert d1.shape() == bout['transformed'].shape()
assert d1.shape() == bout['common_model'].shape()
assert d1.shape()[1] == bout['transformation_matrix'].shape[0]
btransformed = np.dot(d1.data, bout['transformation_matrix'])
np.testing.assert_almost_equal(0, np.sum(bout['transformed'].data-btransformed))
# Test procrustes brain_data
bout = d1.align(out['common_model'], method='procrustes')
assert d1.shape() == bout['transformed'].shape()
assert d1.shape() == bout['common_model'].shape()
assert d1.shape()[1] == bout['transformation_matrix'].shape[0]
centered = d1.data - np.mean(d1.data, 0)
btransformed = (np.dot(centered/np.linalg.norm(centered), bout['transformation_matrix'])*bout['scale'])
np.testing.assert_almost_equal(0, np.sum(bout['transformed'].data-btransformed), decimal=5)
np.testing.assert_almost_equal(0, np.sum(out['transformed'][0].data - bout['transformed'].data))
# Test hyperalignment on Brain_Data over time (axis=1)
sim = Simulator()
y = [0, 1]
n_reps = 10
s1 = create_sphere([0, 0, 0], radius=5)
d1 = sim.create_data(y, 1, reps=n_reps, output_dir=None).apply_mask(s1)
d2 = sim.create_data(y, 2, reps=n_reps, output_dir=None).apply_mask(s1)
d3 = sim.create_data(y, 3, reps=n_reps, output_dir=None).apply_mask(s1)
data = [d1, d2, d3]
out = align(data, method='procrustes', axis=1)
assert len(data) == len(out['transformed'])
assert len(data) == len(out['transformation_matrix'])
assert data[0].shape() == out['common_model'].shape()
centered = data[0].data.T-np.mean(data[0].data.T, 0)
transformed = (np.dot(centered/np.linalg.norm(centered), out['transformation_matrix'][0])*out['scale'][0])
np.testing.assert_almost_equal(0,np.sum(out['transformed'][0].data-transformed.T), decimal=5)
bout = d1.align(out['common_model'], method='deterministic_srm', axis=1)
assert d1.shape() == bout['transformed'].shape()
assert d1.shape() == bout['common_model'].shape()
assert d1.shape()[0] == bout['transformation_matrix'].shape[0]
btransformed = np.dot(d1.data.T, bout['transformation_matrix'])
np.testing.assert_almost_equal(0, np.sum(bout['transformed'].data-btransformed.T))
bout = d1.align(out['common_model'], method='probabilistic_srm', axis=1)
assert d1.shape() == bout['transformed'].shape()
assert d1.shape() == bout['common_model'].shape()
assert d1.shape()[0] == bout['transformation_matrix'].shape[0]
btransformed = np.dot(d1.data.T, bout['transformation_matrix'])
np.testing.assert_almost_equal(0, np.sum(bout['transformed'].data-btransformed.T))
bout = d1.align(out['common_model'], method='procrustes', axis=1)
assert d1.shape() == bout['transformed'].shape()
assert d1.shape() == bout['common_model'].shape()
assert d1.shape()[0] == bout['transformation_matrix'].shape[0]
centered = d1.data.T-np.mean(d1.data.T, 0)
btransformed = (np.dot(centered/np.linalg.norm(centered), bout['transformation_matrix'])*bout['scale'])
np.testing.assert_almost_equal(0, np.sum(bout['transformed'].data-btransformed.T), decimal=5)
np.testing.assert_almost_equal(0, np.sum(out['transformed'][0].data-bout['transformed'].data))
def test_align():
# Test hyperalignment matrix
sim = Simulator()
y = [0, 1]
n_reps = 10
s1 = create_sphere([0, 0, 0], radius=3)
d1 = sim.create_data(y, 1, reps=n_reps, output_dir=None).apply_mask(s1)
d2 = sim.create_data(y, 2, reps=n_reps, output_dir=None).apply_mask(s1)
d3 = sim.create_data(y, 3, reps=n_reps, output_dir=None).apply_mask(s1)
data = [d1.data, d2.data, d3.data]
out = align(data, method="deterministic_srm")
assert len(data) == len(out["transformed"])
assert len(data) == len(out["transformation_matrix"])
assert data[0].shape == out["common_model"].shape
transformed = np.dot(data[0], out["transformation_matrix"][0])
np.testing.assert_almost_equal(np.sum(out["transformed"][0] -
transformed.T),
0,
decimal=3)
assert len(out["isc"]) == out["transformed"][0].shape[0]
out = align(data, method="probabilistic_srm")
assert len(data) == len(out["transformed"])
assert len(data) == len(out["transformation_matrix"])
assert data[0].shape == out["common_model"].shape
transformed = np.dot(data[0], out["transformation_matrix"][0])
np.testing.assert_almost_equal(np.sum(out["transformed"][0] -
transformed.T),
0,
decimal=3)
assert len(out["isc"]) == out["transformed"][0].shape[0]
out2 = align(data, method="procrustes")
assert len(data) == len(out2["transformed"])
assert data[0].shape == out2["common_model"].shape
assert len(data) == len(out2["transformation_matrix"])
assert len(data) == len(out2["disparity"])
centered = data[0] - np.mean(data[0], 0)
transformed = (np.dot(centered / np.linalg.norm(centered),
out2["transformation_matrix"][0]) * out2["scale"][0])
np.testing.assert_almost_equal(np.sum(out2["transformed"][0] -
transformed.T),
0,
decimal=3)
assert out2["transformed"][0].shape == out2["transformed"][0].shape
assert (out2["transformation_matrix"][0].shape ==
out2["transformation_matrix"][0].shape)
assert len(out2["isc"]) == out["transformed"][0].shape[0]
# Test hyperalignment on Brain_Data
data = [d1, d2, d3]
out = align(data, method="deterministic_srm")
assert len(data) == len(out["transformed"])
assert len(data) == len(out["transformation_matrix"])
assert data[0].shape() == out["common_model"].shape
transformed = np.dot(d1.data, out["transformation_matrix"][0].data.T)
np.testing.assert_almost_equal(np.sum(out["transformed"][0].data -
transformed),
0,
decimal=3)
assert len(out["isc"]) == out["transformed"][0].shape[1]
out = align(data, method="probabilistic_srm")
assert len(data) == len(out["transformed"])
assert len(data) == len(out["transformation_matrix"])
assert data[0].shape() == out["common_model"].shape
transformed = np.dot(d1.data, out["transformation_matrix"][0].data.T)
np.testing.assert_almost_equal(np.sum(out["transformed"][0].data -
transformed),
0,
decimal=3)
assert len(out["isc"]) == out["transformed"][0].shape[1]
out2 = align(data, method="procrustes")
assert len(data) == len(out2["transformed"])
assert data[0].shape() == out2["common_model"].shape()
assert len(data) == len(out2["transformation_matrix"])
assert len(data) == len(out2["disparity"])
centered = data[0].data - np.mean(data[0].data, 0)
transformed = (np.dot(centered / np.linalg.norm(centered),
out2["transformation_matrix"][0].data) *
out2["scale"][0])
np.testing.assert_almost_equal(np.sum(out2["transformed"][0].data -
transformed),
0,
decimal=3)
assert out2["transformed"][0].shape() == out2["transformed"][0].shape()
assert (out2["transformation_matrix"][0].shape ==
out2["transformation_matrix"][0].shape)
assert len(out2["isc"]) == out2["transformed"][0].shape()[1]
# Test hyperalignment on matrix over time (axis=1)
sim = Simulator()
y = [0, 1]
n_reps = 10
s1 = create_sphere([0, 0, 0], radius=5)
d1 = sim.create_data(y, 1, reps=n_reps, output_dir=None).apply_mask(s1)
d2 = sim.create_data(y, 2, reps=n_reps, output_dir=None).apply_mask(s1)
d3 = sim.create_data(y, 3, reps=n_reps, output_dir=None).apply_mask(s1)
data = [d1.data, d2.data, d3.data]
out = align(data, method="deterministic_srm", axis=1)
assert len(data) == len(out["transformed"])
assert len(data) == len(out["transformation_matrix"])
assert data[0].shape == out["common_model"].shape
transformed = np.dot(data[0].T, out["transformation_matrix"][0].data)
np.testing.assert_almost_equal(np.sum(out["transformed"][0] - transformed),
0,
decimal=3)
assert len(out["isc"]) == out["transformed"][0].shape[1]
out = align(data, method="probabilistic_srm", axis=1)
assert len(data) == len(out["transformed"])
assert len(data) == len(out["transformation_matrix"])
assert data[0].shape == out["common_model"].shape
transformed = np.dot(data[0].T, out["transformation_matrix"][0])
np.testing.assert_almost_equal(np.sum(out["transformed"][0] - transformed),
0,
decimal=3)
assert len(out["isc"]) == out["transformed"][0].shape[1]
out2 = align(data, method="procrustes", axis=1)
assert len(data) == len(out2["transformed"])
assert data[0].shape == out2["common_model"].shape
assert len(data) == len(out2["transformation_matrix"])
assert len(data) == len(out2["disparity"])
centered = data[0] - np.mean(data[0], 0)
transformed = (np.dot(
(centered / np.linalg.norm(centered)).T,
out2["transformation_matrix"][0].data,
) * out2["scale"][0])
np.testing.assert_almost_equal(np.sum(out2["transformed"][0] -
transformed),
0,
decimal=3)
assert out2["transformed"][0].shape == out2["transformed"][0].shape
assert (out2["transformation_matrix"][0].shape ==
out2["transformation_matrix"][0].shape)
assert len(out2["isc"]) == out2["transformed"][0].shape[0]
# Test hyperalignment on Brain_Data over time (axis=1)
data = [d1, d2, d3]
out = align(data, method="deterministic_srm", axis=1)
assert len(data) == len(out["transformed"])
assert len(data) == len(out["transformation_matrix"])
assert data[0].shape() == out["common_model"].shape
transformed = np.dot(d1.data.T, out["transformation_matrix"][0].data).T
np.testing.assert_almost_equal(np.sum(out["transformed"][0].data -
transformed),
0,
decimal=5)
assert len(out["isc"]) == out["transformed"][0].shape[0]
out = align(data, method="probabilistic_srm", axis=1)
assert len(data) == len(out["transformed"])
assert len(data) == len(out["transformation_matrix"])
assert data[0].shape() == out["common_model"].shape
transformed = np.dot(d1.data.T, out["transformation_matrix"][0].data).T
np.testing.assert_almost_equal(np.sum(out["transformed"][0].data -
transformed),
0,
decimal=5)
assert len(out["isc"]) == out["transformed"][0].shape[0]
out2 = align(data, method="procrustes", axis=1)
assert len(data) == len(out2["transformed"])
assert data[0].shape() == out2["common_model"].shape()
assert len(data) == len(out2["transformation_matrix"])
assert len(data) == len(out2["disparity"])
centered = data[0].data.T - np.mean(data[0].data.T, 0)
transformed = (np.dot(centered / np.linalg.norm(centered),
out2["transformation_matrix"][0].data) *
out2["scale"][0]).T
np.testing.assert_almost_equal(np.sum(out2["transformed"][0].data -
transformed),
0,
decimal=5)
assert out2["transformed"][0].shape() == out2["transformed"][0].shape()
assert (out2["transformation_matrix"][0].shape ==
out2["transformation_matrix"][0].shape)
assert len(out2["isc"]) == out2["transformed"][0].shape()[1]
def test_hyperalignment():
sim = Simulator()
y = [0, 1]
n_reps = 10
s1 = create_sphere([0, 0, 0], radius=3)
d1 = sim.create_data(y, 1, reps=n_reps, output_dir=None).apply_mask(s1)
d2 = sim.create_data(y, 2, reps=n_reps, output_dir=None).apply_mask(s1)
d3 = sim.create_data(y, 3, reps=n_reps, output_dir=None).apply_mask(s1)
data = [d1, d2, d3]
# Test deterministic brain_data
out = align(data, method="deterministic_srm")
bout = d1.align(out["common_model"], method="deterministic_srm")
assert d1.shape() == bout["transformed"].shape
assert d1.shape() == bout["common_model"].shape
assert d1.shape()[1] == bout["transformation_matrix"].shape()[0]
btransformed = np.dot(d1.data, bout["transformation_matrix"].data.T)
np.testing.assert_almost_equal(
0, np.sum(bout["transformed"].data - btransformed))
# Test probabilistic brain_data
bout = d1.align(out["common_model"], method="probabilistic_srm")
assert d1.shape() == bout["transformed"].shape
assert d1.shape() == bout["common_model"].shape
assert d1.shape()[1] == bout["transformation_matrix"].shape()[0]
btransformed = np.dot(d1.data, bout["transformation_matrix"].data.T)
np.testing.assert_almost_equal(
0, np.sum(bout["transformed"].data - btransformed))
# Test procrustes brain_data
out = align(data, method="procrustes")
centered = data[0].data - np.mean(data[0].data, 0)
transformed = (np.dot(centered / np.linalg.norm(centered),
out["transformation_matrix"][0].data) *
out["scale"][0])
bout = d1.align(out["common_model"], method="procrustes")
assert d1.shape() == bout["transformed"].shape()
assert d1.shape() == bout["common_model"].shape()
assert d1.shape()[1] == bout["transformation_matrix"].shape()[0]
centered = d1.data - np.mean(d1.data, 0)
btransformed = (np.dot(centered / np.linalg.norm(centered),
bout["transformation_matrix"].data) * bout["scale"])
np.testing.assert_almost_equal(0,
np.sum(bout["transformed"].data -
btransformed),
decimal=5)
np.testing.assert_almost_equal(
0, np.sum(out["transformed"][0].data - bout["transformed"].data))
# Test over time
sim = Simulator()
y = [0, 1]
n_reps = 10
s1 = create_sphere([0, 0, 0], radius=5)
d1 = sim.create_data(y, 1, reps=n_reps, output_dir=None).apply_mask(s1)
d2 = sim.create_data(y, 2, reps=n_reps, output_dir=None).apply_mask(s1)
d3 = sim.create_data(y, 3, reps=n_reps, output_dir=None).apply_mask(s1)
data = [d1, d2, d3]
out = align(data, method="deterministic_srm", axis=1)
bout = d1.align(out["common_model"], method="deterministic_srm", axis=1)
assert d1.shape() == bout["transformed"].shape
assert d1.shape() == bout["common_model"].shape
assert d1.shape()[0] == bout["transformation_matrix"].shape()[0]
btransformed = np.dot(d1.data.T, bout["transformation_matrix"].data.T)
np.testing.assert_almost_equal(
0, np.sum(bout["transformed"].data - btransformed.T))
out = align(data, method="probabilistic_srm", axis=1)
bout = d1.align(out["common_model"], method="probabilistic_srm", axis=1)
assert d1.shape() == bout["transformed"].shape
assert d1.shape() == bout["common_model"].shape
assert d1.shape()[0] == bout["transformation_matrix"].shape()[0]
btransformed = np.dot(d1.data.T, bout["transformation_matrix"].data.T)
np.testing.assert_almost_equal(
0, np.sum(bout["transformed"].data - btransformed.T))
out = align(data, method="procrustes", axis=1)
bout = d1.align(out["common_model"], method="procrustes", axis=1)
assert d1.shape() == bout["transformed"].shape()
assert d1.shape() == bout["common_model"].shape()
assert d1.shape()[0] == bout["transformation_matrix"].shape()[0]
centered = d1.data.T - np.mean(d1.data.T, 0)
btransformed = (np.dot(centered / np.linalg.norm(centered),
bout["transformation_matrix"].data) * bout["scale"])
np.testing.assert_almost_equal(0,
np.sum(bout["transformed"].data -
btransformed.T),
decimal=5)
np.testing.assert_almost_equal(
0, np.sum(out["transformed"][0].data - bout["transformed"].data))
def test_align():
# Test hyperalignment matrix
sim = Simulator()
y = [0, 1]
n_reps = 10
s1 = create_sphere([0, 0, 0], radius=3)
d1 = sim.create_data(y, 1, reps=n_reps, output_dir=None).apply_mask(s1)
d2 = sim.create_data(y, 2, reps=n_reps, output_dir=None).apply_mask(s1)
d3 = sim.create_data(y, 3, reps=n_reps, output_dir=None).apply_mask(s1)
data = [d1.data.T, d2.data.T, d3.data.T]
out = align(data, method='deterministic_srm')
assert len(data) == len(out['transformed'])
assert len(data) == len(out['transformation_matrix'])
assert data[0].shape == out['common_model'].shape
transformed = np.dot(data[0].T, out['transformation_matrix'][0])
np.testing.assert_almost_equal(
0, np.sum(out['transformed'][0] - transformed.T))
out = align(data, method='probabilistic_srm')
assert len(data) == len(out['transformed'])
assert len(data) == len(out['transformation_matrix'])
assert data[0].shape == out['common_model'].shape
transformed = np.dot(data[0].T, out['transformation_matrix'][0])
np.testing.assert_almost_equal(
0, np.sum(out['transformed'][0] - transformed.T))
out2 = align(data, method='procrustes')
assert len(data) == len(out2['transformed'])
assert data[0].shape == out2['common_model'].shape
assert len(data) == len(out2['transformation_matrix'])
assert len(data) == len(out2['disparity'])
centered = data[0].T - np.mean(data[0].T, 0)
transformed = (np.dot(centered / np.linalg.norm(centered),
out2['transformation_matrix'][0]) * out2['scale'][0])
np.testing.assert_almost_equal(
0, np.sum(out2['transformed'][0] - transformed.T))
assert out['transformed'][0].shape == out2['transformed'][0].shape
assert out['transformation_matrix'][0].shape == out2[
'transformation_matrix'][0].shape
# Test hyperalignment on Brain_Data
data = [d1, d2, d3]
out = align(data, method='deterministic_srm')
assert len(data) == len(out['transformed'])
assert len(data) == len(out['transformation_matrix'])
assert data[0].shape() == out['common_model'].shape()
transformed = np.dot(d1.data, out['transformation_matrix'][0])
np.testing.assert_almost_equal(
0, np.sum(out['transformed'][0].data - transformed))
out = align(data, method='probabilistic_srm')
assert len(data) == len(out['transformed'])
assert len(data) == len(out['transformation_matrix'])
assert data[0].shape() == out['common_model'].shape()
transformed = np.dot(d1.data, out['transformation_matrix'][0])
np.testing.assert_almost_equal(
0, np.sum(out['transformed'][0].data - transformed))
out2 = align(data, method='procrustes')
assert len(data) == len(out2['transformed'])
assert data[0].shape() == out2['common_model'].shape()
assert len(data) == len(out2['transformation_matrix'])
assert len(data) == len(out2['disparity'])
centered = data[0].data - np.mean(data[0].data, 0)
transformed = (np.dot(centered / np.linalg.norm(centered),
out2['transformation_matrix'][0]) * out2['scale'][0])
np.testing.assert_almost_equal(
0, np.sum(out2['transformed'][0].data - transformed))
assert out['transformed'][0].shape() == out2['transformed'][0].shape()
assert out['transformation_matrix'][0].shape == out2[
'transformation_matrix'][0].shape
# Test hyperalignment on matrix over time (axis=1)
sim = Simulator()
y = [0, 1]
n_reps = 10
s1 = create_sphere([0, 0, 0], radius=5)
d1 = sim.create_data(y, 1, reps=n_reps, output_dir=None).apply_mask(s1)
d2 = sim.create_data(y, 2, reps=n_reps, output_dir=None).apply_mask(s1)
d3 = sim.create_data(y, 3, reps=n_reps, output_dir=None).apply_mask(s1)
data = [d1.data.T, d2.data.T, d3.data.T]
out = align(data, method='deterministic_srm', axis=1)
assert len(data) == len(out['transformed'])
assert len(data) == len(out['transformation_matrix'])
assert data[0].shape == out['common_model'].shape
transformed = np.dot(data[0], out['transformation_matrix'][0])
np.testing.assert_almost_equal(0,
np.sum(out['transformed'][0] - transformed))
out = align(data, method='probabilistic_srm', axis=1)
assert len(data) == len(out['transformed'])
assert len(data) == len(out['transformation_matrix'])
assert data[0].shape == out['common_model'].shape
transformed = np.dot(data[0], out['transformation_matrix'][0])
np.testing.assert_almost_equal(0,
np.sum(out['transformed'][0] - transformed))
out2 = align(data, method='procrustes', axis=1)
assert len(data) == len(out2['transformed'])
assert data[0].shape == out2['common_model'].shape
assert len(data) == len(out2['transformation_matrix'])
assert len(data) == len(out2['disparity'])
centered = data[0] - np.mean(data[0], 0)
transformed = (np.dot(centered / np.linalg.norm(centered),
out2['transformation_matrix'][0]) * out2['scale'][0])
np.testing.assert_almost_equal(
0, np.sum(out2['transformed'][0] - transformed))
assert out['transformed'][0].shape == out2['transformed'][0].shape
assert out['transformation_matrix'][0].shape == out2[
'transformation_matrix'][0].shape
# Test hyperalignment on Brain_Data over time (axis=1)
data = [d1, d2, d3]
out = align(data, method='deterministic_srm', axis=1)
assert len(data) == len(out['transformed'])
assert len(data) == len(out['transformation_matrix'])
assert data[0].shape() == out['common_model'].shape()
transformed = np.dot(d1.data.T, out['transformation_matrix'][0])
np.testing.assert_almost_equal(
0, np.sum(out['transformed'][0].data - transformed.T))
out = align(data, method='probabilistic_srm', axis=1)
assert len(data) == len(out['transformed'])
assert len(data) == len(out['transformation_matrix'])
assert data[0].shape() == out['common_model'].shape()
transformed = np.dot(d1.data.T, out['transformation_matrix'][0])
np.testing.assert_almost_equal(
0, np.sum(out['transformed'][0].data - transformed.T))
out2 = align(data, method='procrustes', axis=1)
assert len(data) == len(out2['transformed'])
assert data[0].shape() == out2['common_model'].shape()
assert len(data) == len(out2['transformation_matrix'])
assert len(data) == len(out2['disparity'])
centered = data[0].data.T - np.mean(data[0].data.T, 0)
transformed = (np.dot(centered / np.linalg.norm(centered),
out2['transformation_matrix'][0]) * out2['scale'][0])
np.testing.assert_almost_equal(
0, np.sum(out2['transformed'][0].data - transformed.T))
assert out['transformed'][0].shape() == out2['transformed'][0].shape()
assert out['transformation_matrix'][0].shape == out2[
'transformation_matrix'][0].shape
def test_brain_data(tmpdir):
sim = Simulator()
r = 10
sigma = 1
y = [0, 1]
n_reps = 3
output_dir = str(tmpdir)
sim.create_data(y, sigma, reps=n_reps, output_dir=output_dir)
shape_3d = (91, 109, 91)
shape_2d = (6, 238955)
y=pd.read_csv(os.path.join(str(tmpdir.join('y.csv'))), header=None,index_col=None).T
holdout=pd.read_csv(os.path.join(str(tmpdir.join('rep_id.csv'))),header=None,index_col=None).T
flist = glob.glob(str(tmpdir.join('centered*.nii.gz')))
# Test load list
dat = Brain_Data(data=flist,Y=y)
# Test load file
assert Brain_Data(flist[0])
# Test to_nifti
d = dat.to_nifti()
assert d.shape[0:3] == shape_3d
# Test load nibabel
assert Brain_Data(d)
# Test shape
assert dat.shape() == shape_2d
# Test Mean
assert dat.mean().shape()[0] == shape_2d[1]
# Test Std
assert dat.std().shape()[0] == shape_2d[1]
# Test add
new = dat + dat
assert new.shape() == shape_2d
# Test subtract
new = dat - dat
assert new.shape() == shape_2d
# Test multiply
new = dat * dat
assert new.shape() == shape_2d
# Test Iterator
x = [x for x in dat]
assert len(x) == len(dat)
assert len(x[0].data.shape) == 1
# # Test T-test
out = dat.ttest()
assert out['t'].shape()[0] == shape_2d[1]
# # # Test T-test - permutation method
# out = dat.ttest(threshold_dict={'permutation':'tfce','n_permutations':50,'n_jobs':1})
# assert out['t'].shape()[0]==shape_2d[1]
# Test Regress
dat.X = pd.DataFrame({'Intercept':np.ones(len(dat.Y)), 'X1':np.array(dat.Y).flatten()},index=None)
out = dat.regress()
assert out['beta'].shape() == (2,shape_2d[1])
# Test indexing
assert out['t'][1].shape()[0] == shape_2d[1]
# Test threshold
i=1
tt = threshold(out['t'][i], out['p'][i], .05)
assert isinstance(tt,Brain_Data)
# Test write
dat.write(os.path.join(str(tmpdir.join('test_write.nii'))))
assert Brain_Data(os.path.join(str(tmpdir.join('test_write.nii'))))
# Test append
assert dat.append(dat).shape()[0]==shape_2d[0]*2
# Test distance
distance = dat.distance(method='euclidean')
assert isinstance(distance,Adjacency)
assert distance.square_shape()[0]==shape_2d[0]
# Test predict
stats = dat.predict(algorithm='svm', cv_dict={'type': 'kfolds','n_folds': 2}, plot=False,**{'kernel':"linear"})
# Support Vector Regression, with 5 fold cross-validation with Platt Scaling
# This will output probabilities of each class
stats = dat.predict(algorithm='svm', cv_dict=None, plot=False,**{'kernel':'linear', 'probability':True})
assert isinstance(stats['weight_map'],Brain_Data)
# Logistic classificiation, with 2 fold cross-validation.
stats = dat.predict(algorithm='logistic', cv_dict={'type': 'kfolds', 'n_folds': 2}, plot=False)
assert isinstance(stats['weight_map'],Brain_Data)
# Ridge classificiation,
stats = dat.predict(algorithm='ridgeClassifier', cv_dict=None,plot=False)
assert isinstance(stats['weight_map'],Brain_Data)
# Ridge
stats = dat.predict(algorithm='ridge', cv_dict={'type': 'kfolds', 'n_folds': 2,'subject_id':holdout}, plot=False,**{'alpha':.1})
# Lasso
stats = dat.predict(algorithm='lasso', cv_dict={'type': 'kfolds', 'n_folds': 2,'stratified':dat.Y}, plot=False,**{'alpha':.1})
# PCR
stats = dat.predict(algorithm='pcr', cv_dict=None, plot=False)
# Test Similarity
r = dat.similarity(stats['weight_map'])
assert len(r) == shape_2d[0]
r2 = dat.similarity(stats['weight_map'].to_nifti())
assert len(r2) == shape_2d[0]
# Test apply_mask - might move part of this to test mask suite
s1 = create_sphere([12, 10, -8], radius=10)
assert isinstance(s1, nb.Nifti1Image)
s2 = Brain_Data(s1)
masked_dat = dat.apply_mask(s1)
assert masked_dat.shape()[1] == np.sum(s2.data != 0)
# Test extract_roi
mask = create_sphere([12, 10, -8], radius=10)
assert len(dat.extract_roi(mask)) == shape_2d[0]
# Test r_to_z
z = dat.r_to_z()
assert z.shape() == dat.shape()
# Test copy
d_copy = dat.copy()
assert d_copy.shape() == dat.shape()
# Test detrend
detrend = dat.detrend()
assert detrend.shape() == dat.shape()
# Test standardize
s = dat.standardize()
assert s.shape() == dat.shape()
assert np.isclose(np.sum(s.mean().data), 0, atol=.1)
s = dat.standardize(method='zscore')
assert s.shape() == dat.shape()
assert np.isclose(np.sum(s.mean().data), 0, atol=.1)
# Test Sum
s = dat.sum()
assert s.shape() == dat[1].shape()
# Test Groupby
s1 = create_sphere([12, 10, -8], radius=10)
s2 = create_sphere([22, -2, -22], radius=10)
mask = Brain_Data([s1, s2])
d = dat.groupby(mask)
assert isinstance(d, Groupby)
# Test Aggregate
mn = dat.aggregate(mask, 'mean')
assert isinstance(mn, Brain_Data)
assert len(mn.shape()) == 1
# Test Threshold
s1 = create_sphere([12, 10, -8], radius=10)
s2 = create_sphere([22, -2, -22], radius=10)
mask = Brain_Data(s1)*5
mask = mask + Brain_Data(s2)
m1 = mask.threshold(thresh=.5)
m2 = mask.threshold(thresh=3)
m3 = mask.threshold(thresh='98%')
assert np.sum(m1.data > 0) > np.sum(m2.data > 0)
assert np.sum(m1.data > 0) == np.sum(m3.data > 0)
# Test Regions
r = mask.regions(min_region_size=10)
m1 = Brain_Data(s1)
m2 = r.threshold(1, binarize=True)
# assert len(r)==2
assert len(np.unique(r.to_nifti().get_data())) == 2 # JC edit: I think this is what you were trying to do
diff = m2-m1
assert np.sum(diff.data) == 0
def test_data(tmpdir):
sim = Simulator()
r = 10
sigma = 1
y = [0, 1]
n_reps = 3
output_dir = str(tmpdir)
sim.create_data(y, sigma, reps=n_reps, output_dir=output_dir)
shape_3d = (91, 109, 91)
shape_2d = (6, 238955)
y = pd.read_csv(os.path.join(str(tmpdir.join('y.csv'))),
header=None,
index_col=None).T
flist = glob.glob(str(tmpdir.join('centered*.nii.gz')))
# Test load list
dat = Brain_Data(data=flist, Y=y)
# Test load file
assert Brain_Data(flist[0])
# Test to_nifti
d = dat.to_nifti()
assert d.shape[0:3] == shape_3d
# Test load nibabel
assert Brain_Data(d)
# Test shape
assert dat.shape() == shape_2d
# Test Mean
assert dat.mean().shape()[0] == shape_2d[1]
# Test Std
assert dat.std().shape()[0] == shape_2d[1]
# Test add
new = dat + dat
assert new.shape() == shape_2d
# Test subtract
new = dat - dat
assert new.shape() == shape_2d
# Test multiply
new = dat * dat
assert new.shape() == shape_2d
# Test Iterator
x = [x for x in dat]
assert len(x) == len(dat)
assert len(x[0].data.shape) == 1
# # Test T-test
out = dat.ttest()
assert out['t'].shape()[0] == shape_2d[1]
# # # Test T-test - permutation method
# out = dat.ttest(threshold_dict={'permutation':'tfce','n_permutations':50,'n_jobs':1})
# assert out['t'].shape()[0]==shape_2d[1]
# Test Regress
dat.X = pd.DataFrame(
{
'Intercept': np.ones(len(dat.Y)),
'X1': np.array(dat.Y).flatten()
},
index=None)
out = dat.regress()
assert out['beta'].shape() == (2, shape_2d[1])
# Test indexing
assert out['t'][1].shape()[0] == shape_2d[1]
# Test threshold
i = 1
tt = threshold(out['t'][i], out['p'][i], .05)
assert isinstance(tt, Brain_Data)
# Test write
dat.write(os.path.join(str(tmpdir.join('test_write.nii'))))
assert Brain_Data(os.path.join(str(tmpdir.join('test_write.nii'))))
# Test append
assert dat.append(dat).shape()[0] == shape_2d[0] * 2
# Test distance
distance = dat.distance(method='euclidean')
assert distance.shape == (shape_2d[0], shape_2d[0])
# Test predict
stats = dat.predict(algorithm='svm',
cv_dict={
'type': 'kfolds',
'n_folds': 2,
'n': len(dat.Y)
},
plot=False,
**{'kernel': "linear"})
# Support Vector Regression, with 5 fold cross-validation with Platt Scaling
# This will output probabilities of each class
stats = dat.predict(algorithm='svm',
cv_dict=None,
plot=False,
**{
'kernel': 'linear',
'probability': True
})
assert isinstance(stats['weight_map'], Brain_Data)
# Logistic classificiation, with 5 fold stratified cross-validation.
stats = dat.predict(algorithm='logistic',
cv_dict={
'type': 'kfolds',
'n_folds': 5,
'n': len(dat.Y)
},
plot=False)
assert isinstance(stats['weight_map'], Brain_Data)
# Ridge classificiation, with 5 fold between-subject cross-validation, where data for each subject is held out together.
stats = dat.predict(algorithm='ridgeClassifier', cv_dict=None, plot=False)
assert isinstance(stats['weight_map'], Brain_Data)
# Test Similarity
r = dat.similarity(stats['weight_map'])
assert len(r) == shape_2d[0]
r2 = dat.similarity(stats['weight_map'].to_nifti())
assert len(r2) == shape_2d[0]
# Test apply_mask - might move part of this to test mask suite
s1 = create_sphere([41, 64, 55], radius=10)
assert isinstance(s1, nb.Nifti1Image)
s2 = Brain_Data(s1)
masked_dat = dat.apply_mask(s1)
assert masked_dat.shape()[1] == np.sum(s2.data != 0)
# Test extract_roi
mask = create_sphere([41, 64, 55], radius=10)
assert len(dat.extract_roi(mask)) == shape_2d[0]
# Test r_to_z
z = dat.r_to_z()
assert z.shape() == dat.shape()
# Test copy
d_copy = dat.copy()
assert d_copy.shape() == dat.shape()
# Test detrend
detrend = dat.detrend()
assert detrend.shape() == dat.shape()
from nltools.analysis import Predict, Roc
from nltools.data import Brain_Data
from nltools.stats import threshold
from nltools.mask import create_sphere
import matplotlib.pyplot as plt
import shutil
import tempfile
tmp_dir = os.path.join(tempfile.gettempdir(), str(os.times()[-1]))
###############################################################################
# Create data
tic = time() #Start Timer
sim = Simulator()
r = 10
sigma = .5
cor = .8
cov = .6
n_trials = 10
n_subs = 5
s1 = create_sphere([41, 64, 55], radius=r)
sim.create_cov_data(cor,
cov,
sigma,
mask=s1,
reps=n_trials,
n_sub=n_subs,
output_dir=tmp_dir)
print 'Simulate Data: Elapsed: %.2f seconds' % (time() - tic) #Stop timer
