markovmodel / PyEMMA

🚂 Python API for Emma's Markov Model Algorithms 🚂
http://pyemma.org
GNU Lesser General Public License v3.0
311 stars 119 forks source link

Implied Timescales estimator failure #1470

Closed ghost closed 4 years ago

ghost commented 4 years ago

Hello,

I'm trying to build a Markov State Model for an RNA system and I have issues with the msm estimation and validation. When I run

tica = coor.tica(inp, lag=20, dim=30) tica_concatenated = np.concatenate(tica.get_output()) cl = coor.cluster_kmeans(tica, k=200, max_iter=100) its = msm.its(cl.dtrajs, nits=10, errors='bayes')

I get the error:

Screen Shot 2020-08-18 at 2 46 56 PM

When I plot its, this is what I see:

Figure_3-cluster-its

I don't understand where the issue is.

Your help is greatly appreciated.