Trying to use retime to sum daily values, but I don't want days that only have NaN values to become 0

20 次查看（过去 30 天）
Rita Leavitt
Rita Leavitt 2021-7-9
评论： Peter Perkins 2022-3-21
Basically I have a timetable or rainfall values (a column with half-hourly time increments and a column with the rainfall at each time point). I used retime to get daily sums:
sum_rainfall_2018 = retime(rainfall_2018, 'daily','sum');
This mostly works, but I have a number of days that only have NaN values, and retime is summing those days as 0. I need them to stay NaN because I need to know where my gaps are. Is there a simple way to do this?

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

回答（3 个）

Jonas
Jonas 2021-7-10
weird, the normal behavior of sum is giving NaN if at least one of the summed values is NaN. try @sum instead of 'sum' and if this also does not work, use a self defined function which forces the exact output
Peter Perkins
Peter Perkins 2021-7-27
In the doc for retime (admittedly, not really in flashing bold letters): "All the listed methods omit NaNs, NaTs, and other missing data indicators, except for func. To include missing data indicators, specify func as a function handle to a function that includes them when aggregating data."
Haven't tried it, but I believe that if you pass in @sum, not 'sum', you will get what you want.
Everett Snieder
Everett Snieder 2022-3-17
OK, I figured it out... basically when you pass an empty array into sum(), it returns 0. So to handle big time gaps (larger than your timestep), you need to create a special sum() function and pass the handle to retime:
function y = sum2(x)
if isempty(x)
y = nan();
else
y = sum(x);
end
  1 个评论
Peter Perkins
Peter Perkins 2022-3-21
Ran in:
Right, prod and sum return the "unitary element" when given an empty:
prod([])
ans = 1
sum([])
ans = 0
But I don't think you need to do this:
tt = timetable([1;2;NaN;4;NaN;NaN],'RowTimes',datetime(2022,3,[21;21;22;22;23;23]))
tt = 6×1 timetable
Time Var1 ___________ ____ 21-Mar-2022 1 21-Mar-2022 2 22-Mar-2022 NaN 22-Mar-2022 4 23-Mar-2022 NaN 23-Mar-2022 NaN
retime(tt,'daily','sum') % 'sum' removes NaNs
ans = 3×1 timetable
Time Var1 ___________ ____ 21-Mar-2022 3 22-Mar-2022 4 23-Mar-2022 0
retime(tt,'daily',@sum) % @sum does not
ans = 3×1 timetable
Time Var1 ___________ ____ 21-Mar-2022 3 22-Mar-2022 NaN 23-Mar-2022 NaN
retime(tt,'daily',@(x)sum(x,'omitnan')) % equivalent to 'sum'
ans = 3×1 timetable
Time Var1 ___________ ____ 21-Mar-2022 3 22-Mar-2022 4 23-Mar-2022 0
It may be that your desired behavior is, "remove NaNs before summing, unless all NaN, in which case return NaN." For that, you would need your own function, because the standard behavior of sum on empty is to return 0.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Descriptive Statistics 的更多信息

标签

产品


版本

R2019b

Community Treasure Hunt

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

Start Hunting!

Translated by