биннинг числа с неправильным результатом - PullRequest
Новая
       9

биннинг числа с неправильным результатом

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

У меня проблемы с выводом после того, как я перебрал числовой вектор.

Я пытаюсь ограничить продолжительность пребывания, которая была рассчитана заранее с помощью функции difftime. Предоставлять весь код не имеет смысла, поскольку это только фон. Тем не менее, когда я бен, я не получаю правильный ответ. Вот продолжительность пребывания, назначенная с лос. 

dput(los)
c(61.0416666666667, 61.0416666666667, 61.0416666666667, 2, 2, 3, 3)

Вот мои перерывы. Я использовал na.rm внутри, так как попробовал несколько методов. Я сдал na.rm с TRUE, FALSE и вынул его из своих перерывов. 

breaks <- c(0, 0.8, 0.16, 
           1.0, 1.8, 1.16, 
           2.0, 2.8, 2.16, 
           3.0, 3.8, 3.16, 
           4.0, 4.8, 4.16, 
           5.0, 5.8, 5.16, 
           6.0, 6.8, 6.16, 
           7.0, 14.0, 21.0, 28.0, max(los)) #, , na.rm = FALSE

Тем не менее, следующий код пробовал 

dt_los$losbinned <- cut(dt_los$LOS, 
                         breaks = breaks,
                         labels = c("0hrs", "8hrs", "16hrs", "1 d", 
                                 "1 d 8hrs", "1 d 16hrs", "2 d", 
                                 "2 d 8hrs", "2 d 16hrs", "3 d", 
                                 "3 d 8hrs", "3 d 16hrs", "4 d", 
                                 "4 d 8hrs", "4 d 16hrs", "5 d", 
                                 "5 d 8hrs", "5 d 16hrs", "6 d",  
                                 "6 d 8hrs","6 d 16hrs", "7 - 14 d", 
                                 "14 - 21 d", "21 - 28 d", "> 28 d"), 
                         right = FALSE)#

с различными параметрами, переданными для 'right', дает мне следующее:

когда верно = ЛОЖЬ Я не получаю LOS для бен 61,14 для категории "> 28 дней". BНо все же получите правильные корзины для других 2,00 и 3,00. 

structure(list(IDcol = 101:107, Admissions = structure(c(1539160200, 
1539160200, 1539160200, 1539154800, 1539154800, 1539154800, 1539154800
), class = c("POSIXct", "POSIXt"), tzone = "Europe/London"), 
    Discharges = structure(c(1544434200, 1544434200, 1544434200, 
    1539327600, 1539327600, 1539414000, 1539414000), class = c("POSIXct", 
    "POSIXt"), tzone = "Europe/London"), Admission_type = c("Elective", 
    "Emergency", "Emergency", "Elective", "Emergency", "Elective", 
    "Emergency"), LOS = c(61.0416666666667, 61.0416666666667, 
    61.0416666666667, 2, 2, 3, 3), Ward_code = c("DSN", "DSN", 
    "DNA", "NAS", "BAS", "BAS", "BAS"), Same_day_discharge = c(FALSE, 
    FALSE, FALSE, FALSE, FALSE, FALSE, FALSE), Spell_type = c("Elective", 
    "Emergency", "Emergency", "Elective", "Emergency", "Elective", 
    "Emergency"), Adm_period = c(TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE), losbinned = structure(c(NA, NA, NA, 7L, 7L, 
    10L, 10L), .Label = c("0hrs", "8hrs", "16hrs", "1 d", "1 d 8hrs", 
    "1 d 16hrs", "2 d", "2 d 8hrs", "2 d 16hrs", "3 d", "3 d 8hrs", 
    "3 d 16hrs", "4 d", "4 d 8hrs", "4 d 16hrs", "5 d", "5 d 8hrs", 
    "5 d 16hrs", "6 d", "6 d 8hrs", "6 d 16hrs", "7 - 14 d", 
    "14 - 21 d", "21 - 28 d", "> 28 d"), class = "factor")), row.names = c(NA, 
-7L), class = c("tbl_df", "tbl", "data.frame"))

