This may be caused by the different length of data in the "pl.plot()".
Please add the following line before "for i in modes:" to check the dimension of "SymP_LOE", "SymN_LOE", and so on:
print SymP_LOE.shape()
According to the manual:
http://docs.quantumwise.com/manuals/Types/InelasticTransmissionSpectrum/InelasticTransmissionSpectrum.html#NL.Analysis.InelasticTransmissionSpectrum.InelasticTransmissionSpectrum.inelasticTransmissionAsymmetricNegativeBiasThe shape of "SymP_LOE" should be (number_of_energies, number_of_kpoints, number_of_qpoints, number_of_spins) for mode = integer.
pl.plot(E,SymP_LOE[i,:,0,0],'-',label='mode %d, + (LOE)' %i)
So if you change it to following line, it may work:
pl.plot(E,SymP_LOE[i,:,0,0,0],'-',label='mode %d, + (LOE)' %i)
The same changes are applied to other three lines:
pl.plot(E,AsymP_LOE[i,:,0,0,0],'-',label='mode %d, + (LOE)' %i)
pl.plot(E,SymP_XLOE[i,:,0,0,0],'-',label='mode %d, + (XLOE)' %i)
pl.plot(E,AsymP_XLOE[i,:,0,0,0],'-',label='mode %d, + (XLOE)' %i)
Please try it.