hey大家好久不見,本期帶來一篇繪圖復刻,居然已經出到第11篇了,不知道有朝一日有沒有機會破百,本期繪制的是《PNAS》期刊中pnas.2200057120文章的figure03,文章題目為Intranasal delivery of full-length anti-Nogo-A antibody: A potential alternative route for therapeutic antibodies to central nervous system targets可在以下網址下載該文章pdf原文:
- https://www.pnas.org/doi/epdf/10.1073/pnas.2200057120
本期要復刻的繪圖長這樣:
本人復刻的效果長這樣:
雖然還是有些許不同但是大部分元素已經還原出來了,以下為教程部分:
請盡量使用新版本的MATLAB,比較舊版本的MATLAB無XJitter功能,本人目前較忙近期不會有空開發類似功能,因此在我未有空開發相關功能前,請盡量使用比較新的版本!!!!!!!!~
教程部分
0 數據準備
這里隨便隨機生成了一點數據,文章中的五個顏色都被我提取了出來,我們這次先取用前兩個顏色:
figure('Units','normalized','Position',[.2,.3,.36,.45],'Color','w');rng(24)
% 隨機生成了兩組數據
DataA = rand(7,1)*ones(1,5).*2+rand(7,5)./3.5;
DataB = rand(7,1)*ones(1,5)./2+rand(7,5)./3.5;
meanData = [mean(DataA,2), mean(DataB,2)];% 文章圖片中的顏色數據
CList = [188,188,240; 160,161,166; 237,187,128;177,202,233; 245,185,192]./255;
CList = CList([1,2],:);
% CList = CList([4,5],:);% 橫坐標標簽文本
NameList = {'Cortex';'Hippocampus';'Cerebellum';'Brainstem';'Cervical spinal cord';'Thoracic spinal cord';'Lumbar spinal cord'};
1 基礎繪圖
柱狀圖
hold on
% 繪制柱狀圖 ---------------------------------------------------------------
barHdl = bar(meanData,'EdgeColor','none','FaceAlpha',.5,'BarWidth',.7);
% 修改配色
barHdl(1).FaceColor = CList(1,:);
barHdl(2).FaceColor = CList(2,:);
散點圖
% 繪制散點圖 ---------------------------------------------------------------
XA = barHdl(1).XEndPoints.'*ones(1,size(DataA,2));
scatter(XA(:),DataA(:),55,'filled','o','MarkerEdgeColor','k','LineWidth',.8,...'CData',CList(1,:),'XJitter','rand','XJitterWidth',0.15)
XB = barHdl(2).XEndPoints.'*ones(1,size(DataB,2));
scatter(XB(:),DataB(:),55,'filled','o','MarkerEdgeColor','k','LineWidth',.8,...'CData',CList(2,:),'XJitter','rand','XJitterWidth',0.15)
誤差棒
這里直接用的標準差:
% 繪制誤差棒 ---------------------------------------------------------------
errorbar(barHdl(1).XEndPoints,meanData(:,1),std(DataA,0,2),'vertical',...'LineStyle','none','LineWidth',1,'Color','k')
errorbar(barHdl(2).XEndPoints,meanData(:,2),std(DataB,0,2),'vertical',...'LineStyle','none','LineWidth',1,'Color','k')
3 坐標區域簡單修飾
主要是調整字體,并修改X軸刻度標簽:
% 繪制誤差棒 ---------------------------------------------------------------
errorbar(barHdl(1).XEndPoints,meanData(:,1),std(DataA,0,2),'vertical',...'LineStyle','none','LineWidth',1,'Color','k')
errorbar(barHdl(2).XEndPoints,meanData(:,2),std(DataB,0,2),'vertical',...'LineStyle','none','LineWidth',1,'Color','k')
4 截斷坐標軸
這里用的我之前寫的一個函數,為了怕大家找起來麻煩,放在了文末及文末的gitee倉庫里。
以下是調用文末截斷坐標軸函數實現的:
% 截斷坐標軸
truncAxis('y',[1,1.4]);
5 坐標軸精修
包括隱藏baseline,添加Y軸標簽,修改坐標軸粗細及刻度方向等:
fig = gcf;
ax1 = fig.Children(1);
ax2 = fig.Children(2);% 隱藏基線
ax1.XColor = 'none';
ax1.Children(end).BaseLine.Color = 'none';% 坐標軸修飾
ax1.LineWidth = 1.5;
ax2.LineWidth = 1.5;
ax1.TickDir = 'out';
ax2.TickDir = 'out';% 增添Y軸標簽
ax1.YLabel.String = 'ug of Ab/ g wet weight of tissue';
ax1.YLabel.Position = [0 1.4 -1];
ax1.YLabel.FontSize = 15;
6 添加顯著性標注
就一個柱子一個柱子加吧,改改序號和字符就能得到比較符合想法的效果:
% 隨便加點顯著性標志N = 2; % 第二個柱
S = '***';
% 要修改右側柱請改成barHdl(2)及DataB(N,:)
X = barHdl(1).XEndPoints(N);
Y = max(DataA(N,:))+.1;
errorbar(X,Y,.2,'horizontal','LineStyle','none','LineWidth',1,'Color','k')
text(X,Y,S,'FontSize',15,'FontWeight','bold','FontName','Arial',...'HorizontalAlignment','center','VerticalAlignment','baseline')N = 3; % 第三個柱
S = '****';
% 要修改右側柱請改成barHdl(2)及DataB(N,:)
X = barHdl(1).XEndPoints(N);
Y = max(DataA(N,:))+.1;
errorbar(X,Y,.2,'horizontal','LineStyle','none','LineWidth',1,'Color','k')
text(X,Y,S,'FontSize',15,'FontWeight','bold','FontName','Arial',...'HorizontalAlignment','center','VerticalAlignment','baseline')N = 6; % 第六個柱
S = '****';
% 要修改右側柱請改成barHdl(2)及DataB(N,:)
X = barHdl(1).XEndPoints(N);
Y = max(DataA(N,:))+.1;
errorbar(X,Y,.2,'horizontal','LineStyle','none','LineWidth',1,'Color','k')
text(X,Y,S,'FontSize',15,'FontWeight','bold','FontName','Arial',...'HorizontalAlignment','center','VerticalAlignment','baseline')
更多顏色展示
完整代碼
教程完整代碼
% code for pnas.2200057120fig03
% by slandarer
% bar chart with trunc-axisfigure('Units','normalized','Position',[.2,.3,.36,.45],'Color','w');rng(24)
% 隨機生成了兩組數據
DataA = rand(7,1)*ones(1,5).*2+rand(7,5)./3.5;
DataB = rand(7,1)*ones(1,5)./2+rand(7,5)./3.5;
meanData = [mean(DataA,2), mean(DataB,2)];% 文章圖片中的顏色數據
CList = [188,188,240; 160,161,166; 237,187,128;177,202,233; 245,185,192]./255;
CList = CList([1,2],:);
% CList = CList([4,5],:);% 橫坐標標簽文本
NameList = {'Cortex';'Hippocampus';'Cerebellum';'Brainstem';'Cervical spinal cord';'Thoracic spinal cord';'Lumbar spinal cord'};hold on
% 繪制柱狀圖 ---------------------------------------------------------------
barHdl = bar(meanData,'EdgeColor','none','FaceAlpha',.5,'BarWidth',.7);
% 修改配色
barHdl(1).FaceColor = CList(1,:);
barHdl(2).FaceColor = CList(2,:);
% 繪制散點圖 ---------------------------------------------------------------
XA = barHdl(1).XEndPoints.'*ones(1,size(DataA,2));
scatter(XA(:),DataA(:),55,'filled','o','MarkerEdgeColor','k','LineWidth',.8,...'CData',CList(1,:),'XJitter','rand','XJitterWidth',0.15)
XB = barHdl(2).XEndPoints.'*ones(1,size(DataB,2));
scatter(XB(:),DataB(:),55,'filled','o','MarkerEdgeColor','k','LineWidth',.8,...'CData',CList(2,:),'XJitter','rand','XJitterWidth',0.15)
% 繪制誤差棒 ---------------------------------------------------------------
errorbar(barHdl(1).XEndPoints,meanData(:,1),std(DataA,0,2),'vertical',...'LineStyle','none','LineWidth',1,'Color','k')
errorbar(barHdl(2).XEndPoints,meanData(:,2),std(DataB,0,2),'vertical',...'LineStyle','none','LineWidth',1,'Color','k')% 坐標區域簡單修飾
ax = gca;
ax.XTick = 1:length(barHdl(1).XEndPoints);
ax.XTickLabel = NameList(:);
ax.FontName = 'Arial';
ax.FontWeight = 'bold';
ax.FontSize = 11;
ax.XTickLabelRotation = 35;% 截斷坐標軸
truncAxis('y',[1,1.4]);
fig = gcf;
ax1 = fig.Children(1);
ax2 = fig.Children(2);% 隱藏基線
ax1.XColor = 'none';
ax1.Children(end).BaseLine.Color = 'none';% 坐標軸修飾
ax1.LineWidth = 1.5;
ax2.LineWidth = 1.5;
ax1.TickDir = 'out';
ax2.TickDir = 'out';% 增添Y軸標簽
ax1.YLabel.String = 'ug of Ab/ g wet weight of tissue';
ax1.YLabel.Position = [0 1.4 -1];
ax1.YLabel.FontSize = 15;% -------------------------------------------------------------------------
% 隨便加點顯著性標志N = 2; % 第二個柱
S = '***';
% 要修改右側柱請改成barHdl(2)及DataB(N,:)
X = barHdl(1).XEndPoints(N);
Y = max(DataA(N,:))+.1;
errorbar(X,Y,.2,'horizontal','LineStyle','none','LineWidth',1,'Color','k')
text(X,Y,S,'FontSize',15,'FontWeight','bold','FontName','Arial',...'HorizontalAlignment','center','VerticalAlignment','baseline')N = 3; % 第三個柱
S = '****';
% 要修改右側柱請改成barHdl(2)及DataB(N,:)
X = barHdl(1).XEndPoints(N);
Y = max(DataA(N,:))+.1;
errorbar(X,Y,.2,'horizontal','LineStyle','none','LineWidth',1,'Color','k')
text(X,Y,S,'FontSize',15,'FontWeight','bold','FontName','Arial',...'HorizontalAlignment','center','VerticalAlignment','baseline')N = 6; % 第六個柱
S = '****';
% 要修改右側柱請改成barHdl(2)及DataB(N,:)
X = barHdl(1).XEndPoints(N);
Y = max(DataA(N,:))+.1;
errorbar(X,Y,.2,'horizontal','LineStyle','none','LineWidth',1,'Color','k')
text(X,Y,S,'FontSize',15,'FontWeight','bold','FontName','Arial',...'HorizontalAlignment','center','VerticalAlignment','baseline')
截斷坐標軸工具函數完整代碼
function truncAxis(varargin)
% @author : slandarer
% gzh : slandarer隨筆% 獲取參數
if isa(varargin{1},'matlab.graphics.axis.Axes')ax=varargin{1};varargin(1)=[];
elseax=gca;
end
hold(ax,'on');
% box(ax,'off')
ax.XAxisLocation='bottom';
ax.YAxisLocation='left';
axisPos=ax.Position;
axisXLim=ax.XLim;
axisYLim=ax.YLim;
axisXScale=diff(axisXLim);
axisYScale=diff(axisYLim);
truncRatio=1/20;
Xtrunc=[];Ytrunc=[];
for i=1:length(varargin)-1switch truecase strcmpi('X',varargin{i}),Xtrunc=varargin{i+1};case strcmpi('Y',varargin{i}),Ytrunc=varargin{i+1};end
end
switch truecase isempty(Xtrunc)% 復制坐標區域ax2=copyAxes(ax);% 修改軸基礎屬性ax2.XTickLabels=[];ax2.XColor='none';% 修改坐標區域范圍ax.YLim=[axisYLim(1),Ytrunc(1)];ax2.YLim=[Ytrunc(2),axisYLim(2)];% 坐標區域重定位ax.Position(4)=axisPos(4)*(1-truncRatio)/(axisYScale-diff(Ytrunc))*(Ytrunc(1)-axisYLim(1));ax2.Position(2)=axisPos(2)+ax.Position(4)+axisPos(4)*truncRatio;ax2.Position(4)=axisPos(4)*(1-truncRatio)/(axisYScale-diff(Ytrunc))*(axisYLim(2)-Ytrunc(2));% 鏈接軸范圍變動linkaxes([ax,ax2],'x')% 添加線和標識符if strcmp(ax.Box,'on')ax.Box='off';ax2.Box='off';annotation('line',[1,1].*(ax.Position(1)+ax.Position(3)),[ax.Position(2),ax.Position(2)+ax.Position(4)],'LineStyle','-','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[1,1].*(ax.Position(1)+ax.Position(3)),[ax2.Position(2),ax2.Position(2)+ax2.Position(4)],'LineStyle','-','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[ax.Position(1),ax.Position(1)+ax.Position(3)],[1,1].*(ax2.Position(2)+ax2.Position(4)),'LineStyle','-','LineWidth',ax.LineWidth,'Color',ax.XColor);elseannotation('line',[ax.Position(1),ax.Position(1)+ax.Position(3)],[1,1].*(ax.Position(2)+ax.Position(4)),'LineStyle',':','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[ax.Position(1),ax.Position(1)+ax.Position(3)],[1,1].*(ax2.Position(2)),'LineStyle',':','LineWidth',ax.LineWidth,'Color',ax.XColor);endcreateSlash([ax.Position(1)-.2,ax.Position(2)+ax.Position(4)-.2,.4,.4])createSlash([ax.Position(1)-.2,ax2.Position(2)-.2,.4,.4])createSlash([ax.Position(1)+ax.Position(3)-.2,ax.Position(2)+ax.Position(4)-.2,.4,.4])createSlash([ax.Position(1)+ax.Position(3)-.2,ax2.Position(2)-.2,.4,.4])case isempty(Ytrunc) % 復制坐標區域ax2=copyAxes(ax);% 修改軸基礎屬性ax2.YTickLabels=[];ax2.YColor='none';% 修改坐標區域范圍ax.XLim=[axisXLim(1),Xtrunc(1)];ax2.XLim=[Xtrunc(2),axisXLim(2)];% 坐標區域重定位ax.Position(3)=axisPos(3)*(1-truncRatio)/(axisXScale-diff(Xtrunc))*(Xtrunc(1)-axisXLim(1));ax2.Position(1)=axisPos(1)+ax.Position(3)+axisPos(3)*truncRatio;ax2.Position(3)=axisPos(3)*(1-truncRatio)/(axisXScale-diff(Xtrunc))*(axisXLim(2)-Xtrunc(2));% 鏈接軸范圍變動linkaxes([ax,ax2],'y')% 添加線和標識符if strcmp(ax.Box,'on')ax.Box='off';ax2.Box='off';annotation('line',[ax.Position(1),ax.Position(1)+ax.Position(3)],[1,1].*(ax.Position(2)+ax.Position(4)),'LineStyle','-','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[ax2.Position(1),ax2.Position(1)+ax2.Position(3)],[1,1].*(ax.Position(2)+ax.Position(4)),'LineStyle','-','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[1,1].*(ax2.Position(1)+ax2.Position(3)),[ax2.Position(2),ax2.Position(2)+ax2.Position(4)],'LineStyle','-','LineWidth',ax.LineWidth,'Color',ax.XColor);elseannotation('line',[1,1].*(ax.Position(1)+ax.Position(3)),[ax2.Position(2),ax2.Position(2)+ax2.Position(4)],'LineStyle',':','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[1,1].*(ax2.Position(1)),[ax2.Position(2),ax2.Position(2)+ax2.Position(4)],'LineStyle',':','LineWidth',ax.LineWidth,'Color',ax.XColor);endcreateSlash([ax.Position(1)+ax.Position(3)-.2,ax.Position(2)-.2,.4,.4])createSlash([ax2.Position(1)-.2,ax.Position(2)-.2,.4,.4])createSlash([ax.Position(1)+ax.Position(3)-.2,ax.Position(2)+ax.Position(4)-.2,.4,.4])createSlash([ax2.Position(1)-.2,ax.Position(2)+ax.Position(4)-.2,.4,.4])case (~isempty(Ytrunc))&(~isempty(Ytrunc))% 復制坐標區域ax2=copyAxes(ax);ax3=copyAxes(ax);ax4=copyAxes(ax);% 修改軸基礎屬性ax2.XTickLabels=[];ax2.XColor='none';ax3.XTickLabels=[];ax3.XColor='none';ax3.YTickLabels=[];ax3.YColor='none';ax4.YTickLabels=[];ax4.YColor='none';% 修改坐標區域范圍ax.YLim=[axisYLim(1),Ytrunc(1)];ax.XLim=[axisXLim(1),Xtrunc(1)];ax2.XLim=[axisXLim(1),Xtrunc(1)];ax2.YLim=[Ytrunc(2),axisYLim(2)];ax3.XLim=[Xtrunc(2),axisXLim(2)];ax3.YLim=[Ytrunc(2),axisYLim(2)];ax4.XLim=[Xtrunc(2),axisXLim(2)];ax4.YLim=[axisYLim(1),Ytrunc(1)];% 坐標區域重定位ax.Position(3)=axisPos(3)*(1-truncRatio)/(axisXScale-diff(Xtrunc))*(Xtrunc(1)-axisXLim(1));ax.Position(4)=axisPos(4)*(1-truncRatio)/(axisYScale-diff(Ytrunc))*(Ytrunc(1)-axisYLim(1));ax2.Position(2)=axisPos(2)+ax.Position(4)+axisPos(4)*truncRatio;ax2.Position(3)=axisPos(3)*(1-truncRatio)/(axisXScale-diff(Xtrunc))*(Xtrunc(1)-axisXLim(1));ax2.Position(4)=axisPos(4)*(1-truncRatio)/(axisYScale-diff(Ytrunc))*(axisYLim(2)-Ytrunc(2));ax3.Position(1)=axisPos(1)+ax.Position(3)+axisPos(3)*truncRatio;ax3.Position(2)=axisPos(2)+ax.Position(4)+axisPos(4)*truncRatio;ax3.Position(3)=axisPos(3)*(1-truncRatio)/(axisXScale-diff(Xtrunc))*(axisXLim(2)-Xtrunc(2));ax3.Position(4)=axisPos(4)*(1-truncRatio)/(axisYScale-diff(Ytrunc))*(axisYLim(2)-Ytrunc(2));ax4.Position(1)=axisPos(1)+ax.Position(3)+axisPos(3)*truncRatio;ax4.Position(3)=axisPos(3)*(1-truncRatio)/(axisXScale-diff(Xtrunc))*(axisXLim(2)-Xtrunc(2));ax4.Position(4)=axisPos(4)*(1-truncRatio)/(axisYScale-diff(Ytrunc))*(Ytrunc(1)-axisYLim(1));% 鏈接軸范圍變動linkaxes([ax3,ax2],'y')linkaxes([ax4,ax3],'x')linkaxes([ax,ax2],'x')linkaxes([ax,ax4],'y')% 添加線和標識符if strcmp(ax.Box,'on')ax.Box='off';ax2.Box='off';ax3.Box='off';ax4.Box='off';annotation('line',[ax.Position(1),ax.Position(1)+ax.Position(3)],[1,1].*(ax2.Position(2)+ax2.Position(4)),'LineStyle','-','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[ax3.Position(1),ax3.Position(1)+ax3.Position(3)],[1,1].*(ax2.Position(2)+ax2.Position(4)),'LineStyle','-','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[1,1].*(ax4.Position(1)+ax4.Position(3)),[ax3.Position(2),ax3.Position(2)+ax3.Position(4)],'LineStyle','-','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[1,1].*(ax4.Position(1)+ax4.Position(3)),[ax4.Position(2),ax4.Position(2)+ax4.Position(4)],'LineStyle','-','LineWidth',ax.LineWidth,'Color',ax.XColor);elseannotation('line',[1,1].*(ax.Position(1)+ax.Position(3)),[ax2.Position(2),ax2.Position(2)+ax2.Position(4)],'LineStyle',':','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[1,1].*(ax3.Position(1)),[ax2.Position(2),ax2.Position(2)+ax2.Position(4)],'LineStyle',':','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[1,1].*(ax.Position(1)+ax.Position(3)),[ax.Position(2),ax.Position(2)+ax.Position(4)],'LineStyle',':','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[1,1].*(ax3.Position(1)),[ax.Position(2),ax.Position(2)+ax.Position(4)],'LineStyle',':','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[ax.Position(1),ax.Position(1)+ax.Position(3)],[1,1].*(ax.Position(2)+ax.Position(4)),'LineStyle',':','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[ax.Position(1),ax.Position(1)+ax.Position(3)],[1,1].*(ax2.Position(2)),'LineStyle',':','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[ax4.Position(1),ax4.Position(1)+ax4.Position(3)],[1,1].*(ax.Position(2)+ax.Position(4)),'LineStyle',':','LineWidth',ax.LineWidth,'Color',ax.XColor);annotation('line',[ax4.Position(1),ax4.Position(1)+ax4.Position(3)],[1,1].*(ax2.Position(2)),'LineStyle',':','LineWidth',ax.LineWidth,'Color',ax.XColor);endcreateSlash([ax.Position(1)-.2,ax.Position(2)+ax.Position(4)-.2,.4,.4])createSlash([ax.Position(1)-.2,ax2.Position(2)-.2,.4,.4])createSlash([ax4.Position(1)+ax4.Position(3)-.2,ax.Position(2)+ax.Position(4)-.2,.4,.4])createSlash([ax4.Position(1)+ax4.Position(3)-.2,ax2.Position(2)-.2,.4,.4])createSlash([ax.Position(1)+ax.Position(3)-.2,ax.Position(2)-.2,.4,.4])createSlash([ax.Position(1)+ax.Position(3)-.2,ax2.Position(2)+ax2.Position(4)-.2,.4,.4])createSlash([ax4.Position(1)-.2,ax.Position(2)-.2,.4,.4])createSlash([ax4.Position(1)-.2,ax2.Position(2)+ax2.Position(4)-.2,.4,.4])% 修改當前坐標區域,方便legend添加set(gcf,'currentAxes',ax3)
end
% 復制原坐標區域全部可復制屬性function newAX=copyAxes(ax)axStruct=get(ax);fNames=fieldnames(axStruct);newAX=axes('Parent',ax.Parent);coeList={'CurrentPoint','XAxis','YAxis','ZAxis','BeingDeleted',...'TightInset','NextSeriesIndex','Children','Type','Legend'};for n=1:length(coeList)coePos=strcmp(fNames,coeList{n});fNames(coePos)=[];endfor n=1:length(fNames)newAX.(fNames{n})=ax.(fNames{n});endcopyobj(ax.Children,newAX)end
% 添加截斷標識符函數function createSlash(pos)anno=annotation('textbox');anno.String='/';anno.LineStyle='none';anno.FontSize=15;anno.Position=pos;anno.FitBoxToText='on';anno.VerticalAlignment='middle';anno.HorizontalAlignment='center';end
end
完
以上已經是本文全部內容,若懶得復制代碼,可以去以下gitee倉庫獲取全部代碼:
https://gitee.com/slandarer/PLTreprint/