How can I use loop to produce multiple graphs?

Hello guys! I have a simple problem but I cannot seem to be able to solve it.

I want to produce 11 graphs within a loop, based on the distinction of my Daily Data according to year (11 years 2006-2016, 11 graphs). Here is what I have done so far, which it doesn't work though:

DTindex  = datenum(Daily_Data.year);
lim   = min(DTindex);
Hindd = DTindex - lim +1;
YEAR = accumarray(Hindd, datenum(Daily_Data.year), [], @mean );
Deaths_allyears = accumarray(Hindd, Daily_Data.Daily_Deaths,   [], @sum );

for YEAR=1:11
DTindex  = datenum(Daily_Data.month); %Daily_Data.month should focus on each year separately within every loop
lim   = min(DTindex);
Hindd = DTindex - lim +1;
MONTH = accumarray(Hindd, datenum(Daily_Data.month), [], @mean );%MONTH should focus on each year separately within every loop
Deaths_permonth = accumarray(Hindd, Daily_Data.Daily_Deaths,   [], @sum ); %Deaths_permonth should focus on each year separately within every loop
m = (1:12).';
Monthly_Temp = accumarray(Hindd, Daily_Data.Daily_T,   [], @nanmean );%Monthly_Temp should focus on each year separately within every loop

yyaxis left
plot(m,Deaths_permonth)
xlabel('Month'); xlim([1 12]); xticks([1 2 3 4 5 6 7 8 9 10 11 12]);
xticklabels({'Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov','Dec'});
ylabel('Monthly Deaths')
yyaxis right
plot(m,Monthly_Temp)
ylabel('Temperature (oC)'); 
grid on;
title({'City-specific distribution of deaths by month and temperature', YEAR});
end
    

Hopefully I want to produce 1 graph for each year (like the attached, which is accummulative for all years).

I would really appreciate your help.

PS. I'm on R2019a.

Answer 1

[Mathieu NOE]

hello Daphne

I modified a bit your code - see the comments

only issue for me in line :

Monthly_Temp = accumarray(Hindd_ck, Daily_Data_this_year.Daily_T,   [], @nanmean );

i could not make this work with "nanmean" so I changed to regular "mean"

hope it helps

Daily_Data = readtable('Daily_Data_Thessaloniki.xlsx');
DTindex  = datenum(Daily_Data.year);
lim   = min(DTindex);
Hindd = DTindex - lim +1;
YEAR = accumarray(Hindd, datenum(Daily_Data.year), [], @mean );
Deaths_allyears = accumarray(Hindd, Daily_Data.Daily_Deaths,   [], @sum );

for ck=1:numel(YEAR)
    ind = find(Hindd == ck); % this will give which rows to select in Daily_Data table for the "ck" year
    Daily_Data_this_year = Daily_Data(ind,:); % extracted table for one (ck) single year
    DTindex  = datenum(Daily_Data_this_year.month); %Daily_Data.month should focus on each year separately within every loop
    lim_ck   = min(DTindex); % changed name to make sure not to overwritte a variable used before the for loop
    Hindd_ck = DTindex - lim_ck +1; % changed name to make sure not to overwritte a variable used before the for loop
    MONTH = accumarray(Hindd_ck, datenum(Daily_Data_this_year.month), [], @mean );%MONTH should focus on each year separately within every loop
    Deaths_permonth = accumarray(Hindd_ck, Daily_Data_this_year.Daily_Deaths,   [], @sum ); %Deaths_permonth should focus on each year separately within every loop
    m = (1:12).';
%     Monthly_Temp = accumarray(Hindd_ck, Daily_Data_this_year.Daily_T,   [], @nanmean );%Monthly_Temp should focus on each year separately within every loop
    Monthly_Temp = accumarray(Hindd_ck, Daily_Data_this_year.Daily_T,   [], @mean );%Monthly_Temp should focus on each year separately within every loop
    yyaxis left
    figure(ck),plot(m,Deaths_permonth)
    xlabel('Month'); xlim([1 12]); xticks([1 2 3 4 5 6 7 8 9 10 11 12]);
    xticklabels({'Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov','Dec'});
    ylabel('Monthly Deaths')
    yyaxis right
    plot(m,Monthly_Temp)
    ylabel('Temperature (oC)'); 
    grid on;
    title({'City-specific distribution of deaths by month and temperature', YEAR(ck)});
