Использование padr с сгущением приводит к ошибке "отсутствует значение там, где требуется ИСТИНА / ЛОЖЬ"
Я пытался заставить padr работать с моим набором данных без особого успеха, хотя я могу заставить примеры работать:
# I have a few datetime columns so I convert all to POSIXct with UTC.
> df <- mutate_at(DATABASE, vars(ends_with("time")), funs(ymd_hms(., tz = "UTC", locale = Sys.getlocale("LC_TIME"))))
> df <- as_tibble(df)
> head(df, 20)
# A tibble: 20 x 2
charttime sbp
<dttm> <dbl>
1 2101-10-20 22:30:01 NA
2 2101-10-20 18:45:00 62
3 2101-10-20 19:00:00 66
4 2101-10-20 19:12:00 NA
5 2101-10-20 19:14:00 NA
6 2101-10-20 19:15:00 217
7 2101-10-20 19:26:00 NA
8 2101-10-20 19:30:00 102
9 2101-10-20 19:45:00 94
10 2101-10-20 19:59:00 NA
11 2101-10-20 20:00:00 80
12 2101-10-20 20:04:00 NA
13 2101-10-20 20:15:00 91
14 2101-10-20 20:30:00 86
15 2101-10-20 20:45:00 96
16 2101-10-20 21:00:00 73
17 2101-10-20 21:15:00 84
18 2101-10-20 21:30:00 96
19 2101-10-20 21:45:00 100
20 2101-10-20 21:51:00 NA
> df$charttime %>% get_interval # should say 'sec'
[1] "sec"
> df %>% thicken(interval='hour')
Error in if (to_date) x <- as.Date(x, tz = attr(x, "tzone")) :
missing value where TRUE/FALSE needed
Но на примере падра это работает:
> coffee %>% thicken(interval='day')
time_stamp amount time_stamp_day
1 2016-07-07 03:11:21 3.14 2016-07-07
2 2016-07-07 03:46:48 2.98 2016-07-07
3 2016-07-09 07:25:17 4.11 2016-07-09
4 2016-07-10 04:45:11 3.14 2016-07-10
> coffee$time_stamp %>% get_interval # should say 'sec'
[1] "sec"
Я не смог понять, почему мой набор данных не работает и как интерпретировать ошибку.
Обновление 1
Вот еще один, более полный пример того, что я пытаюсь сделать. Я также включил CSV с небольшим фрагментом данных, с которыми я работаю, чтобы эта проблема была более воспроизводимой. Я пробовал это на двух машинах, и я получил тот же результат.
Вы заметите, что в примере выше и примере ниже первое значение диаграммы времени отличается. (2101-10-20 22:30:01 изменяется на 2101-10-20 22:30:00). Я хотел, чтобы интервал был "сек" вместо "мин", поэтому я вручную изменил значение. В любом случае это приводит к одной и той же проблеме.
> packageVersion("tidyverse")
[1] ‘1.1.1’
> packageVersion("lubridate")
[1] ‘1.6.0’
> packageVersion("padr")
[1] ‘0.3.0’
> library(tidyverse)
> library(lubridate)
> library(padr)
>
> df <- read.csv("padr_data.csv")
> df <- mutate_at(df, vars(ends_with("time")), funs(ymd_hms(., tz = "UTC", locale = Sys.getlocale("LC_TIME"))))
> df$sbp <- as.numeric(df$sbp)
> summary(df)
charttime sbp
Min. :2101-10-20 18:30:00 Min. : 62.0
1st Qu.:2101-10-20 19:33:45 1st Qu.: 84.5
Median :2101-10-20 20:52:30 Median : 95.0
Mean :2101-10-20 21:08:22 Mean :100.9
3rd Qu.:2101-10-20 22:26:15 3rd Qu.:102.0
Max. :2101-10-21 00:42:00 Max. :217.0
NA's :12
> lapply(df, class)
$charttime
[1] "POSIXct" "POSIXt"
$sbp
[1] "numeric"
> df$charttime %>% get_interval
[1] "min"
>
> # this does not work
> df[!is.na(df$charttime),] %>%
+ thicken(interval = 'hour')
Error in if (to_date) x <- as.Date(x, tz = attr(x, "tzone")) :
missing value where TRUE/FALSE needed
>
> # this does not work
> df %>%
+ thicken(interval = 'hour')
Error in if (to_date) x <- as.Date(x, tz = attr(x, "tzone")) :
missing value where TRUE/FALSE needed
2 ответа
Решение
Похоже, что падр плохо работает с датами, установленными в будущем! В частности, дата, которая в будущем будет длиться более 20 лет, не будет работать. Я открою проблему с разработчиком падра, чтобы посмотреть, как можно улучшить код.
> packageVersion("tidyverse")
[1] ‘1.1.1’
> packageVersion("lubridate")
[1] ‘1.6.0’
> packageVersion("padr")
[1] ‘0.3.0’
> library(tidyverse)
> library(lubridate)
> library(padr)
>
> df <- read.csv("padr_data.csv")
> df <- mutate_at(df, vars(ends_with("time")), funs(ymd_hms(., tz = "UTC", locale = Sys.getlocale("LC_TIME"))- dyears(63)))
>
> df$sbp <- as.numeric(df$sbp)
> #df <- na.omit(df)
>
> summary(df)
charttime sbp
Min. :2038-11-04 18:30:00 Min. : 62.0
1st Qu.:2038-11-04 19:33:45 1st Qu.: 84.5
Median :2038-11-04 20:52:30 Median : 95.0
Mean :2038-11-04 21:08:22 Mean :100.9
3rd Qu.:2038-11-04 22:26:15 3rd Qu.:102.0
Max. :2038-11-05 00:42:00 Max. :217.0
NA's :12
> lapply(df, class)
$charttime
[1] "POSIXct" "POSIXt"
$sbp
[1] "numeric"
> df$charttime %>% get_interval
[1] "min"
>
> # this does not work
> df[!is.na(df$charttime),] %>%
+ thicken(interval = 'hour')
Error in if (to_date) x <- as.Date(x, tz = attr(x, "tzone")) :
missing value where TRUE/FALSE needed
In addition: Warning messages:
1: In round_down_core(a, b) : NAs introduced by coercion to integer range
2: In round_down_core(a, b) : NAs introduced by coercion to integer range
+ Изменить dyears(63) в dyears(64)
> df <- mutate_at(df, vars(ends_with("time")), funs(ymd_hms(., tz = "UTC", locale = Sys.getlocale("LC_TIME"))- dyears(64)))
>
> df$sbp <- as.numeric(df$sbp)
> #df <- na.omit(df)
>
> summary(df)
charttime sbp
Min. :2037-11-04 18:30:00 Min. : 62.0
1st Qu.:2037-11-04 19:33:45 1st Qu.: 84.5
Median :2037-11-04 20:52:30 Median : 95.0
Mean :2037-11-04 21:08:22 Mean :100.9
3rd Qu.:2037-11-04 22:26:15 3rd Qu.:102.0
Max. :2037-11-05 00:42:00 Max. :217.0
NA's :12
> lapply(df, class)
$charttime
[1] "POSIXct" "POSIXt"
$sbp
[1] "numeric"
> df$charttime %>% get_interval
[1] "min"
>
> # this does work
> df[!is.na(df$charttime),] %>%
+ thicken(interval = 'hour')
charttime sbp charttime_hour
1 2037-11-04 18:30:00 NA 2037-11-04 18:00:00
2 2037-11-04 18:45:00 62 2037-11-04 18:00:00
3 2037-11-04 19:00:00 66 2037-11-04 19:00:00
4 2037-11-04 19:12:00 NA 2037-11-04 19:00:00
5 2037-11-04 19:14:00 NA 2037-11-04 19:00:00
6 2037-11-04 19:15:00 217 2037-11-04 19:00:00
7 2037-11-04 19:26:00 NA 2037-11-04 19:00:00
8 2037-11-04 19:30:00 102 2037-11-04 19:00:00
9 2037-11-04 19:45:00 94 2037-11-04 19:00:00
10 2037-11-04 19:59:00 NA 2037-11-04 19:00:00
11 2037-11-04 20:00:00 80 2037-11-04 20:00:00
12 2037-11-04 20:04:00 NA 2037-11-04 20:00:00
13 2037-11-04 20:15:00 91 2037-11-04 20:00:00
14 2037-11-04 20:30:00 86 2037-11-04 20:00:00
15 2037-11-04 20:45:00 96 2037-11-04 20:00:00
16 2037-11-04 21:00:00 73 2037-11-04 21:00:00
17 2037-11-04 21:15:00 84 2037-11-04 21:00:00
18 2037-11-04 21:30:00 96 2037-11-04 21:00:00
19 2037-11-04 21:45:00 100 2037-11-04 21:00:00
20 2037-11-04 21:51:00 NA 2037-11-04 21:00:00
21 2037-11-04 22:00:00 NA 2037-11-04 22:00:00
22 2037-11-04 22:15:00 123 2037-11-04 22:00:00
23 2037-11-04 22:30:00 125 2037-11-04 22:00:00
24 2037-11-04 22:45:00 132 2037-11-04 22:00:00
25 2037-11-04 23:00:00 88 2037-11-04 23:00:00
26 2037-11-04 23:15:00 NA 2037-11-04 23:00:00
27 2037-11-04 23:45:00 NA 2037-11-04 23:00:00
28 2037-11-05 00:00:00 102 2037-11-05 00:00:00
29 2037-11-05 00:28:00 NA 2037-11-05 00:00:00
30 2037-11-05 00:42:00 NA 2037-11-05 00:00:00
РЕШЕНИЕ 1 - НЕ РАБОТАЕТ
Не очень хорошо знаю пакет, но я бы попробовал две вещи:
- фильтрация значений NA
- декларирование
byаргумент
Попробуй это
df[!is.na(df$sbp),] %>% thicken(interval='hour', by = 'charttime')
РЕШЕНИЕ 2 - НЕ РАБОТАЕТ
попробуйте принуждение df в фрейм данных вместо таблицы, также попробуйте принуждение charttime в дату позже:
df <- data.frame(df)
df$charttime <- as.POSIXct(df$charttime)
РЕШЕНИЕ 3 - НЕ РАБОТАЕТ
Вы можете иметь некоторые NAна вашем charttime, попробуй это:
df[!is.na(df$charttime),] %>% thicken(interval = 'hour')
Я попытался переименовать переменную, но это не проблема. Извините, но я пока не могу комментировать. Пожалуйста, скажите мне, если это сработало.