首页 > 解决方案 > 找到R中数字总和恰好为70的概率?

问题描述

问题是从掷出的 20 个骰子中找出数字总和正好等于 70 的概率。

我创建了一个矩阵

A = matrix(NA, nrow=20, ncol=1e4)
for (i in 1:20) {
  A[i,]=sample(1:6,1e4,replace=TRUE)
}

我被困在下一步该怎么做。我应该创建一个长度为 10^4 的向量 s,使得 s[i] = 位置 i 中的列的总和。

标签: rsimulationprobability

解决方案


您可以创建一个 20 x 10,000 的样本骰子矩阵,然后计算列总和,然后计算恰好是 70 的那些。概率估计就是这个结果除以 10,000。所以这个单行就足够了:

length(which(colSums(replicate(1e4, sample(1:6, 20, replace = TRUE))) == 70))/1e4
#> [1] 0.0526

推荐阅读