-1

我有一个如下的数据集，它有很多列。有一些列的标题是：

baked_hamburgur,spinach,mashed_potato,cabbages,jello,rolls,brown,milk,coffee,water,cakes,vanilla,chocolate,fruitsalad

还有其他专栏，但我现在只对上述专栏感兴趣。

这些列的每一行中的值是：是或否。

此数据的屏幕截图如下，因为我无法在问题本身中附加/共享此文件。

我的数据集有很多列，用于对病人进行数据分析

dput(head()) 输出如下：

> dput(head(illness_data))

structure(list(Age = structure(c(18L, 26L, 22L, 25L, 29L, 13L
), .Label = c("10", "106", "11", "12", "14", "15", "16", "17", 
"18", "19", "2", "20", "22", "23", "24", "25", "26", "27", "28", 
"30", "31", "32", "33", "34", "36", "38", "39", "4", "42", "43", 
"44", "45", "46", "48", "5", "7", "8", "9", "seven"), class = "factor"), 
    sex = structure(c(3L, 2L, 3L, 3L, 2L, 3L), .Label = c("-1", 
    "Female", "Male"), class = "factor"), timesupper = c(2000L, 
    1830L, 1830L, 1930L, 1930L, 1930L), ill = structure(c(1L, 
    1L, 1L, 1L, 1L, 1L), .Label = "yes", class = "factor"), onsetdate = structure(c(4L, 
    4L, 4L, 1L, 1L, 4L), .Label = c("18-Apr", "18-Jun", "18/4", 
    "19-Apr"), class = "factor"), onsettime = c(30L, 30L, 30L, 
    2230L, 2230L, 200L), baked_hamburgur = structure(c(2L, 2L, 
    2L, 2L, 2L, 1L), .Label = c("no", "yes"), class = "factor"), 
    spinach = structure(c(2L, 2L, 2L, 2L, 2L, 1L), .Label = c("no", 
    "yes"), class = "factor"), mashed_potato = structure(c(2L, 
    2L, 1L, 1L, 2L, 1L), .Label = c("no", "yes"), class = "factor"), 
    cabbages = structure(c(1L, 2L, 1L, 2L, 1L, 1L), .Label = c("no", 
    "yes"), class = "factor"), jello = structure(c(1L, 1L, 1L, 
    2L, 2L, 1L), .Label = c("no", "yes"), class = "factor"), 
    rolls = structure(c(2L, 1L, 1L, 1L, 2L, 1L), .Label = c("no", 
    "yes"), class = "factor"), brown = structure(c(1L, 1L, 1L, 
    1L, 2L, 1L), .Label = c("no", "yes"), class = "factor"), 
    milk = structure(c(1L, 1L, 1L, 1L, 1L, 1L), .Label = c("no", 
    "yes"), class = "factor"), coffee = structure(c(2L, 2L, 2L, 
    1L, 2L, 1L), .Label = c("no", "yes"), class = "factor"), 
    water = structure(c(1L, 1L, 1L, 2L, 2L, 1L), .Label = c("no", 
    "yes"), class = "factor"), cakes = structure(c(1L, 1L, 2L, 
    1L, 1L, 1L), .Label = c("no", "yes"), class = "factor"), 
    vanilla = structure(c(2L, 2L, 2L, 2L, 2L, 2L), .Label = c("no", 
    "yes"), class = "factor"), chocolate = structure(c(1L, 2L, 
    2L, 1L, 1L, 2L), .Label = c("no", "yes"), class = "factor"), 
    fruitsalad = structure(c(1L, 1L, 1L, 1L, 1L, 1L), .Label = c("no", 
    "yes", "yes</pre></body></html>Ztext/plain\b\v\035(F]l~Ó_Ý\026R\002\001"
    ), class = "factor")), .Names = c("Age", "sex", "timesupper", 
"ill", "onsetdate", "onsettime", "baked_hamburgur", "spinach", 
"mashed_potato", "cabbages", "jello", "rolls", "brown", "milk", 
"coffee", "water", "cakes", "vanilla", "chocolate", "fruitsalad"
), row.names = c(NA, 6L), class = "data.frame")

