Uncertainties

[1]:
import smpl
from smpl import plot
from smpl import io
from smpl import functions as f
import numpy as np
np.random.seed(1337)
print(smpl.__version__)
1.3.0.11
[2]:
import uncertainties.unumpy as unp
data = np.loadtxt(io.find_file('test_linear_data.txt',3))
xdata = data[:,0]
xerr = data[:,2]
ydata = data[:,1]
yerr = data[:,3]
x = unp.uarray(xdata,xerr)
y = unp.uarray(ydata,yerr)
[3]:
plot.data(x,y,label="data",fmt=None)
plot.data(x,y,label="data",fmt="step",init=True)
plot.data(x,y,label="data",fmt="hist",init=True)
../../../_images/example_plot_uncertainties_plot_uncertainties_3_0.png
../../../_images/example_plot_uncertainties_plot_uncertainties_3_1.png
../../../_images/example_plot_uncertainties_plot_uncertainties_3_2.png
[4]:
plot.data(plot.unv(x),plot.unv(y),label="data",fmt="hist")
../../../_images/example_plot_uncertainties_plot_uncertainties_4_0.png
[5]:
x= np.random.randn(1000)
print(isinstance(x, (list, tuple, np.ndarray)))
plot.data(x,bins=20,label="data",fmt=None)
plot.data(x,bins=20,label="data",fmt="step",init=True)

True
../../../_images/example_plot_uncertainties_plot_uncertainties_5_1.png
../../../_images/example_plot_uncertainties_plot_uncertainties_5_2.png

Visually increased uncertainties by the scaling given by data_sigmas.

[6]:
plot.data(x,bins=20,data_sigmas=5,label="data",fmt=None)
plot.data(x,bins=20,data_sigmas=5,label="data",fmt="step",init=True)
../../../_images/example_plot_uncertainties_plot_uncertainties_7_0.png
../../../_images/example_plot_uncertainties_plot_uncertainties_7_1.png
[ ]: