How to modify this code for calculation seasonal values for each year (not all years)

Question

BN 2020-4-2

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/514962-how-to-modify-this-code-for-calculation-seasonal-values-for-each-year-not-all-years

评论： Ameer Hamza 2020-4-4

CELL.mat

I have a 1x3 cell contains three 360x3 tables. For every table, I want to calculate the average of rrr24 for each season and write them on the 5th column. For example:

I have a code as you can see below, but it calculates the average of each season for all 30 years (not yearly), and didn't write averages on the 5th column too.

I would like just a new 1x3 cell with tables that have the 5th column as I described above.

meanTableDataArray = cell(length(CELL),4);
for j = 1:length(CELL)
    sampleTableData = CELL{j};
    sampleTableData.month = month(sampleTableData.dates);
    sampleTableData.seasons = floor(sampleTableData.month ./3);
    sampleTableData.seasons(sampleTableData.seasons ==4 ) = 0;
    sampleTableData.seasons = categorical(sampleTableData.seasons, [0 1 2 3], ["Spr", "Sum", "Aut", "Win"]);
    [group, mean_table] = findgroups(sampleTableData(:, 'seasons'));
    mean_table.rrr24 = splitapply(@mean, sampleTableData.rrr24, group);
    for k = 1:4
        meanTableDataArray{j,k} = mean_table(k,:); % extracts as a table
    end
end

Thank you.

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

darova 2020-4-3

I like your drawing

请先登录，再进行评论。

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

Answer 1

Ameer Hamza 2020-4-3

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/514962-how-to-modify-this-code-for-calculation-seasonal-values-for-each-year-not-all-years#answer_423743

在 MATLAB Online 中打开

Try this

data = load('CELL.mat');
CELL = data.CELL;
month_names = month(datetime(1,1,1):calmonths(1):datetime(1,12,1), 'name');
seasons = ["Spr", "Sum", "Aut", "Win"];
for i=1:numel(CELL)
    CELL{i}.grid_name = [];
    CELL{i}.month = month_names(month(CELL{i}.dates))';
    CELL{i}.seasons = seasons(floor((month(CELL{i}.dates)+2)/3))';
    grps = findgroups(findgroups(year(CELL{i}.dates), CELL{i}.seasons));
    avg_value = accumarray(grps, CELL{i}.rrr24, [], @mean);
    CELL{i}.Average = avg_value(grps);
end

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

BN 2020-4-3

在 MATLAB Online 中打开

Dear Ameer Hamza,

Thank you.

Just a little problem. I saw that the code classify January and February (and other months) in the wrong season, for example January in the Spring. I know that this is my problem because in the first code in my question I typed Spring first, so I changed your code to this:

month_names = month(datetime(1,1,1):calmonths(1):datetime(1,12,1), 'name');
seasons = ["Win", "Spr", "Sum", "Aut"];

But it is not accomplished and this time for example code considers March in the winter.

I thought about this code today about 8 hours; thought I could solve this problem if I add an extra row above all rows with NaN value for rrr24 and 12/01/1988 for dates column but even with this method, It is not fixed.

Furthermore, since it is not possible to have a table with the different number of rows for each column I want to ask you if you please guide me how I can store the output of the code in another cell in order to have only one number for each season:

Thank you so much again.