Re: [AMBER] Percent Occupancey with PCA

From: Daniel Roe <daniel.r.roe.gmail.com>
Date: Sat, 3 May 2014 01:11:50 -0600

Hi,

You can generate histograms of the principal component projections (I
think this may be what you want) using CPPTRAJ, e.g.:

parm mytopology.parm7
# Average the trajectory
rms first $MASK
average Avg.rst7 restart
createcrd crd1
run
# Fit to average structure
reference Avg.rst7.1
crdaction crd1 rms reference $MASK
# Calculate and diagonalize covariance matrix
crdaction crd1 matrix covar $MASK name Covar
runanalysis diagmatrix Covar out evecs.dat vecs $NVECS name Evecs
# Now calculate projection for trajectory
crdaction crd1 projection T1 modes Evecs beg 1 end $NVECS $MASK
# Histogram projection
runanalysis multihist T1 bins 100 normint out hist.agr

Replace $MASK with the appropriate atom mask and $NVECS with the
desired # of vectors.

Hope this helps,

-Dan

On Fri, May 2, 2014 at 9:17 AM, Nicole Ippolito
<ippolito.nicole89.gmail.com> wrote:
> Dear Amber Users,
>
> Is there a way to generate the percent occupancy with principal component
> analysis using ptraj?
>
> Thank you,
>
> Nicole
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber



-- 
-------------------------
Daniel R. Roe, PhD
Department of Medicinal Chemistry
University of Utah
30 South 2000 East, Room 201
Salt Lake City, UT 84112-5820
http://home.chpc.utah.edu/~cheatham/
(801) 587-9652
(801) 585-6208 (Fax)
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Sat May 03 2014 - 00:30:02 PDT
Custom Search