在R中计算投掷20枚骰子点数之和恰好等于70的概率

• 生成列和向量s
  你可以直接调用R内置的按列求和函数colSums()，返回结果就是长度为1e4的向量，对应矩阵每一列的元素和：

s <- colSums(A)

• 计算点数和等于70的概率
  统计向量s中取值等于70的样本占总模拟次数的比例，就是你需要的模拟概率：

dice_prob <- sum(s == 70) / ncol(A)

如果需要更高的计算精度，可以把你初始化矩阵时的模拟次数1e4调整为更大的数值（例如1e5），样本量越大模拟结果越接近理论概率。

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