How to plot confident interval of ecdf() as shade?

Hello all,

I used ecdf() to plot my data, it calculates upper and lower band, but I don't know how to plot them like this:

a = rand(30,1);
[f1,x1,up,low] = ecdf(a);

I used some of file exchange function such as ciplot but it not work.

Anyone knows how I can achive such a plot? Thanks

Answer 1

[Star Strider]

Try this:

a = rand(30,1);
[f1,x1,up,low] = ecdf(a);
idx = ~isnan(up) & ~isnan(low);

figure
plot(x1(idx), f1(idx), '-r')
hold on
patch([x1(idx); flipud(x1(idx))], [up(idx); flipud(low(idx))], 'b', 'FaceAlpha',0.25, 'EdgeColor','none')
hold off
grid
legend('f1', '95% Confidence Interval', 'Location','SE')

The data are random, so for example for one run:

.

[BN]

Thank you so much. It's awesome.

[Star Strider]

As always, my pleasure!