<html><head></head><body><div style="font-family: Verdana;font-size: 12.0px;"><div>
<div>Trying to help I would say I am trying to remember some things I did some time ago. That said, your missing the time column as your first input. You are making a 2D energy plot is what your telling it to do, you can also tell it not to expect a time which would be correct for your case, however this is a landscape your generating. It is already a total energy landscape in 1 nm by 2 nm 2D volumn, what your doing is trying to use the boltzman inversion equation in the g_sham to break it down into gibbs, enth and entro as .xpm plots. Your bins are the "pixels" or points in the landscape I believe, thus the binwidth is the xmin and max, but I amnot 100% sure how the numbering goes. You can check it against the image file generated (look at it with gedit, the pixel portion will tell you).</div>
<div> </div>
<div>One way to get it to work, would be to supply a free energy single column (or time, then single energy column), however someone might correct me on this (I hope) if I am wrong.</div>
<div> </div>
<div>However your origional question was to look up or coordinate a line graph gibbs energy plot, with a frame in your trajectory. Thus I believe you want to look at a linear plot of your energy change over time. If you generated the eigen vectors already, there are methods in these techniques to look at frames versus energy (total). You can also if you designed your experiment as a pulled/steered free energy experiment, use the g_wham, with no bootstraping or autocorrelation, etc...with one trajectory, and it is a sham. Otherwise, your using the wrong things to get what you want (g_ scripts in gromacs), and need to read the manual again. You can, with an index for only what you want, extract total energy and plot it as well versus frame, which might give you what you want also...using just g_energy.</div>
<div> </div>
<div>I hope that helps,</div>
<div>
<div> </div>
<div>Stephan Watkins</div>
<div name="quote" style="margin:10px 5px 5px 10px; padding: 10px 0 10px 10px; border-left:2px solid #C3D9E5; word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">
<div style="margin:0 0 10px 0;"><b>Gesendet:</b> Donnerstag, 19. Dezember 2013 um 00:09 Uhr<br/>
<b>Von:</b> jkrieger@mrc-lmb.cam.ac.uk<br/>
<b>An:</b> "gmx-users@gromacs.org" <gmx-users@gromacs.org><br/>
<b>Betreff:</b> Re: [gmx-users] Bin files of g_sham</div>
<div name="quoted-content">You're not asking it to read the bin index file in your command. Maybe that explains what you're seeing?<br/>
<br/>
On 18 Dec 2013, at 13:51, Ankita Naithani <ankitanaithani@gmail.com> wrote:<br/>
<br/>
> Hi Lloyd,<br/>
><br/>
> So, my input to g_sham is a file with following lines:<br/>
><br/>
> @ title "2D projection of trajectory"<br/>
> @ xaxis label "projection on eigenvector 1 (nm)"<br/>
> @ yaxis label "projection on eigenvector 2 (nm)"<br/>
> @TYPE xy<br/>
> 0.16750 0.18278<br/>
> -0.01448 -0.14834<br/>
> -0.32757 0.22259<br/>
> -0.33832 0.63817<br/>
><br/>
> This file is named AM_12.xvg<br/>
><br/>
> and my command line is g_sham -f AM_12.xvg -ls AM_FEL_12.xpm<br/>
><br/>
> My bin index file has the following lines:<br/>
><br/>
> [ 18 ]<br/>
> 7105<br/>
> [ 19 ]<br/>
> 7104<br/>
> [ 46 ]<br/>
> 7124<br/>
> [ 48 ]<br/>
> 7123<br/>
> [ 49 ]<br/>
> 7106<br/>
> [ 50 ]<br/>
> 6497<br/>
> 6498<br/>
> [ 51 ]<br/>
> 6499<br/>
> 6500<br/>
> [ 52 ]<br/>
> 6496<br/>
> 7103<br/>
> [ 53 ]<br/>
> 7102<br/>
> [ 55 ]<br/>
> 6672<br/>
> [ 75 ]<br/>
> 5745<br/>
> 5749<br/>
> and so on.<br/>
><br/>
> I am not sure how can I look at the default binwidth.<br/>
><br/>
> Kind regards,<br/>
><br/>
> Ankita<br/>
><br/>
><br/>
> On Tue, Dec 17, 2013 at 6:03 PM, lloyd riggs <lloyd.riggs@gmx.ch> wrote:<br/>
><br/>
>><br/>
>> You need to explain what you are doing better, what you want to do, and<br/>
>> show your input (command line, .xvg lines of input and output as example)<br/>
>> and you might get more help on here. If you read the -h from g_sham, it is<br/>
>> making a single plot from your input, and in this case your .xpm would be a<br/>
>> matrix with bins on a 2D surface with points taken across the entire run,<br/>
>> as an example of what I mean. I believe you can also use a single line<br/>
>> (time, energy as column 2), and get a single line out as well, with a<br/>
>> meaningless "landscape" still put out as .xpm, etc...look at your binwidth<br/>
>> default, input and what you want out is all I could say at this point. You<br/>
>> might be trying to do something else however...<br/>
>><br/>
>> Stephan<br/>
>> *Gesendet:* Dienstag, 17. Dezember 2013 um 08:38 Uhr<br/>
>><br/>
>> *Von:* "Ankita Naithani" <ankitanaithani@gmail.com><br/>
>> *An:* "Discussion list for GROMACS users" <gmx-users@gromacs.org><br/>
>> *Betreff:* Re: [gmx-users] Bin files of g_sham<br/>
>> I would like to add that maybe the bins store Frame numbers. Because I seem<br/>
>> to have around 21000 frames and I guess I am unable to see a number higher<br/>
>> than that. Apologies for the confusion about earlier number. So, in this<br/>
>> case if I want a confirmation at a particular frame number, how should I go<br/>
>> about that? Could you please help me with that?<br/>
>><br/>
>> Kind regards,<br/>
>><br/>
>> Ankita<br/>
>><br/>
>><br/>
>> On Tue, Dec 17, 2013 at 7:20 AM, Ankita Naithani<br/>
>> <ankitanaithani@gmail.com>wrote:<br/>
>><br/>
>>> Hi Llyod,<br/>
>>><br/>
>>> Thank you for your reply. I actually wanted to extract frames<br/>
>>> corresponding to the g_sham output. The input for g_sham in my case was<br/>
>>> column 1 had time points, column 2 and column 3 had Principal component 1<br/>
>>> and Principal component 2 respectively. So, the binning which has been<br/>
>>> done, i.e. if I say that the bin number 786 has the following few numbers<br/>
>>> like : 3876, 9870, 22002 etc. These should technically correspond to the<br/>
>>> time points right which I can then dump using trjconv too maybe.<br/>
>>><br/>
>>> Also, I am a little confused because if the bins have all the time<br/>
>> points,<br/>
>>> I don't seem to find the time points uptil the end of simulation in that<br/>
>>> case i.e. I am unable to see higher time scale. My simulation was for<br/>
>> 48ns<br/>
>>> and had 3fs time step. I am unable to trace anything beyond 22000 or<br/>
>> 23000<br/>
>>> number in my bins.<br/>
>>><br/>
>>> I hope you could help me solve this confusion.<br/>
>>><br/>
>>> Kind regards,<br/>
>>><br/>
>>> Ankita<br/>
>>><br/>
>>><br/>
>>> On Sun, Dec 15, 2013 at 11:17 PM, lloyd riggs <lloyd.riggs@gmx.ch><br/>
>> wrote:<br/>
>>><br/>
>>>><br/>
>>>> I will attempt an answer, however there might be a better response. It<br/>
>>>> depends on your input .xvg (what the first column is) and what you are<br/>
>>>> telling it to do (defualt is time), and if these correspond to your time<br/>
>>>> frames (can be timepoint or siply frame # I believe). g_sham seems (to<br/>
>> me)<br/>
>>>> to work like a comand line spread sheet manipulator with some built in<br/>
>>>> functions. If you want a particular frame from a trajectory time point,<br/>
>>>> then use g_traj with the first/last set to the frame you want but<br/>
>> there's<br/>
>>>> also 2 other tools that do this the same way...<br/>
>>>><br/>
>>>> Stephan Watkins<br/>
>>>> *Gesendet:* Sonntag, 15. Dezember 2013 um 14:49 Uhr<br/>
>>>> *Von:* "Ankita Naithani" <ankitanaithani@gmail.com><br/>
>>>> *An:* "Discussion list for GROMACS users" <gmx-users@gromacs.org><br/>
>>>> *Betreff:* [gmx-users] Bin files of g_sham<br/>
>><br/>
>>>> Hi,<br/>
>>>><br/>
>>>> I had a query regarding the bin files generated by g_sham.<br/>
>>>> Does the bin index contain timepoints? i.e. there are 1024 bins so each<br/>
>>>> bin<br/>
>>>> has a number of time points? And supposes bin 756 has numbers like 3876,<br/>
>>>> 20020, 7896 so these would all correspond to time and we could then dump<br/>
>>>> the representative frame corresponding to these time points?<br/>
>>>><br/>
>>>> I am a little confused as to what is stored in the bins? Timepoints or<br/>
>>>> frame number and then alternatively how to get the representative<br/>
>>>> frame/snapshot stored at that specific bin number.<br/>
>>>><br/>
>>>> Kind regards,<br/>
>>>><br/>
>>>> --<br/>
>>>> Ankita Naithani<br/>
>>>> --<br/>
>>>> Gromacs Users mailing list<br/>
>>>><br/>
>>>> * Please search the archive at<br/>
>>>> <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List</a> before<br/>
>>>> posting!<br/>
>>>><br/>
>>>> * Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br/>
>>>><br/>
>>>> * For (un)subscribe requests visit<br/>
>>>> <a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users" target="_blank">https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users</a> or<br/>
>>>> send a mail to gmx-users-request@gromacs.org.<br/>
>>>><br/>
>>>> --<br/>
>>>> Gromacs Users mailing list<br/>
>>>><br/>
>>>> * Please search the archive at<br/>
>>>> <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List</a> before<br/>
>>>> posting!<br/>
>>>><br/>
>>>> * Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br/>
>>>><br/>
>>>> * For (un)subscribe requests visit<br/>
>>>> <a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users" target="_blank">https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users</a> or<br/>
>>>> send a mail to gmx-users-request@gromacs.org.<br/>
>>><br/>
>>><br/>
>>> --<br/>
>>> Ankita Naithani<br/>
>><br/>
>><br/>
>><br/>
>> --<br/>
>> Ankita Naithani<br/>
>> --<br/>
>> Gromacs Users mailing list<br/>
>><br/>
>> * Please search the archive at<br/>
>> <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List</a> before<br/>
>> posting!<br/>
>><br/>
>> * Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br/>
>><br/>
>> * For (un)subscribe requests visit<br/>
>> <a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users" target="_blank">https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users</a> or<br/>
>> send a mail to gmx-users-request@gromacs.org.<br/>
>><br/>
>> --<br/>
>> Gromacs Users mailing list<br/>
>><br/>
>> * Please search the archive at<br/>
>> <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List</a> before<br/>
>> posting!<br/>
>><br/>
>> * Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br/>
>><br/>
>> * For (un)subscribe requests visit<br/>
>> <a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users" target="_blank">https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users</a> or<br/>
>> send a mail to gmx-users-request@gromacs.org.<br/>
><br/>
><br/>
> --<br/>
> Ankita Naithani<br/>
> --<br/>
> Gromacs Users mailing list<br/>
><br/>
> * Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List</a> before posting!<br/>
><br/>
> * Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br/>
><br/>
> * For (un)subscribe requests visit<br/>
> <a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users" target="_blank">https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users</a> or send a mail to gmx-users-request@gromacs.org.<br/>
--<br/>
Gromacs Users mailing list<br/>
<br/>
* Please search the archive at <a href="http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List" target="_blank">http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List</a> before posting!<br/>
<br/>
* Can't post? Read <a href="http://www.gromacs.org/Support/Mailing_Lists" target="_blank">http://www.gromacs.org/Support/Mailing_Lists</a><br/>
<br/>
* For (un)subscribe requests visit<br/>
<a href="https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users" target="_blank">https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users</a> or send a mail to gmx-users-request@gromacs.org.</div>
</div>
</div>
</div></div></body></html>