diff --git a/pyemma/msm/tests/test_bayesian_hmsm.py b/pyemma/msm/tests/test_bayesian_hmsm.py
index <HASH>..<HASH> 100644
--- a/pyemma/msm/tests/test_bayesian_hmsm.py
+++ b/pyemma/msm/tests/test_bayesian_hmsm.py
@@ -320,5 +320,19 @@ class TestBHMMSpecialCases(unittest.TestCase):
             assert strajs[0][0] == 2
             assert strajs[0][6] == 2
 
+    def test_initialized_bhmm(self):
+        import pyemma.datasets
+        import pyemma.msm
+
+        obs = pyemma.datasets.load_2well_discrete().dtraj_T100K_dt10
+
+        init_hmm = pyemma.msm.estimate_hidden_markov_model(obs, 2, 10)
+        bay_hmm = pyemma.msm.estimators.BayesianHMSM(nstates=init_hmm.nstates, lag=init_hmm.lag,
+                                                     init_hmsm=init_hmm)
+        bay_hmm.estimate(obs)
+
+        assert np.isclose(bay_hmm.stationary_distribution.sum(), 1)
+
+
 if __name__ == "__main__":
     unittest.main()
[TestBHMMSpecialCases] added test for initialized bhmm