How to shade area between standard deviations in a plot?

22 views (last 30 days)
Hi all,
I would like to shade the area between standard deviations.
The data is attached. My code is as follows:
load 'mean'
load 'std';
bottom_sd=(mean-sd);
top_sd=(mean+sd);
figure(1)
plot(mean,(bottom_sd),'r');
hold on;
plot(mean,(top_sd),'r');
patch([mean(:); flipud(mean(:))], [mean(:)-sd(:); flipud(mean(:)+sd(:))], [0.6 0.7 0.8])
hold off
But it gives me this. Cold you please help?

Accepted Answer

Jan
Jan on 28 Mar 2022
Edited: Jan on 18 Aug 2022
data1 = load('mean.mat');
m = data1.mean;
data2 = loadt('std.mat');
s = data2.sd; % Do not use "mean" and "std" as variables!
x = (1:numel(m)).';
m = m(:); % [EDITED] Be sure that the data are column vectors
s = s(:); % [EDITED]
figure;
axes('NextPlot', 'add'); % as: hold on
size(x)
size(low)
patch([x; flip(x)], [m + s; flip(m - s)], [0.6 0.7 0.8])
plot(x, m);
Your code mixes x and y coordinates.
  2 Comments
MANDRAKE
MANDRAKE on 18 Aug 2022
Edited: MANDRAKE on 18 Aug 2022
Thanks a lot Jan.
The line with the patch function should be something like
patch([x; flip(x)], [m + s, flip(m - s)], [0.6 0.7 0.8])
to avoid the possible error about having two vectors of different lengths.
Jan
Jan on 18 Aug 2022
@MANDRAKE: Thanks for this comment. I've inserted some code to ensure, that the data are column vectors.

Sign in to comment.

More Answers (0)

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!