Daily average precip for many years

5 次查看（过去 30 天）
Thomas Burbey
Thomas Burbey 2018-5-17
评论： dpb 2018-5-18
I have a large matrix of two columns containing dates (eg 1/1/1991) and precip values. I want to create the average daily precip for all years from 1991-2016 such that the day of the year is in column 1 and the ave precip value is in column 2. I'm getting bogged down I believe because of leap years (not sure how to eliminate that date, which is fine). Any help would be greatly appreciated.
  2 个评论
the cyclist
the cyclist 2018-5-17
Can you (ideally) upload your data in a MAT file, or (less ideally) give an example -- not just a description -- of how your input data are stored?
Thomas Burbey
Thomas Burbey 2018-5-18
Here is the precip file you requested showing you the format.

请先登录，再进行评论。

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

采纳的回答

dpb
dpb 2018-5-18
编辑：dpb 2018-5-18
Just taking a stab presuming the file is something recognizable--
precip=table2timetable(readtable('yourfile.ext')); % read the file; convert to timetable
daily=retime(precip,'daily','mean');
Done!
"Here is the precip file..."
Oh...only daily data and already in table--OK, then create the grouping variable
>> precip.DOM=day(precip.DATE); % day of month each entry
>>daily=varfun(@mean,precip,'groupingvariables','DOM','inputvariables','Precipitation');
...
ADDENDUM OK, answer the question asked... :)
>> precip.DOY=day(precip.DATE,'dayofyear');
>> daily=varfun(@mean,precip,'groupingvariables','DOY','inputvariables','Precipitation');
>> whos daily
Name Size Bytes Class Attributes
daily 366x3 10149 table
>> daily(1:10,:)
ans =
10×3 table
DOY GroupCount mean_Precipitation
___ __________ __________________
1 26 3.27692307692308
2 26 3.8
3 26 2.5
4 26 2.84230769230769
5 26 3.94615384615385
6 26 5.21923076923077
7 26 2.13461538461539
8 26 2.53461538461538
9 26 3.65
10 26 3.33461538461538
>>
should work; I believe the 'dayofyear' option came about R2104 or so...there's just too much stuff to try to recall introduced too rapidly.
  11 个评论
显示 9更早的评论
Thomas Burbey
Thomas Burbey 2018-5-18
Everything works fine. I gave you the matrix after I had eliminated NaN, which is what the raw precip file has when no precip for a day is measured. Sorry for the confusion. Thanks again for your help.
dpb
dpb 2018-5-18
Ah, so! That 'splains things, indeed! :)
The time-related stuff has gotten pretty deep as well as facilities to use the various varieties of tables. Seems to me like there's an excess of types that are similar but not quite the same that makes for very difficult learning task; I wish TMW would slow down some and work more on consistency of interfaces and methods.

请先登录，再进行评论。

更多回答（1 个）

Razvan Carbunescu
Razvan Carbunescu 2018-5-18
An alternative if using R2018a is groupsummary, it has a few different ways of selecting days for each week, month, year from a date:
>> GT = groupsummary(precip,'DATE','dayofyear','mean','Precipitation');
>> head(GT)
ans =
8×3 table
dayofyear_DATE GroupCount mean_Precipitation
______________ __________ __________________
1 26 3.2769
2 26 3.8
3 26 2.5
4 26 2.8423
5 26 3.9462
6 26 5.2192
7 26 2.1346
8 26 2.5346
  2 个评论
Thomas Burbey
Thomas Burbey 2018-5-18
That is cool! I wish I had that latest version. SO helpful
dpb
dpb 2018-5-18
There's only a one-liner difference to generate the grouping variable, though, so it's only a relatively minor syntactic sugar refinement.
Note same results as above which is another independent confirmation of correct result.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Dates and Time 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by