MATLAB Answers

0

Timetable Monthly Average over Many Years

Asked by Nathan Gjersvik on 2 Oct 2019
Latest activity Commented on by Nathan Gjersvik on 3 Oct 2019
I have a daily precipitation data file with two columns.The first column is "date", with date format yyyy-mm-dd and the second column is "precip" with double format. My data ranges from 1980-01-01 to 2010-12-31, and I would like to determine the average precipitation values for each month over that interval.
Example of desired result:
Month avgPrecip
Jan 0.0056
Feb 0.0034
...
for each month.
I have attached my csv file with the data and my current progess. I can get the monthly average of the data within each year using a timetable and retime, but I am stuck trying to average the monthly average data for the entire data range. Thank you!
precip = readtable('CR_precip.csv');
precip = table2timetable(precip);
monthlyAvg = retime(precip,'monthly','mean');

  0 Comments

Sign in to comment.

Products


Release

R2018a

3 Answers

Answer by Guillaume
on 2 Oct 2019
 Accepted Answer

The easiest:
precip = table2timetable(readtable('CR_precip.csv'));
monthlyavg = groupsummary(precip, 'day', 'monthofyear', 'mean')

  1 Comment

This works and is very simple! Thank you very much!

Sign in to comment.


Answer by QuanCCC on 2 Oct 2019

you can simply make another column of year-month-day-time, merge it to your data. Then transform it to timetable, time column will be your new column. Then use the same 'monthly','mean' on all months.
t1 = datetime(1980,01,01,12,0,0);
t2 = datetime(2010,12,31,12,0,0);
t = t1:t2;
t = t';
clear t1 t2
NewData = timetable(t, YourOldData); % merge data
MonthlyAve = retime(NewData, 'monthly', 'mean');

  0 Comments

Sign in to comment.


Answer by Akira Agata
on 3 Oct 2019

Looking at your csv data, some additional options will be needed.
(1) To specify the delimiter in your csv data, 'Delimiter' option should be added in readtabe function
(2) 'Format' option should be added in datetime function
The following is an example:
T = readtable('CR_precip.csv','Delimiter',',');
T.day = datetime(T.day,'Format','MM/dd/yyyy');
TT = table2timetable(T);
TT = retime(TT,'monthly','mean');

  0 Comments

Sign in to comment.