Colour bar isn't correct for subplot contour

Question

Alistair McQueen 2022-2-23

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1656605-colour-bar-isn-t-correct-for-subplot-contour

评论： Alistair McQueen 2022-2-23

采纳的回答： Simon Chan

在 MATLAB Online 中打开

Hello,

Been looking through the forums for help with this, and found some solutions but did not fix my issue.

(Code attached at bottom). The z-values on the left plot do not quite match up with my globally defined colourbar.

I know contours, in this instance, work with discrete points, so unsure if the issue is with that.

I could plot two separate bars which address the issue (see commented out) but I would rather just have the one bar.

Is it possible I need to simulate more combinations such that the plot is more 'exact'?

Thanks!

%The action of drug on logsitic growth model - sensitivity analysis%
%Preamble
clc
clear all
clear workspace
%Parameters - Fixed:
%All the different initial conditions
MPC1 = [10000 000];
MPC2 = [7500 2500];
MPC3 = [5000 5000];
MPC4 = [2500 7500];
MPC5 = [000 10000];
g1 = 1.0392; %Phase 1 transiton
g2 = 1.0392; %Phase 2 transition
x0 = 10000; %Initial cell count
K = 100000; %Maximal cell count
time = linspace(0,15,1000); %time of simulation
%Parameters that vary:
cd = 1e-3;
kmax = [1 0.95 0.9 0.85 0.8 0.75 0.7 0.65 0.6 0.55 0.5];
k50 = [1e-3 7.5e-4 2.5e-4 1e-4 5e-5 1e-5 5e-6 1e-6 7.5e-7 5e-7 1e-7];
%Initial condition 1
for ii = 1:length(kmax)
    for kk = 1:length(k50)
        M = (cd*kmax(ii))/(cd+k50(kk));
        ode = @(t,x) [2*g2*x(2)*(1 - ((x(1)+x(2))/K)) - g1*x(1)*(1-M);...
            -g2*x(2)*(1 - ((x(1)+x(2))/K)) + g1*x(1)*(1-M)];
        [t, x(:,:,ii,kk)] = ode45(ode, time, MPC1);
    end
end
%Initial condition 5
for ii = 1:length(kmax)
    for kk = 1:length(k50)
        M = (cd*kmax(ii))/(cd+k50(kk));
        ode = @(t,x) [2*g2*x(2)*(1 - ((x(1)+x(2))/K)) - g1*x(1)*(1-M);...
            -g2*x(2)*(1 - ((x(1)+x(2))/K)) + g1*x(1)*(1-M)];
        [t, x5(:,:,ii,kk)] = ode45(ode, time, MPC5);
    end
end
%Data Extraction:
%This informs us how long our data set is
CellS = x(:,:);
CellS5 = x5(:,:);
%Loop length
LL = length(CellS(1,:))/2;
%Adding the first and second phase numbers together
counter = 0; %Initialise counter for initial grouping
for ii = 1:LL
    T_1(:,ii) = CellS(:,ii+counter) + CellS(:,(ii+1)+counter);
    T_5(:,ii) = CellS5(:,ii+counter) + CellS5(:,(ii+1)+counter);
    %Skips to next group
    counter = counter + 1;
end
%Specific times - Halfway point
for ii = 1:LL
    T1(:,ii) = T_1(length(t)/2,ii);
    T5(:,ii) = T_5(length(t)/2,ii);
end
x_contour = k50';
y_contour = kmax';
[Y,X] = meshgrid(y_contour,x_contour);
%Collects data by k50 - each entry is the different kmax values for a given
%k50
%MPC1
Z1 = [T1(:,1:11);T1(:,12:22);T1(:,23:33);T1(:,34:44);...
    T1(:,45:55);T1(:,56:66);T1(:,67:77);T1(:,78:88);...
    T1(:,89:99);T1(:,100:110);T1(:,111:121)];
%MPC5
Z5 = [T5(:,1:11);T5(:,12:22);T5(:,23:33);T5(:,34:44);...
    T5(:,45:55);T5(:,56:66);T5(:,67:77);T5(:,78:88);...
    T5(:,89:99);T5(:,100:110);T5(:,111:121)];
fig=figure(1)
%Tile 1
subplot(1,2,1)
contourf(Y,X,Z1,'Linewidth',2) %G1D1
colormap(parula)
% cb = colorbar;
% cb.LineWidth = 1.5;
set(gca,'YScale','log')
xlabel('k_{max}'), ylabel('k_{50} [M]')
xticks(0.5:0.1:1)
get(gca,'fontname')  % shows you what you are using.
set(gca,'linewidth', 2,'fontsize',24,'fontname','Helectiva') % Sets the width of the axis lines, font size, font
set(gca,'TickDir','out')
set(gca,'YScale','log')
%Tile 2
subplot(1,2,2)
contourf(Y,X,Z5,'Linewidth',2)
% cb = colorbar;
% cb.LineWidth = 1.5;
set(gca,'YScale','log')
xlabel('k_{max}'), ylabel('k_{50} [M]')
xticks(0.5:0.1:1)
get(gca,'fontname')  % shows you what you are using.
set(gca,'linewidth', 2,'fontsize',24,'fontname','Helectiva') % Sets the width of the axis lines, font size, font
set(gca,'TickDir','out')
set(gca,'YScale','log')
h = axes(fig,'visible','off'); 
c = colorbar(h,'Position',[0.93 0.168 0.022 0.7]);  % attach colorbar to h
c.LineWidth = 1.5;
c.FontSize = 24;
c.Label.String = 'Cell Number';
bottom = min(min(min(Z1)),min(min(Z5)));
top  = max(max(max(Z1)),max(max(Z5)));
caxis(h,[bottom, top]); 

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Alistair McQueen 2022-2-23

