dplyr でグループごとの quantile をカラムに追加したくなったんですが、複数カラムを追加する方法がよくわかりませんでした。
例えば、iris データに対して、Species ごとの Sepal.Length の quantile が知りたかったとします。
それはこのようにすれば一応は結果が見れるんですが、
R
iris %>% group_by(Species) %>%
summarize(sepal_length_quantile = quantile(Sepal.Length) %>% list %>% as.character)
結果
Species sepal_length_quantile
1 setosa c(4.3, 4.8, 5, 5.2, 5.8)
2 versicolor c(4.9, 5.6, 5.9, 6.3, 7)
3 virginica c(4.9, 6.225, 6.5, 6.9, 7.9)
やはりこの sepal_length_quantile を複数列に分けて見たいですよね。
色々やり方を試してみて、下記に落ち着きました。
R
iris %>% group_by(Species) %>%
summarize(sepal_length_quantile = quantile(Sepal.Length) %>% list) %>%
cbind(do.call(rbind, .$sepal_length_quantile)) %>%
select(-sepal_length_quantile)
結果
Species 0% 25% 50% 75% 100%
1 setosa 4.3 4.800 5.0 5.2 5.8
2 versicolor 4.9 5.600 5.9 6.3 7.0
3 virginica 4.9 6.225 6.5 6.9 7.9
一旦リストとして保存して、do.call で matrix を作り、cbind で結合するという方法です。
ちなみに、リストとして保存するところは dplyr::do を使った方が正式かもしれません。
R
iris %>% group_by(Species) %>%
do(sepal_length_quantile = quantile(.$Sepal.Length)) %>%
cbind(do.call(rbind, .$sepal_length_quantile)) %>%
select(-sepal_length_quantile)
結果
Species 0% 25% 50% 75% 100%
1 setosa 4.3 4.800 5.0 5.2 5.8
2 versicolor 4.9 5.600 5.9 6.3 7.0
3 virginica 4.9 6.225 6.5 6.9 7.9
enjoy!