2个骰子结果的不对称密度图滚动

时间:2014-10-02 07:54:05

标签: r plot distribution probability

我在https://stats.stackexchange.com/questions/117578/density-distribution-of-outcomes-of-2-dice-rolled上发布了以下问题,但没有得到任何回复。

此问题与:A histogram with a bar for each frequency value

有关

滚动两个骰子并绘制总和。直方图如预期,但密度图显示不同密度2和12,图是不对称的。为什么会这样? 在2个骰子滚动的结果中,2的概率等于12的概率。为什么密度图具有不相等的值?

num.dices <- 2L
num.rolls <- 100000L
outcomes <- matrix(sample(1:6, num.dices * num.rolls, replace = TRUE),
                   nrow = num.rolls, ncol = num.dices)
sums <- rowSums(outcomes)

直方图:

ggplot(data.frame(sums), aes(x=factor(sums)))+geom_histogram()

enter image description here

密度图:

ggplot(data.frame(sums), aes(x=factor(sums), fill=factor(sums)))+geom_density()

enter image description here

我也尝试过:

ggplot(data.frame(sums),aes(x = factor(sums),fill = factor(sums)))+ geom_density(aes(y = ..count ..))

enter image description here

1 个答案:

答案 0 :(得分:0)

在我看来,你正在使用11 diff进行11种不同的内核估计。 bandwitdths,改为:

ggplot(data.frame(sums), aes(x=sums, fill=2))+geom_density()

或者你可以添加group = 1,如果你坚持使用额外的参数:

ggplot(data.frame(sums), aes(x=sums, fill=factor(sums)))+geom_density(aes(group=1))