from spectrum import parma, marple_data
p = parma(marple_data, 15, 15, 30, NFFT=4096)
p.plot(norm=True, color='red', linewidth=2) # same options as pylab.plot