Azzera filtri
Azzera filtri

How to get monthly average from daily data?

44 visualizzazioni (ultimi 30 giorni)
frankovaT
frankovaT il 27 Set 2018
Modificato: jonas il 27 Set 2018
So i have this Matrix, consisting of 4 columns: Column 1: Year Column 2: Month Column 3: day (1 to 30) Column 4: Some Value
I want to make monthly average of "some value" without breaking the matrix into some new matrix and calculate mean for each moths matrix seperately then rejoin them. I want to make an average matrix using for loop and find function for each myear seperately. So I dont like to have an average from all the first months of all the years together. Does anyone have an Idea?
I have sth like this, but i donno why it is not working:
k=1
%Calculating Monthly Average
for i=1961:1990 %Years Looping
for j=1:12 %Months Looping
idx=find((date(:,1)==i) & (date(:,2)==j)); %Find all indeces for the specific month
meandata1=mean(date(idx,4),1);
Monthly_Ave(k,:)=[i,j,meandata1;]; %Store in a new matrix the mean values per month
k=k+1;
end
end
  5 Commenti
frankovaT
frankovaT il 27 Set 2018
Dear jonas I truly appreciate your help i attached an Excel file of my data.
jonas
jonas il 27 Set 2018
and answer updated as promised!

Accedi per commentare.

Risposta accettata

jonas
jonas il 27 Set 2018
Modificato: jonas il 27 Set 2018
Not sure what you mean by "without breaking the matrix into some new matrix", but here we go:
x is column with months,
y is data
[G,name] = findgroups(x)
out = splitapply(@mean,y,G)
This will give you max 12 values, one for each month of the year. If you have several years, and you want 12 values per year, then you can do the following trick.
x1 is months (double)
x2 is years (double)
x=cellstr(num2str([x1, x2]))
[G,name] = findgroups(x)
out = splitapply(@mean,y,G)
Better yet, put your data in a timetable and use retime to calculate the monthly average. Here's an example using your data:
data=xlsread('data.xlsx')
t=datetime(data(:,1:3))
T=timetable(t,data(:,4))
T2=retime(T,'monthly','mean')
  6 Commenti
jonas
jonas il 27 Set 2018
Modificato: jonas il 27 Set 2018
There's a typo here:
Monthly_Ave(k,:)=[i,j,meandate1;];
should be meandata1. Also date is a built in function, so don't name your variable that! Both methods give the same results, so that's reassuring :)
frankovaT
frankovaT il 27 Set 2018
Oh god! you are an angle! I was so unsure about my code that i didnt even tried to look for typos:))) My first for loop!!! thanks :)

Accedi per commentare.

Più risposte (1)

KSSV
KSSV il 27 Set 2018
Let M be your month column, extracted and A be your data.
B = zeros(12,1) ; % means
for i = 1:12
B(i) = mean(A(M==i)) ;
end
  1 Commento
frankovaT
frankovaT il 27 Set 2018
maybe I should have mentioned that i have 50 years of data, so if I ll have like 50, Marches, and this is going to make an average from all the Marches for instance. but i want to make an average form each daily vaule in month of a year seperately.

Accedi per commentare.

Community Treasure Hunt

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

Start Hunting!

Translated by