error bar on the graph

Question

ahmad albngali 2020-4-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/519674-error-bar-on-the-graph

评论： Star Strider 2020-6-22

采纳的回答： Star Strider

hi

after i have plot between x y and z ( scatter3(x, y, z))

x=[2.6,5.2,14,23.3,28.3,2.6,5.2,14,23.3,28.3,2.6,5.2,14,23.3,28.3];

y=[3.5,6.9,17.7,29.1,35.2,3.5,6.8,17.5,28.7,34.8,3.4,6.6,16.9,27.6,33.5];

z=[450,450,450,450,450,350,350,350,350,350,250,250,250,250,250];

i want to make error bar for y in the same graph

and i have the standard deviation for y

which is ( 0.047140452, 0.047140452, 0.047140452, 0.124721913, 0.047140452, 0.047140452, 0.047140452, 0.047140452, 0.124721913, 0.124721913, 0.047140452, 0.047140452, 0.047140452, 0.124721913, 0.124721913)

so could you tell me what is the code for that

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Star Strider 2020-4-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/519674-error-bar-on-the-graph#answer_427464

在 MATLAB Online 中打开

The error bars are difficult to see, so I multiplied them by 10 here:

x=[2.6,5.2,14,23.3,28.3,2.6,5.2,14,23.3,28.3,2.6,5.2,14,23.3,28.3];	
y=[3.5,6.9,17.7,29.1,35.2,3.5,6.8,17.5,28.7,34.8,3.4,6.6,16.9,27.6,33.5];
z=[450,450,450,450,450,350,350,350,350,350,250,250,250,250,250];
ysd = [0.047140452, 0.047140452, 0.047140452, 0.124721913, 0.047140452, 0.047140452, 0.047140452, 0.047140452, 0.124721913, 0.124721913, 0.047140452, 0.047140452, 0.047140452, 0.124721913, 0.124721913];
xr = reshape(x, [], 3);
yr = reshape(y, [], 3);
zr = reshape(z, [], 3);
ysdr = reshape(ysd, [], 3)*10;
figure
hold on
for k1 = 1:size(xr,2)
    for k2 = 1:size(xr,1)
    plot3((xr(k2,k1)*[1 1]), (yr(k2,k1)+ysdr(k2,k1)*[-1 1]), (yr(2,k1)*[1 1]), '-r')
    plot3((xr(k2,k1)), (yr(k2,k1)), (yr(2,k1)), '.r')
    end
end
hold off
view(-60,30)
grid on
xlabel('X')
ylabel('Y')
zlabel('Z')
% legend('Y\pm10\sigma')

Experiment to get the result you want.

16 个评论
显示 14更早的评论隐藏 14更早的评论

Star Strider 2020-4-23

在 MATLAB Online 中打开

With the extrapolated values, reshape will not work (added points must correspond to the existing (5x3) format), so the code needs to be rearranged slightly and the extrapolated values plotted separately, after the plotting loop:

x=[2.6,5.2,14,23.3,28.3,2.6,5.2,14,23.3,28.3,2.6,5.2,14,23.3,28.3];	
y=[3.5,6.9,17.7,29.1,35.2,3.5,6.8,17.5,28.7,34.8,3.4,6.6,16.9,27.6,33.5];
z=[450,450,450,450,450,350,350,350,350,350,250,250,250,250,250];
ysd = [0.047140452, 0.047140452, 0.047140452, 0.124721913, 0.047140452, 0.047140452, 0.047140452, 0.047140452, 0.124721913, 0.124721913, 0.047140452, 0.047140452, 0.047140452, 0.124721913, 0.124721913];
xr = reshape(x, [], 3);
yr = reshape(y, [], 3);
zr = reshape(z, [], 3);
ysdr = reshape(ysd, [], 3)*10;
YfromXZ = scatteredInterpolant(x(:), z(:), y(:));
queryx = 20;    
queryz = 300;  
YfromXZ.ExtrapolationMethod = 'linear';
matchingy = YfromXZ(queryx, queryz);
figure;scatter3(x, y, z); hold('on'); plot3(queryx, matchingy, queryz, 'r*');
errYfromXZ = scatteredInterpolant(x(:), z(:), y(:)+ysd(:));
errYfromXZ.ExtrapolationMethod = 'linear';
matchingyerr = errYfromXZ(queryx, queryz);
errY = (matchingyerr - matchingy)
figure
hold on
for k1 = 1:size(xr,2)
    for k2 = 1:size(xr,1)
    plot3((xr(k2,k1)*[1 1]), (yr(k2,k1)+ysdr(k2,k1)*[-1 1]), (yr(2,k1)*[1 1]), '-r')
    plot3((xr(k2,k1)), (yr(k2,k1)), (yr(2,k1)), '.r')
    end
