Dear Hannes,
Thank you. The problem has been solved.
Best regards,
Cat
> Subject: RE: [AMBER] PCA analysis: the first vector is the movement of whole molecule.
> From: hannes.loeffler.stfc.ac.uk
> To: amber.ambermd.org
> Date: Fri, 31 Jul 2009 08:27:54 +0100
>
> Ok, I see that you don't really understand what you are doing.
>
> The '*' mask means ALL atoms. Unless this selection is not just the
> molecule(s) you are interested in you are obviously trying to fix the
> wrong part of your system. So make sure that you use here the same mask
> as in the 'matrix' command.
>
> Why simulation length matters? PCA is a statistical method as you are
> calculating the variance-covariance matrix as an _average over time_.
> This should make it clear to you that you will need sufficient sampling
> before you can rely on your results.
>
> Here are two sample input files for computing the
> eigenvectors/eigenmodes and the projection of the eigenvectors onto the
> trajectory for a protein with 100 residues. The examples are just meant
> to guide you. So adapt it to your problem.
>
> -------------------------- PCA.in ------------------------
> rms first :1-100.CA
>
> # PCA analysis
> matrix covar name PCA :1-100.CA
> analyze matrix PCA vecs 25 out evec.pev
> -------------------------- PCA.in ------------------------
>
> -------------------------- proj.in ------------------------
> # centering and imaging not necessary if selection is just one molecule
> center :1-100 mass origin
> image origin center
>
> rms first :1-100.CA
>
> projection :1-100.CA modes evec.pev out proj.ppj beg 1 end 25
> -------------------------- proj.in ------------------------
>
>
> Hope that helps,
> Hannes.
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
_________________________________________________________________
With Windows Live, you can organize, edit, and share your photos.
http://www.microsoft.com/middleeast/windows/windowslive/products/photo-gallery-edit.aspx
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Aug 19 2009 - 20:13:29 PDT
