Using 'plotResiduals' in a for loop.. is it possible to using indexing to store the results and plot then together?

I have a series of data fits I've performed using 'fitnlm' which I've acquired using a for loop like so.

    fittables = cell(1, numel(files));
    modelfun = @(sigma,x) exp(-x.^2 / (2 * sigma));

for k = 1:length(files)
    table{k} = table(x{k}, y{k}, 'VariableNames', {'x', 'y'});
    beta0 = 1; 
    mdl{k} = fitnlm(table{k}, modelfun, beta0);
end

I want to use the 'plotResiduals()' function to plot the residuals of each of these fits onto one graph. I can get it to work on a single fit but I can't work out how to index the function. I've tried the following two ways..

residuals{k} = plotResiduals(mdl{k},'caseorder');

which will run but it just overwrites the plot with each iteration,

and

residuals(k) = plotResiduals(mdl(k),'caseorder');

which results in the error "Undefined function 'plotResiduals' for input arguments of type 'cell' ".

Is there any way to plot these residuals on the same plot within the for loop? Indexing doesn't seem to be working.

[Jakob B. Nielsen]

Do you have

hold on

active? It retains plots in the current axes, whereas if hold is off, plots are overwritten.

[mel1708]

Here is the whole plotting section of the code which is inside the for loop. I did include hold on and hold off as you can see. Should it matter that the other plots are in between hold on and hold off? Or should I be starting a new hold on and hold off section for plotting the residuals?

     figure(11);
     hold on
     plot(dose{k},sf{k},'Color',col(k,:),'Marker','x','LineStyle','none','DisplayName',txt);
     h = errorbar(dose{k},sf{k},errs{k},'Marker','none','LineStyle','none','Color',[col(k,1), col(k,2), col(k,3)],'HandleVisibility','off');
     legend('Location','southoutside');
     xlabel('Dose (Gy)'); ylabel('Surviving Fraction'); grid on; legend show;
     title('Rayleigh Distribution Fit - Nonlinear Regression Model');
     plot(dose{k},yfitRayl{k},'Color',col(k,:),'HandleVisibility','off');
     set(gcf,'color','w');
     
     figure(12);
     resid(k) = plotResiduals(mdlRayl(k),'caseorder');
     hold off

Answer 1

[Chidvi Modala]

It appears that plotResidual respects "hold on" for other than probability plots. You can create the probability plot directly using normplot, and not using probplot or plotResiduals. For example:

load carsmall
tbl = table(MPG,Weight);
tbl.Year = categorical(Model_Year);
mdl = fitlm(tbl,'MPG ~ Year + Weight^2');
normplot(mdl.Residuals.Raw); 
hold on; 
h1= plotResiduals(mdl)