code: def _indicator_multitraj(ss, i, j):
      
      iseqs = ss[i]
      res = []
      for iseq in iseqs:
          x = np.zeros(len(iseq))
          I = np.where(iseq == j)
          x[I] = 1.0
          res.append(x)
      return res
docs: Returns conditional sequence for transition i -> j given all conditional sequences