You are in luck in not using categorical variables, because this approach will not work with them (although I’ve not had the opportunity to read through the complete R2019a Release Notes yet, so there may be options I’m not aware of).
Try this:
y=[0.095 0.037 0.132; 0.1 0.058 0.158; 0 0 0.02; 0 0 0]
StdErr = rand(3,4)*0.01; % Create ‘Standard Error’ Matrix
figure
hBar = bar(y, 0.8); % Return ‘bar’ Handle
for k1 = 1:size(y,2)
ctr(k1,:) = bsxfun(@plus, hBar(k1).XData, hBar(k1).XOffset'); % Note: ‘XOffset’ Is An Undocumented Feature, This Selects The ‘bar’ Centres
ydt(k1,:) = hBar(k1).YData; % Individual Bar Heights
end
hold on
errorbar(ctr, ydt, StdErr, '.r') % Plot Error Bars
set(gca, 'XTickLabel', {'Rural' 'Exurban' 'Suburban' 'Urban'});
You did not supply your ‘Standard Error’ matrix, so I created one. Note that it must be the same size as the one I created, in order for this code to work.
