How to assign individual colors to grouped and stacked elements in bar chart?

Hi all!

%Öffnen des benötigten workspace:
load('sensitivity_analysis.mat');

Data = zeros(9,3,3);

Data(:,:,1) = production_wo_elec;
Data(:,:,2) = production_variable;
Data(:,:,3) = usephase;

groupLabels = { 'hc small', 'hc medium', 'hc large', 'mix small', 'mix medium', 'mix large','green small','green medium','green large'};     % set labels
plotBarStackGroups(Data, groupLabels); % plot groups of stacked bars

function [] = plotBarStackGroups(Data, groupLabels)
%%Plot a set of stacked bars, but group them according to labels provided.
%%Params:

stackData is a 3D matrix (i.e., stackData(i, j, k) => (Group, Stack, StackElement))

groupLabels is a CELL type (i.e., { 'a', 1 , 20, 'because' };)

%%Copyright 2011 Evan Bollig (bollig at scs DOT fsu ANOTHERDOT edu

NumGroupsPerAxis = size(Data, 1);
NumStacksPerGroup = size(Data, 2);

% Count off the number of bins
groupBins = 1:NumGroupsPerAxis;
MaxGroupWidth = 0.65; % Fraction of 1. If 1, then we have all bars in groups touching
groupOffset = MaxGroupWidth/NumStacksPerGroup;
figure
    hold on; 
for i=1:NumStacksPerGroup

      Y = squeeze(Data(:,i,:));

      % Center the bars:

      internalPosCount = i - ((NumStacksPerGroup+1) / 2);

      % Offset the group draw positions:
      groupDrawPos = (internalPosCount)* groupOffset + groupBins;

      h(i,:) = bar(Y, 'stacked','Facecolor','flat');
      set(h(i,:),'BarWidth',groupOffset);
      set(h(i,:),'XData',groupDrawPos);

     if i==1
      h(1,1).CData=[189/255 183/255 107/255];
      h(1,2).CData=[ 77/255  77/255 1];
      h(1,3).CData=[173/255 1 47/255];
     elseif i==2
          h(2,1).CData=[189/255 183/255 107/255];
          h(2,2).CData=[ 77/255  77/255 1];
     elseif i==3  
          h(3,1).CData=[189/255 183/255 107/255];
          h(3,2).CData=[ 77/255  77/255 1];

          h(3,3).CData=[0 0 0];
     end
  end
  %h=gca;
  %h.FaceColor='flat';
  %h.CData(2,:)=[.5 0 .5];
  hold off;
  set(gca,'XTickMode','manual');
  set(gca,'XTick',1:NumGroupsPerAxis);
  set(gca,'XTickLabelMode','manual');
  set(gca,'XTickLabel',groupLabels);
  end

What I want to do is: in the first three Groups the middle Elements should be black instead of blue (as they represent small/medium/large amounts of electricity generated with hard coal) while in the groups 7 to 9 the middle Elements should be green instead of blue cause they represent electricity from renewable energy.

Cheers, Sonja

[jonas]

Can you provide the data so that we can run the code?

[Sonja Glas]

Hi jonas, Data(:,:,1) has 0.0241 for all values Data(:,:,2) has these values [0.0131 0.0131 0.0131; 0.0183 0.0183 0.0183; 0.0366 0.0366 0.0366; 0.0080 0.0080 0.0080; 0.0113 0.0113 0.0113; 0.0225 0.0225 0.0225; 0.0017 0.0017 0.0017; 0.0024 0.0024 0.0024; 0.0048 0.0048 0.0048] Data(:,:,3) has these values [0.0265 0.1087 0.2034; 0.0265 0.1087 0.2034; 0.0265 0.1087 0.2034; 0.0265 0.1087 0.2034; 0.0265 0.1087 0.2034; 0.0265 0.1087 0.2034; 0.0265 0.1087 0.2034; 0.0265 0.1087 0.2034; 0.0265 0.1087 0.2034] Do you have an idea?

Thanks in advance, Sonja

Answer 1

[jonas]

Note that I have changed the solution completely after the discussion below. You may now change the color of any single bar-segments. This solution is also independent of third-party functions. Load the data (given in the comments) and ticklabels before running.

%Fix data structure

Data=permute(Data,[2 1 3])
Data=reshape(Data,27,3)

%set some input
Groups=9;
Stacks=3;
NumInGroup=length(Data)/Groups;
NumBars=Groups*NumInGroup;

%Plot
figure;hold on;
for i=1:Groups*Stacks;
    h(i,1:3)=bar([Data(i,:);nan(1,3)],'stacked');
end

%Group and set xdata
x1=1:Groups
x0=x1-0.25;
x2=x1+0.25;
pos=[x0;x1;x2];
xpos=pos(:)';
for i=1:Stacks*Groups
    set(h(i,:),'xdata',xpos(i))
end
set(h,'barwidth',0.25)

%Set labels
set(gca,'xtick',[1:Groups],...
    'xticklabels',groupLabels)

%Set colors
set(h(:,1),'facecolor',[.7 .7 .7])
set(h(1:9,2),'facecolor',[0 0 1])
set(h(10:end-9,2),'facecolor',[0.7 0.3 0.4])
set(h(end-8:end,2),'facecolor',[0 1 0])
set(h(1:3:end,3),'facecolor',[1 1 1])
set(h(2:3:end,3),'facecolor',[0.5 0.1 0.9])
set(h(3:3:end,3),'facecolor',[0 0 0])

[Sonja Glas]

Hi Jonas, thanks a lot for your solution! Your answer works somehow, but to me it's mysterious. The numbers for assigning the colors to the bar elements make no sense at all to me. I am trying now to assign the [.5 .5 .5] color only to the "hc s/m/l" middle elements (first three groups)

set(h(2:3:8),'facecolor',neongreen)
    set(h(10:2:2),'facecolor',grey)
%set(h(11:3:end),'facecolor',grey)

set(h(10:2:2),'facecolor',grey) e.g. doesn't do anything. I tried to make sense out of the line with the neongreen changing numbers gradually...?????

If you could just tell which line makes the middle elements of groups one to three grey instead of groups 1 to 6. I'd appreciate! Then I only have to add some boxes to the legend and it'll be perfect! Thanks in advance. Sonja

[Sonja Glas]

Is it that I have to plot in three parts?

plotBarStackGroups(Data(1:end-6,:,:),{'',''});
plotBarStackGroups(Data(end-5:end-3,:,:),{'',''});
plotBarStackGroups(Data(end-2:end,:,:),groupLabels);

??? but then also the xdata stuff changes and I don't know how... Matlab should come around with some direct Solutions!

[jonas]

You are correct. If you want the first three groups to have unique properties then you need three plots. Problem is that the handles refer to one bar (or bar segment) of each group, so you cannot single out one group.

I can write something a bit later today if you havnt solved it. On my phone atm.

[jonas]

I've updated the solution. This one is more general and has no third-party functions, so it should be more help for anyone wanting to use it in the future.

[Sonja Glas]

That's cool - I'll try to implement is tomorrow. Thanks a lot!

[jonas]

Happy to help. Let me know if it works!

[Sonja Glas]

It worked perfectly! Great! Thanks again!