# How to shade area between standard deviations in a plot?

24 views (last 30 days)
Tomaszzz on 28 Mar 2022
Edited: MANDRAKE on 18 Aug 2022
Hi all,
I would like to shade the area between standard deviations.
The data is attached. My code is as follows:
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

Jan on 28 Mar 2022
Edited: Jan on 18 Aug 2022
m = data1.mean;
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.
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 on 18 Aug 2022
@MANDRAKE: Thanks for this comment. I've inserted some code to ensure, that the data are column vectors.