【发布时间】:2015-06-22 06:44:07
【问题描述】:
有人可以解释这里出了什么问题吗?我想模拟 10,000 个 20 面掷骰子。我使用了这段代码:
x <- sample(1:20,10000,replace=T)
但这给了我这个:
hist(x)
好像是12以上的问题:
我在这里不明白什么? 谢谢
【问题讨论】:
标签: r
有人可以解释这里出了什么问题吗?我想模拟 10,000 个 20 面掷骰子。我使用了这段代码:
x <- sample(1:20,10000,replace=T)
但这给了我这个:
hist(x)
好像是12以上的问题:
我在这里不明白什么? 谢谢
【问题讨论】:
标签: r
这实际上与您的sample 无关,而是hist。
如果你这样做
set.seed(1)
x <- sample(1:20,10000,replace=T)
table(x)
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
513 522 482 495 459 549 506 505 518 498 495 492 440 490 459 509 496 528 511 533
你会发现它是随机的。但是hist 复制了您的图表。如果你数一下柱子,你会发现有 19 个而不是 20 个。
试试这个:
bins <- seq(0, 20, by=1)
hist(x, breaks=bins)
给出一个条形高度均匀的图表,因为所有 20 个条形都显示了(即 1 和 2 没有折叠在一起)。
【讨论】:
barplot(table(x))