Dense [k-1]x4 transition frequency matrix

def transitions(self):
        """Dense [k-1]x4 transition frequency matrix"""
        if self._transitions is not None:
            return self._transitions
        transitions = self.array.astype(np.float)
        transitions /= transitions.sum(1)[:, np.newaxis]
        self._transitions = transitions
        return transitions