gplotmatrix2

[haxes,hax_big] = gplotmatrix2(X,Y,group1,group2) plots the given X and
Y data, each Nx1 vectors, according to grouping variables group1 and
group2. gplotmatrix2 splits the X and Y into g1 x g2 subplots, where g1
is the number of groups in group1, and g2 is the number of groups in
group2.
[haxes,hax_big] = gplotmatrix2(X,Y,group1,group2,xlbl,ylbl,legendlbl)
accepts optional parameters xlbl, ylbl, and legendlbl. xlbl and ylbl
describe the x- and y- labels, respectively. legendlbl is the title
text object for the legend.
Example:
idx1 = randi([1 5], 1000, 1);
idx2 = randi([1 15], 1000, 1);
catnames1 = {'0''', '5''', '30''', '60''', '180'''};
catnames2 = {'a', 'b', 'c', 'd', 'e', ...
'f', 'g', 'h', 'i', 'j', ...
'k', 'l', 'm', 'n', 'o'};
T = categorical(idx1, 1:length(catnames1), catnames1); % group 1
M = categorical(idx2, 1:length(catnames2), catnames2); % group 2
X = 100*randn(1000,1);
Y = 20*randn(1000,1);
gplotmatrix2(X,Y,M,T,'Cell Area','Cytoplasm EGFR','Time min')

See also GPLOTMATRIX, GSCATTER, GSCATTER2.