How can I use loop to produce multiple graphs?

Question

Daphne PARLIARI 2021-2-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/755884-how-can-i-use-loop-to-produce-multiple-graphs

评论： Mathieu NOE 2021-2-26

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.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Mathieu NOE 2021-2-25

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/755884-how-can-i-use-loop-to-produce-multiple-graphs#answer_633184

在 MATLAB Online 中打开

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

4 个评论
显示 2更早的评论隐藏 2更早的评论

Mathieu NOE 2021-2-26

在 MATLAB Online 中打开

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