r - 这种类型的图是什么，你如何在 R 中绘制它？

Question

我认为这将是一种叫做累积和发现累积频率图和累积流量图的东西。但是，我认为图像中的图表也不是，因为累积图表从 0 开始，但我的变量不是。此外，密度图听起来最接近，但它是 1 区域上的分布，但我想显示频率。

基本上，变量是主变量的子部分，我想显示这些子变量何时收敛以创建峰值。本质上，这些变量相加以显示累积界限。

Answer 1

使用ggplot2您可以使用该geom_area()功能

library(ggplot2)
library(gcookbook) # For the data set

ggplot(uspopage, aes(x=Year, y=Thousands, fill=AgeGroup)) + geom_area()

Answer 2

感谢您分享更多关于您的数据是什么样子的信息。

让我们以休斯顿警察局公开的犯罪统计数据为例。在本例中，我们使用 2015 年 1 月的数据集。

library(ggplot2)

crime <- gdata::read.xls('http://www.houstontx.gov/police/cs/xls/jan15.xls')

# There's a single case in there where the offense type is called '1',
# that doesn't make sense to us so we'll remove it.
crime <- crime[!crime$Offense.Type == '1', ]
crime$Offense.Type <- droplevels(crime$Offense.Type)

有 10 列，但我们感兴趣的列如下所示：

# Hour Offense.Type
# 8   Auto Theft
# 13  Theft
# 5   Auto Theft
# 13  Theft
# 18  Theft
# 18  Theft

正如您所提到的，问题在于每一行都是一个事件。我们需要一种方法来获取每小时的频率以传递给geom_area().

第一种方式是让ggplot2处理，不需要预先格式化数据。

p <- ggplot(crime, aes(x=Hour, fill=Offense.Type)) 
p + geom_area(aes(y = ..count..), stat='density')

另一种方法是使用 R'stable()和 reshape2's预先格式化频率表melt()：

library(reshape2)
crime.counts <- table(crime$Hour, crime$Offense.Type)
crime.counts.l <- melt(crime.counts,
                        id.vars = c('Hour'),
                        value.name = "NumberofCrimes")

names(crime.counts.l) <- c("Hour", "Offense.Type", "numberOfCrimes")
p <- ggplot(crime.counts.l, aes(x = Hour,
                                 y = numberOfCrimes,
                                 fill = Offense.Type))
p + geom_area()