How can I draw 95% CI plot in my data?

Question

Bluemin 2021-3-26

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/783991-how-can-i-draw-95-ci-plot-in-my-data

评论： Muhammad Sajjad 2022-10-4

采纳的回答： Star Strider

I have 1 data(100x1 matrix).

In Matlab, I want to draw 95% ci plot in my data

so, I found good code

like this

x = 1:100; % Create Independent Variable

y = randn(50,100); % Create Dependent Variable ‘Experiments’ Data

N = size(y,1); % Number of ‘Experiments’ In Data Set

yMean = mean(y); % Mean Of All Experiments At Each Value Of ‘x’

ySEM = std(y)/sqrt(N); % Compute ‘Standard Error Of The Mean’ Of All Experiments At Each Value Of ‘x’

CI95 = tinv([0.025 0.975], N-1); % Calculate 95% Probability Intervals Of t-Distribution

yCI95 = bsxfun(@times, ySEM, CI95(:)); % Calculate 95% Confidence Intervals Of All Experiments At Each Value Of ‘x’

figure

plot(x, yMean,'p') % Plot Mean Of All Experiments

hold on

plot(x, yCI95+yMean,'-r') % Plot 95% Confidence Intervals Of All Experiments

hold off

grid

but this graph is not what I wanted

How can i draw graph? like this..

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Muhammad Sajjad 2022-10-4

how to plot median confidence interval

请先登录，再进行评论。

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

Answer 1

Star Strider 2021-3-26

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/783991-how-can-i-draw-95-ci-plot-in-my-data#answer_658611

在 MATLAB Online 中打开

Replace the second plot call with:

patch([x, fliplr(x)], [yCI95(1,:) fliplr(yCI95(2,:))], 'b', 'EdgeColor','none', 'FaceAlpha',0.25)

so the full code is now:

x = 1:100;                                          % Create Independent Variable
y = randn(50,100);                                  % Create Dependent Variable ‘Experiments’ Data
N = size(y,1);                                      % Number of ‘Experiments’ In Data Set
yMean = mean(y);                                    % Mean Of All Experiments At Each Value Of ‘x’
ySEM = std(y)/sqrt(N);                              % Compute ‘Standard Error Of The Mean’ Of All Experiments At Each Value Of ‘x’
CI95 = tinv([0.025 0.975], N-1);                    % Calculate 95% Probability Intervals Of t-Distribution
yCI95 = bsxfun(@times, ySEM, CI95(:));              % Calculate 95% Confidence Intervals Of All Experiments At Each Value Of ‘x’
figure
plot(x, yMean,'p')                                  % Plot Mean Of All Experiments
hold on
% plot(x, yCI95+yMean,'-r')                           % Plot 95% Confidence Intervals Of All Experiments
patch([x, fliplr(x)], [yCI95(1,:) fliplr(yCI95(2,:))], 'b', 'EdgeColor','none', 'FaceAlpha',0.25)
hold off
grid