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
# 16 Strains:
strains = {
0: 'C57BL6J',
1: 'BALB',
2: 'A',
3: '129S1',
4: 'DBA',
5: 'C3H',
6: 'AKR',
7: 'SWR',
8: 'SJL',
9: 'FVB',
