MATLAB: Confidence band around linear least-squares line

Question

I am using lsline in Matlab 2016b to add a linear least-squares line to a scatter plot. I would like to add a 95% confidence band around that fit line, such that it looks like this (plot is made with the python function seaborn):

However, lsline returns no fit parameters from which to construct the 95% confidence band, and the only Matlab function I could find that does return these, is nlpredci, but that function is used for something else (nonlinear regression prediction)

Answer 1

One option is the Statistics and Machine Learning Toolbox fitlm (link) function. Then use the predict function (linked at the end of that page).

Another option is to sue polyfit and polyval with the File Exchange polypredci function, or the Statistics and Machine Learning Toolbox polyconf (link) function.