R 如何在数据集中找到超过6个值的最大增量
我有一个核心温度值的数据集。每10秒钟，我收集一个持续数小时的温度值。我试图在6个温度值范围内找到温度的最高升高（Tslope）。我尝试了滚动平均的方法，但我甚至没有接近解决方案。 数据集如下所示：R 如何在数据集中找到超过6个值的最大增量,r,moving-average,highest,R,Moving Average,Highest,我有一个核心温度值的数据集。每10秒钟，我收集一个持续数小时的温度值。我试图在6个温度值范围内找到温度的最高升高（Tslope）。我尝试了滚动平均的方法，但我甚至没有接近解决方案。 数据集如下所示： R001 <- R001_G2_Q3_T Test Subject Datum Zeit Temperatur Timepoint 1 G2 R001 2018-10-27 1899-12-31 16:51:45 38.4
R001 <- R001_G2_Q3_T
Test Subject Datum Zeit Temperatur Timepoint
1 G2 R001 2018-10-27 1899-12-31 16:51:45 38.4 Q3
2 G2 R001 2018-10-27 1899-12-31 16:51:55 38.4 Q3
3 G2 R001 2018-10-27 1899-12-31 16:52:05 38.4 Q3
4 G2 R001 2018-10-27 1899-12-31 16:52:15 38.4 Q3
5 G2 R001 2018-10-27 1899-12-31 16:52:25 38.4 Q3
6 G2 R001 2018-10-27 1899-12-31 16:52:35 38.3 Q3
7 G2 R001 2018-10-27 1899-12-31 16:52:45 38.4 Q3
8 G2 R001 2018-10-27 1899-12-31 16:52:55 38.5 Q3
9 G2 R001 2018-10-27 1899-12-31 16:53:06 38.3 Q3
10 G2 R001 2018-10-27 1899-12-31 16:53:16 38.4 Q3
11 G2 R001 2018-10-27 1899-12-31 16:53:26 38.4 Q3
12 G2 R001 2018-10-27 1899-12-31 16:53:36 38.4 Q3
13 G2 R001 2018-10-27 1899-12-31 16:53:46 38.4 Q3
14 G2 R001 2018-10-27 1899-12-31 16:53:56 38.4 Q3
15 G2 R001 2018-10-27 1899-12-31 16:54:06 38.3 Q3
16 G2 R001 2018-10-27 1899-12-31 16:54:16 38.3 Q3
17 G2 R001 2018-10-27 1899-12-31 16:54:26 38.3 Q3
18 G2 R001 2018-10-27 1899-12-31 16:54:36 38.4 Q3
19 G2 R001 2018-10-27 1899-12-31 16:54:47 38.4 Q3
20 G2 R001 2018-10-27 1899-12-31 16:54:57 38.3 Q3
21 G2 R001 2018-10-27 1899-12-31 16:55:07 38.4 Q3
22 G2 R001 2018-10-27 1899-12-31 16:55:17 38.4 Q3
23 G2 R001 2018-10-27 1899-12-31 16:55:27 38.4 Q3
24 G2 R001 2018-10-27 1899-12-31 16:55:37 38.4 Q3
25 G2 R001 2018-10-27 1899-12-31 16:55:47 38.5 Q3
26 G2 R001 2018-10-27 1899-12-31 16:55:57 38.4 Q3
27 G2 R001 2018-10-27 1899-12-31 16:56:07 38.4 Q3
28 G2 R001 2018-10-27 1899-12-31 16:56:17 38.4 Q3
29 G2 R001 2018-10-27 1899-12-31 16:56:28 38.5 Q3
30 G2 R001 2018-10-27 1899-12-31 16:56:38 38.4 Q3
31 G2 R001 2018-10-27 1899-12-31 16:56:48 38.4 Q3
32 G2 R001 2018-10-27 1899-12-31 16:56:58 38.5 Q3
33 G2 R001 2018-10-27 1899-12-31 16:57:08 38.5 Q3
[...]
R001以下是使用data.table
对您的问题的尝试
样本数据
library( data.table )
DT <- fread("ID Test Subject Datum Zeit Temperatur Timepoint
1 G2 R001 2018-10-27 '1899-12-31 16:51:45' 38.4 Q3
2 G2 R001 2018-10-27 '1899-12-31 16:51:55' 38.4 Q3
3 G2 R001 2018-10-27 '1899-12-31 16:52:05' 38.4 Q3
4 G2 R001 2018-10-27 '1899-12-31 16:52:15' 38.4 Q3
5 G2 R001 2018-10-27 '1899-12-31 16:52:25' 38.4 Q3
6 G2 R001 2018-10-27 '1899-12-31 16:52:35' 38.3 Q3
7 G2 R001 2018-10-27 '1899-12-31 16:52:45' 38.4 Q3
8 G2 R001 2018-10-27 '1899-12-31 16:52:55' 38.5 Q3
9 G2 R001 2018-10-27 '1899-12-31 16:53:06' 38.3 Q3
10 G2 R001 2018-10-27 '1899-12-31 16:53:16' 38.4 Q3
11 G2 R001 2018-10-27 '1899-12-31 16:53:26' 38.4 Q3
12 G2 R001 2018-10-27 '1899-12-31 16:53:36' 38.4 Q3
13 G2 R001 2018-10-27 '1899-12-31 16:53:46' 38.4 Q3
14 G2 R001 2018-10-27 '1899-12-31 16:53:56' 38.4 Q3
15 G2 R001 2018-10-27 '1899-12-31 16:54:06' 38.3 Q3
16 G2 R001 2018-10-27 '1899-12-31 16:54:16' 38.3 Q3
17 G2 R001 2018-10-27 '1899-12-31 16:54:26' 38.3 Q3
18 G2 R001 2018-10-27 '1899-12-31 16:54:36' 38.4 Q3
19 G2 R001 2018-10-27 '1899-12-31 16:54:47' 38.4 Q3
20 G2 R001 2018-10-27 '1899-12-31 16:54:57' 38.3 Q3
21 G2 R001 2018-10-27 '1899-12-31 16:55:07' 38.4 Q3
22 G2 R001 2018-10-27 '1899-12-31 16:55:17' 38.4 Q3
23 G2 R001 2018-10-27 '1899-12-31 16:55:27' 38.4 Q3
24 G2 R001 2018-10-27 '1899-12-31 16:55:37' 38.4 Q3
25 G2 R001 2018-10-27 '1899-12-31 16:55:47' 38.5 Q3
26 G2 R001 2018-10-27 '1899-12-31 16:55:57' 38.4 Q3
27 G2 R001 2018-10-27 '1899-12-31 16:56:07' 38.4 Q3
28 G2 R001 2018-10-27 '1899-12-31 16:56:17' 38.4 Q3
29 G2 R001 2018-10-27 '1899-12-31 16:56:28' 38.5 Q3
30 G2 R001 2018-10-27 '1899-12-31 16:56:38' 38.4 Q3
31 G2 R001 2018-10-27 '1899-12-31 16:56:48' 38.4 Q3
32 G2 R001 2018-10-27 '1899-12-31 16:56:58' 38.5 Q3
33 G2 R001 2018-10-27 '1899-12-31 16:57:08' 38.5 Q3", quote = "'")
输出
# period max_deltaT
# 1: 2018-10-27T16:51 0.0
# 2: 2018-10-27T16:52 0.1
# 3: 2018-10-27T16:53 0.1
# 4: 2018-10-27T16:54 0.1
# 5: 2018-10-27T16:55 0.1
# 6: 2018-10-27T16:56 0.1
# 7: 2018-10-27T16:57 0.0
我猜这就是你想要的：
DT[, maxTslope_last6 := frollapply(Temperatur, n = 6, function(x) max(diff(x)))]
> DT
ID Test Subject Datum Zeit Temperatur Timepoint maxTslope_last6
1: 1 G2 R001 2018-10-27 1899-12-31 16:51:45 38.4 Q3 NA
2: 2 G2 R001 2018-10-27 1899-12-31 16:51:55 38.4 Q3 NA
3: 3 G2 R001 2018-10-27 1899-12-31 16:52:05 38.4 Q3 NA
4: 4 G2 R001 2018-10-27 1899-12-31 16:52:15 38.4 Q3 NA
5: 5 G2 R001 2018-10-27 1899-12-31 16:52:25 38.4 Q3 NA
6: 6 G2 R001 2018-10-27 1899-12-31 16:52:35 38.3 Q3 0.0
7: 7 G2 R001 2018-10-27 1899-12-31 16:52:45 38.4 Q3 0.1
8: 8 G2 R001 2018-10-27 1899-12-31 16:52:55 38.5 Q3 0.1
9: 9 G2 R001 2018-10-27 1899-12-31 16:53:06 38.3 Q3 0.1
10: 10 G2 R001 2018-10-27 1899-12-31 16:53:16 38.4 Q3 0.1
11: 11 G2 R001 2018-10-27 1899-12-31 16:53:26 38.4 Q3 0.1
12: 12 G2 R001 2018-10-27 1899-12-31 16:53:36 38.4 Q3 0.1
13: 13 G2 R001 2018-10-27 1899-12-31 16:53:46 38.4 Q3 0.1
14: 14 G2 R001 2018-10-27 1899-12-31 16:53:56 38.4 Q3 0.1
15: 15 G2 R001 2018-10-27 1899-12-31 16:54:06 38.3 Q3 0.0
16: 16 G2 R001 2018-10-27 1899-12-31 16:54:16 38.3 Q3 0.0
17: 17 G2 R001 2018-10-27 1899-12-31 16:54:26 38.3 Q3 0.0
18: 18 G2 R001 2018-10-27 1899-12-31 16:54:36 38.4 Q3 0.1
19: 19 G2 R001 2018-10-27 1899-12-31 16:54:47 38.4 Q3 0.1
20: 20 G2 R001 2018-10-27 1899-12-31 16:54:57 38.3 Q3 0.1
21: 21 G2 R001 2018-10-27 1899-12-31 16:55:07 38.4 Q3 0.1
22: 22 G2 R001 2018-10-27 1899-12-31 16:55:17 38.4 Q3 0.1
23: 23 G2 R001 2018-10-27 1899-12-31 16:55:27 38.4 Q3 0.1
24: 24 G2 R001 2018-10-27 1899-12-31 16:55:37 38.4 Q3 0.1
25: 25 G2 R001 2018-10-27 1899-12-31 16:55:47 38.5 Q3 0.1
26: 26 G2 R001 2018-10-27 1899-12-31 16:55:57 38.4 Q3 0.1
27: 27 G2 R001 2018-10-27 1899-12-31 16:56:07 38.4 Q3 0.1
28: 28 G2 R001 2018-10-27 1899-12-31 16:56:17 38.4 Q3 0.1
29: 29 G2 R001 2018-10-27 1899-12-31 16:56:28 38.5 Q3 0.1
30: 30 G2 R001 2018-10-27 1899-12-31 16:56:38 38.4 Q3 0.1
31: 31 G2 R001 2018-10-27 1899-12-31 16:56:48 38.4 Q3 0.1
32: 32 G2 R001 2018-10-27 1899-12-31 16:56:58 38.5 Q3 0.1
33: 33 G2 R001 2018-10-27 1899-12-31 16:57:08 38.5 Q3 0.1
ID Test Subject Datum Zeit Temperatur Timepoint maxTslope_last6
请添加一些预期结果。您好，请查看此链接。希望这会有所帮助。谢谢，我已经添加了我喜欢的输出，分别是我想尝试的作为计算结果的输出。在16:55和16:56之间，温度始终在38.4和38.5之间，我看不到任何0.3°C的变化。你能更好地解释一下你预期结果背后的计算吗？这只是一个例子，与数据无关。我想找出一分钟内温度的最大上升量。我认为一种方法是取六个值并计算平均值，而不是向前移动一个值并再次计算平均值等等。然后计算每个平均值之间的差值，然后找到最大值。我希望这有助于理解我的请求。感谢这个好的解决方案。它使我更接近我所需要的。但我不需要一分钟内最高值和最低值之间的差值，而是一分钟内值的增加或减少。例如，在前6个值期间，平均增加0.1°，然后向前移动一个值，然后取下6个值，平均增加0.2°，依此类推。然后，找到这个增量的最大值应该很容易。我建议您创建一些更具代表性的样本数据和所需的输出..谢谢。我添加了另一个图像解释。这可能有助于理解我的愿望。
DT[, maxTslope_last6 := frollapply(Temperatur, n = 6, function(x) max(diff(x)))]
> DT
ID Test Subject Datum Zeit Temperatur Timepoint maxTslope_last6
1: 1 G2 R001 2018-10-27 1899-12-31 16:51:45 38.4 Q3 NA
2: 2 G2 R001 2018-10-27 1899-12-31 16:51:55 38.4 Q3 NA
3: 3 G2 R001 2018-10-27 1899-12-31 16:52:05 38.4 Q3 NA
4: 4 G2 R001 2018-10-27 1899-12-31 16:52:15 38.4 Q3 NA
5: 5 G2 R001 2018-10-27 1899-12-31 16:52:25 38.4 Q3 NA
6: 6 G2 R001 2018-10-27 1899-12-31 16:52:35 38.3 Q3 0.0
7: 7 G2 R001 2018-10-27 1899-12-31 16:52:45 38.4 Q3 0.1
8: 8 G2 R001 2018-10-27 1899-12-31 16:52:55 38.5 Q3 0.1
9: 9 G2 R001 2018-10-27 1899-12-31 16:53:06 38.3 Q3 0.1
10: 10 G2 R001 2018-10-27 1899-12-31 16:53:16 38.4 Q3 0.1
11: 11 G2 R001 2018-10-27 1899-12-31 16:53:26 38.4 Q3 0.1
12: 12 G2 R001 2018-10-27 1899-12-31 16:53:36 38.4 Q3 0.1
13: 13 G2 R001 2018-10-27 1899-12-31 16:53:46 38.4 Q3 0.1
14: 14 G2 R001 2018-10-27 1899-12-31 16:53:56 38.4 Q3 0.1
15: 15 G2 R001 2018-10-27 1899-12-31 16:54:06 38.3 Q3 0.0
16: 16 G2 R001 2018-10-27 1899-12-31 16:54:16 38.3 Q3 0.0
17: 17 G2 R001 2018-10-27 1899-12-31 16:54:26 38.3 Q3 0.0
18: 18 G2 R001 2018-10-27 1899-12-31 16:54:36 38.4 Q3 0.1
19: 19 G2 R001 2018-10-27 1899-12-31 16:54:47 38.4 Q3 0.1
20: 20 G2 R001 2018-10-27 1899-12-31 16:54:57 38.3 Q3 0.1
21: 21 G2 R001 2018-10-27 1899-12-31 16:55:07 38.4 Q3 0.1
22: 22 G2 R001 2018-10-27 1899-12-31 16:55:17 38.4 Q3 0.1
23: 23 G2 R001 2018-10-27 1899-12-31 16:55:27 38.4 Q3 0.1
24: 24 G2 R001 2018-10-27 1899-12-31 16:55:37 38.4 Q3 0.1
25: 25 G2 R001 2018-10-27 1899-12-31 16:55:47 38.5 Q3 0.1
26: 26 G2 R001 2018-10-27 1899-12-31 16:55:57 38.4 Q3 0.1
27: 27 G2 R001 2018-10-27 1899-12-31 16:56:07 38.4 Q3 0.1
28: 28 G2 R001 2018-10-27 1899-12-31 16:56:17 38.4 Q3 0.1
29: 29 G2 R001 2018-10-27 1899-12-31 16:56:28 38.5 Q3 0.1
30: 30 G2 R001 2018-10-27 1899-12-31 16:56:38 38.4 Q3 0.1
31: 31 G2 R001 2018-10-27 1899-12-31 16:56:48 38.4 Q3 0.1
32: 32 G2 R001 2018-10-27 1899-12-31 16:56:58 38.5 Q3 0.1
33: 33 G2 R001 2018-10-27 1899-12-31 16:57:08 38.5 Q3 0.1
ID Test Subject Datum Zeit Temperatur Timepoint maxTslope_last6