How to draw the below type of graph in MATLAB?

% import bar plot data from spreadsheet
opts = detectImportOptions('data.xlsx','Range',[1 1 18 18]);
data = readtable('data.xlsx',opts);
% import line plot data
ave = readmatrix('data.xlsx',"Range","B20");
% Create x vector of 13 months to give room for Dec bar plot
d = datetime(2019,1,1,'Format',"MMM") + calmonths(0:12);
% Plot ave data
plot(d(1:end-1),ave,'k--','LineWidth',1.5)
ylim([0 16])
% Plot 13th month separately to give room for Dec bar plot
hold on
plot(d(end),0)
hold off
% Turn off x-label for Jan 2020
ax = gca;
ax.XTickLabel{13} = "";
% Use position of current plot to determine placement of bar plots
p = get(gca,'position'); % [left bottom width height]
dx = [0 p(3).*days(diff(d))/365]; % width of each month
% Create bin data for bar plots
edges = 1:16;
% plot bar plots by positioning a new axes on top of the current axis
for i = 1:length(d)-1
    a = axes('Position',[p(1)+sum(dx(1:i)) p(2) .75*dx(i+1) p(4)]);
    b = barh(a,edges,data{:,i+1});
    axis off
end

Not perfect, but it should be able to get your started.

8 个评论
显示 6更早的评论隐藏 6更早的评论

Cris LaPierre 2019-9-11

编辑：Cris LaPierre 2019-9-11

在 MATLAB Online 中打开

The y axis on the bar plots still starts at 0 even though the edges start at 1, but it's a valid critique. Best to ensure X and Y lims are consistant. That can be accomplished by adding an explicit xlim and ylim when creating the bar plots.

I hardcoded the limits using the values contained in the data and looking at the original image. This was just a quick attempt. They can be replaced with code by someone with more interest in creating a dynamic solution.

Also noticed that barh centers the bar around the x value while the original places the bar between the edges. That can be fixed by adding an offset of 0.5 when plotting. I've left the axis on here to show that x is the same in all bar plots, and y aligns with y axis of line plot for all.

That doesn't look particularly good, so here it is w/o the axes.

I noticed the original also had the line in front of the bar plots. It's possible to do that as well (use uistack).

Here's the updated code:

% import bar plot data from spreadsheet
opts = detectImportOptions('ChanghoLEE_data.xlsx','Range',[1 1 18 18]);
data = readtable('ChanghoLEE_data.xlsx',opts);
% import line plot data
ave = readmatrix('ChanghoLEE_data.xlsx',"Range","B20");
% Create x vector of 13 months to give room for Dec bar plot
d = datetime(2019,1,1,'Format',"MMM") + calmonths(0:12);
% Plot ave data
plot(d(1:end-1),ave,'k--','LineWidth',1.5)
ylim([0 16])
% Plot 13th month separately to give room for Dec bar plot
hold on
plot(d(end),0)
hold off
% Turn off x-label for Jan 2020
ax = gca;
ax.XTickLabel{13} = "";
% Use position of current plot to determine placement of bar plots
p = get(gca,'position'); % [left bottom width height]
dx = [0 p(3).*days(diff(d))/365]; % width of each month
% Create bin data for bar plots
edges = 1:16;
% plot bar plots by positioning a new axes on top of the current axis
for i = 1:length(d)-1
    a = axes('Position',[p(1)+sum(dx(1:i)) p(2) .75*dx(i+1) p(4)]);
    barh(a,edges+.5,data{:,i+1});
    ylim([0 16])
    xlim([0 14])
    axis off
end
% Bring line in front of bar plots
uistack(ax,'top')
set(ax,'color','none')

Cris LaPierre 2019-9-11

Good point. I thought about doing it that way, but went with the summary axis first so that it could set where I needed to place the histograms. Had to hunt for the uistack fxn because of that. I'd never used it before.

Changho LEE 2019-9-11

This graph is part of my teacher's research in 1975.

His research was so beautiful that I wanted to make it MATLAB and use it for other research.

You made it so perfect and beautiful.

It looks really elegant.

This is a wonderful thing I never thought of.

Thank you very much for creating such a wonderful work of art.

Thank you very much.

请先登录，再进行评论。

Answer 2

Adam Danz 2019-9-10

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/479767-how-to-draw-the-below-type-of-graph-in-matlab#answer_391291

编辑：Adam Danz 2019-9-11

在 MATLAB Online 中打开

"Feb" in your excel file has a $ character which needs fixed. If you have many excel sheets, you can fix it from within the code but I assumed it's fixed prior to running this.

The first section loads the data and prepares it for histogram().

The 2nd section defines the plot layout. You can change the margin size etc.

The 3rd section does the plotting and saves each axis handle.

The 4th section cleans up the axes and does all of the cosmetics.

The 5th section overlays an additonal axis at the same scale as the underlying subplots (as dpb suggested) and plots the mean curve.

See comments for details.

data = readtable('data.xlsx'); 
edgesStr = data{:,1};  %save edges to an array
edges = str2double(cellfun(@(x)regexp(x,'^\d+','match'),edgesStr(1:end-2))); %edges, numeric
data(:,1) = [];     %remove edges from data table
data = fillmissing(data,'constant', 0); %replace NaNs with 0
% Determine subplot dimensions (all in normalized units)
nSubplots = size(data,2);                           %number of subplots needed
LRmargins = [.1,.08];                              %left and right figure margins
TBmargins = [.1, .1];                               %top and bottom figure margins
subplotWidth = (1-sum(LRmargins))/nSubplots;        %width of each subplot
subplotHeight = 1-sum(TBmargins);                   %height of each subplot
subplotPos = linspace(LRmargins(1),1-LRmargins(2)-subplotWidth, nSubplots); % lateral position of each subplot
% Loop through each subplot, create subplot and plot the data
figure();
binEdges = [edges(:).',edges(end)+1]; 
sph = gobjects(1,nSubplots); %store subplot handles
for i = 1:nSubplots
    sph(i) = axes('Units', 'Normalize', 'Position', [subplotPos(i),TBmargins(1),subplotWidth,subplotHeight]); 
    histogram(sph(i),'BinEdges', binEdges, 'BinCounts', data{1:17,i},'Orientation','Horizontal'); 
end
% Clean up
maxBarHeight = max(data{1:17,:},[],'all'); %max bar height
% remove the x axis ticks for all subplots
set(sph,'XTick',[])
% Set the y tick for the left most subplot as the bin lables
set(sph,'YTick', binEdges(1:end-1))
set(sph(2:end), 'YTickLabel',[])
% Make sure all axes have the same ylim
set(sph,'YLim',binEdges([1,end-1]))
% Make sure all axes have same scale
set(sph,'XLim',[0,maxBarHeight])
% Set the month labels along the horizontal, bottom axis
arrayfun(@(x,y)set(x.XLabel,'String',y),sph,data.Properties.VariableNames);
% Label bins
ylabel(sph(1),'Days')
% overlay summary axis
sph2 = axes('Position',[subplotPos(1),TBmargins(1),1-sum(LRmargins),subplotHeight]); 
% plot mean curve
plot(data{end,:}, 'r-s','Linewidth',1.5)
% Scale the axis to match the range of subplots
xlim([.5,12.5])
ylim(ylim(sph(1)))
% Turn off axis 
axis(sph2, 'off')