Plotyy and error bars

Question

Clement Wong 2011-8-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/12946-plotyy-and-error-bars

I'd like to do a plotyy function, and I have confidence intervals for each point of my data. I can't seem to figure out a way to do both of these at once. My x data is stored in an array RC, and the two y data are in arrays L and D. The confidence intervals are in ciL and ciD, and are not necessarily symmetric. Is there an easy way to add these as error bars to my plot?

2 个评论
显示无隐藏无

Jan 2011-8-3

If you post the code you have currently, it would be easier to insert the changes.

Clement Wong 2011-8-3

The code looks something like this:

x = 1;

for rc = 1000:1000:10000

[beta,r,j,covb,mse] = nlfin(parameters including rc);

ci = nlparci(beta,r,'covar',covb);

L(x) = beta(2);

D(x) = beta(3);

ciL(x,:) = ci(2,:);

ciD(x,:) = ci(3,:);

x = x+1;

end

L,D,ciL,ciD are all preallocated for speed. But the idea is that now, I have the confidence intervals for all 10 data points as I sweep rc. RC goes from 1,000 to 10,000 in 1,000 step size. I'm plotting the values for L and D against RC using plotyy(RC,L,RC,D,semilogx).

ciL and ciD contain the confidence intervals for each point, in order. That is, if L(1) is 5, then ciL(1) could be something like [3 8]. I hope this makes sense.

请先登录，再进行评论。

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

Answer 1

Patrick Kalita 2011-8-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/12946-plotyy-and-error-bars#answer_17768

在 MATLAB Online 中打开

When customizing graphs made with plotyy, it is usually useful to store axes handles that it provides as its return value. With those axes handles you can (1) turn hold mode on for each one and (2) add errorbars to each.

Here's an example:

% Make up some data to plot
x = 1:5;
y1 = rand(1,5);
y2 = rand(1,5) + 100;
L1 = 0.01 * ones(1,5);
U1 = 0.02 * ones(1,5);
L2 = 0.03 * ones(1,5);
U2 = 0.04 * ones(1,5);
% Store the axes handles produced by PLOTYY
ax = plotyy(x, y1, x, y2);
% Use HOLD and ERRORBAR, passing axes handles to the functions.
hold(ax(1), 'on');
errorbar(ax(1), x, y1, L1, U1);
hold(ax(2), 'on');
errorbar(ax(2), x, y2, L2, U2);

2 个评论
显示无隐藏无

Clement Wong 2011-8-3

I have one question. For the errorbars, do I have to convert my confidence interval (ci) into amount below and amount above the given value? That is, say my data point is 5, and my ci is [3 9]. Then would I need to make L = 2 (5 - 3 = 2), and U = 4 (9 - 5 = 4)? That's how it seems to work. Is there any way to just input a range for the errobars? I want to be able to input the code errorbar(x,y,ci(1,:),ci(2,:)), because the nlparci method that I'm using to create my ci's give them in that form. It's not much trouble to change it to the other form, but it would be easier this way, if there's a way to do it.

Thanks again!