def test_filter_series():
series = Query(proj_name).filter_series(series_name,
subjects=test_subject_name)
assert_equal(len(series), 1)
assert_equal(int(series[0]['serieno']), series_number)
def test_get_studies():
studies = Query(proj_name).get_studies(test_subject_name)
assert_equal(len(studies), len(study_modalities))
assert_equal(studies[study_name_tuple[0]], study_name_tuple[1])
if crop_plot_time is not None:
epochs.crop(*crop_plot_time)
epochs.plot_image(pick, fig=imgfig)
# mnefig[0].get_axes()[1].set_title('')
stats = dict()
stats['mean'] = np.mean(delays)
stats['std'] = np.std(delays)
stats['median'] = np.median(delays)
stats['q10'] = np.percentile(delays, 10.)
stats['q90'] = np.percentile(delays, 90.)
stats['max_amp'] = np.max(epochs._data[:, pick, :]) # ovr epochs & times
stats['min_amp'] = np.min(epochs._data[:, pick, :]) # ovr epochs & times
return (delays, stats)
if __name__ == '__main__':
from stormdb.access import Query
proj_name = 'MEG_service'
subj_id = '0032'
qy = Query(proj_name)
series = qy.filter_series('audvis*', subjects=subj_id)
cur_series = series[0]
delays = extract_delays(cur_series,
stim_chan='STI101',
misc_chan='MISC001',
trig_codes=[2],
plot_figures=True)
def test_get_subjects():
subs = Query(proj_name).get_subjects()
assert_true(isinstance(subs, list))
assert_equal(subs[0], test_subject_name) # 0001 is excluded in DB!
autobad="on",
st=True,
st_buflen=30,
st_corr=0.95,
movecomp=True,
cal=cal,
ctc=ctc,
mx_args='',
maxfilter_bin='maxfilter',
force=True)
# path to submit_to_isis
cmd = "/usr/local/common/meeg-cfin/configurations/bin/submit_to_isis"
proj_code = "MINDLAB2011_24-MEG-readiness"
db = Query(proj_code)
proj_folder = os.path.join('/projects', proj_code)
script_dir = proj_folder + '/scripts/'
conditions = {"classic": "2", "plan": "4", "interupt": "3"}
included_subjects = db.get_subjects()
# just test with first one!
included_subjects = included_subjects[6:]
for sub in included_subjects:
MEG_study = db.get_studies(sub, modality='MEG')
for condition in conditions:
in_name = db.get_files(sub, MEG_study[0], 'MEG', conditions[condition])
