MATLAB | 全网最全边际图绘制模板(直方图、小提琴图、箱线图、雨云图、散点图... ...)

如标题所言，这应该是全网最全的边际图绘制模板，中心图有8种格式，边际图有11种格式，共计88种组合，另外模板中给了8款配色，我愿称其为888组合，只需要更换一下数据就能绘制出各种类型的边际图:

甚至可以随意拖动改变X,Y轴范围：

---

教程部分

以下代码展示顺序仅为讲解方便，实际顺序请见完整代码部分。

0 数据准备

就每组数据都是两列的数值矩阵，之后放在元胞数组中即可，如下例：

% 构造三个符合高斯分布的点集
PntSet1=mvnrnd([2 3],[1 0;0 2],300);
PntSet2=mvnrnd([6 7],[1 0;0 2],300);
PntSet3=mvnrnd([14 9],[1 0;0 1],300);
% 将数据放进元胞数组
PntSet={PntSet1,PntSet2,PntSet3};

1 坐标区域绘制

就简单的构造了三个坐标区域，X,Y轴标签啥的都可自行修改：

% figure图床及主要区域axes坐标区域创建及基础修饰
fig=figure('Units','normalized','Position',[.3,.2,.5,.63]);
axM=axes('Parent',fig);hold on;
set(axM,'Position',[.08,.08,.65,.65],'LineWidth',1.1,'Box','on','TickDir','in',...
    'XMinorTick','on','YMinorTick','on','XGrid','on','YGrid','on','GridLineStyle','--',...
    'FontName','Times New Roman','FontSize',12,'GridAlpha',.09)
axM.XLabel.String='Main XXXXX';
axM.YLabel.String='Main YYYYY';
axM.XLabel.FontSize=14;
axM.YLabel.FontSize=14;
% -------------------------------------------------------------------------
% 右侧axes坐标区域创建及基础修饰
axR=axes('Parent',fig);hold on;
set(axR,'Position',[.75,.08,.23,.65],'LineWidth',1.1,'Box','off','TickDir','out',...
    'XMinorTick','on','YMinorTick','on','XGrid','on','YGrid','on','GridLineStyle','--',...
    'FontName','Times New Roman','FontSize',12,'GridAlpha',.09,'TickLength',[.006 .015],'YTickLabel','')
axR.XLabel.String='Y-axis data statistics';
axR.XLabel.FontSize=14;
% -------------------------------------------------------------------------
% 上方axes坐标区域创建及基础修饰
axU=axes('Parent',fig);hold on;
set(axU,'Position',[.08,.76,.65,.2],'LineWidth',1.1,'Box','off','TickDir','out',...
    'XMinorTick','on','YMinorTick','on','XGrid','on','YGrid','on','GridLineStyle','--',...
    'FontName','Times New Roman','FontSize',12,'GridAlpha',.09,'TickLength',[.006 .015],'XTickLabel','')
axU.YLabel.String='X-axis data statistics';
axU.YLabel.FontSize=14;

2 配色

可以自行替换RGB三元组，同时本模板也贴心的提供了八款配色，只需修改:

• colorList=Cn
  中n的数值即可：

% 配色预设(here to change the default color list)
C1=[211 43 43;61 96 137;249 206 61;76 103 86;80 80 80]./255;
C2=[102,173,194;36,59,66;232,69,69;194,148,102;54,43,33]./255;
C3=[223,122,94;60 64 91;130 178 154;244,241,222;240 201 134]./255;
C4=[122,117,119;255,163,25;135,146,73;126,15,4;30,93,134]./255;
C5=[198,199,201;38,74,96;209,80,51;241,174,44;12,13,15]./255;
C6=[235,75,55;77,186,216;58,84,141;2,162,136;245,155,122]./255;
C7=[23,23,23;121,17,36;31,80,91;44,9,75;61,36,42]./255;
C8=[47,62,66;203,129,70;0 64 115;152,58,58;20 72 83]./255;
colorList=C2;

以下是配色展示：

C1

C2

C3

C4

C5

C6

C7

C8

3 绘图格式修改

中心图及边缘图格式可以通过修改：

• mainType
• marginType

俩参数进行修改，mainType为中心图类型，可选1-8，marginType为边际图类型，可选1-11,

% 在这改边缘图及中心图种类
% (here to change the type of marginal plot and main plot)
%    |
%    |
%    |
% \  |  /
%  \ | /
%   \|/
%     
mainType=5;
marginType=7;
%  
%   /|\
%  / | \
% /  |  \
%    |
%    |
%    |

4 中心图

中心图绘制代码：