untitled.png

To emphasise this, plotting with two separate colourbars gives the results attached.

This is more prominent when changing the cd value to 1e-5.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Simon Chan 2022-2-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1656605-colour-bar-isn-t-correct-for-subplot-contour#answer_902740

在 MATLAB Online 中打开

Below is only a workaround to use each colorbar for each subplot and finally hidden them.

%Preamble
clc
clear all
clear workspace
%Parameters - Fixed:
%All the different initial conditions
MPC1 = [10000 000];
MPC2 = [7500 2500];
MPC3 = [5000 5000];
MPC4 = [2500 7500];
MPC5 = [000 10000];
g1 = 1.0392; %Phase 1 transiton
g2 = 1.0392; %Phase 2 transition
x0 = 10000; %Initial cell count
K = 100000; %Maximal cell count
time = linspace(0,15,1000); %time of simulation
%Parameters that vary:
cd = 1e-3;
kmax = [1 0.95 0.9 0.85 0.8 0.75 0.7 0.65 0.6 0.55 0.5];
k50 = [1e-3 7.5e-4 2.5e-4 1e-4 5e-5 1e-5 5e-6 1e-6 7.5e-7 5e-7 1e-7];
%Initial condition 1
for ii = 1:length(kmax)
    for kk = 1:length(k50)
        M = (cd*kmax(ii))/(cd+k50(kk));
        ode = @(t,x) [2*g2*x(2)*(1 - ((x(1)+x(2))/K)) - g1*x(1)*(1-M);...
            -g2*x(2)*(1 - ((x(1)+x(2))/K)) + g1*x(1)*(1-M)];
        [t, x(:,:,ii,kk)] = ode45(ode, time, MPC1);
    end
end
%Initial condition 5
for ii = 1:length(kmax)
    for kk = 1:length(k50)
        M = (cd*kmax(ii))/(cd+k50(kk));
        ode = @(t,x) [2*g2*x(2)*(1 - ((x(1)+x(2))/K)) - g1*x(1)*(1-M);...
            -g2*x(2)*(1 - ((x(1)+x(2))/K)) + g1*x(1)*(1-M)];
        [t, x5(:,:,ii,kk)] = ode45(ode, time, MPC5);
    end
end
%Data Extraction:
%This informs us how long our data set is
CellS = x(:,:);
CellS5 = x5(:,:);
%Loop length
LL = length(CellS(1,:))/2;
%Adding the first and second phase numbers together
counter = 0; %Initialise counter for initial grouping
for ii = 1:LL
    T_1(:,ii) = CellS(:,ii+counter) + CellS(:,(ii+1)+counter);
    T_5(:,ii) = CellS5(:,ii+counter) + CellS5(:,(ii+1)+counter);
    %Skips to next group
    counter = counter + 1;
end
%Specific times - Halfway point
for ii = 1:LL
    T1(:,ii) = T_1(length(t)/2,ii);
    T5(:,ii) = T_5(length(t)/2,ii);
end
x_contour = k50';
y_contour = kmax';
[Y,X] = meshgrid(y_contour,x_contour);
%Collects data by k50 - each entry is the different kmax values for a given
%k50
%MPC1
Z1 = [T1(:,1:11);T1(:,12:22);T1(:,23:33);T1(:,34:44);...
    T1(:,45:55);T1(:,56:66);T1(:,67:77);T1(:,78:88);...
    T1(:,89:99);T1(:,100:110);T1(:,111:121)];
%MPC5
Z5 = [T5(:,1:11);T5(:,12:22);T5(:,23:33);T5(:,34:44);...
    T5(:,45:55);T5(:,56:66);T5(:,67:77);T5(:,78:88);...
    T5(:,89:99);T5(:,100:110);T5(:,111:121)];
fig=figure(1);
%Tile 1
subplot(1,2,1)
contourf(Y,X,Z1,'Linewidth',2); %G1D1
colormap(parula)
cb1 = colorbar;
cb1.LineWidth = 1.5;
ax1=gca;
set(gca,'YScale','log')
xlabel('k_{max}'), ylabel('k_{50} [M]');
xticks(0.5:0.1:1);
get(gca,'fontname');  % shows you what you are using.
set(gca,'linewidth', 2,'fontsize',24,'fontname','Helectiva'); % Sets the width of the axis lines, font size, font
set(gca,'TickDir','out');
set(gca,'YScale','log');
%Tile 2
subplot(1,2,2)
contourf(Y,X,Z5,'Linewidth',2);
cb2 = colorbar;
cb2.LineWidth = 1.5;
ax2=gca;
set(gca,'YScale','log');
xlabel('k_{max}'), ylabel('k_{50} [M]');
xticks(0.5:0.1:1);
get(gca,'fontname')  % shows you what you are using.
set(gca,'linewidth', 2,'fontsize',24,'fontname','Helectiva'); % Sets the width of the axis lines, font size, font
set(gca,'TickDir','out')
set(gca,'YScale','log')
h = axes(fig,'visible','off'); 
bottom = min(min(min(Z1)),min(min(Z5)));
top  = max(max(max(Z1)),max(max(Z5)));
caxis(ax1,[bottom, top]);                   % Set the same limit to first subplot
caxis(ax2,[bottom, top]);                   % Set the same limit to second subplots
cb1.Position = [0 0 0 0];                   % Reduce the position for cb1
cb2.Position = [0 0 0 0];                   % Reduce the position for cb2
c = colorbar(h,'Position',[0.93 0.168 0.022 0.7]);  % attach colorbar to h
c.LineWidth = 1.5;
c.FontSize = 24;
c.Label.String = 'Cell Number';
caxis(cb1.Limits);                          % Set the new colorbar limits