完整的 dput 命令输出如下：

> dput(illness_data)

structure(list(Age = structure(c(18L, 26L, 22L, 25L, 29L, 13L, 
36L, 8L, 11L, 7L, 24L, 10L, 8L, 35L, 34L, 6L, 22L, 39L, 12L, 
9L, 36L, 17L, 9L, 20L, 37L, 27L, 32L, 30L, 21L, 24L, 3L, 18L, 
33L, 16L, 5L, 31L, 28L, 14L, 19L, 38L, 2L, 4L, 23L, 1L, 18L, 
15L), .Label = c("10", "106", "11", "12", "14", "15", "16", "17", 
"18", "19", "2", "20", "22", "23", "24", "25", "26", "27", "28", 
"30", "31", "32", "33", "34", "36", "38", "39", "4", "42", "43", 
"44", "45", "46", "48", "5", "7", "8", "9", "seven"), class = "factor"), 
    sex = structure(c(3L, 2L, 3L, 3L, 2L, 3L, 3L, 3L, 2L, 3L, 
    3L, 2L, 3L, 3L, 2L, 2L, 3L, 3L, 3L, 2L, 2L, 1L, 3L, 3L, 3L, 
    2L, 2L, 3L, 2L, 3L, 3L, 3L, 2L, 3L, 2L, 2L, 3L, 3L, 3L, 3L, 
    3L, 3L, 2L, 3L, 2L, 3L), .Label = c("-1", "Female", "Male"
    ), class = "factor"), timesupper = c(2000L, 1830L, 1830L, 
    1930L, 1930L, 1930L, 2200L, 1900L, 1930L, NA, NA, NA, NA, 
    2200L, NA, NA, NA, 2200L, NA, NA, 2200L, 2200L, NA, NA, 2200L, 
    NA, NA, NA, NA, NA, 1900L, NA, 1100L, NA, NA, NA, 2200L, 
    1930L, 1930L, 2200L, NA, NA, 1930L, 1930L, NA, NA), ill = structure(c(1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L
    ), .Label = "yes", class = "factor"), onsetdate = structure(c(4L, 
    4L, 4L, 1L, 1L, 4L, 4L, 2L, 4L, 4L, 4L, 1L, 1L, 4L, 1L, 3L, 
    1L, 4L, 1L, 1L, 4L, 4L, 1L, 1L, 4L, 1L, 4L, 4L, 1L, 4L, 4L, 
    1L, 1L, 1L, 1L, 1L, 4L, 4L, 1L, 4L, 4L, 4L, 4L, 1L, 4L, 1L
    ), .Label = c("18-Apr", "18-Jun", "18/4", "19-Apr"), class = "factor"), 
    onsettime = c(30L, 30L, 30L, 2230L, 2230L, 200L, 100L, 2300L, 
    200L, 1030L, 30L, 2215L, 2200L, 100L, 2300L, 2145L, 2145L, 
    100L, 2300L, 2100L, 100L, 100L, 2115L, 2330L, 100L, 2130L, 
    230L, 200L, 2130L, 30L, 100L, 2230L, 1500L, 2400L, 2300L, 
    2230L, 100L, 230L, 2330L, 100L, 30L, 30L, 100L, 2400L, 215L, 
    2300L), baked_hamburgur = structure(c(2L, 2L, 2L, 2L, 2L, 
    1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 
    1L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 
    2L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 2L), .Label = c("no", 
    "yes"), class = "factor"), spinach = structure(c(2L, 2L, 
    2L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 
    1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 
    1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 2L), .Label = c("no", 
    "yes"), class = "factor"), mashed_potato = structure(c(2L, 
    2L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 
    2L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 
    2L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 2L
    ), .Label = c("no", "yes"), class = "factor"), cabbages = structure(c(1L, 
    2L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 
    2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 
    2L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L
    ), .Label = c("no", "yes"), class = "factor"), jello = structure(c(1L, 
    1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 
    1L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 
    1L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 2L
    ), .Label = c("no", "yes"), class = "factor"), rolls = structure(c(2L, 
    1L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 
    2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 
    2L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L
    ), .Label = c("no", "yes"), class = "factor"), brown = structure(c(1L, 
    1L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 
    2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 1L, 
    1L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L
    ), .Label = c("no", "yes"), class = "factor"), milk = structure(c(1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 
    1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L
    ), .Label = c("no", "yes"), class = "factor"), coffee = structure(c(2L, 
    2L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 
    1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 1L, 
    2L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L
    ), .Label = c("no", "yes"), class = "factor"), water = structure(c(1L, 
    1L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 
    2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 
    1L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 1L
    ), .Label = c("no", "yes"), class = "factor"), cakes = structure(c(1L, 
    1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 
    1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L
    ), .Label = c("no", "yes"), class = "factor"), vanilla = structure(c(2L, 
    2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 
    2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L
    ), .Label = c("no", "yes"), class = "factor"), chocolate = structure(c(1L, 
    2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 
    1L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, NA, 1L, 1L, 
    2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 1L
    ), .Label = c("no", "yes"), class = "factor"), fruitsalad = structure(c(1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
    1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 
    1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 3L
    ), .Label = c("no", "yes", "yes</pre></body></html>Ztext/plain\b\v\035(F]l~Ó_Ý\026R\002\001"
    ), class = "factor")), .Names = c("Age", "sex", "timesupper", 
"ill", "onsetdate", "onsettime", "baked_hamburgur", "spinach", 
"mashed_potato", "cabbages", "jello", "rolls", "brown", "milk", 
"coffee", "water", "cakes", "vanilla", "chocolate", "fruitsalad"
), class = "data.frame", row.names = c(NA, -46L))

