Residue Plot

from smpl import plot
from smpl import io
from smpl import functions as f
import numpy as np
import uncertainties.unumpy as unp
import matplotlib.pyplot as plt
import smpl
smpl.__version__

dat = np.loadtxt(io.find_file("dat.dat",3), skiprows=0, delimiter=" ")


xdata = dat[:,0]
ydata = dat[:,1]
ymodel = dat[:,-2]

plot.data(xdata,ydata,logy=True,ss=False)
plot.data(xdata,ymodel,fmt="-",logy=True,init=False,label="model")

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)

plot.fit(x,y,function=f.line,residue=True,xaxis="t",yaxis="s",units=["l","b"],sigmas=1)

plot.fit(x,y,function=f.line,residue=True,xaxis="t",yaxis="s",units=["l","b"],sigmas=1,residue_err=False)