我在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()
密度图:
ggplot(data.frame(sums), aes(x=factor(sums), fill=factor(sums)))+geom_density()
我也尝试过:
ggplot(data.frame(sums),aes(x = factor(sums),fill = factor(sums)))+ geom_density(aes(y = ..count ..))
答案 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))