R 已将这些列正确读取为因子类型变量。

现在，这些列中的每一个都对应于医院中每个病人的消费量。

我想知道病人吃的最多的食物，使用 R。

请建议一个好的方法。谢谢！

请注意，除了下面此 URL 中提到的选项外，我没有尝试过任何其他选项。但是，我无法让它工作。

使用 R 计算因子列

4

2 回答 2

1

我不完全确定您在寻找什么，但这将计算食物的食用频率（使用tidyverse包装）：

library(tidyverse)
illness_data_summed <- illness_data %>%
  mutate_at(vars(-Age, -sex, -timesupper,-onsetdate,-onsettime), ~ifelse(. == "yes", 1,0)) %>% 
  summarise_at(vars(-Age, -sex, -timesupper,-onsetdate,-onsettime, -ill), ~sum(., na.rm = TRUE)) 

illness_data_summed[which(illness_data_summed == max(illness_data_summed))]

所以首先我将 yes 转换为 1，no 转换为 0，这使得总和表示特定食物被吃掉的次数。我对所有列执行此操作，除了那些您不感兴趣的列（由 - in vars 表示），但如果需要，您也可以反转它（例如，当您想要转换的 var 数量低于您不想要的想转换）。

最后一部分将导致：

  vanilla
1      43
于 2019-12-13T11:44:39.007 回答
1

由于所有行都有ill = 'yes'，我们可以在每一列中计算其中包含的值的数量yes。基本 R 方法可以是：

head(sort(colSums(illness_data[7:20] == "yes"), decreasing = TRUE), 5)

#        vanilla baked_hamburgur           cakes         spinach   mashed_potato 
#             43              29              27              26              23

我选择了第 7 到 20 列，因为这些是唯一存在食物的列。另外我只选择了前 5 个值，您可以通过更改head命令中的数字 5 来选择任何值。

于 2019-12-18T10:56:57.867 回答