когда я передаю право = ИСТИНА, вывод для 61.04 складывается в «> 28 d», что является желаемым ответом, но я не получаю правильные ячейки для 2.0 и 3.0, которые складываются за 1 д 16 часов для 2,0 и 2 д 16 часов в течение 3. И снова, они должны быть объединены в 2, соответственно 3. 

 structure(list(IDcol = 101:107, Admissions = structure(c(1539160200, 
1539160200, 1539160200, 1539154800, 1539154800, 1539154800, 1539154800
), class = c("POSIXct", "POSIXt"), tzone = "Europe/London"), 
    Discharges = structure(c(1544434200, 1544434200, 1544434200, 
    1539327600, 1539327600, 1539414000, 1539414000), class = c("POSIXct", 
    "POSIXt"), tzone = "Europe/London"), Admission_type = c("Elective", 
    "Emergency", "Emergency", "Elective", "Emergency", "Elective", 
    "Emergency"), LOS = c(61.0416666666667, 61.0416666666667, 
    61.0416666666667, 2, 2, 3, 3), Ward_code = c("DSN", "DSN", 
    "DNA", "NAS", "BAS", "BAS", "BAS"), Same_day_discharge = c(FALSE, 
    FALSE, FALSE, FALSE, FALSE, FALSE, FALSE), Spell_type = c("Elective", 
    "Emergency", "Emergency", "Elective", "Emergency", "Elective", 
    "Emergency"), Adm_period = c(TRUE, TRUE, TRUE, TRUE, TRUE, 
    TRUE, TRUE), losbinned = structure(c(25L, 25L, 25L, 6L, 6L, 
    9L, 9L), .Label = c("0hrs", "8hrs", "16hrs", "1 d", "1 d 8hrs", 
    "1 d 16hrs", "2 d", "2 d 8hrs", "2 d 16hrs", "3 d", "3 d 8hrs", 
    "3 d 16hrs", "4 d", "4 d 8hrs", "4 d 16hrs", "5 d", "5 d 8hrs", 
    "5 d 16hrs", "6 d", "6 d 8hrs", "6 d 16hrs", "7 - 14 d", 
    "14 - 21 d", "21 - 28 d", "> 28 d"), class = "factor")), row.names = c(NA, 
-7L), class = c("tbl_df", "tbl", "data.frame"))

Фактические и ожидаемые результаты должны быть указаны в правильных ячейках для моего срока пребывания. Для 61.04 -> «> 28d», для 2 -> «2 d», для 3 -> «3 d».

Если это можно сделать с помощью Tidyverse, это было бы замечательно. Но уважая ящики, которые я назначил. Тем не менее, я знаю, что это еще не сделано. Поэтому хорошо с исправленным кодом, который я придумал, но исправил.

1 Ответ

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

Бункеры функции вырезки являются эксклюзивными и включительно.

Из справки функции вырезания: Метки уровня фактора создаются как "(b1, b2]", "(b2, b3]" и т. Д. Для right = TRUE и как "[b1, b2)"

Чтобы включить самое низкое значение (или самое высокое значение в этом случае), требуется опция include.lowest=TRUE. Это сделает первый контейнер эксклюзивным для "[b1, b2]".

Попробуйте: 

labels<-c("0hrs", "8hrs", "16hrs", "1 d", 
          "1 d 8hrs", "1 d 16hrs", "2 d", 
          "2 d 8hrs", "2 d 16hrs", "3 d", 
          "3 d 8hrs", "3 d 16hrs", "4 d", 
          "4 d 8hrs", "4 d 16hrs", "5 d", 
          "5 d 8hrs", "5 d 16hrs", "6 d",  
          "6 d 8hrs","6 d 16hrs", "7 - 14 d", 
          "14 - 21 d", "21 - 28 d", "> 28 d")

dt_los$losbinned <- cut(los, breaks=breaks, labels=labels, right=FALSE, include.lowest = TRUE)
...