Hi,
On Wed, Jun 18, 2008 at 1:07 PM, <fatima.chami.durham.ac.uk> wrote:
>
> Hi folks,
>
> I am trying to calculate the cross correlation time function according to
> the
> ired approach. The ired eigen vector file shows undefined values:
>
>
> Eigenvector file: IRED
> 2 2
> 1.00000 1.00000
> ****
> 1 1.10699
> -0.70711 0.70711
> ****
> 2 0.89301
> -0.70711 -0.70711
>
>
> my ptraj script is as follow:
> -------------------------------------------------------
> trajin stack8w_md8_no_wat.mdcrd
> rms first
>> If you do "rms fit", you will remove the overall-rotational motion.
>> Don't superimpose your trajectories for iRED.
> vector v1 :4.C11 ired :4.O2
> vector v2 :5.C11 ired :5.O2
> matrix ired name residue_anal
> analyze matrix residue_anal vecs 6 out ired.vec thermo
>> vecs 2 (you have only two vectors)
>
> vector v1 :4.C11 corrired :4.O2 order 2 modes ired.vec beg 1 end 6 npair
> 1
> vector v2 :5.C11 corrired :5.O2 order 2 modes ired.vec beg 1 end 6 npair
> 2
>> beg 1 end 2
>
> analyze timecorr vec1 v1 vec2 v2 out vired.out
> ----------------------------------------------------------
>
> The logfile of ptraj
>
> Checking coordinates: stack8w_md8_no_wat.mdcrd
> ptraj(), analyzeMatrix: parameter "thermo" only
> works for MATRIX_MWCOVAR, returning.
> FYI: Number of read evecs is 2, number of requested evecs is 6
> WARNING in ptraj(), analyze timecorr: only calculating IRED corr for vec1,
> ignoring vec2
>
>
> best wishes
> Fatima
> -----------------------------------------------------------------------
> The AMBER Mail Reflector
> To post, send mail to amber.scripps.edu
> To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
> to majordomo.scripps.edu
>
--
Best wishes,
Myunggi Yi PhD
==================================
KLB 419
Institute of Molecular Biophysics
Florida State University
Tallahassee, FL 32306
Office: (850) 645-1334
http://www.scs.fsu.edu/~myunggi
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber.scripps.edu
To unsubscribe, send "unsubscribe amber" (in the *body* of the email)
to majordomo.scripps.edu
Received on Sun Jun 22 2008 - 06:07:15 PDT
