Hi Mark,<div><br></div><div>That is true. But what does the manual refer to when it says </div><div><br></div><div><span class="Apple-style-span" style="font-family: arial, sans-serif; font-size: 13px; border-collapse: collapse; ">the first graph is the average, the rest are the individual angles</span></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></font></div><div><span class="Apple-style-span" style="font-size: 13px; "></span><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">what and where are the rest? I am only curious because there might be a hidden feature that could be useful for me.</span></font></div>
<div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;"><br></span></font></div><div><font class="Apple-style-span" face="arial, sans-serif"><span class="Apple-style-span" style="border-collapse: collapse;">Amit<br>
</span></font><br><div class="gmail_quote">On Wed, Dec 2, 2009 at 12:21 AM, Mark Abraham <span dir="ltr"><<a href="mailto:Mark.Abraham@anu.edu.au">Mark.Abraham@anu.edu.au</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Amit Choubey wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi Mark,<br>
<br><div class="im">
Ok here is the command line i used :-<br>
<br>
g_angle -f after_md.trr -type dihedral -n chain1.ndx -b 800 -e 804 -all <br>
output is <br>
# This file was created Tue Dec 1 20:20:41 2009<br>
# by the following command:<br>
# g_angle -f after_md.trr -type dihedral -n chain1.ndx -b 800 -e 804 -all<br>
#<br>
# g_angle is part of G R O M A C S:<br>
#<br>
# Go Rough, Oppose Many Angry Chinese Serial killers<br>
#<br>
@ title "Dihedral Distribution: C34_C36_C37_C38"<br>
</div></blockquote>
<br>
This title makes it look like there is only one angle in your index file.<br>
<br>
Mark<br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div></div><div class="h5">
@ xaxis label "Degrees"<br>
@ yaxis label ""<br>
@TYPE xy<br>
@ subtitle "average angle: 9.87432\So\N"<br>
@with g0<br>
@ world xmin -180<br>
@ world xmax 180<br>
@ world ymin 0<br>
@ world ymax 0.0451172<br>
@ xaxis tick major 60<br>
@ xaxis tick minor 30<br>
@ yaxis tick major 0.005<br>
@ yaxis tick minor 0.0025<br>
-180 0.015625<br>
-179 0.011719<br>
-178 0.015625<br>
-177 0.000000<br>
-176 0.035156<br>
-175 0.031250<br>
-174 0.019531<br>
-173 0.019531<br>
-172 0.023438<br>
-171 0.015625<br>
-170 0.015625<br>
<br>
I did not include the distribution for rest of the angles.<br>
The second column is the fractional distribution of dihedrals with the corresponding angle in the 1st column.<br>
<br>
The distribution is averaged over time t = 800ps and 802 ps snapshots. I checked this part.<br>
<br>
Amit<br>
<br>
<br>
<br>
<br></div></div><div><div></div><div class="h5">
On Tue, Dec 1, 2009 at 8:15 PM, Mark Abraham <<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a> <mailto:<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a>>> wrote:<br>
<br>
Amit Choubey wrote:<br>
<br>
Hi Mark,<br>
<br>
It does seem that the -all option averages the distributions.<br>
The thing is i was confused with "the first graph is the<br>
average, the rest are the individual angles." because when i<br>
tried it previously i only found 2 columns which were angle and<br>
distribution.<br>
<br>
But now i checked again explicitly by doing few more g_angle<br>
commands and it seems that -all does do what i wanted to do. I<br>
still dont know what the "rest are the individual angles" means.<br>
<br>
<br>
I haven't ever used it, but would have expected the output to be<br>
along the lines of<br>
<br>
Time "Average over all angles" Angle-1 Angle-2 ...<br>
1 234 23 345<br>
2 123 234 232<br>
2 223 ...<br>
3 323<br>
<br>
What output do you get and what seems wrong to you?<br>
<br>
Mark<br>
<br>
On Tue, Dec 1, 2009 at 6:52 PM, Mark Abraham<br>
<<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a> <mailto:<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a>><br>
<mailto:<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a><br>
<mailto:<a href="mailto:Mark.Abraham@anu.edu.au" target="_blank">Mark.Abraham@anu.edu.au</a>>>> wrote:<br>
<br>
Amit Choubey wrote:<br>
<br>
Hi everyone,<br>
<br>
I am trying to analyze dihedrals of molecules after a<br>
simulation. I am able to calculate dihedral distribution<br>
at ant<br>
particular time using g<br>
-angle and proper group using<br>
<br>
g_angle -f after_md.trr -b 800 -type dihedral -n<br>
chain1.ndx -od<br>
dihed12<br>
<br>
This gives me the distribution at t=800 ps . What i<br>
really want<br>
is to see the dihedral distribution at all time steps<br>
(after 800<br>
ps) in one go, and may be do some averaging later. Is there a<br>
straight forward way to do this. I can imagine that i could<br>
write a script file that could do this but then i dont<br>
want to<br>
deal with a lot of files, i rather want everything to be<br>
saved<br>
in a single file with multiple columns. Could someone<br>
suggest a<br>
trick for this?<br>
<br>
<br>
The first paragraph of g_angle -h suggests there's an option<br>
where<br>
"the first graph is the average, the rest are the individual<br>
angles." Does this work?<br>
<br>
<br>
Mark<br>
-- gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br></div></div>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a> <mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>>><div class="im"><br>
<br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before<br>
posting!<br>
Please don't post (un)subscribe requests to the list. Use the www<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>>.<br>
<br>
Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
<br>
<br>
-- gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a>><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before<br>
posting!<br>
Please don't post (un)subscribe requests to the list. Use the www<br>
interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a><br>
<mailto:<a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>>.<br>
Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
<br>
<br>
</div></blockquote><div><div></div><div class="h5">
-- <br>
gmx-users mailing list <a href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>
<a href="http://lists.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://lists.gromacs.org/mailman/listinfo/gmx-users</a><br>
Please search the archive at <a href="http://www.gromacs.org/search" target="_blank">http://www.gromacs.org/search</a> before posting!<br>
Please don't post (un)subscribe requests to the list. Use the www interface or send it to <a href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org</a>.<br>
Can't post? Read <a href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php</a><br>
</div></div></blockquote></div><br></div>