Events = Intervals([[-2, -1], [0, 1], [3, 5], [8, 11]])
ISDT = 2 # IST Threshold
ECG = Events.copy().connect_gaps(ISDT) # connect gaps <= ISDT
ETm = Events.copy().trim(ISDT) # trim inteversl <= ISDT
print("Events (%d intervals):" % Events.num(), Events)
# print ("Events connected gaps <= %1.3f:" % ISDT,
# Events.copy().connect_gaps(ISDT))
print("Events trim intervals <= %1.3f:" % ISDT, Events.copy().trim(ISDT))
print("IS/AS computation version 1")
IS1 = Events.complement().trim(ISDT) # complement is new object
AS1 = IS1.complement()
print("IS (%d):" % IS1.num(), IS1)
print("AS (%d):" % AS1.num(), AS1)
print("IS/AS computation version 2 (mathematically equivalent)")
AS2 = Events.copy().connect_gaps(ISDT) # trim alters object
IS2 = AS2.complement()
print("AS2 (%d):" % AS2.num(), AS2)
print("IS2 (%d):" % IS2.num(), IS2)
print("IS/AS built into Intervals")
print("AS1:", Events.ASs(ISDT))
print("IS1:", Events.ISs(ISDT))
# Data set consists of 1921 Mouse days (22 hours each) from 170 Mice and
from intervals import Intervals
Events = Intervals([[-2, -1], [0, 1], [3, 5], [8, 11]])
ISDT = 2 # IST Threshold
ECG = Events.copy().connect_gaps(ISDT) # connect gaps <= ISDT
ETm = Events.copy().trim(ISDT) # trim inteversl <= ISDT
print("Events (%d intervals):" % Events.num(), Events)
# print ("Events connected gaps <= %1.3f:" % ISDT,
# Events.copy().connect_gaps(ISDT))
print("Events trim intervals <= %1.3f:" % ISDT, Events.copy().trim(ISDT))
print("IS/AS computation version 1")
IS1 = Events.complement().trim(ISDT) # complement is new object
AS1 = IS1.complement()
print("IS (%d):" % IS1.num(), IS1)
print("AS (%d):" % AS1.num(), AS1)
print("IS/AS computation version 2 (mathematically equivalent)")
AS2 = Events.copy().connect_gaps(ISDT) # trim alters object
IS2 = AS2.complement()
print("AS2 (%d):" % AS2.num(), AS2)
print("IS2 (%d):" % IS2.num(), IS2)
print("IS/AS built into Intervals")
print("AS1:", Events.ASs(ISDT))
print("IS1:", Events.ISs(ISDT))
# Data set consists of 1921 Mouse days (22 hours each) from 170 Mice and
