gscatter2

[haxes,hax_big] = gscatter2(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 KxJ subplots, where the product of K and J is more than or equal to the number of groups in group1. K and J is automatically determined, and their product will be as close to the number of groups in group1 as possible (with the condition that K and J differ by at most 1).

[haxes,hax_big] = gscatter2(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);
gscatter2(X,Y,M,T,'Cell EFGR','Cell PIP2','Time min')

See also gplotmatrix, gscatter, gplotmatrix2.