I want to plot the time serie data of number of hydrogen bond (between protein-ligand) from a MD simulation with Gromacs as in this picture:

http://igem.org/wiki/images/3/3e/Zju_model1_10.png

Instead of the classic Grace style.

https://www.researchgate.net/profile/Trevor_Marshall2/publication/305254443/figure/fig2/AS:391130341494799@1470264148666/A-plot-from-the-Gromacs-g-hbond-software-of-the-instantaneous-number-of-hydrogen-bonds.png

The first plot seems to come from R. Trie it, but all I could get is that:

R_hbond.PNG (see attached)

Using:

qplot(x, y, xlab = 'Time (ns)', ylab = 'Number of Hhonds', geom = c("path", "smooth"))

x : time

y: # h bonds (see attached file Lig_H_bond.txt)

Tried different geom options without success.

It would really nice if someone can help me out.

Thanks a lot!

Similar questions and discussions