【问题标题】:Histogram for probability distribution in RR中概率分布的直方图
【发布时间】:2021-12-03 05:30:34
【问题描述】:

我试图弄清楚如何在掷五个骰子时绘制 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)

【问题讨论】:

  • 请问您这个问题的问题/目的是什么?
  • 我只是被要求编写几行代码来绘制掷 5 个骰子时在 6 秒数上的概率分布的直方图。
  • 谢谢!这样我就可以了解,为什么这是一个更好的方法?

标签: r histogram probability-distribution


【解决方案1】:

你可以试试

library(stats)
pd <- function(n){
  df <- data.frame(
    prob = dbinom(0:n, n, 1/6),
    m = 0:n
  )
  ggplot(df, aes(m, prob)) +
    geom_col(fill = "steelblue") +
    geom_line(color = "red", lwd = 1)
}
pd(5)

【讨论】:

    猜你喜欢
    • 2018-04-16
    • 2014-04-24
    • 2021-08-19
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    相关资源
    最近更新 更多