# Tecott Lab HCM Data
# working with AS/IS, IS Duration Thresholds (ISDT)
# using class intervals objects
# C. Hillar, Jan 2016
from __future__ import print_function, absolute_import, division
import numpy as np
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)
# Tecott Lab HCM Data
# working with AS/IS, IS Duration Thresholds (ISDT)
# using class intervals objects
# C. Hillar, Jan 2016
from __future__ import print_function, absolute_import, division
import numpy as np
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)
