Warm tip: This article is reproduced from stackoverflow.com, please click
data-visualization matlab plot

Shaded plot in Matlab

发布于 2021-02-07 15:55:43

I would like to plot a function in Matlab with a shaded area indicating the uncertainty over it (e.g., confidence interval). This can be achieved by using the fill function to create a color patch. For example

x = linspace(0, 2*pi, 100);
f = cos(x);
fUp = cos(x) + 1;
fLow = cos(x) - 1;

x2 = [x, fliplr(x)];

plot(x, f, 'k')
hold on
fill(x2, [f, fliplr(fUp)], 0.7 * ones(1, 3), 'linestyle', 'none', 'facealpha', 0.4);
fill(x2, [fLow, fliplr(f)], 0.7 * ones(1, 3), 'linestyle', 'none', 'facealpha', 0.4);

This creates a shaded gray area between the functions fLow and fUp, with f in the middle represented as a solid black line, as in the picture below.

enter image description here

I would like now to have the shaded area degrade its color when we approach the lower (resp. upper) bound of the confidence interval. In particular, I would like that while approaching its boundaries, the shaded area gets brighter and brighter. Is there a way to do it? I'm doing two separate patches because I think it may be necessary for my purpose.

Questioner
G. Gare
Viewed
0
分享
obchardon 2020-10-13 17:24

You can split your CI into n subarea:

x = linspace(0, 2*pi, 100);
f = cos(x);
n = 20;  % step number
g = 0.3; % grayscale intensity
fUp = cos(x) + linspace(0,1,n).';
fLow = cos(x) - linspace(0,1,n).';

x2 = [x, fliplr(x)];

plot(x, f, 'k')
hold on
fill(x2, [repmat(f,n,1), fliplr(fUp)], g * ones(1, 3), 'linestyle', 'none', 'facealpha', [1/n]);
fill(x2, [fLow, repmat(fliplr(f),n,1)], g * ones(1, 3), 'linestyle', 'none', 'facealpha', [1/n]);

Which produce:

enter image description here

The subarea are overlapping and produce a maximum facealpha of n*(1/n) * g = g

Noticed that this method is not really memory efficient (since it produce n subarea on each side) and will only works with a linear shading.

If your CI is non linear then you should adjust this part:

% Prediction              Linear CI
%     ↓                       ↓
    cos(x)       +     linspace(0,1,n).';
    cos(x)       -     linspace(0,1,n).';

to

% Prediction                   Non linear CI
%     ↓                              ↓
    cos(x)       +     your_non_linear_CI_distribution;
    cos(x)       -     your_non_linear_CI_distribution;
G. Gare 1970-01-01 08:00:00

Cool! Perfect! As you said, It required a bit of extra work if the width of the confidence interval depends on x but otherwise it's great.

Related issues

1
How to find differential of input function while writing in MATLAB for Newton Raphson method
2
Does Matlab's new function argument validation allow for unmatched name-value arguments?
3
Defining a matrix by avoiding the use of for loops
4
How to set the DPI of Java Swing apps on Windows/Linux?
5
How is stereo rectification implemented?
6
How to unfold a Matrix on Matlab?
7
Inconsistent row/column dimensions error in Scilab
8
Matlab Structures: How can I pass a concatenated string as a field to a structure?
9
Matlab Structures: How to compare the contents of structure? Getting the following error
10
adaptive elliptical structuring element in MATLAB