带竖线散点的核密度图绘制模板
今天带来一个带竖线散点的核密度图绘制模板，作为模板，想要画出属于自己的图，需要修改的就只有数据，格式，配色三样，以下给出模板使用详细讲解：
part1 数据
这里的数据需要通过Data(n).X=...的格式设置，例如，三组数据：
|
1 2 3 4 5 |
% 导入数据，更多的数据也请使用 % Data(n).X=......的格式 Data(1).X=mvnrnd(40,60,300); Data(2).X=mvnrnd(60,60,600); Data(3).X=mvnrnd(80,60,900); |
五组数据(数据组数超过五组则后期需要设置更长的配色)：
|
1 2 3 4 5 6 7 |
% 导入数据，更多的数据也请使用 % Data(n).X=......的格式 Data(1).X=mvnrnd(40,60,300); Data(2).X=mvnrnd(60,60,600); Data(3).X=mvnrnd(80,60,900); Data(4).X=mvnrnd(100,60,1200); Data(5).X=mvnrnd(120,60,1200); |
这部分代码里用的是随机生成的正态分布，可以将其更改为自己的数据。
part2 格式
有俩基础格式：
|
1 2 3 |
% 一些基础设置 scatterSep='off'; % 是否分开绘制竖线散点 totalRatio='off'; % 是否各组按比例绘制 |
totalRatio属性能调节各组图像是否各组按比例绘制，如果将其设置为on：
scatterSep属性能调节是否分开绘制竖线散点，如果将其设置为on：
part3 配色
这里准备了八组数据，只需要修改colorList=Cn即可：
|
1 2 3 4 5 6 7 8 9 10 |
% 配色列表 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=[244,241,222;223,122,94;60 64 91;130 178 154;240 201 134]./255; C4=[126,15,4;122,117,119;255,163,25;135,146,73;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;2,162,136;58,84,141;245,155,122]./255; C7=[23,23,23;121,17,36;44,9,75;31,80,91;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:
part4 模板完整代码
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 |
% line area % @author : slandarer % @公众号 : slandarer随笔
clc;clear
% 导入数据，更多的数据也请使用 % Data(n).X=......的格式 Data(1).X=mvnrnd(40,60,300); Data(2).X=mvnrnd(60,60,600); Data(3).X=mvnrnd(80,60,900); Data(4).X=mvnrnd(100,60,1200); Data(5).X=mvnrnd(120,60,1200);
% 一些基础设置 scatterSep='off'; % 是否分开绘制竖线散点 totalRatio='on'; % 是否各组按比例绘制
% 配色列表 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=[244,241,222;223,122,94;60 64 91;130 178 154;240 201 134]./255; C4=[126,15,4;122,117,119;255,163,25;135,146,73;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;2,162,136;58,84,141;245,155,122]./255; C7=[23,23,23;121,17,36;44,9,75;31,80,91;61,36,42]./255; C8=[47,62,66;203,129,70;0 64 115;152,58,58;20 72 83]./255; colorList=C2;
% =========================================================================
% 图像绘制 ax=gca;hold on N=length(Data); areaHdl(N)=nan; lgdStrs{N}='';
% 计算各类数据量 K=arrayfun(@(x) length(x.X),Data); % 循环绘图 for n=1:N [f,xi]=ksdensity(Data(n).X); if strcmp(totalRatio,'on') f=f.*K(n)./sum(K); end areaHdl(n)=area(xi,f,'FaceColor',colorList(n,:),... 'EdgeColor',colorList(n,:),'FaceAlpha',.5,'LineWidth',1.5); lgdStrs{n}=['Group ',num2str(n)]; end
% 绘制图例 lgd=legend(areaHdl,lgdStrs{:}); lgd.AutoUpdate='off'; lgd.Location='best';
% 调整轴范围 posSep=ax.YLim(2)-0; if strcmp(scatterSep,'on') ax.YLim(1)=-posSep/6*N; else ax.YLim(1)=-posSep/6; end ax.XLim=ax.XLim; totalSep=diff(ax.YLim);
for n=1:N dy=strcmp(scatterSep,'on'); LY=ones(1,K(n)).*[(-posSep/6).*(.1+dy.*(n-1));(-posSep/6.)*(.9+dy.*(n-1));nan]; LX=[Data(n).X(:)';Data(n).X(:)';ones(1,K(n)).*nan]; line(LX(:),LY(:),'Color',[colorList(n,:),.4],'lineWidth',1) end
% 坐标区域修饰 ax.Box='on'; ax.BoxStyle='full'; ax.LineWidth=1; ax.FontSize=11; ax.FontName='Arial'; ax.TickDir='out'; ax.TickLength=[.005,.1]; ax.YTick(ax.YTick<-eps)=[]; ax.Title.String='area plot with | scatter'; ax.Title.FontSize=14; ax.XLabel.String='XXXXX'; ax.YLabel.String='YYYYY';
% 绘制基准线及框线 fplot(@(t)t.*0,'Color',ax.XColor,'LineWidth',ax.LineWidth); |
