I don't understand what you want exactly.
Your g_energy command does exponential averaging, that happens
on the printed data points. So there g_analyze or any program will do fine.
For the original data g_energy gives the exact standard deviation
over all MD steps, called RMSD.
But do you really want the standard deviation, or do you want
a standard error estimate?
Berk
Date: Fri, 28 Aug 2009 09:30:37 -0300
From: [email protected]
To: [email protected]
Subject: [gmx-users] Re: g_energy and g_analyze give different averages
Is there a correct way to obtain the standard deviation for these data sets?
>
> Ragnarok sdf wrote:
> > When analysing FEP simulations. After running g_energy -f fep000.edr -f2
> > fep_000-005.edr for obtaining the dF = -kT ln < exp(-(EB-EA)/kT) >A I
> tried
> > to obtain the standard deviation for this ensemble average using
> g_analyze
> > -f runavg.xvg, but I've noticed that the average values are quite
> different
> >
> >
> g_energy takes all values during your simulation into account, g_analyze
> only the printed datapoints.
>
>
_________________________________________________________________
See all the ways you can stay connected to friends and family
http://www.microsoft.com/windows/windowslive/default.aspx
_______________________________________________
gmx-users mailing list [email protected]
http://lists.gromacs.org/mailman/listinfo/gmx-users
Please search the archive at http://www.gromacs.org/search before posting!
Please don't post (un)subscribe requests to the list. Use the
www interface or send it to [email protected].
Can't post? Read http://www.gromacs.org/mailing_lists/users.php
