MmPbSaStat.py error

[Sajad Ahrari]

Dear Users,

I am trying to run "MmPbSaStat.py" after running the "g_mmpbsa -s ../../5es1_100ns_nowater.tpr -f  5es1_5RC_last5ns_dt80.xtc  -i pbsa.mdp -n ../../../MD_preparation/index.ndx -pdie 2 -pbsa -decomp" command. However, I was faced with the following error. Could you please help me know what I am missing? my input files are attached.

Thank you


D:\MD-Lab\Ahrari\MARK4\MARK4\analysis\inactive_drug_interaction\docking_AutoDock \5es1\1\MD-out\analysis\MMPBSA>python MmPbSaStat.py -m energy_MM.xvg -p polar.xv g -a apolar.xvg
Traceback (most recent call last):
  File "MmPbSaStat.py", line 303, in <module>
    main()
  File "MmPbSaStat.py", line 62, in main
    cTmp.CalcEnergy(args,frame_wise,0)
  File "MmPbSaStat.py", line 81, in CalcEnergy
    polEn = ReadData(self.PolFile,n=4)
  File "MmPbSaStat.py", line 245, in ReadData
    x_temp.append(float(data[i][j]))
IndexError: index 2 is out of bounds for axis 0 with size 2

[g_mmpbsa Mailing List]

polar solvation energy is incomplete at last frame. If you open polar.xvg in any text editor, you can see that polar energy is not calculated for ligand and complex at time 95000 ps. Also, for apolar, no calculation is performed for last frame (95000 ps). Number of frames should be same in all three input files.

[Khaled Darwish]

I am now trying to calculate the average and get the summary of free binding energy but this message appears

python MmPbSaStat.py -m energy_MM.xvg -p polar.xvg -a sasa.xvg

Traceback (most recent call last):

  File "MmPbSaStat.py", line 44, in <module>
    import numpy as np
ImportError: No module named numpy

Can you help me with this?