Add analysis for symmetry corrected RMSD#92
Conversation
|
The difference between the old ligand RMSD and the symmetry corrected ligand RMSD on the skipped dataset is quite big and the symmetry corrected ligand RMSD is greater (esp. at frame 5) which seemed weird. Questions:
|
IAlibay
left a comment
There was a problem hiding this comment.
Overally this looks great, just a few comments.
Co-authored-by: Irfan Alibay <IAlibay@users.noreply.github.com>
Co-authored-by: Josh Horton <joshua.horton@openforcefield.org>
|
pre-commit.ci autofix |
Codecov Report✅ All modified and coverable lines are covered by tests. Additional details and impacted files@@ Coverage Diff @@
## main #92 +/- ##
==========================================
+ Coverage 97.84% 98.05% +0.21%
==========================================
Files 9 9
Lines 417 462 +45
==========================================
+ Hits 408 453 +45
Misses 9 9 ☔ View full report in Codecov by Sentry. 🚀 New features to boost your workflow:
|
IAlibay
left a comment
There was a problem hiding this comment.
Couple of small things, mostly the rdkit addition to pyproject.toml. Otherwise this lgtm!
| if ligand: | ||
| # For now, leave it at the normal RMSD | ||
| lig_rmsd = RMSDAnalysis(ligand, mass_weighted=True).run(step=skip) | ||
| # state_lig = select_state_atoms(u, end_state="A").select_atoms("resname UNK") |
There was a problem hiding this comment.
Could you migrate this to an issue and remove the comments here?
| "pass an rdmol directly." | ||
| ) | ||
| self._mol = rdmol if rdmol is not None else atomgroup.convert_to("RDKIT") | ||
| self._aprops = np.array([atom.GetAtomicNum() for atom in self._mol.GetAtoms()]) |
There was a problem hiding this comment.
It might be good to add a check here that ensures that the atomgroup has the same number of atoms as the rdmol.
No description provided.