Skip to content
This repository has been archived by the owner on Jan 26, 2022. It is now read-only.

What are the units of distance in dihedral clustering? #443

Open
vvoelz opened this issue Oct 8, 2014 · 2 comments
Open

What are the units of distance in dihedral clustering? #443

vvoelz opened this issue Oct 8, 2014 · 2 comments
Labels

Comments

@vvoelz
Copy link
Contributor

vvoelz commented Oct 8, 2014

Suppose I perform dihedral clustering and MSMBuilder reports that the cluster radius is 1.8. What are the units here? I know that the the angles phi_i get projected 2-fold onto x1_i = cos(phi_I) and x2_i = sin(phi_i). What is the "rmsd" that gets reported?

@vvoelz vvoelz added the docs label Oct 8, 2014
@schwancr
Copy link
Contributor

schwancr commented Oct 8, 2014

I don't really think it has a meaningful unit, unfortunately. If you convert each dihedral angle into a complex number (via the Euler form exp(i \phi)), then the distance corresponds to the L2 norm on a vector of these complex numbers.

That, unfortunately, doesn't really have a nice unit associated with it.

@vvoelz
Copy link
Contributor Author

vvoelz commented Oct 9, 2014

Actually that’s exactly the info I need. I have N=15 angles, and you’re telling me that distances are computed as L2 norms. So, if the reported distance is d = ( \sum_{i=1}^{N} ||\delta z_i ||^2 )^{1/2}, then the "average" length of a complex difference vector is z = d N^{-1/2}. Thus, we can use the law of cosines,
z^2 = (1)^2 + (1)^2 - 2(1)(1)\cos(\phi) = 2(1 - \cos(\phi) )
to back-calculate \phi, the "average" dihedral angle that separates any two dihedral clusters.

In my specific case; d = 1.8, so \phi = \cos^{-1}( 1 - d^2/(2N) ) = 26.8 degrees

Is this reasoning sound, or am I missing something?

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

No branches or pull requests

2 participants