我在R中有一个包含43个变量和80多行的数据帧。我想基于一个变量 - 地理区域对数据进行分组,然后计算变量的不同值(多少0,1s,2s,3s和NA等)。
我知道group_by
中的summarize
和tidyverse
函数,我知道我可以使用像“sum”和“mean”这样的函数,但我想要计算
我试过了
est1 <- df %>%
group_by(region) %>%
summarize(count)
数据如下所示:
iso3 Country WHOregion WBIncomeGroup UrbanSanPol UrbanSanWom UrbanSanExt RuralSanPol RuralSanWom
<chr> <chr> <chr> <chr> <chr> <chr> <chr> <chr> <chr>
1 AFG Afghanistan EMRO Low income 0 <NA> <NA> 1 1
2 ALB Albania EURO Upper middle income 1 0 0 1 0
3 ARG Argentina PAHO Upper middle income 1 0 0.5 1 0
4 AZE Azerbaijan EURO Upper middle income 1 1 0.5 1 1
5 BDI Burundi AFRO Low income 1 1 0.5 1 1
6 BFA Burkina Faso AFRO Low income 1 1 1 1 1
但这不是我想要的。有人可以帮忙吗?
structure(list(iso3 = c("AFG", "ALB", "ARG", "AZE", "BDI", "BFA",
"BGD", "BIH", "BLR", "BOL"), Country = c("Afghanistan", "Albania",
"Argentina", "Azerbaijan", "Burundi", "Burkina Faso", "Bangladesh",
"Bosnia and Herzegovina", "Belarus", "Bolivia (Plurinational State of)"
), WHOregion = c("EMRO", "EURO", "PAHO", "EURO", "AFRO", "AFRO",
"SEARO", "EURO", "EURO", "PAHO"), WBIncomeGroup = c("Low income",
"Upper middle income", "Upper middle income", "Upper middle income",
"Low income", "Low income", "Lower middle income", "Upper middle income",
"Upper middle income", "Lower middle income"), UrbanSanPol = c("0",
"1", "1", "1", "1", "1", "1", "1", "1", "1"), UrbanSanWom = c(NA,
"0", "0", "1", "1", "1", "1", "0", NA, "0"), UrbanSanExt = c(NA,
"0", "0.5", "0.5", "0.5", "1", "0.5", "0", "0.5", "0"), RuralSanPol = c("1",
"1", "1", "1", "1", "1", "1", "1", "1", "1"), RuralSanWom = c("1",
"0", "0", "1", "1", "1", "1", "0", NA, "0"), RuralSanExt = c("0.5",
"0", "0", "0.5", "0.5", "1", "0.5", "0", "0.5", "0.5"), UrbanDWPol = c("0",
"1", "1", "1", "1", "1", "1", "1", "1", "1"), UrbanDWWom = c(NA,
"0", "0", "1", "1", "1", "1", "0", NA, "0"), UrbanDWExt = c(NA,
"0", "0.5", "1", "0", "0.5", "0.5", "0.5", "0.5", "0"), RuralDWPol = c("1",
"1", "1", "1", "1", "1", "1", "1", "1", "1"), RuralDWWom = c("1",
"0", "0", "1", "1", "1", "1", "0", NA, "0"), RuralDWExt = c("0.5",
"0", "0", "1", "0.5", "1", "0.5", "0.5", "0.5", "0.5"), HygienePol = c("1",
"1", "0", "1", "1", "1", "1", "1", "1", "0"), HygieneWom = c("1",
NA, NA, "1", "1", "1", "1", "0", NA, "0"), HygieneExt = c("0.5",
NA, NA, "0", "0.5", "0", "0.5", "0", "0.5", "0"), WASHHealthPol = c("1",
"1", "0", "1", "1", "1", "1", "1", "0", "0"), WASHHealthWom = c("0",
NA, NA, "1", "1", "1", "1", "0", NA, "0"), WASHHealthExt = c("0",
NA, "0.5", "1", "0", "0.5", "0", "0", NA, "0"), WpollutionPol = c("1",
"1", "1", "1", "1", "1", "1", "1", "1", "0"), WpollutionWom = c("1",
NA, "0", "1", "1", "1", "1", "0", NA, "0"), WpollutionExt = c("0",
NA, "0", "1", "0", "0.5", "0", "0", "0.5", "0"), WQMPol = c("1",
"1", "1", "1", "1", "1", "1", "1", "1", "0"), WQMWom = c("1",
NA, "0", "1", "1", "1", "1", "0", NA, "0"), WQMExt = c("0", NA,
"0", "1", "0", "0.5", "0", "0", "0.5", "0"), WatRightPol = c("0",
"1", "1", "1", NA, "1", "1", "1", "1", "1"), WatRightWom = c("0",
NA, "0", "1", NA, "1", "1", "0", NA, "0"), WatRightExt = c("0",
NA, "0.5", "1", NA, "1", "0", "0", "0.5", "0.5"), WRMPol = c("1",
"1", "1", "1", "1", "1", "1", "1", "1", "1"), WRMWom = c("0",
NA, "0", "1", "1", "1", "1", "0", NA, "0"), WRMExt = c("0", NA,
"0.5", "1", "0.5", "1", "0", "0", "0.5", "0"), EnvProtPol = c("1",
"1", "1", "1", "1", "1", "1", "1", "1", "1"), EnvProtWom = c("0",
NA, "0", "1", "1", "1", "1", "0", NA, "0"), EnvProtExt = c("0",
NA, "0", "1", "0", "1", "0", "0", "0.5", "0"), `SDG regions` = c("Central Asia (M49) and Southern Asia (MDG=M49)",
"Northern America (M49) and Europe (M49)", "Latin America and the Caribbean (MDG=M49)",
"Western Asia (M49) and Northern Africa (M49)", "Sub-Saharan Africa (M49)",
"Sub-Saharan Africa (M49)", "Central Asia (M49) and Southern Asia (MDG=M49)",
"Northern America (M49) and Europe (M49)", "Northern America (M49) and Europe (M49)",
"Latin America and the Caribbean (MDG=M49)"), M49_level1 = c("Asia (M49)",
"Europe (M49)", "Latin America and the Caribbean (MDG=M49)",
"Asia (M49)", "Sub-Saharan Africa (M49)", "Sub-Saharan Africa (M49)",
"Asia (M49)", "Europe (M49)", "Europe (M49)", "Latin America and the Caribbean (MDG=M49)"
), M49_level2 = c("Southern Asia (MDG=M49)", "Southern Europe (M49)",
"South America (M49)", "Western Asia (M49)", "Eastern Africa (M49)",
"Western Africa (M49)", "Southern Asia (MDG=M49)", "Southern Europe (M49)",
"Eastern Europe (M49)", "South America (M49)"), LDCs = c("Least Developed Countries (LDCs)",
NA, NA, NA, "Least Developed Countries (LDCs)", "Least Developed Countries (LDCs)",
"Least Developed Countries (LDCs)", NA, NA, NA), LLDCS_SIDS = c("Landlocked developing countries (LLDCs)",
NA, NA, "Landlocked developing countries (LLDCs)", "Landlocked developing countries (LLDCs)",
"Landlocked developing countries (LLDCs)", NA, NA, NA, "Landlocked developing countries (LLDCs)"
), `Income group` = c("Low income", "Upper middle income", "Upper middle income",
"Upper middle income", "Low income", "Low income", "Lower middle income",
"Upper middle income", "Upper middle income", "Lower middle income"
)), .Names = c("iso3", "Country", "WHOregion", "WBIncomeGroup",
"UrbanSanPol", "UrbanSanWom", "UrbanSanExt", "RuralSanPol", "RuralSanWom",
"RuralSanExt", "UrbanDWPol", "UrbanDWWom", "UrbanDWExt", "RuralDWPol",
"RuralDWWom", "RuralDWExt", "HygienePol", "HygieneWom", "HygieneExt",
"WASHHealthPol", "WASHHealthWom", "WASHHealthExt", "WpollutionPol",
"WpollutionWom", "WpollutionExt", "WQMPol", "WQMWom", "WQMExt",
"WatRightPol", "WatRightWom", "WatRightExt", "WRMPol", "WRMWom",
"WRMExt", "EnvProtPol", "EnvProtWom", "EnvProtExt", "SDG regions",
"M49_level1", "M49_level2", "LDCs", "LLDCS_SIDS", "Income group"
), row.names = c(NA, -10L), class = c("tbl_df", "tbl", "data.frame"
在这里输入代码
答案 0 :(得分:0)
想象一下,我有三列。第一个是国家列表(法国,德国等)。第二个是地区列表(亚洲,欧洲),第三个是每个国家的离散值(奥运金牌数量)。我想按区域对所有数据进行分组,并计算每个区域的次数,0次发生,1次发生,2次发生。
根据你在评论中所说的内容,以及我所理解的内容:
解释:
df %>% select(continent,countries,medals) %>% group_by(continent,countries) %>% summarize(count =n())
另一种解释:你想要的是每个大陆以及由它赢得的不同数量的奖牌。
这些数字是该国赢得的奥运奖牌。
df <- as.data.frame(matrix(c("Asia","Asia","Asia","Asia","Europe","Europe","India","China","Bangladesh","Japan","Spain", "Italy",6,3,4,4,3,3),ncol = 3))
df %>% group_by(V1,V3) %>% summarise(count= n()) %>% spread(V3,count)
给我一个输出
# A tibble: 2 x 4
# Groups: V1 [2]
V1 `3` `4` `6`
* <fctr> <int> <int> <int>
1 Asia 1 2 1
2 Europe 2 NA NA
答案 1 :(得分:-1)
试试这个;需要dplyr和tidyverse
distinct_cnt <- input_df %>%
gather(variable, value) %>%
group_by(variable) %>%
summarise(n_distinct(value))