% 主要区域散点图绘制
switch mainType
    case 1 % 散点图
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            scatter(axM,tPntSet(:,1),tPntSet(:,2),70,'filled','CData',colorList(i,:),'MarkerFaceAlpha',.6)
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,legendStr,'FontSize',14,'Box','off','Location','best');
    case 2 % 等高线图
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            gridx1=linspace(min(tPntSet(:,1)),max(tPntSet(:,1)),round(sqrt(size(tPntSet,1))));
            gridx2=linspace(min(tPntSet(:,2)),max(tPntSet(:,2)),round(sqrt(size(tPntSet,1))));
            [x1,x2]=meshgrid(gridx1, gridx2);
            [m,n]=size(x1);
            x1=x1(:);
            x2=x2(:);
            xi=[x1,x2];
            h=ksdensity(tPntSet,xi);
            contour(axM,reshape(x1,[m,n]),reshape(x2,[m,n]),reshape(h,[m,n]),...
                'LineWidth',1.2,'EdgeColor',colorList(i,:))
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,legendStr,'FontSize',14,'Box','off','Location','best');
    case 3 % 回归曲线
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            [p,S]=polyfit(tPntSet(:,1),tPntSet(:,2),1); 
            x1=linspace(min(tPntSet(:,1)),max(tPntSet(:,1)),100);
            [y_fit,delta]=polyval(p,x1,S);
            % 绘制原始数据、线性拟合和 95% 预测区间 y±2Δ。
            uy=y_fit+2.*delta;
            dy=y_fit-2.*delta;
            % 绘制原始数据
            sHdl(i)=scatter(axM,tPntSet(:,1),tPntSet(:,2),30,'filled','CData',colorList(i,:),'MarkerFaceAlpha',.9);
            % 绘制拟合曲线
            plot(axM,x1,y_fit,'Color',colorList(i,:),'LineWidth',2.5)
            % 绘制置信区间
            fill(axM,[x1,x1(end:-1:1)],[uy,dy(end:-1:1)],colorList(i,:),'EdgeColor','none','FaceAlpha',.15)
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,sHdl,legendStr,'FontSize',14,'Box','off','Location','best');
    case 4 % 散点+凸包
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            [k,av]=convhull(tPntSet);
            sHdl(i)=scatter(axM,tPntSet(:,1),tPntSet(:,2),30,'filled','CData',colorList(i,:),'MarkerFaceAlpha',.9);
            fill(axM,tPntSet(k,1),tPntSet(k,2),colorList(i,:),'FaceAlpha',0.3,'EdgeColor',colorList(i,:),'LineWidth',2)
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,sHdl,legendStr,'FontSize',14,'Box','off','Location','best');
    case 5 % 散点+不规则图形
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            sHdl(i)=scatter(axM,tPntSet(:,1),tPntSet(:,2),30,'filled','CData',colorList(i,:),'MarkerFaceAlpha',.9);
            R=max(max(tPntSet(:,1))-min(tPntSet(:,1)),max(tPntSet(:,2))-min(tPntSet(:,2)))./round(sqrt(size(tPntSet,1)));
            t=linspace(0,2*pi,100);
            [k,~]=convhull(tPntSet);
            kPntSet=tPntSet(k,:);
            kX=repmat(kPntSet(:,1),[1,100])+repmat(cos(t).*R,[length(k),1]);
            kY=repmat(kPntSet(:,2),[1,100])+repmat(sin(t).*R,[length(k),1]);
            kX=kX(:);kY=kY(:);
            [k,av]=convhull([kX,kY]);
            fill(axM,kX(k),kY(k),colorList(i,:),'FaceAlpha',0.3,'EdgeColor',colorList(i,:),'LineWidth',2)
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,sHdl,legendStr,'FontSize',14,'Box','off','Location','best');
    case 6 % 椭圆形
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            sHdl(i)=scatter(axM,tPntSet(:,1),tPntSet(:,2),30,'filled','CData',colorList(i,:),'MarkerFaceAlpha',.9);
            Mu=mean(tPntSet);
            Y=tPntSet-repmat(Mu,size(tPntSet,1),1);
            Sigma=cov(Y);% Sigma=(Y'*(ones(size(pntSet,1),1).*Y))./(size(pntSet,1)-1);
            [X,Y]=getEllipse(Mu,Sigma,9.21,100);
            fill(axM,X,Y,colorList(i,:),'FaceAlpha',0.3,'EdgeColor',colorList(i,:),'LineWidth',2);
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,sHdl,legendStr,'FontSize',14,'Box','off','Location','best');
    case 7 % 簇状
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            sHdl(i)=scatter(axM,tPntSet(:,1),tPntSet(:,2),30,'filled','CData',colorList(i,:),'MarkerFaceAlpha',.9);
            Mu=mean(tPntSet);
            LX=[tPntSet(:,1),Mu(1).*ones(size(tPntSet,1),1),nan(size(tPntSet,1),1)]';
            LY=[tPntSet(:,2),Mu(2).*ones(size(tPntSet,1),1),nan(size(tPntSet,1),1)]';
            plot(axM,LX(:),LY(:),'Color',colorList(i,:),'LineWidth',1.2);
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,sHdl,legendStr,'FontSize',14,'Box','off','Location','best');
    case 8 % 误差棒
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            Mu=mean(tPntSet);
            SX=std(tPntSet(:,1));
            SY=std(tPntSet(:,2));
            errorbar(axM,Mu(1),Mu(2),SY,SY,SX,SX,'o','LineWidth',1.5,'Color',colorList(i,:),'MarkerSize',10);
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,legendStr,'FontSize',14,'Box','off','Location','best');
end

中心图格式展示：

1 散点图

2 等高线图

3 回归曲线

由于数据不是很线形，因此置信区间比较大：

4 多边形边界

5 不规则图形边界

6 椭圆置信区间

7 与质心连线

点少点更好看：

8 误差棒

5 边际图

边际图绘制代码：

