def lagcor_is_trimcor(self, h): a = self.a b = self.b assert_almost_equal(lagcor(a, b, h), corsubs(trim(a, h=h), trim(b, h=-h)))
def lagcor_equal_subcor(self, key, standardized=False): x, y = self.ts[key][0:2] assert_almost_equal(lagcor(x, y, 0), corsubs(x, y))
def lagcor_is_shiftcor(self, h): a = self.a b = self.b outlen = a.shape[-1] - h assert_almost_equal(lagcor(a, b, h), corsubs(shift(a, h), shift(b, 0, outlen)))
def lagcor_is_trimcor(self, h): a = self.a b = self.b assert_almost_equal(lagcor(a,b,h), corsubs(trim(a, h=h), trim(b, h=-h)))
def lagcor_is_shiftcor(self, h): a = self.a b = self.b outlen = a.shape[-1]-h assert_almost_equal(lagcor(a, b, h), corsubs(shift(a, h), shift(b, 0, outlen)))
def lagcor_equal_subcor(self, key, standardized=False): x, y = self.ts[key][0:2] assert_almost_equal(lagcor(x,y,0), corsubs(x,y))
for condname in config['conds']: cond = E.get_cond(condname) for run in E.iter_runs(condname): run.threshold(cond.attrs['threshold'], save=True) # Alignment ------------------------------------------- a1 = E.f['rois/audenv_3mm_gdboth_wordscram_mean_thr0.2'][ ...] #TODO remove hardcoded roi? aud_env = cond['audio_env'] # Calculate lags over range, then set offset lags = xrange(-5, 5) for run in E.iter_runs(condname): data = run.load(roi=a1, standardized=True, threshold=True) #slice(None,None,None)) corrs = [lagcor(data, aud_env, h=h) for h in lags] print corrs offset = lags[np.argmax(corrs)] run.attrs['offset'] = offset # With offsets down, can generate composite if condname != "dummy": E.gen_composite(condname) # ISC Calculations ---------------------------------------- cond1 = E.get_cond('wordscram') # conditions for calculations cond2 = E.get_cond('wordscram2') # Run calculations isc_within(E, cond1) isc_between(E, cond1, cond2, method='inter-subject')
# Iterate over conditions --------------------------------- for condname in config['conds']: cond = E.get_cond(condname) for run in E.iter_runs(condname): run.threshold(cond.attrs['threshold'], save=True) # Alignment ------------------------------------------- a1 = E.f['rois/audenv_3mm_gdboth_wordscram_mean_thr0.2'][...] #TODO remove hardcoded roi? aud_env = cond['audio_env'] # Calculate lags over range, then set offset lags = xrange(-5,5) for run in E.iter_runs(condname): data = run.load(roi=a1, standardized=True, threshold=True)#slice(None,None,None)) corrs = [lagcor(data, aud_env, h=h) for h in lags] print corrs offset = lags[np.argmax(corrs)] run.attrs['offset'] = offset # With offsets down, can generate composite if condname != "dummy":E.gen_composite(condname) # ISC Calculations ---------------------------------------- cond1 = E.get_cond('wordscram') # conditions for calculations cond2 = E.get_cond('wordscram2') # Run calculations isc_within(E, cond1) isc_between(E, cond1, cond2, method='inter-subject')