test_timedelta.py 6.8 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206
  1. from datetime import timedelta
  2. import numpy as np
  3. import pytest
  4. import pandas as pd
  5. from pandas import (
  6. DataFrame,
  7. Series,
  8. )
  9. import pandas._testing as tm
  10. from pandas.core.indexes.timedeltas import timedelta_range
  11. def test_asfreq_bug():
  12. df = DataFrame(data=[1, 3], index=[timedelta(), timedelta(minutes=3)])
  13. result = df.resample("1T").asfreq()
  14. expected = DataFrame(
  15. data=[1, np.nan, np.nan, 3],
  16. index=timedelta_range("0 day", periods=4, freq="1T"),
  17. )
  18. tm.assert_frame_equal(result, expected)
  19. def test_resample_with_nat():
  20. # GH 13223
  21. index = pd.to_timedelta(["0s", pd.NaT, "2s"])
  22. result = DataFrame({"value": [2, 3, 5]}, index).resample("1s").mean()
  23. expected = DataFrame(
  24. {"value": [2.5, np.nan, 5.0]},
  25. index=timedelta_range("0 day", periods=3, freq="1S"),
  26. )
  27. tm.assert_frame_equal(result, expected)
  28. def test_resample_as_freq_with_subperiod():
  29. # GH 13022
  30. index = timedelta_range("00:00:00", "00:10:00", freq="5T")
  31. df = DataFrame(data={"value": [1, 5, 10]}, index=index)
  32. result = df.resample("2T").asfreq()
  33. expected_data = {"value": [1, np.nan, np.nan, np.nan, np.nan, 10]}
  34. expected = DataFrame(
  35. data=expected_data, index=timedelta_range("00:00:00", "00:10:00", freq="2T")
  36. )
  37. tm.assert_frame_equal(result, expected)
  38. def test_resample_with_timedeltas():
  39. expected = DataFrame({"A": np.arange(1480)})
  40. expected = expected.groupby(expected.index // 30).sum()
  41. expected.index = timedelta_range("0 days", freq="30T", periods=50)
  42. df = DataFrame(
  43. {"A": np.arange(1480)}, index=pd.to_timedelta(np.arange(1480), unit="T")
  44. )
  45. result = df.resample("30T").sum()
  46. tm.assert_frame_equal(result, expected)
  47. s = df["A"]
  48. result = s.resample("30T").sum()
  49. tm.assert_series_equal(result, expected["A"])
  50. def test_resample_single_period_timedelta():
  51. s = Series(list(range(5)), index=timedelta_range("1 day", freq="s", periods=5))
  52. result = s.resample("2s").sum()
  53. expected = Series([1, 5, 4], index=timedelta_range("1 day", freq="2s", periods=3))
  54. tm.assert_series_equal(result, expected)
  55. def test_resample_timedelta_idempotency():
  56. # GH 12072
  57. index = timedelta_range("0", periods=9, freq="10L")
  58. series = Series(range(9), index=index)
  59. result = series.resample("10L").mean()
  60. expected = series.astype(float)
  61. tm.assert_series_equal(result, expected)
  62. def test_resample_offset_with_timedeltaindex():
  63. # GH 10530 & 31809
  64. rng = timedelta_range(start="0s", periods=25, freq="s")
  65. ts = Series(np.random.randn(len(rng)), index=rng)
  66. with_base = ts.resample("2s", offset="5s").mean()
  67. without_base = ts.resample("2s").mean()
  68. exp_without_base = timedelta_range(start="0s", end="25s", freq="2s")
  69. exp_with_base = timedelta_range(start="5s", end="29s", freq="2s")
  70. tm.assert_index_equal(without_base.index, exp_without_base)
  71. tm.assert_index_equal(with_base.index, exp_with_base)
  72. def test_resample_categorical_data_with_timedeltaindex():
  73. # GH #12169
  74. df = DataFrame({"Group_obj": "A"}, index=pd.to_timedelta(list(range(20)), unit="s"))
  75. df["Group"] = df["Group_obj"].astype("category")
  76. result = df.resample("10s").agg(lambda x: (x.value_counts().index[0]))
  77. expected = DataFrame(
  78. {"Group_obj": ["A", "A"], "Group": ["A", "A"]},
  79. index=pd.TimedeltaIndex([0, 10], unit="s", freq="10s"),
  80. )
  81. expected = expected.reindex(["Group_obj", "Group"], axis=1)
  82. expected["Group"] = expected["Group_obj"]
  83. tm.assert_frame_equal(result, expected)
  84. def test_resample_timedelta_values():
  85. # GH 13119
  86. # check that timedelta dtype is preserved when NaT values are
  87. # introduced by the resampling
  88. times = timedelta_range("1 day", "6 day", freq="4D")
  89. df = DataFrame({"time": times}, index=times)
  90. times2 = timedelta_range("1 day", "6 day", freq="2D")
  91. exp = Series(times2, index=times2, name="time")
  92. exp.iloc[1] = pd.NaT
  93. res = df.resample("2D").first()["time"]
  94. tm.assert_series_equal(res, exp)
  95. res = df["time"].resample("2D").first()
  96. tm.assert_series_equal(res, exp)
  97. @pytest.mark.parametrize(
  98. "start, end, freq, resample_freq",
  99. [
  100. ("8H", "21h59min50s", "10S", "3H"), # GH 30353 example
  101. ("3H", "22H", "1H", "5H"),
  102. ("527D", "5006D", "3D", "10D"),
  103. ("1D", "10D", "1D", "2D"), # GH 13022 example
  104. # tests that worked before GH 33498:
  105. ("8H", "21h59min50s", "10S", "2H"),
  106. ("0H", "21h59min50s", "10S", "3H"),
  107. ("10D", "85D", "D", "2D"),
  108. ],
  109. )
  110. def test_resample_timedelta_edge_case(start, end, freq, resample_freq):
  111. # GH 33498
  112. # check that the timedelta bins does not contains an extra bin
  113. idx = timedelta_range(start=start, end=end, freq=freq)
  114. s = Series(np.arange(len(idx)), index=idx)
  115. result = s.resample(resample_freq).min()
  116. expected_index = timedelta_range(freq=resample_freq, start=start, end=end)
  117. tm.assert_index_equal(result.index, expected_index)
  118. assert result.index.freq == expected_index.freq
  119. assert not np.isnan(result[-1])
  120. @pytest.mark.parametrize("duplicates", [True, False])
  121. def test_resample_with_timedelta_yields_no_empty_groups(duplicates):
  122. # GH 10603
  123. df = DataFrame(
  124. np.random.normal(size=(10000, 4)),
  125. index=timedelta_range(start="0s", periods=10000, freq="3906250n"),
  126. )
  127. if duplicates:
  128. # case with non-unique columns
  129. df.columns = ["A", "B", "A", "C"]
  130. result = df.loc["1s":, :].resample("3s").apply(lambda x: len(x))
  131. expected = DataFrame(
  132. [[768] * 4] * 12 + [[528] * 4],
  133. index=timedelta_range(start="1s", periods=13, freq="3s"),
  134. )
  135. expected.columns = df.columns
  136. tm.assert_frame_equal(result, expected)
  137. @pytest.mark.parametrize("unit", ["s", "ms", "us", "ns"])
  138. def test_resample_quantile_timedelta(unit):
  139. # GH: 29485
  140. dtype = np.dtype(f"m8[{unit}]")
  141. df = DataFrame(
  142. {"value": pd.to_timedelta(np.arange(4), unit="s").astype(dtype)},
  143. index=pd.date_range("20200101", periods=4, tz="UTC"),
  144. )
  145. result = df.resample("2D").quantile(0.99)
  146. expected = DataFrame(
  147. {
  148. "value": [
  149. pd.Timedelta("0 days 00:00:00.990000"),
  150. pd.Timedelta("0 days 00:00:02.990000"),
  151. ]
  152. },
  153. index=pd.date_range("20200101", periods=2, tz="UTC", freq="2D"),
  154. ).astype(dtype)
  155. tm.assert_frame_equal(result, expected)
  156. def test_resample_closed_right():
  157. # GH#45414
  158. idx = pd.Index([pd.Timedelta(seconds=120 + i * 30) for i in range(10)])
  159. ser = Series(range(10), index=idx)
  160. result = ser.resample("T", closed="right", label="right").sum()
  161. expected = Series(
  162. [0, 3, 7, 11, 15, 9],
  163. index=pd.TimedeltaIndex(
  164. [pd.Timedelta(seconds=120 + i * 60) for i in range(6)], freq="T"
  165. ),
  166. )
  167. tm.assert_series_equal(result, expected)