In [1]:
# %load Plot.py
import pandas as pd
import matplotlib.pyplot as plt
import matplotlib 
%matplotlib inline

f, (ax1, ax2) = plt.subplots(1, 2)
#plt.suptitle('OPLS-AA/M in NAMD vs GMX')
df = pd.read_csv('ANALYSIS_ZZZ.csv',delim_whitespace=True)
df[['RES','BOND','ANGLE','DIHEDRAL','IMPROPER','VDWL']].plot(kind='box',ax=ax1,rot=45)
ax1.set_ylabel('GMX - NAMD')
df[['RES','ELEC','TOTAL']].plot(kind='box',ax=ax2,rot=41)
ax2.set_ylabel('GMX - NAMD')
#plt.show()
plt.tight_layout()
plt.savefig('RES_ZZZ.pdf')



In [6]:
df.sort_values(['BOND']).tail(1)


Out[6]:
RES BOND ANGLE DIHEDRAL IMPROPER VDWL ELEC TOTAL
8 HP 0.0843 -0.00797 -0.00323 0.00623 -0.02116 1.37902 1.4372

In [7]:
df.sort_values(['ANGLE']).tail(1)


Out[7]:
RES BOND ANGLE DIHEDRAL IMPROPER VDWL ELEC TOTAL
9 I -0.03197 0.04611 0.01671 0.00193 -0.00596 -0.68315 -0.65624

In [8]:
df.sort_values(['DIHEDRAL']).tail(1)


Out[8]:
RES BOND ANGLE DIHEDRAL IMPROPER VDWL ELEC TOTAL
9 I -0.03197 0.04611 0.01671 0.00193 -0.00596 -0.68315 -0.65624

In [9]:
df.sort_values(['IMPROPER']).tail(1)


Out[9]:
RES BOND ANGLE DIHEDRAL IMPROPER VDWL ELEC TOTAL
12 M 0.02584 -0.00613 0.01191 0.00849 0.01041 -0.80623 -0.75573

In [10]:
df.sort_values(['VDWL']).tail(1)


Out[10]:
RES BOND ANGLE DIHEDRAL IMPROPER VDWL ELEC TOTAL
14 Q -0.03312 0.00527 0.00728 0.0041 0.02398 -1.00903 -1.00152

In [ ]: