R语言绘制掷5枚骰子出现6的次数的概率分布直方图

要绘制掷5枚骰子时，6出现次数对应的概率分布直方图，可直接运行以下R代码实现：

library(polynom)
library(ggplot2)
library(gridExtra)

pd <- function(n){
  s = 6
  p<-polynomial(c(0,rep(1/s,s)))
  k <- data.frame(as.vector(p^n)[-(1:n)])
  k$no <- n:(s*n)
  names(k)[1] <- "p"
  ggplot(data.frame(k)) +
    geom_bar(stat = "identity",aes(no, p), size = 1, fill = "orange", col = "black") +
    labs(title = paste(n, ifelse(n == 1,"die roll", "dice rolls"), sep = " "),
         x = "Outcome", y = "Probability")
}

grid.arrange(pd(5), ncol =2)

代码运行后生成的概率分布效果如下：

内容的提问来源于stack exchange，提问作者nodes4codes