Corrected typo in plotting of figure 2 (mh/k/min->mg/kg/min)

plotFigures.m

@@ -84,7 +84,7 @@ plotArea(colorMA,simM1,sim_time);
 
 yticks([0:1:5])
 xlabel('time (min)')
-ylabel('mh/kg/min')
+ylabel('mg/kg/min')
 set(gca,'FontSize',20);
 sgtitle('Glucose uptake in muscle and adipose tissue combined','FontSize',20)
 legend('M0','M1','data')
@@ -120,7 +120,7 @@ plotErrorbar(colorM,DATA.time,DATA.Muscle,DATA.SEM1); hold on
 plotErrorbar(colorA,DATA.time,DATA.Adipose,DATA.SEM2);
 
 xlabel('time (min)')
-ylabel('mh/kg/min')
+ylabel('mg/kg/min')
 legend([m,a],'muscle tissue','adipose tissue')
 set(gca,'FontSize',20);
 xticks([0:200:400]);
@@ -165,13 +165,12 @@ sim_muscleB = sim.variablevalues(:,ismember(IQMvariables(model),'U_idm'))';
 [sampledParams] = sample_params(100,allParams);
 boundries=get_maxmin(sim_time, model, squeeze(sampledParams),paramsFixed);
 
+% adipose
 figure()
 set(gcf, 'Units', 'Normalized', 'OuterPosition', [0, 0.04, 0.4, 0.7]);
 figSets()
 sgtitle('Glucose uptake in adipose tissue','FontSize',20)
-yticks([0:1:3]);
 
-%M2a
 subplot(1,2,1);
 plotSim(colorA,sim_adiposeA,sim_time); hold on
 plotErrorbar(colorA,DATA.time,DATA.Adipose,DATA.SEM2); 
@@ -180,7 +179,6 @@ yticks([0:1:3])
 set(gca,'FontSize',15);
 set(gca,'linewidth',1.4);
 
-% M2b
 subplot(1,2,2);
 plotUncertainty(colorA,boundries,sim_time,'a'); hold on
 plotSim(colorA,sim_adiposeB,sim_time); 
@@ -196,13 +194,13 @@ figure()
 set(gcf, 'Units', 'Normalized', 'OuterPosition', [0, 0.04, 0.4, 0.7]);
 figSets()
 sgtitle('Glucose uptake in muscle tissue','FontSize',20)
-yticks([0:1:3]);
 
 subplot(1,2,1);
 plotSim(colorM,sim_muscleA,sim_time); hold on
 plotErrorbar(colorM,DATA.time,DATA.Muscle,DATA.SEM1); 
 title('M2a','FontSize',12);
 yticks([0:1:3])
+ylim([0 3.5])
 set(gca,'FontSize',15);
 set(gca,'linewidth',1.4);
 
@@ -212,9 +210,10 @@ plotErrorbar(colorM,DATA.time,DATA.Muscle,DATA.SEM1);
 plotSim(colorM,sim_muscleB,sim_time); 
 title('M2b','FontSize',12);
 yticks([0:1:3])
+ylim([0 3.5])
 set(gca,'FontSize',15);
-legend('uncertainty','data','simulation','Location', 'northeast')
 set(gca,'linewidth',1.4)
+legend('uncertainty','data','simulation','Location', 'northeast')
  
 
 %% fig5
@@ -230,8 +229,8 @@ iozzo.clearance=100*iozzo.clearance/iozzo.clearance(1);
 [~,paramsAll]=IQMparameters(model);
 [paramsBI,paramsI,paramsB,params]=deal(paramsAll);
 
-params(90) = 0; % removing insulin infusion, parameter INS_inf
-paramsB(90) = 0; % removing insulin infusion, parameter INS_inf
+params(90) = 0; 
+paramsB(90) = 0; 
 paramsBI(99) = 30*paramsBI(99); % adding bradykinin, increases blood foow two-folds
 
 [initsI] = simulateSteadyState(model, paramsI);