Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

PMI2 RMF hierarchies of input and output RMFs different #224

Open
shruthivis opened this issue Nov 15, 2016 · 0 comments
Open

PMI2 RMF hierarchies of input and output RMFs different #224

shruthivis opened this issue Nov 15, 2016 · 0 comments

Comments

@shruthivis
Copy link

There is a difference in the hierarchies between output and input RMFs of the PMI version I am using.
For example, after sampling:
/salilab/park1/shruthi/sampcon/pmibugs/1AVX/1AVX_4.1
output/initial.0.rmf3 for example has a hierarchy like System->State 0-> Molecule A, Molecule B.

But after getting the top 500 models from the run using AnalysisReplicaExchangeMacro (Script: /salilab/park1/shruthi/sampcon/get_top_models_each_run.py)
the hierarchy becomes State 0->Molecule A, Molecule B (missing System).
e.g. all_models.499/0.0.rmf3 in the same example directory.

Apparently Chimera can read both kinds of RMFs but PMI itself cannot.
When I try to run the Analysis macro again to do clustering etc. it fails (script: /salilab/park1/shruthi/sampcon/cluster_k1.py)
PMI does not throw up an error but simply outputs an empty distance matrix and AVERAGE_RMSD in the stat files in the cluster directories is 10000000000.0, because it does not read any coordinates from the RMFs.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant