I have:
Temperature profiles array for 5 days. such that, temp(heigth,day). however, there are multiple profile in some days which makes the number of profile equal 11.
size of temp = (3650 by 11)
size of day = (1 by 11)
day=[1,1,1,2,2,3,3,3,4,4,5].
I want to:
average temperature for every similar days
such that:
dailyday=[1,2,3,4,5]
temp=(3650 by 5)

 Réponse acceptée

madhan ravi
madhan ravi le 17 Mai 2019

1 vote

splitapply(@(x) mean(x,2),Temp,findgroups(Day))

5 commentaires

madhan ravi
madhan ravi le 17 Mai 2019
Modifié(e) : madhan ravi le 17 Mai 2019
Prior to 2015b:
dailyday=unique(Day)
M=mat2cell(Temp,size(Temp,1),histc(Day,dailyday));
temp=cell2mat(cellfun(@(x) mean(x,2), M, 'un',0))
Kafayat Olayinka
Kafayat Olayinka le 17 Mai 2019
I recieved this message while running the second line.
Error using mat2cell (line 89)
Input arguments, D1 through D2, must sum to each dimension of the input matrix size
Kafayat Olayinka
Kafayat Olayinka le 17 Mai 2019
Recieved this error while compiling the first code with splitapply*
Error using splitapply (line 132)
The function '@(x)mean(x,2)' returned a non-scalar value when applied to the 5th group of data.
To compute nonscalar values for each group, create an anonymous function to return each value in a scalar
cell:
@(x){mean(x,2)}
madhan ravi
madhan ravi le 17 Mai 2019
Modifié(e) : madhan ravi le 17 Mai 2019
C=splitapply(@(x) {mean(x,2)},Temp,findgroups(Day));
[C{:}]
%or
[~,~,c]=unique(Day);
C = accumarray(c,1:size(Temp,2),[],@(x){mean(Temp(:,x),2)});
[C{:}]
Kafayat Olayinka
Kafayat Olayinka le 17 Mai 2019
Worked perfectly. Thank you.

Connectez-vous pour commenter.

Plus de réponses (1)

Walter Roberson
Walter Roberson le 17 Mai 2019

0 votes

[mean(temp(:,1:3),2), mean(temp(:,4:5)), mean(temp(:,6:8)), mean(temp(:,9:10)), temp(:,11)]

1 commentaire

Kafayat Olayinka
Kafayat Olayinka le 17 Mai 2019
Thanks. How should I write the loop assuming Days are more than 5 days and each day may have multiple temperature profiles. How should I average each temp profile per day?

Connectez-vous pour commenter.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by