end

[Daphne PARLIARI]

It works like a charm, thank you so much!

What should I change, if I want Monthly_Temp to accumulate not only each year's temperatures, but the temperatures from all timeseries (2006-2016)?

This way, the temperature curve will be the same for all 11 graphs.

[Mathieu NOE]

hello again

I had to do 2 for loops - maybe there is a better way to code it , I leave it to the experts...

so first loop will do the maths for each month , then I do the average for the temperature curve , then the plots .

Daily_Data = readtable('Daily_Data_Thessaloniki.xlsx');
DTindex  = datenum(Daily_Data.year);
lim   = min(DTindex);
Hindd = DTindex - lim +1;
YEAR = accumarray(Hindd, datenum(Daily_Data.year), [], @mean );
Deaths_allyears = accumarray(Hindd, Daily_Data.Daily_Deaths,   [], @sum );

for ck=1:numel(YEAR)
    ind = find(Hindd == ck); % this will give which rows to select in Daily_Data table for the "ck" year
    Daily_Data_this_year = Daily_Data(ind,:); % extracted table for one (ck) single year
    DTindex  = datenum(Daily_Data_this_year.month); %Daily_Data.month should focus on each year separately within every loop
    lim_ck   = min(DTindex); % changed name to make sure not to overwritte a variable used before the for loop
    Hindd_ck = DTindex - lim_ck +1; % changed name to make sure not to overwritte a variable used before the for loop
    MONTH = accumarray(Hindd_ck, datenum(Daily_Data_this_year.month), [], @mean );%MONTH should focus on each year separately within every loop
    Deaths_permonth(:,ck) = accumarray(Hindd_ck, Daily_Data_this_year.Daily_Deaths,   [], @sum ); %Deaths_permonth should focus on each year separately within every loop
%     Monthly_Temp = accumarray(Hindd_ck, Daily_Data_this_year.Daily_T,   [], @nanmean );%Monthly_Temp should focus on each year separately within every loop
    Monthly_Temp(:,ck) = accumarray(Hindd_ck, Daily_Data_this_year.Daily_T,   [], @mean );%Monthly_Temp should focus on each year separately within every loop
end

% averaged temperature accross all years
Monthly_Temp_all_avg = mean(Monthly_Temp,2,'omitnan');
    m = (1:12).';

for ck=1:numel(YEAR)
%     Monthly_Temp = accumarray(Hindd_ck, Daily_Data_this_year.Daily_T,   [], @nanmean );%Monthly_Temp should focus on each year separately within every loop
    yyaxis left
    figure(ck),plot(m,Deaths_permonth(:,ck))
    xlabel('Month'); xlim([1 12]); xticks([1 2 3 4 5 6 7 8 9 10 11 12]);
    xticklabels({'Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov','Dec'});
    ylabel('Monthly Deaths')
    yyaxis right
    plot(m,Monthly_Temp_all_avg)
    ylabel('Temperature (oC)'); 
    grid on;
    title({'City-specific distribution of deaths by month and temperature', YEAR(ck)});
end

[Daphne PARLIARI]

It works better than fine!

Mathieu I am grateful for your help.

[Mathieu NOE]

My pleasure !

Answer 2

[Bob Thompson]

If you want individual graphs, add the following inside the loop, before the yyaxis command:

...
figure(i)
yyaxis left
...

If you want all of the plots on the same graph, put a 'hold on' command either inside or before the loop that generates the plots.

If you want all of the plots in the same figure, but different subplots, then use the subplots command, and indexing to select the appropriate pane.