pip install pymvpa2 python import numpy as np import mvpa2.suite as mvpa python dataset = mvpa.Dataset.from_hdf5('data.hdf5') python clf = mvpa.LinearCSVMC() fsel = mvpa.RFE(clf, 5) fsel.train(dataset) selected_features = fsel.get_mask() python tst_data = dataset[selected_features].copy() tst_data.sa['targets'] = dataset.sa['targets'] tst_data.sa['chunks'] = dataset.sa['chunks'] cv = mvpa.CrossValidation(clf, mvpa.NFoldPartitioner()) error = mvpa.mean(cv(tst_data)) python sensitivity = mvpa.map2nifti( cv.get_sensitivity_analyzer()(dataset)) sensitivity.to_filename('sensitivity_map.nii.gz') python import numpy as np import mvpa2.suite as mvpa dataset = mvpa.Dataset.from_hdf5('data.hdf5') clf = mvpa.LinearCSVMC() fsel = mvpa.RFE(clf, 5) fsel.train(dataset) selected_features = fsel.get_mask() tst_data = dataset[selected_features].copy() tst_data.sa['targets'] = dataset.sa['targets'] tst_data.sa['chunks'] = dataset.sa['chunks'] cv = mvpa.CrossValidation(clf, mvpa.NFoldPartitioner()) error = mvpa.mean(cv(tst_data)) sensitivity = mvpa.map2nifti( cv.get_sensitivity_analyzer()(dataset)) sensitivity.to_filename('sensitivity_map.nii.gz')