<div>I think g_covar -ascii is what you are looking for.</div>
<div>I have a Maple script that does this kind of graph. If you want you can contact me off list.<br> </div>
<div>Regards.</div>
<div> </div>
<div>Pedro.<br> </div>
<div><span class="gmail_quote">2007/9/18, Dhananjay <<a href="mailto:dhananjay.c.joshi@gmail.com">dhananjay.c.joshi@gmail.com</a>>:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid">Well, I have done 30ns MD simulation for a protein. The protein consists of 4 loop regions. <br>Using g_cover programme the covariance matrix have been generated and using g_anaeig programme with option -rmsf , along first 8 eigen vectors the rms fluction of c-alpha atoms were plotted.
<br>But this is giving me just the information that particular atom is fluctuating along particular direction.<br>I want to know over a 30 ns simulation whether the motions of two atoms or group of two atoms are correlated or anti correlated. For this I want to form a dynamics cross correlation map in which I could get the precise information of group of atoms. The function used for this kind of analysis is
<br><br>C(i,j) = < delta <span style="FONT-WEIGHT: bold">r</span>(i) * delta <span style="FONT-WEIGHT: bold">r</span>(j) > / sqrt < sqr(delta <span style="FONT-WEIGHT: bold">r</span>(i) ) > . sqrt < sqr(delta
<span style="FONT-WEIGHT: bold">r</span>(j) ) ><br><br>the positive C(i,j) -> motions are correlated<br>the negative C(i,j) -> motions are anti-correlated<br><br>I want to plot a 2-D map indicated correlated and anti-correlated motions.
<br><br>My question is whether it is possible in GROMACS to plot this king of map or could you please suggest any other free software which could read the trajectories generated by GROMACS and plot the map.<br><br>Again thanking you in advance ......
<br><br><br>-- Dhananjay<br><br><br><br><br>
<div><span class="gmail_quote">On 9/18/07, <b class="gmail_sendername">David van der Spoel</b> <<a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:spoel@xray.bmc.uu.se" target="_blank">spoel@xray.bmc.uu.se
</a>> wrote:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0pt 0pt 0pt 0.8ex; BORDER-LEFT: rgb(204,204,204) 1px solid">Dhananjay wrote:<br>> This may be bit different question from main theme.<br>><br>> I want to form dynamic cross correlation map as I have trajectories
<br>> generated by GROMACS 3.3.<br>><br>> Please suggest me free software which can read the .trr file for<br>> plotting DCCM (dynamic cross correlation map )<br>try g_covar or otherwise explain in more detail (equations) what you
<br>want to do.<br><br>><br>> Thanking you in advance....<br>><br>><br>><br>> -- Dhananjay<br>><br>><br>> ------------------------------------------------------------------------<br>><br>> _______________________________________________
<br>> gmx-users mailing list <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org</a><br>> <a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">
http://www.gromacs.org/mailman/listinfo/gmx-users</a> <br>> Please search the archive at <a onclick="return top.js.OpenExtLink(window,event,this)" 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<br>> www interface or send it to <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users-request@gromacs.org" target="_blank">
gmx-users-request@gromacs.org</a>.<br>> Can't post? Read <a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">http://www.gromacs.org/mailing_lists/users.php
</a><br><br><br>--<br>David van der Spoel, Ph.D.<br>Molec. Biophys. group, Dept. of Cell & Molec. Biol., Uppsala University.<br>Box 596, 75124 Uppsala, Sweden. Phone: +46184714205. Fax: +4618511755.<br><a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:spoel@xray.bmc.uu.se" target="_blank">
spoel@xray.bmc.uu.se</a> <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:spoel@gromacs.org" target="_blank">spoel@gromacs.org</a> <a onclick="return top.js.OpenExtLink(window,event,this)" href="http://folding.bmc.uu.se/" target="_blank">
http://folding.bmc.uu.se</a><br>_______________________________________________<br>gmx-users mailing list <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users@gromacs.org" target="_blank">gmx-users@gromacs.org
</a><br><a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please search the archive at
<a onclick="return top.js.OpenExtLink(window,event,this)" 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
<br>www interface or send it to <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users-request@gromacs.org" target="_blank">gmx-users-request@gromacs.org </a>.<br>Can't post? Read <a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">
http://www.gromacs.org/mailing_lists/users.php</a><br></blockquote></div><span class="sg"><br><br clear="all"><br>-- <br>Dhananjay C Joshi<br>Project Assistant <br>LT & LSB, C D F D<br>ECIL Road, Nacharam<br>Hyderabad-500 076, INDIA
<br>Tel : +91-40-27151344<br>Fax : +91-40-27155610 </span><br>_______________________________________________<br>gmx-users mailing list <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users@gromacs.org">
gmx-users@gromacs.org</a><br><a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailman/listinfo/gmx-users" target="_blank">http://www.gromacs.org/mailman/listinfo/gmx-users</a><br>Please search the archive at
<a onclick="return top.js.OpenExtLink(window,event,this)" 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
<br>www interface or send it to <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:gmx-users-request@gromacs.org">gmx-users-request@gromacs.org</a>.<br>Can't post? Read <a onclick="return top.js.OpenExtLink(window,event,this)" href="http://www.gromacs.org/mailing_lists/users.php" target="_blank">
http://www.gromacs.org/mailing_lists/users.php</a><br></blockquote></div><br>