In [3]:

    

%matplotlib inline
import numpy as np
import spectraplotpy as spp
import matplotlib.pyplot as plt
from glob import glob
#for easier development do autoreload
%load_ext autoreload
%autoreload 1
%aimport spectraplotpy


    

In [4]:

    

#load files
files = glob('../sampledata/01-CD-Aviv62DS/PEP-average/4RNSX.*')
spectra = [ spp.Spectrum(spp.AvivImporter(f).dataset) for f in files]

# Two subplots, unpack the axes array immediately
f, (ax1, ax2) = plt.subplots(2, 1, sharex=True)

for s in spectra:
    s.plot(axes=ax1)
avg = spp.average_spectra(*spectra, error_type='st_dev')
avg.errorbar(axes=ax2)


    

    
Out[4]:





<Container object of 3 artists>






    

In [7]:

    

f, a = plt.subplots()
avg.errorbar(axes=a)

f, a = plt.subplots()
avg.errorfill(axes=a, alpha_fill=0.2)