2

我想得到一个结合了观察密度和 cdf 的图。

通常的问题是两者的比例相差甚远。如何解决这个问题,即使用两个比例,或者,重新调整其中一个数据系列(最好在 ggplot 中,因为我想将数据的计算和显示分开)。

这是到目前为止的代码:

>dput(tmp) 产量

structure(list(drivenkm = c(8, 11, 21, 4, 594, 179, 19, 7, 10, 36)), .Names = "drivenkm", class = c("data.table", "data.frame" ), row.names = c(NA, -10L), .internal.selfref = <pointer: 0x223cb78>)

然后我做

p = ggplot(data = tmp, aes(x = drivenkm)) + geom_histogram(aes(y = ..density..), alpha = 0.2, binwidth = 3) + stat_ecdf(aes(x = drivenkm)); print(p)

我得到的是以下内容:

在此处输入图像描述

显然,天秤已经很远了。如何解决这个问题,以便可以以合理的方式解释直方图和 cdf?

谢谢!

4

1 回答 1

5

密度由 binwidth 缩放,因此面积总和为 1。因此y,您的直方图也应该乘以这个:

p = ggplot(data = tmp, aes(x = drivenkm)) +
   geom_histogram(aes(y = 3*..density..), alpha = 0.2, binwidth = 3) +
   stat_ecdf(aes(x = drivenkm))

在此处输入图像描述

于 2014-01-14T09:54:45.120 回答