[gmx-users] dPCA analysis
dong at pampas.chem.purdue.edu
Thu Feb 16 21:20:29 CET 2006
According to David's instruction at
gromax 3.3 can be used to do dPCA (dihedrals PCA) analysis by
1. g_angle -or to get angle.trr
2. g_covar -f angle.trr -s *.tpr to get eigenvec.trr
3. g_anaeig -f angle.trr -v eigenvec.trr -proj to get the projection of
angle.trr on eigenvec.trr
However, when I tried it in the way I mentioned above, I did not get
what I expected. Therefore, I download Dr. Mu's code from
and did the analysis again following Mu's recipe. I got a totally
After my investigation, I have a question regarding to the dPCA method
in gromacs 3.3. In the conventional PCA, we need to remove the
translational and rotational motion before generating the covariance
matrix. That is what I think g_covar is actually doing.
In dPCA, we should not do any modification of coordinates in angle.trr.
I mean, no translation and rotation removal. This should be so because
the dihedrals (or cos and sin) are internal coordinates.
Then, I think we can not use g_covar directly to analyze angle.trr,
rather we need to modify g_covar.
Has anyone used the gmx3.3 to do dPCA analysis? If so please give me
your comments on my thoughts. Thank you in advance!
More information about the gromacs.org_gmx-users