end
plot3(queryx*[1 1], matchingy+errY*[-1 1]*10, queryz*[1 1], '-g')                       % <- Added
plot3(queryx, matchingy, queryz, '.g')                                                  % <- Added
hold off
view(-60,30)
grid on
xlabel('X')
ylabel('Y')
zlabel('Z')

The extrapolated point and errorbar are plotted in green to make them more visible. (All the errorbars are multiplied by 10 because they would not otherwise be visible.)

.

ahmad albngali 2020-6-16

在 MATLAB Online 中打开

hi as you remember i used this code

This code was used to find the Relationship of DEq, CTDIvol and scanning length, also the code was used to show the error bar for DEqvalues (The error bars for all DEqvalues are difficult to see in the graph , so I multiplied them by 5)

% x= CTDIv
x=[2.6,5.2,14,23.3,28.3,2.6,5.2,14,23.3,28.3,2.6,5.2,14,23.3,28.3 ];    
% y= planar average equilibrium dose (DEq)
y=[3.5,6.9,17.7,29.1,35.2,3.5,6.8,17.5,28.7,34.8,3.4,6.6,16.9,27.6,33.5 ];
% z= scanning length 
z=[450,450,450,450,450,350,350,350,350,350,250,250,250,250,250];
% ysd= DEq standard deviation 
ysd = [0.047140452, 0.047140452, 0.047140452, 0.124721913, 0.047140452, 0.047140452, 0.047140452, 0.047140452, 0.124721913, 0.124721913, 0.047140452, 0.047140452, 0.047140452, 0.124721913, 0.124721913];
xr = reshape(x, [], 3);
yr = reshape(y, [], 3);
zr = reshape(z, [], 3);
ysdr = reshape(ysd, [], 3)*5;
 
%for visualisation
figure
hold on
for k1 = 1:size(xr,2)
    for k2 = 1:size(xr,1)
    plot3((xr(k2,k1)*[1 1]), (yr(k2,k1)+ysdr(k2,k1)*[-1 1]), (zr(2,k1)*[1 1]), '-r')
    plot3((xr(k2,k1)), (yr(k2,k1)), (zr(2,k1)), '.r')
    end
end
hold off
view(-60,30)
grid on

so my question

the Error bars are too small. and i want to Increase the dot size and fit a grid or lines to the dots ??

ahmad albngali 2020-6-22

hi i try to do it but it is confusing

could you please show me hoe i do it in my examble ( in my code )

Star Strider 2020-6-22

在 MATLAB Online 中打开

Add them to the plot3 calls:

plot3((xr(k2,k1)*[1 1]), (yr(k2,k1)+ysdr(k2,k1)*[-1 1]), (zr(2,k1)*[1 1]), '-r', 'LineWidth',1.5)
plot3((xr(k2,k1)), (yr(k2,k1)), (zr(2,k1)), '.r', 'MarkerSize',10)

so the complete code is now:

% x= CTDIv
x=[2.6,5.2,14,23.3,28.3,2.6,5.2,14,23.3,28.3,2.6,5.2,14,23.3,28.3 ];    
% y= planar average equilibrium dose (DEq)
y=[3.5,6.9,17.7,29.1,35.2,3.5,6.8,17.5,28.7,34.8,3.4,6.6,16.9,27.6,33.5 ];
% z= scanning length 
z=[450,450,450,450,450,350,350,350,350,350,250,250,250,250,250];
% ysd= DEq standard deviation 
ysd = [0.047140452, 0.047140452, 0.047140452, 0.124721913, 0.047140452, 0.047140452, 0.047140452, 0.047140452, 0.124721913, 0.124721913, 0.047140452, 0.047140452, 0.047140452, 0.124721913, 0.124721913];
xr = reshape(x, [], 3);
yr = reshape(y, [], 3);
zr = reshape(z, [], 3);
ysdr = reshape(ysd, [], 3)*5;
 
%for visualisation
figure
hold on
for k1 = 1:size(xr,2)
    for k2 = 1:size(xr,1)
    plot3((xr(k2,k1)*[1 1]), (yr(k2,k1)+ysdr(k2,k1)*[-1 1]), (zr(2,k1)*[1 1]), '-r', 'LineWidth',1.5)
    plot3((xr(k2,k1)), (yr(k2,k1)), (zr(2,k1)), '.r', 'MarkerSize',10)
    end
end
hold off
view(-60,30)
grid on

Change the values to create the result you want.

.

请先登录，再进行评论。

error bar on the graph

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

16 个评论
显示 14更早的评论隐藏 14更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

error bar on the graph

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

采纳的回答

16 个评论 显示 14更早的评论隐藏 14更早的评论

更多回答（0 个）

另请参阅

类别

标签

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

16 个评论
显示 14更早的评论隐藏 14更早的评论