Confidence band around linear least-squares line

Question

z8080 2018-4-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/392636-confidence-band-around-linear-least-squares-line

评论： z8080 2018-4-5

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)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2018-4-4

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/392636-confidence-band-around-linear-least-squares-line#answer_313489

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.

7 个评论
显示 5更早的评论隐藏 5更早的评论

Star Strider 2018-4-4

在 MATLAB Online 中打开

First, this will plot your data and produce a plot similar to the one you posted:

x = [62  97  85  78  79  68  91  98  89  54  61  93  76  51  105  24  62  94  92  91  94  81  73  101  63  95  79  81  81  68 ];
y = [1.53  1.35  1.20  1.06  1.23  0.55  0.58  1.38  1.63  0.98  1.63  1.21  1.10  1.03  0.85  0.47  0.85  1.78  2.23  1.55  0.80  1.35  1.25  1.45  1.10  1.38  1.00  1.00  0.79  1.63 ];
[p,S] = polyfit(x,y,1);
xv = linspace(min(x), max(x), 150);
[y_ext,delta] = polyconf(p,xv,S);
figure
plot(x, y, 'pg', 'MarkerFaceColor','g')
hold on
plot(xv, y_ext, '-g')
patch([xv fliplr(xv)], [(y_ext+delta) fliplr((y_ext-delta))], 'g', 'FaceAlpha',0.1, 'EdgeColor','none')
hold off

Experiment with the colours and other arguments to get the result you want.

Second, it should be symmetrical, since that is the way the confidence interval with symmetric distributions (such as the normal and t-distributions) are defined. Also, with respect to symmetry, it will be symmetrical about the (mean(x),mean(y)) point. It diverges as the distance from the means increases.

Star Strider 2018-4-5

As always my pleasure.

The width isn’t constant. It may look that way if you use only the original ‘x’ values. However if you use my code (with the ‘xv’ vector) you can easily see that it is not constant, and diverges appropriately from the mean values to the extremes.

It’s not your inexperience. When I plot a nonlinear regression, confidence intervals, or anything else that needs nigher-than normal resolution to depict correctly, I use something like ‘xv’ instead of only using the original *‘x’*values, unless they are sufficiently high resolution. The image you posted appears to have used the same sort of approach.

z8080 2018-4-5

Makes perfect sense now.

请先登录，再进行评论。

Confidence band around linear least-squares line

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

7 个评论
显示 5更早的评论隐藏 5更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

Community Treasure Hunt

Confidence band around linear least-squares line

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

7 个评论 显示 5更早的评论隐藏 5更早的评论

更多回答（0 个）

另请参阅

类别

标签

产品

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

7 个评论
显示 5更早的评论隐藏 5更早的评论