% 边际图绘制
axR.YLim=axM.YLim;
axU.XLim=axM.XLim;
for i=1:length(PntSet)
tPntSet=PntSet{i};
switch marginType
    case 1 % 柱状图
        histogram(axR,tPntSet(:,2),'FaceColor',colorList(i,:),'Orientation','horizontal',...
            'EdgeColor',[1 1 1]*.1,'FaceAlpha',0.6,'LineWidth',.8)
        histogram(axU,tPntSet(:,1),'FaceColor',colorList(i,:),...
            'EdgeColor',[1 1 1]*.1,'FaceAlpha',0.6,'LineWidth',.8)
    case 2 % 核密度面
        [f,yi]=ksdensity(tPntSet(:,2));
        fill(axR,[f,0],[yi,yi(1)],colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor','none')
        [f,xi]=ksdensity(tPntSet(:,1));
        fill(axU,[xi,xi(1)],[f,0],colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor','none')
    case 3 % 核密度线
        [f,yi]=ksdensity(tPntSet(:,2));
        plot(axR,f,yi,'Color',colorList(i,:),'LineWidth',2)
        [f,xi]=ksdensity(tPntSet(:,1));
        plot(axU,xi,f,'Color',colorList(i,:),'LineWidth',2)
    case 4 % 核密度面+线
        [f,yi]=ksdensity(tPntSet(:,2));
        fill(axR,[f,0],[yi,yi(1)],colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor','none');
        plot(axR,f,yi,'Color',colorList(i,:),'LineWidth',1.2)
        [f,xi]=ksdensity(tPntSet(:,1));
        fill(axU,[xi,xi(1)],[f,0],colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor','none');
        plot(axU,xi,f,'Color',colorList(i,:),'LineWidth',1.2)
    case 5 % 直方图+密度线
        HR=histogram(axR,tPntSet(:,2),'FaceColor',colorList(i,:),'Orientation','horizontal',...
            'EdgeColor',[1 1 1]*.1,'FaceAlpha',0.6,'LineWidth',.8);
        HU=histogram(axU,tPntSet(:,1),'FaceColor',colorList(i,:),...
            'EdgeColor',[1 1 1]*.1,'FaceAlpha',0.6,'LineWidth',.8);
        [f,yi]=ksdensity(tPntSet(:,2));
        plot(axR,f.*size(tPntSet,1).*HR.BinWidth,yi,'Color',colorList(i,:),'LineWidth',2)
        [f,xi]=ksdensity(tPntSet(:,1));
        plot(axU,xi,f.*size(tPntSet,1).*HU.BinWidth,'Color',colorList(i,:),'LineWidth',2)
    case 6 % 箱线图
        tjDataY=tPntSet(:,2);fullDataY=tjDataY;
        outliBool=isoutlier(tjDataY,'quartiles');
        outli=tjDataY(outliBool);
        scatter(axR,repmat(i,[sum(outliBool),1]),outli,45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.8)
        tjDataY(outliBool)=[];
        qt25=quantile(fullDataY,0.25);
        qt75=quantile(fullDataY,0.75);
        med=median(fullDataY);
        plot(axR,[i,i],[max(tjDataY),qt75],'LineWidth',1.2,'Color',[0,0,0])
        plot(axR,[i,i],[min(tjDataY),qt25],'LineWidth',1.2,'Color',[0,0,0])
        plot(axR,[-.15,.15]+i,[max(tjDataY),max(tjDataY)],'LineWidth',1.2,'Color',[0,0,0])
        plot(axR,[-.15,.15]+i,[min(tjDataY),min(tjDataY)],'LineWidth',1.2,'Color',[0,0,0])
        plot(axR,[-.25,.25]+i,[med,med],'LineWidth',1.2,'Color',colorList(i,:))
        fill(axR,i+.25.*1.*[-1 1 1 -1],[qt25,qt25,qt75,qt75],colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        axR.XLim=[0,length(PntSet)+1];
        % -----------------------------------------------------------------
        tjDataX=tPntSet(:,1);fullDataX=tjDataX;
        outliBool=isoutlier(tjDataX,'quartiles');
        outli=tjDataX(outliBool);
        scatter(axU,outli,repmat(i,[sum(outliBool),1]),45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.8)
        tjDataX(outliBool)=[];
        qt25=quantile(fullDataX,0.25);
        qt75=quantile(fullDataX,0.75);
        med=median(fullDataX);
        plot(axU,[max(tjDataX),qt75],[i,i],'LineWidth',1.2,'Color',[0,0,0])
        plot(axU,[min(tjDataX),qt25],[i,i],'LineWidth',1.2,'Color',[0,0,0])
        plot(axU,[max(tjDataX),max(tjDataX)],[-.15,.15]+i,'LineWidth',1.2,'Color',[0,0,0])
        plot(axU,[min(tjDataX),min(tjDataX)],[-.15,.15]+i,'LineWidth',1.2,'Color',[0,0,0])
        plot(axU,[med,med],[-.25,.25]+i,'LineWidth',1.2,'Color',colorList(i,:))
        fill(axU,[qt25,qt25,qt75,qt75],i+.25.*1.*[-1 1 1 -1],colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        axU.YLim=[0,length(PntSet)+1];
    case 7 % 小提琴图
        tjDataY=tPntSet(:,2);fullDataY=tjDataY;
        [f,yi]=ksdensity(tjDataY);
        f(yi>max(tjDataY))=[];yi(yi>max(tjDataY))=[];
        f(yi<min(tjDataY))=[];yi(yi<min(tjDataY))=[];
        fill(axR,[f,-f(end:-1:1)].*1.2+i,[yi,yi(end:-1:1)],colorList(i,:),'FaceAlpha',0.3,...
            'EdgeColor',colorList(i,:),'LineWidth',1.2);
        outliBool=isoutlier(tjDataY,'quartiles');
        outli=tjDataY(outliBool);
        scatter(axR,repmat(i,[sum(outliBool),1]),outli,45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.8)
        tjDataY(outliBool)=[];
        qt25=quantile(fullDataY,0.25);
        qt75=quantile(fullDataY,0.75);
        med=median(fullDataY);
        plot(axR,[i,i],[max(tjDataY),qt75],'LineWidth',1.2,'Color',[0,0,0])
        plot(axR,[i,i],[min(tjDataY),qt25],'LineWidth',1.2,'Color',[0,0,0])
        fill(axR,i+.2.*1.*[-1 1 1 -1],[qt25,qt25,qt75,qt75],[1,1,1],...
            'FaceAlpha',0.95,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        plot(axR,[-.2,.2]+i,[med,med],'LineWidth',1.2,'Color',colorList(i,:))
        axR.XLim=[0,length(PntSet)+1];
        % -----------------------------------------------------------------
        tjDataX=tPntSet(:,1);fullDataX=tjDataX;
        [f,xi]=ksdensity(tjDataX);
        f(xi>max(tjDataX))=[];xi(xi>max(tjDataX))=[];
        f(xi<min(tjDataX))=[];xi(xi<min(tjDataX))=[];
        fill(axU,[xi,xi(end:-1:1)],[f,-f(end:-1:1)].*1.2+i,colorList(i,:),'FaceAlpha',0.3,...
            'EdgeColor',colorList(i,:),'LineWidth',1.2);
        outliBool=isoutlier(tjDataX,'quartiles');
        outli=tjDataX(outliBool);
        scatter(axU,outli,repmat(i,[sum(outliBool),1]),45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.8)
        tjDataX(outliBool)=[];
        qt25=quantile(fullDataX,0.25);
        qt75=quantile(fullDataX,0.75);
        med=median(fullDataX);
        plot(axU,[max(tjDataX),qt75],[i,i],'LineWidth',1.2,'Color',[0,0,0])
        plot(axU,[min(tjDataX),qt25],[i,i],'LineWidth',1.2,'Color',[0,0,0])
        fill(axU,[qt25,qt25,qt75,qt75],i+.25.*1.*[-1 1 1 -1],[1,1,1],...
            'FaceAlpha',0.95,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        plot(axU,[med,med],[-.25,.25]+i,'LineWidth',1.2,'Color',colorList(i,:))
        axU.YLim=[0,length(PntSet)+1];
    case 8 % 线条散点
        LY=[tPntSet(:,2),tPntSet(:,2),tPntSet(:,2).*nan]';
        LX=repmat([i-.3,i+.3,nan],[size(tPntSet,1),1])';
        line(axR,LX(:),LY(:),'Color',[colorList(i,:),.4],'lineWidth',1)
        axR.XLim=[0,length(PntSet)+1];
        LX=[tPntSet(:,1),tPntSet(:,1),tPntSet(:,1).*nan]';
        LY=repmat([i-.3,i+.3,nan],[size(tPntSet,1),1])';
        line(axU,LX(:),LY(:),'Color',[colorList(i,:),.4],'lineWidth',1)
        axU.YLim=[0,length(PntSet)+1];
    case 9 % 山脊图
        [f,yi]=ksdensity(tPntSet(:,2));
        fill(axR,[f,0]+(i-1).*.5,[yi,yi(1)],colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor','none');
        plot(axR,f+(i-1).*.5,yi,'Color',colorList(i,:),'LineWidth',1.2)
        [f,xi]=ksdensity(tPntSet(:,1));
        fill(axU,[xi,xi(1)],[f,0]+(i-1).*.5,colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor','none');
        plot(axU,xi,f+(i-1).*.5,'Color',colorList(i,:),'LineWidth',1.2)
    case 10 % 类小提琴图
        tjDataY=tPntSet(:,2);fullDataY=tjDataY;
        [f,yi]=ksdensity(tjDataY);
        f(yi>max(tjDataY))=[];yi(yi>max(tjDataY))=[];
        f(yi<min(tjDataY))=[];yi(yi<min(tjDataY))=[];
        fill(axR,[0,f,0].*1.2+i+.15,[yi(1),yi,yi(end)],colorList(i,:),'FaceAlpha',0.3,...
            'EdgeColor',colorList(i,:),'LineWidth',1.2);
        outliBool=isoutlier(tjDataY,'quartiles');
        outli=tjDataY(outliBool);
        scatter(axR,repmat(i-.15,[sum(outliBool),1]),outli,45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.8)
        tjDataY(outliBool)=[];
        qt25=quantile(fullDataY,0.25);
        qt75=quantile(fullDataY,0.75);
        med=median(fullDataY);
        plot(axR,[i,i]-.2,[max(tjDataY),qt75],'LineWidth',1.2,'Color',[0,0,0])
        plot(axR,[i,i]-.2,[min(tjDataY),qt25],'LineWidth',1.2,'Color',[0,0,0])
        fill(axR,i-.2+.2.*1.*[-1 1 1 -1],[qt25,qt25,qt75,qt75],[1,1,1],...
            'FaceAlpha',0.95,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        plot(axR,[-.2,.2]+i-.2,[med,med],'LineWidth',1.2,'Color',colorList(i,:))
        axR.XLim=[0,length(PntSet)+1];
        % -----------------------------------------------------------------
        tjDataX=tPntSet(:,1);fullDataX=tjDataX;
        [f,xi]=ksdensity(tjDataX);
        f(xi>max(tjDataX))=[];xi(xi>max(tjDataX))=[];
        f(xi<min(tjDataX))=[];xi(xi<min(tjDataX))=[];
        fill(axU,[xi(1),xi,xi(end)],[0,f,0].*1.2+i+.15,colorList(i,:),'FaceAlpha',0.3,...
            'EdgeColor',colorList(i,:),'LineWidth',1.2);
        outliBool=isoutlier(tjDataX,'quartiles');
        outli=tjDataX(outliBool);
        scatter(axU,outli,repmat(i-.2,[sum(outliBool),1]),45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.8)
        tjDataX(outliBool)=[];
        qt25=quantile(fullDataX,0.25);
        qt75=quantile(fullDataX,0.75);
        med=median(fullDataX);
        plot(axU,[max(tjDataX),qt75],[i,i]-.2,'LineWidth',1.2,'Color',[0,0,0])
        plot(axU,[min(tjDataX),qt25],[i,i]-.2,'LineWidth',1.2,'Color',[0,0,0])
        fill(axU,[qt25,qt25,qt75,qt75],i-.2+.25.*1.*[-1 1 1 -1],[1,1,1],...
            'FaceAlpha',0.95,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        plot(axU,[med,med],[-.25,.25]+i-.2,'LineWidth',1.2,'Color',colorList(i,:))
        axU.YLim=[0,length(PntSet)+1];
    case 11 % 雨云图
        tjDataY=tPntSet(:,2);fullDataY=tjDataY;
        [f,yi]=ksdensity(tjDataY);
        f(yi>max(tjDataY))=[];yi(yi>max(tjDataY))=[];
        f(yi<min(tjDataY))=[];yi(yi<min(tjDataY))=[];
        fill(axR,[0,f,0].*1.2+i+.15,[yi(1),yi,yi(end)],colorList(i,:),'FaceAlpha',0.3,...
            'EdgeColor',colorList(i,:),'LineWidth',1.2);
        scatter(axR,tPntSet(:,2).*0+i-.2+.22.*(rand(size(tPntSet,1),1)-.5).*2,tPntSet(:,2),45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.3)
        tjDataY(outliBool)=[];
        qt25=quantile(fullDataY,0.25);
        qt75=quantile(fullDataY,0.75);
        med=median(fullDataY);
        plot(axR,[i,i]-.2,[max(tjDataY),qt75],'LineWidth',1.2,'Color',[0,0,0])
        plot(axR,[i,i]-.2,[min(tjDataY),qt25],'LineWidth',1.2,'Color',[0,0,0])
        fill(axR,i-.2+.2.*1.*[-1 1 1 -1],[qt25,qt25,qt75,qt75],[1,1,1],...
            'FaceAlpha',0.95,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        plot(axR,[-.2,.2]+i-.2,[med,med],'LineWidth',1.2,'Color',colorList(i,:))
        axR.XLim=[0,length(PntSet)+1];
        % -----------------------------------------------------------------
        tjDataX=tPntSet(:,1);fullDataX=tjDataX;
        [f,xi]=ksdensity(tjDataX);
        f(xi>max(tjDataX))=[];xi(xi>max(tjDataX))=[];
        f(xi<min(tjDataX))=[];xi(xi<min(tjDataX))=[];
        fill(axU,[xi(1),xi,xi(end)],[0,f,0].*1.2+i+.15,colorList(i,:),'FaceAlpha',0.3,...
            'EdgeColor',colorList(i,:),'LineWidth',1.2);
        scatter(axU,tPntSet(:,1),tPntSet(:,1).*0+i-.2+.22.*(rand(size(tPntSet,1),1)-.5),45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.3)
        tjDataX(outliBool)=[];
        qt25=quantile(fullDataX,0.25);
        qt75=quantile(fullDataX,0.75);
        med=median(fullDataX);
        plot(axU,[max(tjDataX),qt75],[i,i]-.2,'LineWidth',1.2,'Color',[0,0,0])
        plot(axU,[min(tjDataX),qt25],[i,i]-.2,'LineWidth',1.2,'Color',[0,0,0])
        fill(axU,[qt25,qt25,qt75,qt75],i-.2+.25.*1.*[-1 1 1 -1],[1,1,1],...
            'FaceAlpha',0.95,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        plot(axU,[med,med],[-.25,.25]+i-.2,'LineWidth',1.2,'Color',colorList(i,:))
        axU.YLim=[0,length(PntSet)+1];
end
end
linkaxes([axM,axR],'y')
linkaxes([axM,axU],'x')
%% ========================================================================
% 所使用到的函数
% 置信椭圆定位函数
function [X,Y]=getEllipse(Mu,Sigma,S,pntNum)
% 置信区间 | 95%:5.991  99%:9.21  90%:4.605
% (X-Mu)*inv(Sigma)*(X-Mu)=S

invSig=inv(Sigma);

[V,D]=eig(invSig);
aa=sqrt(S/D(1));
bb=sqrt(S/D(4));

t=linspace(0,2*pi,pntNum);
XY=V*[aa*cos(t);bb*sin(t)];
X=(XY(1,:)+Mu(1))';
Y=(XY(2,:)+Mu(2))';
end

边际图格式展示：

1 直方图

2 核密度填充

3 核密度线

4 核密度线+填充

5 直方图+密度线

6 箱线图

7 小提琴图

8 散点图

9 山脊图

10 类小提琴图

11 雨云图

---

完整代码

% 边际图示例
% marginalPlotDemo
% @author : slandarer
% Zhaoxu Liu / slandarer (2023). marginal plot 
% (https://www.mathworks.com/matlabcentral/fileexchange/123470-marginal-plot), 
% MATLAB Central File Exchange. 检索来源 2023/1/19.


% 构造三个符合高斯分布的点集
PntSet1=mvnrnd([2 3],[1 0;0 2],300);
PntSet2=mvnrnd([6 7],[1 0;0 2],300);
PntSet3=mvnrnd([14 9],[1 0;0 1],300);
% 将数据放进元胞数组
PntSet={PntSet1,PntSet2,PntSet3};
% -------------------------------------------------------------------------
% 配色预设(here to change the default color list)
C1=[211 43 43;61 96 137;249 206 61;76 103 86;80 80 80]./255;
C2=[102,173,194;36,59,66;232,69,69;194,148,102;54,43,33]./255;
C3=[223,122,94;60 64 91;130 178 154;244,241,222;240 201 134]./255;
C4=[122,117,119;255,163,25;135,146,73;126,15,4;30,93,134]./255;
C5=[198,199,201;38,74,96;209,80,51;241,174,44;12,13,15]./255;
C6=[235,75,55;77,186,216;58,84,141;2,162,136;245,155,122]./255;
C7=[23,23,23;121,17,36;31,80,91;44,9,75;61,36,42]./255;
C8=[47,62,66;203,129,70;0 64 115;152,58,58;20 72 83]./255;
colorList=C2;
% -------------------------------------------------------------------------
% 在这改边缘图及中心图种类
% (here to change the type of marginal plot and main plot)
%    |
%    |
%    |
% \  |  /
%  \ | /
%   \|/
%     
mainType=1;
marginType=1;
%  
%   /|\
%  / | \
% /  |  \
%    |
%    |
%    |
%% =========================================================================
% figure图床及主要区域axes坐标区域创建及基础修饰
fig=figure('Units','normalized','Position',[.3,.2,.5,.63]);
axM=axes('Parent',fig);hold on;
set(axM,'Position',[.08,.08,.65,.65],'LineWidth',1.1,'Box','on','TickDir','in',...
    'XMinorTick','on','YMinorTick','on','XGrid','on','YGrid','on','GridLineStyle','--',...
    'FontName','Times New Roman','FontSize',12,'GridAlpha',.09)
axM.XLabel.String='Main XXXXX';
axM.YLabel.String='Main YYYYY';
axM.XLabel.FontSize=14;
axM.YLabel.FontSize=14;
% -------------------------------------------------------------------------
% 右侧axes坐标区域创建及基础修饰
axR=axes('Parent',fig);hold on;
set(axR,'Position',[.75,.08,.23,.65],'LineWidth',1.1,'Box','off','TickDir','out',...
    'XMinorTick','on','YMinorTick','on','XGrid','on','YGrid','on','GridLineStyle','--',...
    'FontName','Times New Roman','FontSize',12,'GridAlpha',.09,'TickLength',[.006 .015],'YTickLabel','')
axR.XLabel.String='Y-axis data statistics';
axR.XLabel.FontSize=14;
% -------------------------------------------------------------------------
% 上方axes坐标区域创建及基础修饰
axU=axes('Parent',fig);hold on;
set(axU,'Position',[.08,.76,.65,.2],'LineWidth',1.1,'Box','off','TickDir','out',...
    'XMinorTick','on','YMinorTick','on','XGrid','on','YGrid','on','GridLineStyle','--',...
    'FontName','Times New Roman','FontSize',12,'GridAlpha',.09,'TickLength',[.006 .015],'XTickLabel','')
axU.YLabel.String='X-axis data statistics';
axU.YLabel.FontSize=14;
%% =========================================================================
% 主要区域散点图绘制
switch mainType
    case 1 % 散点图
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            scatter(axM,tPntSet(:,1),tPntSet(:,2),70,'filled','CData',colorList(i,:),'MarkerFaceAlpha',.6)
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,legendStr,'FontSize',14,'Box','off','Location','best');
    case 2 % 等高线图
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            gridx1=linspace(min(tPntSet(:,1)),max(tPntSet(:,1)),round(sqrt(size(tPntSet,1))));
            gridx2=linspace(min(tPntSet(:,2)),max(tPntSet(:,2)),round(sqrt(size(tPntSet,1))));
            [x1,x2]=meshgrid(gridx1, gridx2);
            [m,n]=size(x1);
            x1=x1(:);
            x2=x2(:);
            xi=[x1,x2];
            h=ksdensity(tPntSet,xi);
            contour(axM,reshape(x1,[m,n]),reshape(x2,[m,n]),reshape(h,[m,n]),...
                'LineWidth',1.2,'EdgeColor',colorList(i,:))
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,legendStr,'FontSize',14,'Box','off','Location','best');
    case 3 % 回归曲线
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            [p,S]=polyfit(tPntSet(:,1),tPntSet(:,2),1); 
            x1=linspace(min(tPntSet(:,1)),max(tPntSet(:,1)),100);
            [y_fit,delta]=polyval(p,x1,S);
            % 绘制原始数据、线性拟合和 95% 预测区间 y±2Δ。
            uy=y_fit+2.*delta;
            dy=y_fit-2.*delta;
            % 绘制原始数据
            sHdl(i)=scatter(axM,tPntSet(:,1),tPntSet(:,2),30,'filled','CData',colorList(i,:),'MarkerFaceAlpha',.9);
            % 绘制拟合曲线
            plot(axM,x1,y_fit,'Color',colorList(i,:),'LineWidth',2.5)
            % 绘制置信区间
            fill(axM,[x1,x1(end:-1:1)],[uy,dy(end:-1:1)],colorList(i,:),'EdgeColor','none','FaceAlpha',.15)
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,sHdl,legendStr,'FontSize',14,'Box','off','Location','best');
    case 4 % 散点+凸包
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            [k,av]=convhull(tPntSet);
            sHdl(i)=scatter(axM,tPntSet(:,1),tPntSet(:,2),30,'filled','CData',colorList(i,:),'MarkerFaceAlpha',.9);
            fill(axM,tPntSet(k,1),tPntSet(k,2),colorList(i,:),'FaceAlpha',0.3,'EdgeColor',colorList(i,:),'LineWidth',2)
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,sHdl,legendStr,'FontSize',14,'Box','off','Location','best');
    case 5 % 散点+不规则图形
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            sHdl(i)=scatter(axM,tPntSet(:,1),tPntSet(:,2),30,'filled','CData',colorList(i,:),'MarkerFaceAlpha',.9);
            R=max(max(tPntSet(:,1))-min(tPntSet(:,1)),max(tPntSet(:,2))-min(tPntSet(:,2)))./round(sqrt(size(tPntSet,1)));
            t=linspace(0,2*pi,100);
            [k,~]=convhull(tPntSet);
            kPntSet=tPntSet(k,:);
            kX=repmat(kPntSet(:,1),[1,100])+repmat(cos(t).*R,[length(k),1]);
            kY=repmat(kPntSet(:,2),[1,100])+repmat(sin(t).*R,[length(k),1]);
            kX=kX(:);kY=kY(:);
            [k,av]=convhull([kX,kY]);
            fill(axM,kX(k),kY(k),colorList(i,:),'FaceAlpha',0.3,'EdgeColor',colorList(i,:),'LineWidth',2)
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,sHdl,legendStr,'FontSize',14,'Box','off','Location','best');
    case 6 % 椭圆形
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            sHdl(i)=scatter(axM,tPntSet(:,1),tPntSet(:,2),30,'filled','CData',colorList(i,:),'MarkerFaceAlpha',.9);
            Mu=mean(tPntSet);
            Y=tPntSet-repmat(Mu,size(tPntSet,1),1);
            Sigma=cov(Y);% Sigma=(Y'*(ones(size(pntSet,1),1).*Y))./(size(pntSet,1)-1);
            [X,Y]=getEllipse(Mu,Sigma,9.21,100);
            fill(axM,X,Y,colorList(i,:),'FaceAlpha',0.3,'EdgeColor',colorList(i,:),'LineWidth',2);
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,sHdl,legendStr,'FontSize',14,'Box','off','Location','best');
    case 7 % 簇状
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            sHdl(i)=scatter(axM,tPntSet(:,1),tPntSet(:,2),30,'filled','CData',colorList(i,:),'MarkerFaceAlpha',.9);
            Mu=mean(tPntSet);
            LX=[tPntSet(:,1),Mu(1).*ones(size(tPntSet,1),1),nan(size(tPntSet,1),1)]';
            LY=[tPntSet(:,2),Mu(2).*ones(size(tPntSet,1),1),nan(size(tPntSet,1),1)]';
            plot(axM,LX(:),LY(:),'Color',colorList(i,:),'LineWidth',1.2);
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,sHdl,legendStr,'FontSize',14,'Box','off','Location','best');
    case 8 % 误差棒
        for i=1:length(PntSet)
            tPntSet=PntSet{i};
            Mu=mean(tPntSet);
            SX=std(tPntSet(:,1));
            SY=std(tPntSet(:,2));
            errorbar(axM,Mu(1),Mu(2),SY,SY,SX,SX,'o','LineWidth',1.5,'Color',colorList(i,:),'MarkerSize',10);
        end
        legendStr{length(PntSet)}='';
        for i=1:length(PntSet)
            legendStr{i}=['Class-',num2str(i)];
        end
        legend(axM,legendStr,'FontSize',14,'Box','off','Location','best');
end
% -------------------------------------------------------------------------
% 边际图绘制
axR.YLim=axM.YLim;
axU.XLim=axM.XLim;
for i=1:length(PntSet)
tPntSet=PntSet{i};
switch marginType
    case 1 % 柱状图
        histogram(axR,tPntSet(:,2),'FaceColor',colorList(i,:),'Orientation','horizontal',...
            'EdgeColor',[1 1 1]*.1,'FaceAlpha',0.6,'LineWidth',.8)
        histogram(axU,tPntSet(:,1),'FaceColor',colorList(i,:),...
            'EdgeColor',[1 1 1]*.1,'FaceAlpha',0.6,'LineWidth',.8)
    case 2 % 核密度面
        [f,yi]=ksdensity(tPntSet(:,2));
        fill(axR,[f,0],[yi,yi(1)],colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor','none')
        [f,xi]=ksdensity(tPntSet(:,1));
        fill(axU,[xi,xi(1)],[f,0],colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor','none')
    case 3 % 核密度线
        [f,yi]=ksdensity(tPntSet(:,2));
        plot(axR,f,yi,'Color',colorList(i,:),'LineWidth',2)
        [f,xi]=ksdensity(tPntSet(:,1));
        plot(axU,xi,f,'Color',colorList(i,:),'LineWidth',2)
    case 4 % 核密度面+线
        [f,yi]=ksdensity(tPntSet(:,2));
        fill(axR,[f,0],[yi,yi(1)],colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor','none');
        plot(axR,f,yi,'Color',colorList(i,:),'LineWidth',1.2)
        [f,xi]=ksdensity(tPntSet(:,1));
        fill(axU,[xi,xi(1)],[f,0],colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor','none');
        plot(axU,xi,f,'Color',colorList(i,:),'LineWidth',1.2)
    case 5 % 直方图+密度线
        HR=histogram(axR,tPntSet(:,2),'FaceColor',colorList(i,:),'Orientation','horizontal',...
            'EdgeColor',[1 1 1]*.1,'FaceAlpha',0.6,'LineWidth',.8);
        HU=histogram(axU,tPntSet(:,1),'FaceColor',colorList(i,:),...
            'EdgeColor',[1 1 1]*.1,'FaceAlpha',0.6,'LineWidth',.8);
        [f,yi]=ksdensity(tPntSet(:,2));
        plot(axR,f.*size(tPntSet,1).*HR.BinWidth,yi,'Color',colorList(i,:),'LineWidth',2)
        [f,xi]=ksdensity(tPntSet(:,1));
        plot(axU,xi,f.*size(tPntSet,1).*HU.BinWidth,'Color',colorList(i,:),'LineWidth',2)
    case 6 % 箱线图
        tjDataY=tPntSet(:,2);fullDataY=tjDataY;
        outliBool=isoutlier(tjDataY,'quartiles');
        outli=tjDataY(outliBool);
        scatter(axR,repmat(i,[sum(outliBool),1]),outli,45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.8)
        tjDataY(outliBool)=[];
        qt25=quantile(fullDataY,0.25);
        qt75=quantile(fullDataY,0.75);
        med=median(fullDataY);
        plot(axR,[i,i],[max(tjDataY),qt75],'LineWidth',1.2,'Color',[0,0,0])
        plot(axR,[i,i],[min(tjDataY),qt25],'LineWidth',1.2,'Color',[0,0,0])
        plot(axR,[-.15,.15]+i,[max(tjDataY),max(tjDataY)],'LineWidth',1.2,'Color',[0,0,0])
        plot(axR,[-.15,.15]+i,[min(tjDataY),min(tjDataY)],'LineWidth',1.2,'Color',[0,0,0])
        plot(axR,[-.25,.25]+i,[med,med],'LineWidth',1.2,'Color',colorList(i,:))
        fill(axR,i+.25.*1.*[-1 1 1 -1],[qt25,qt25,qt75,qt75],colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        axR.XLim=[0,length(PntSet)+1];
        % -----------------------------------------------------------------
        tjDataX=tPntSet(:,1);fullDataX=tjDataX;
        outliBool=isoutlier(tjDataX,'quartiles');
        outli=tjDataX(outliBool);
        scatter(axU,outli,repmat(i,[sum(outliBool),1]),45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.8)
        tjDataX(outliBool)=[];
        qt25=quantile(fullDataX,0.25);
        qt75=quantile(fullDataX,0.75);
        med=median(fullDataX);
        plot(axU,[max(tjDataX),qt75],[i,i],'LineWidth',1.2,'Color',[0,0,0])
        plot(axU,[min(tjDataX),qt25],[i,i],'LineWidth',1.2,'Color',[0,0,0])
        plot(axU,[max(tjDataX),max(tjDataX)],[-.15,.15]+i,'LineWidth',1.2,'Color',[0,0,0])
        plot(axU,[min(tjDataX),min(tjDataX)],[-.15,.15]+i,'LineWidth',1.2,'Color',[0,0,0])
        plot(axU,[med,med],[-.25,.25]+i,'LineWidth',1.2,'Color',colorList(i,:))
        fill(axU,[qt25,qt25,qt75,qt75],i+.25.*1.*[-1 1 1 -1],colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        axU.YLim=[0,length(PntSet)+1];
    case 7 % 小提琴图
        tjDataY=tPntSet(:,2);fullDataY=tjDataY;
        [f,yi]=ksdensity(tjDataY);
        f(yi>max(tjDataY))=[];yi(yi>max(tjDataY))=[];
        f(yi<min(tjDataY))=[];yi(yi<min(tjDataY))=[];
        fill(axR,[f,-f(end:-1:1)].*1.2+i,[yi,yi(end:-1:1)],colorList(i,:),'FaceAlpha',0.3,...
            'EdgeColor',colorList(i,:),'LineWidth',1.2);
        outliBool=isoutlier(tjDataY,'quartiles');
        outli=tjDataY(outliBool);
        scatter(axR,repmat(i,[sum(outliBool),1]),outli,45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.8)
        tjDataY(outliBool)=[];
        qt25=quantile(fullDataY,0.25);
        qt75=quantile(fullDataY,0.75);
        med=median(fullDataY);
        plot(axR,[i,i],[max(tjDataY),qt75],'LineWidth',1.2,'Color',[0,0,0])
        plot(axR,[i,i],[min(tjDataY),qt25],'LineWidth',1.2,'Color',[0,0,0])
        fill(axR,i+.2.*1.*[-1 1 1 -1],[qt25,qt25,qt75,qt75],[1,1,1],...
            'FaceAlpha',0.95,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        plot(axR,[-.2,.2]+i,[med,med],'LineWidth',1.2,'Color',colorList(i,:))
        axR.XLim=[0,length(PntSet)+1];
        % -----------------------------------------------------------------
        tjDataX=tPntSet(:,1);fullDataX=tjDataX;
        [f,xi]=ksdensity(tjDataX);
        f(xi>max(tjDataX))=[];xi(xi>max(tjDataX))=[];
        f(xi<min(tjDataX))=[];xi(xi<min(tjDataX))=[];
        fill(axU,[xi,xi(end:-1:1)],[f,-f(end:-1:1)].*1.2+i,colorList(i,:),'FaceAlpha',0.3,...
            'EdgeColor',colorList(i,:),'LineWidth',1.2);
        outliBool=isoutlier(tjDataX,'quartiles');
        outli=tjDataX(outliBool);
        scatter(axU,outli,repmat(i,[sum(outliBool),1]),45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.8)
        tjDataX(outliBool)=[];
        qt25=quantile(fullDataX,0.25);
        qt75=quantile(fullDataX,0.75);
        med=median(fullDataX);
        plot(axU,[max(tjDataX),qt75],[i,i],'LineWidth',1.2,'Color',[0,0,0])
        plot(axU,[min(tjDataX),qt25],[i,i],'LineWidth',1.2,'Color',[0,0,0])
        fill(axU,[qt25,qt25,qt75,qt75],i+.25.*1.*[-1 1 1 -1],[1,1,1],...
            'FaceAlpha',0.95,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        plot(axU,[med,med],[-.25,.25]+i,'LineWidth',1.2,'Color',colorList(i,:))
        axU.YLim=[0,length(PntSet)+1];
    case 8 % 线条散点
        LY=[tPntSet(:,2),tPntSet(:,2),tPntSet(:,2).*nan]';
        LX=repmat([i-.3,i+.3,nan],[size(tPntSet,1),1])';
        line(axR,LX(:),LY(:),'Color',[colorList(i,:),.4],'lineWidth',1)
        axR.XLim=[0,length(PntSet)+1];
        LX=[tPntSet(:,1),tPntSet(:,1),tPntSet(:,1).*nan]';
        LY=repmat([i-.3,i+.3,nan],[size(tPntSet,1),1])';
        line(axU,LX(:),LY(:),'Color',[colorList(i,:),.4],'lineWidth',1)
        axU.YLim=[0,length(PntSet)+1];
    case 9 % 山脊图
        [f,yi]=ksdensity(tPntSet(:,2));
        fill(axR,[f,0]+(i-1).*.5,[yi,yi(1)],colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor','none');
        plot(axR,f+(i-1).*.5,yi,'Color',colorList(i,:),'LineWidth',1.2)
        [f,xi]=ksdensity(tPntSet(:,1));
        fill(axU,[xi,xi(1)],[f,0]+(i-1).*.5,colorList(i,:),...
            'FaceAlpha',0.3,'EdgeColor','none');
        plot(axU,xi,f+(i-1).*.5,'Color',colorList(i,:),'LineWidth',1.2)
    case 10 % 类小提琴图
        tjDataY=tPntSet(:,2);fullDataY=tjDataY;
        [f,yi]=ksdensity(tjDataY);
        f(yi>max(tjDataY))=[];yi(yi>max(tjDataY))=[];
        f(yi<min(tjDataY))=[];yi(yi<min(tjDataY))=[];
        fill(axR,[0,f,0].*1.2+i+.15,[yi(1),yi,yi(end)],colorList(i,:),'FaceAlpha',0.3,...
            'EdgeColor',colorList(i,:),'LineWidth',1.2);
        outliBool=isoutlier(tjDataY,'quartiles');
        outli=tjDataY(outliBool);
        scatter(axR,repmat(i-.15,[sum(outliBool),1]),outli,45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.8)
        tjDataY(outliBool)=[];
        qt25=quantile(fullDataY,0.25);
        qt75=quantile(fullDataY,0.75);
        med=median(fullDataY);
        plot(axR,[i,i]-.2,[max(tjDataY),qt75],'LineWidth',1.2,'Color',[0,0,0])
        plot(axR,[i,i]-.2,[min(tjDataY),qt25],'LineWidth',1.2,'Color',[0,0,0])
        fill(axR,i-.2+.2.*1.*[-1 1 1 -1],[qt25,qt25,qt75,qt75],[1,1,1],...
            'FaceAlpha',0.95,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        plot(axR,[-.2,.2]+i-.2,[med,med],'LineWidth',1.2,'Color',colorList(i,:))
        axR.XLim=[0,length(PntSet)+1];
        % -----------------------------------------------------------------
        tjDataX=tPntSet(:,1);fullDataX=tjDataX;
        [f,xi]=ksdensity(tjDataX);
        f(xi>max(tjDataX))=[];xi(xi>max(tjDataX))=[];
        f(xi<min(tjDataX))=[];xi(xi<min(tjDataX))=[];
        fill(axU,[xi(1),xi,xi(end)],[0,f,0].*1.2+i+.15,colorList(i,:),'FaceAlpha',0.3,...
            'EdgeColor',colorList(i,:),'LineWidth',1.2);
        outliBool=isoutlier(tjDataX,'quartiles');
        outli=tjDataX(outliBool);
        scatter(axU,outli,repmat(i-.2,[sum(outliBool),1]),45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.8)
        tjDataX(outliBool)=[];
        qt25=quantile(fullDataX,0.25);
        qt75=quantile(fullDataX,0.75);
        med=median(fullDataX);
        plot(axU,[max(tjDataX),qt75],[i,i]-.2,'LineWidth',1.2,'Color',[0,0,0])
        plot(axU,[min(tjDataX),qt25],[i,i]-.2,'LineWidth',1.2,'Color',[0,0,0])
        fill(axU,[qt25,qt25,qt75,qt75],i-.2+.25.*1.*[-1 1 1 -1],[1,1,1],...
            'FaceAlpha',0.95,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        plot(axU,[med,med],[-.25,.25]+i-.2,'LineWidth',1.2,'Color',colorList(i,:))
        axU.YLim=[0,length(PntSet)+1];
    case 11 % 雨云图
        tjDataY=tPntSet(:,2);fullDataY=tjDataY;
        [f,yi]=ksdensity(tjDataY);
        f(yi>max(tjDataY))=[];yi(yi>max(tjDataY))=[];
        f(yi<min(tjDataY))=[];yi(yi<min(tjDataY))=[];
        fill(axR,[0,f,0].*1.2+i+.15,[yi(1),yi,yi(end)],colorList(i,:),'FaceAlpha',0.3,...
            'EdgeColor',colorList(i,:),'LineWidth',1.2);
        scatter(axR,tPntSet(:,2).*0+i-.2+.22.*(rand(size(tPntSet,1),1)-.5).*2,tPntSet(:,2),45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.3)
        tjDataY(outliBool)=[];
        qt25=quantile(fullDataY,0.25);
        qt75=quantile(fullDataY,0.75);
        med=median(fullDataY);
        plot(axR,[i,i]-.2,[max(tjDataY),qt75],'LineWidth',1.2,'Color',[0,0,0])
        plot(axR,[i,i]-.2,[min(tjDataY),qt25],'LineWidth',1.2,'Color',[0,0,0])
        fill(axR,i-.2+.2.*1.*[-1 1 1 -1],[qt25,qt25,qt75,qt75],[1,1,1],...
            'FaceAlpha',0.95,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        plot(axR,[-.2,.2]+i-.2,[med,med],'LineWidth',1.2,'Color',colorList(i,:))
        axR.XLim=[0,length(PntSet)+1];
        % -----------------------------------------------------------------
        tjDataX=tPntSet(:,1);fullDataX=tjDataX;
        [f,xi]=ksdensity(tjDataX);
        f(xi>max(tjDataX))=[];xi(xi>max(tjDataX))=[];
        f(xi<min(tjDataX))=[];xi(xi<min(tjDataX))=[];
        fill(axU,[xi(1),xi,xi(end)],[0,f,0].*1.2+i+.15,colorList(i,:),'FaceAlpha',0.3,...
            'EdgeColor',colorList(i,:),'LineWidth',1.2);
        scatter(axU,tPntSet(:,1),tPntSet(:,1).*0+i-.2+.22.*(rand(size(tPntSet,1),1)-.5),45,'filled',...
            'CData',colorList(i,:),'LineWidth',1.2,'MarkerFaceAlpha',0.3)
        tjDataX(outliBool)=[];
        qt25=quantile(fullDataX,0.25);
        qt75=quantile(fullDataX,0.75);
        med=median(fullDataX);
        plot(axU,[max(tjDataX),qt75],[i,i]-.2,'LineWidth',1.2,'Color',[0,0,0])
        plot(axU,[min(tjDataX),qt25],[i,i]-.2,'LineWidth',1.2,'Color',[0,0,0])
        fill(axU,[qt25,qt25,qt75,qt75],i-.2+.25.*1.*[-1 1 1 -1],[1,1,1],...
            'FaceAlpha',0.95,'EdgeColor',colorList(i,:),'LineWidth',1.2);
        plot(axU,[med,med],[-.25,.25]+i-.2,'LineWidth',1.2,'Color',colorList(i,:))
        axU.YLim=[0,length(PntSet)+1];
end
end
linkaxes([axM,axR],'y')
linkaxes([axM,axU],'x')
%% ========================================================================
% 所使用到的函数
% 置信椭圆定位函数
function [X,Y]=getEllipse(Mu,Sigma,S,pntNum)
% 置信区间 | 95%:5.991  99%:9.21  90%:4.605
% (X-Mu)*inv(Sigma)*(X-Mu)=S

invSig=inv(Sigma);

[V,D]=eig(invSig);
aa=sqrt(S/D(1));
bb=sqrt(S/D(4));

t=linspace(0,2*pi,pntNum);
XY=V*[aa*cos(t);bb*sin(t)];
X=(XY(1,:)+Mu(1))';
Y=(XY(2,:)+Mu(2))';
end
% Zhaoxu Liu / slandarer (2023). marginal plot 
% (https://www.mathworks.com/matlabcentral/fileexchange/123470-marginal-plot), 
% MATLAB Central File Exchange. 检索来源 2023/1/19.