Beispiel #1
0
 def sample_forwards(self,betal,betastarl):
     from pyhsmm.internals.hsmm_messages_interface import sample_forwards_log
     if self.left_censoring:
         raise NotImplementedError
     caBl = np.vstack((np.zeros(betal.shape[1]), np.cumsum(self.aBl[:-1],axis=0)))
     self.stateseq = sample_forwards_log(
             self.trans_matrix, caBl, self.aDl, self.pi_0, betal, betastarl,
             np.empty(betal.shape[0],dtype='int32'))
Beispiel #2
0
 def sample_forwards(self,betal,betastarl):
     from pyhsmm.internals.hsmm_messages_interface import sample_forwards_log
     if self.left_censoring:
         raise NotImplementedError
     caBl = np.vstack((np.zeros(betal.shape[1]),np.cumsum(self.aBl[:-1],axis=0)))
     self.stateseq = sample_forwards_log(
             self.trans_matrix,caBl,self.aDl,self.pi_0,betal,betastarl,
             np.empty(betal.shape[0],dtype='int32'))
     assert not (0 == self.stateseq).all()