首頁 > 軟體

Matlab實現帶豎線散點的核密度圖的繪製

2022-08-21 18:00:27

帶豎線散點的核密度圖繪製模板

今天帶來一個帶豎線散點的核密度圖繪製模板,作為模板,想要畫出屬於自己的圖,需要修改的就只有資料,格式,配色三樣,以下給出模板使用詳細講解:

part1 資料

這裡的資料需要通過Data(n).X=...的格式設定,例如,三組資料:

% 匯入資料,更多的資料也請使用
% 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(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 格式

有倆基礎格式:

% 一些基礎設定
scatterSep='off'; % 是否分開繪製豎線散點
totalRatio='off'; % 是否各組按比例繪製

totalRatio屬效能調節各組影象是否各組按比例繪製,如果將其設定為on

scatterSep屬效能調節是否分開繪製豎線散點,如果將其設定為on

part3 配色

這裡準備了八組資料,只需要修改colorList=Cn即可:

% 配色列表
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 模板完整程式碼

% 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);

到此這篇關於Matlab實現帶豎線散點的核密度圖的繪製的文章就介紹到這了,更多相關Matlab核密度圖內容請搜尋it145.com以前的文章或繼續瀏覽下面的相關文章希望大家以後多多支援it145.com!


IT145.com E-mail:sddin#qq.com