For use with null_lm
null_lm
resid_sigma(model, data, sigma = 1)
model
to extract residuals from
data
used to fit model
sigma
a specific sigma to model