Использование dplyr для выбора данных за определенный час и использования их для расчета вычислений - PullRequest
Новая
       5

Использование dplyr для выбора данных за определенный час и использования их для расчета вычислений

0 голосов
/ 16 января 2019

У меня есть время в Канаде / горное время. Данные за час. За день я хочу взять значение F в 0:29:05 (Fn) (Fn = значение F сразу после полуночи по местному времени), и для каждого дня я хочу рассчитать Z = (Fn-F) / Fn

Тем не менее, ему нужно выбрать Fn для 0:29:05 в Cadadian / Mountain, и для каждого часа в день z нужно вычислять, используя Fn того дня.

Фиктивные данные: 

 datetime <- seq(
 from=as.POSIXct("2012-1-1 0:29:05", tz="Canada/Mountain"),
 to=as.POSIXct("2012-2-1 0:29:05", tz="Canada/Mountain"),
 by="hour")

 #variable F
 F <- runif(745, min = 0, max =2)

df <- as.data.frame(cbind(datetime,F))
library(lubridate)
#make sure its in  "POSIXct" "POSIXt" format
df$datetime <- as_datetime(df$datetime)

Теперь у меня уже была некоторая помощь в dplyr при использовании минутного набора данных - но, очевидно, мое понимание все еще довольно слабое, поскольку я не могу перевести это на мой пример данных часа-1. Ниже приведена моя попытка ... и я думаю, что, возможно, мутация является правильным вариантом в этом случае? 

df2 <- df %>%
group_by(Date = as.Date(datetime)) %>%
mutate(Fn = F[hour(datetime) == 0]), 
z = (Fn - F)/Fn) %>%
ungroup() %>%
select(-Date)

Спасибо.

1 Ответ

0 голосов
/ 16 января 2019 
  library(lubridate)
  library(tidyverse)

datetime <- seq(
   from = as.POSIXct("2012-1-1 0:29:05", tz = "Canada/Mountain"),
   to = as.POSIXct("2012-2-1 0:29:05", tz = "Canada/Mountain"),
   by = "hour"
   )

f <- runif(745, min = 0, max =2) #variable f 
df <- data.frame(datetime, f)

# method using fill function from tidyr package
df2 <- df %>%
   mutate(Date = as.Date(datetime, tz = "Canada/Mountain")) %>% 
   left_join( #this will grab the f value at 0:29:05 of each day
     df %>% filter(hour(datetime) == 0) %>% select(datetime, Fn = f),
     by = 'datetime'
   ) %>% 
   group_by(Date) %>% 
   fill(Fn, Fn, .direction = 'down') %>% #fills in NA values with f values of the following day
   mutate(
     Z = ( Fn - f ) / Fn
   ) %>% 
   ungroup() %>% 
   select(-Date)

# method not using fill
df3 <- df %>%
   mutate(Date = as.Date(datetime, tz = "Canada/Mountain")) %>% 
   left_join( #this will grab the f value at 0:29:05 of each day
     df %>% filter(hour(datetime) == 0) %>% select(datetime, Fn = f),
     by = 'datetime'
   ) %>% 
   group_by(Date) %>% 
   mutate(
     Fn = na.omit(Fn),
     Z = ( Fn - f ) / Fn
   ) %>% 
   ungroup() %>% 
   select(-Date)

# both methods result in the same result, as shown below
# A tibble: 745 x 4
   datetime                f    Fn       Z
   <dttm>              <dbl> <dbl>   <dbl>
 1 2012-01-01 00:29:05 0.590 0.590  0     
 2 2012-01-01 01:29:05 1.57  0.590 -1.66  
 3 2012-01-01 02:29:05 0.537 0.590  0.0900
 4 2012-01-01 03:29:05 0.691 0.590 -0.171 
 5 2012-01-01 04:29:05 0.719 0.590 -0.218 
 6 2012-01-01 05:29:05 0.811 0.590 -0.374 
 7 2012-01-01 06:29:05 0.248 0.590  0.581 
 8 2012-01-01 07:29:05 1.98  0.590 -2.35  
 9 2012-01-01 08:29:05 0.839 0.590 -0.422 
10 2012-01-01 09:29:05 0.733 0.590 -0.242 
# ... with 735 more rows

Просто мысль, но вы не должны называть объекты R F, поскольку она обычно зарезервирована для FALSE значений.

...