def interactiontest(dm):
dm.bl = baseline.baseline(dm)
dm.test = srs.reduce_(srs.window(dm.y, start=-50, end=-1))
lmdata = DataMatrix(length=len(dm) * 2)
lmdata.pupil = -1
lmdata.time = -1
lmdata.condition = -1
lmdata.trialid = -1
for i, row in enumerate(dm):
lmdata.pupil[2 * i] = row.bl
lmdata.time[2 * i] = 0
lmdata.condition[2 * i] = row.c
lmdata.trialid[2 * i] = i
lmdata.pupil[2 * i + 1] = row.test
lmdata.time[2 * i + 1] = 1
lmdata.condition[2 * i + 1] = row.c
lmdata.trialid[2 * i + 1] = i
lm = lme4.lmer(lmdata, 'pupil ~ condition*time + (1|trialid)')
print(lm.p[3])
return lm.p[3] < ALPHA
#print(ef.report)
TNTdm = DataMatrix(length=4)
TNTdm.efficiency = 'Fe','Fd','Ff','Fc'
TNTdm.score = (TNTdES.Fe,TNTdES.Fd,TNTdES.Ff,TNTdES.Fc)
print('TNT Run-1 Power & Efficiency Matrix:')
print(TNTdm)
TNTeffTOT = list(TNTdm.score)
TNTeffSUM = sum(TNTeffTOT)
print("TNT Run-1 Power & Efficiency Sum: "+str(TNTeffSUM))
TNTconditions = list(TNTpseudM.condition)
TNTtrialDm = DataMatrix(length=124)
TNTtrialDm.condition = (TNTconditions)
TNTtrialDm.order = (TNTdES.order)
TNTtrialDm.ITI = (TNTdES.ITI)
TNTtrialDm.onset = (TNTdES.onsets)
print('TNT Run-1 Trial Order Matrix:')
print(TNTtrialDm)
#print("Condition oder:")
#print(TNTorder)
#print("Interstimulus Intervals:")
#print(TNTdES.ITI)#[:120])
#print("Condition Onsets")
#print(TNTdES.onsets)#[:120])
