作为学习dplyr及其同类的问题。
我正在计算一个因子的相对频率和一个以df中的另外两个变量为条件的相对频率。例如:
library(dplyr)
library(tidyr)
set.seed(3457)
pct <- function(x) {x/sum(x)}
foo <- data.frame(x = rep(seq(1:3),20),
y = rep(rep(c("a","b"),each=3),10),
z = LETTERS[floor(runif(60, 1,5))])
bar <- foo %>%
group_by(x, y, z) %>%
tally %>%
mutate(freq = (n / sum(n)) * 100)
head(bar)
我希望输出bar
包含foo$z
的所有级别。即,这里没有C
的案例:
subset(bar, x==2 & y=="a")
我怎样才能bar
计算缺失的等级,所以我得到:
subset(bar, x==2 & y=="a",select = n)
返回4,5,0,1(和select = freq
给40,50,0,10)?
非常感谢。
编辑:使用种子集进行游戏!
答案 0 :(得分:1)
我们可以使用complete
tidyr
bar1 <- bar %>%
complete(z, nesting(x, y), fill = list(n = 0, freq = 0))%>%
select_(.dots = names(bar))
filter(bar1, x==2 & y=="a")
# x y z n freq
# <int> <fctr> <fctr> <dbl> <dbl>
#1 2 a A 4 40
#2 2 a B 5 50
#3 2 a C 0 0
#4 2 a D 1 10