pyslds icon indicating copy to clipboard operation
pyslds copied to clipboard

error on simple example

Open andrewcz opened this issue 6 years ago • 3 comments

hey champion just getting some errors on the simple example

......................... [ 825/1000, 0.04sec avg, ETA 6.19 ] ......................... [ 850/1000, 0.04sec avg, ETA 5.31 ] ......................... [ 875/1000, 0.04sec avg, ETA 4.42 ] ......................... [ 900/1000, 0.04sec avg, ETA 3.54 ] ......................... [ 925/1000, 0.04sec avg, ETA 2.65 ] ......................... [ 950/1000, 0.04sec avg, ETA 1.77 ] ......................... [ 975/1000, 0.04sec avg, ETA 0.88 ] ......................... [ 1000/1000, 0.04sec avg, ETA 0.00 ]

0.04sec avg, 35.35 total

Fitting with VBEM ................/home/andrewcz/anaconda3/lib/python3.6/site-packages/pyhsmm/internals/hmm_states.py:675: RuntimeWarning: divide by zero encountered in log np.log(trans_potential),likelihood_log_potential,alphal,betal, /home/andrewcz/pyslds/pyslds/states.py:594: RuntimeWarning: divide by zero encountered in log params = (np.log(self.trans_matrix), np.log(self.pi_0), aBl, self._normalizer) Traceback (most recent call last): File "simple_demo.py", line 95, in vbem_lls = [update(test_model) for _ in progprint_xrange(N_vbem_iters)] File "simple_demo.py", line 95, in vbem_lls = [update(test_model) for _ in progprint_xrange(N_vbem_iters)] File "simple_demo.py", line 91, in update model.VBEM_step() File "/home/andrewcz/pyslds/pyslds/models.py", line 217, in VBEM_step self._vb_E_step() File "/home/andrewcz/pyslds/pyslds/models.py", line 170, in _vb_E_step state.vb_E_step() File "/home/andrewcz/pyslds/pyslds/states.py", line 586, in vb_E_step H_z = self.vb_E_step_discrete_states() File "/home/andrewcz/pyslds/pyslds/states.py", line 595, in vb_E_step_discrete_states return hmm_entropy(params, self.all_expected_stats) File "/home/andrewcz/pyslds/pyslds/util.py", line 63, in hmm_entropy neg_entropy += np.sum(sum_E_ztztp1T * log_transmatrix) FloatingPointError: invalid value encountered in multiply

best andrew

andrewcz avatar Jun 04 '18 10:06 andrewcz