fail to compute implied timescale after hybrid clustering

Provide extensible software for building Markov State Models.
POST REPLY
User avatar
Lizhe Zhu
Posts: 1
Joined: Tue Dec 18, 2012 3:12 am

fail to compute implied timescale after hybrid clustering

Post by Lizhe Zhu » Sun Jan 06, 2013 2:23 am

Dear developers,

I am testing the lastest version MSMbuilder2.5.1 on my own MD data.

I went through the tutorial completely and could repeat each step on my own computer. This should validate the installation.

But when I try to build models with my own data, I was stuck at when computing the implied Timescale.

My trajectories were converted and clustered into microstates without problem using
ConvertDataToHDF.py -s protein.pdb -i trj_data
Cluster.py rmsd hybrid -d 0.025
But when I use CalculateImpliedTimescales.py as follows,
CalculateImpliedTimescales.py -l 1,25 -i 1 -o Data/ImpliedTimescales.dat
I get the following error
Getting 10 eigenvalues (timescales) for each lagtime...
Building MSMs at the following lag times: [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25]
Selected component 0 with population 0.000084
Traceback (most recent call last):
File "MSMBuilder2/bin/CalculateImpliedTimescales.py", line 87, in <module>
args.symmetrize, args.procs, args.output)
File "MSMBuilder2/bin/CalculateImpliedTimescales.py", line 48, in run
Slide=True, Symmetrize=symmetrize, nProc=nProc)
File "MSMBuilder2/lib/python2.7/site-packages/msmbuilder/MSMLib.py", line 504, in GetImpliedTimescales
lags = result.get(999999)
File "epd-7.3-64bit/lib/python2.7/multiprocessing/pool.py", line 528, in get
raise self._value
ValueError: cannot infer dimensions from zero sized index arrays
Anyone has a clue what the problem is?
Thanks in advance.

Lizhe

User avatar
Robert McGibbon
Posts: 20
Joined: Tue Jul 19, 2011 9:25 am

Re: fail to compute implied timescale after hybrid clusterin

Post by Robert McGibbon » Wed Feb 27, 2013 3:10 pm

Lizhe,

This looks like a shape error. Perhaps your assignments array 1 dimensional instead of the required 2 dimensional? We (the developers) have been trying to improve the quality of the error messages in the 2.6 release.

Also, sorry for not checking this forum for a very long time. Please sign up for the msmbuilder-user email list, which will hopefully be a little more active than these forums. https://mailman.stanford.edu/mailman/li ... ilder-user

-Robert

POST REPLY