def run_perm_test(row):
row_array = row.split(';')
perm_ind = row_array[0]
perm_y = row_array[1]
perm_y = np.asarray(perm_y.split(',')).astype(int)
networks = [
'visual', 'somatomotor', 'dorsal_attention', 'ventral_attention',
'limbic', 'default_mode', 'fronto_parietal'
]
for network in networks:
analysis_name = 'schaefer_{}_noTiv_perm_'.format(network) + str(
perm_ind)
data_folder = '/scratch/tmp/wintern/iq_frankfurt/'
project_folder = '/scratch/tmp/wintern/iq_frankfurt/results/noTivRescaling/perm/schaefer_{}'.format(
network)
os.makedirs(project_folder, exist_ok=True)
cache_dir = '/scratch/tmp/wintern/cache'
# get data
data = IQData(data_folder=data_folder, tiv_rescaled=False)
covariates = np.asarray([data.age, data.gender, data.handedness]).T
data.load_schaefer_network(network, use_cached=True)
# run analysis
pipe = construct_hyperpipe_schaefer(analysis_name, project_folder,
cache_dir)
pipe.groups = data.fsiq
pipe.fit(data.schaefer_network, perm_y, **{'covariates': covariates})
os.remove(pipe.output_settings.pretrained_model_filename)
def run_perm_test(row):
row_array = row.split(';')
perm_ind = row_array[0]
perm_y = row_array[1]
perm_y = np.asarray(perm_y.split(',')).astype(int)
analysis_name = 'schaefer_perm_' + str(perm_ind)
data_folder = '/scratch/tmp/wintern/iq_frankfurt/'
project_folder = '/scratch/tmp/wintern/iq_frankfurt/results/perm/schaefer/'
cache_dir = '/scratch/tmp/wintern/cache'
# get data
data = IQData(data_folder=data_folder)
covariates = np.asarray([data.age, data.gender, data.handedness]).T
data.load_schaefer(use_cached=True)
# run analysis
pipe = construct_hyperpipe_schaefer(analysis_name, project_folder, cache_dir)
pipe.groups = data.fsiq
pipe.fit(data.schaefer, perm_y, **{'covariates': covariates})
os.remove(pipe.output_settings.pretrained_model_filename)
Last updated: 18-02-2019 Author ------ Nils R. Winter [email protected] Translationale Psychiatrie Universitaetsklinikum Muenster """ import sys sys.path.append('/scratch/tmp/wintern/iq_frankfurt/photonai') sys.path.append('/scratch/tmp/wintern/iq_frankfurt/') from analyses.analysis_base import construct_hyperpipe_schaefer from data.data import IQData import numpy as np import os analysis_name = 'shen_no_tiv' data_folder = '/scratch/tmp/wintern/iq_frankfurt/' project_folder = '/scratch/tmp/wintern/iq_frankfurt/results/noTivRescaling/shen' cache_dir = '/scratch/tmp/wintern/cache' # get data data = IQData(data_folder=data_folder, tiv_rescaled=False) covariates = np.asarray([data.age, data.gender, data.handedness]).T data.load_shen(use_cached=False) # run analysis pipe = construct_hyperpipe_schaefer(analysis_name, project_folder, cache_dir) pipe.fit(data.shen, data.fsiq, **{'covariates': covariates}) os.remove(pipe.output_settings.pretrained_model_filename)
