rmoving-averagehighest
How can I find the maximal increase over 6 values in a dataset
I have a data set of core temperature values. Each 10th second I have collected a temperature value, over a duration of several hours. I try to find the highest increase (Tslope) of the temperature within 6 temperature values. I tried the approach of rolling average but I was not even close to a solution. The data set looks like this:
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
[...]
What I expect is that I can say for example: "Between 16:55 and 16:56 was with 0.3°C the highest increase in core temperature"
Maybe it is easier to understand my problem, when I show you additionally an image:
What I like to find is Tslope: Tslope will be calculated iteratively over a one minute period, including 6 successive measurement time points. Of course my data-set contains much more values that there will be a similar curve like in the image.
It would be great if someone could help me in this case. Thanks for your help.
Solution
I'm guessing this is what you are looking for:
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
- How to hide NAs when using xlsx::saveWorkbook?
- How do I retrieve a simple numeric value from a named numeric vector in R?
- Matching pair-wise columns from left to right across rows in one dataframe to another dataframe and adding new columns with matching values
- Income to outcome flow chart in Sankey plotly R
- color mapping in geom_conn_bundle not showing correctly
- Print R package startup message AFTER automatic package conflict messages instead of before
- Summing a set of R dataframe rows (column-wise), while retaining the first n columns
- Added variable / partial regression plots for groups in an interaction?
- how to make a topoplot in R with coordinates variable distribution
- List of all functions in base R?
- Plotting multiple plots for different initial conditions in one graph
- Printing repetitively on the same line in R
- Generating UI/Server based on initial selection
- Subset dataframe based on pickerInput
- How to let user pick the data in R-shiny?
- Couldn't show my simple bar charts separately on Shiny R dashboardBody
- How to programmatically filter contents of a second shiny app displayed via iframe
- How to select specific interesting groups for the boxplot in R Shiny app?
- Crosstable and Plot grouping with reactive values
- Is there a way to make multiple Shiny picker inputs where the selections must be disjoint?
- Delay/avoid duplication of shiny server side functions until after credentials
- Predictions only returns value "1"
- How to display a busy indicator in a shiny app?
- Append doesn't work when writing to CSV in R
- Changing the start date of a gantt chart in DiagrammeR
- Check for installed packages before running install.packages()
- Compare two columns element-wise
- For loop is several times faster in R than in Python using the rpy2 library
- Setting a flag based on two samples' dates
- How to bring pop-ups for long text string selection items forward in front of modal they are rendered from when using virtualSelectInput()?