def pearson_correlation(datadir,
outdir,
usemask=True,
excludef='',
exclude_idx=-1):
slidir = datadir + os.path.sep + au.slices_str()
subjsfile = datadir + os.path.sep + au.subjects_str()
labelsfile = datadir + os.path.sep + au.labels_str()
lst = os.listdir(slidir)
n = au.count_match(lst, au.data_str() + '_' + au.slice_regex())
exclude_log = ''
if exclude_idx > -1:
exclude_log = ' excluding subject ' + str(exclude_idx)
au.log.info('Calculating correlation of ' + slidir + os.path.sep +
au.data_str() + '_' + au.slice_regex() + exclude_log)
for i in range(n):
slino = au.zeropad(i)
dataf = slidir + os.path.sep + au.data_str() + '_' + au.slice_str(
) + '_' + slino + au.ext_str()
maskf = slidir + os.path.sep + au.mask_str() + '_' + au.slice_str(
) + '_' + slino + au.ext_str()
outf = outdir + os.path.sep + au.pearson_str() + '_' + au.slice_str(
) + '_' + slino
if exclude_idx > -1:
outf += '_' + au.excluded_str() + str(exclude_idx) + au.ext_str()
else:
outf += au.ext_str()
if not os.path.isfile(dataf):
au.log.error('Could not find ' + dataf)
continue
if not usemask:
maskf = ''
try:
measure_pearson(dataf, labelsfile, outf, maskf, excludef,
exclude_idx)
except:
au.log.error(
'pearson_correlation: Error measuring correlation on ' + dataf)
au.log.error("Unexpected error: ", sys.exc_info()[0])
exit(1)
def group_stats (datadir, groupname, groupsize, outdir=''):
lst = os.listdir(datadir)
n = au.count_match(lst, groupname + 's_' + au.slice_regex() + au.ext_str())
if not outdir:
outdir = datadir
au.log.info ('Calculating stats from group ' + groupname + ' in ' + outdir)
for i in range(n):
slino = au.zeropad(i)
dataf = datadir + os.path.sep + groupname + 's_' + au.slice_str() + '_' + slino + au.ext_str()
volstats (dataf, groupname, groupsize, outdir)
def group_stats(datadir, groupname, groupsize, outdir=''):
lst = os.listdir(datadir)
n = au.count_match(lst, groupname + 's_' + au.slice_regex() + au.ext_str())
if not outdir:
outdir = datadir
au.log.info('Calculating stats from group ' + groupname + ' in ' + outdir)
for i in range(n):
slino = au.zeropad(i)
dataf = datadir + os.path.sep + groupname + 's_' + au.slice_str(
) + '_' + slino + au.ext_str()
volstats(dataf, groupname, groupsize, outdir)
def pearson_correlation (datadir, outdir, usemask=True, excludef='', exclude_idx=-1):
slidir = datadir + os.path.sep + au.slices_str()
subjsfile = datadir + os.path.sep + au.subjects_str()
labelsfile = datadir + os.path.sep + au.labels_str()
lst = os.listdir(slidir)
n = au.count_match(lst, au.data_str() + '_' + au.slice_regex())
exclude_log = ''
if exclude_idx > -1:
exclude_log = ' excluding subject ' + str(exclude_idx)
au.log.info ('Calculating correlation of ' + slidir + os.path.sep + au.data_str() + '_' + au.slice_regex() + exclude_log)
for i in range(n):
slino = au.zeropad(i)
dataf = slidir + os.path.sep + au.data_str() + '_' + au.slice_str() + '_' + slino + au.ext_str()
maskf = slidir + os.path.sep + au.mask_str() + '_' + au.slice_str() + '_' + slino + au.ext_str()
outf = outdir + os.path.sep + au.pearson_str() + '_' + au.slice_str() + '_' + slino
if exclude_idx > -1:
outf += '_' + au.excluded_str() + str(exclude_idx) + au.ext_str()
else:
outf += au.ext_str()
if not os.path.isfile(dataf):
au.log.error('Could not find ' + dataf)
continue
if not usemask:
maskf = ''
try:
measure_pearson(dataf, labelsfile, outf, maskf, excludef, exclude_idx)
except:
au.log.error('pearson_correlation: Error measuring correlation on ' + dataf)
au.log.error("Unexpected error: ", sys.exc_info()[0] )
exit(1)
