67

主要问题

在尝试使用 ggplot2 制作直方图时,我无法理解为什么日期、标签和中断的处理不像我在 R 中所期望的那样工作。

我在找:

  • 我的约会频率的直方图
  • 在匹配栏下方居中的刻度线
  • 日期标签%Y-b格式
  • 适当的限制;最小化网格空间边缘和最外层条之间的空白空间

我已将我的数据上传到 pastebin以使其可重现。我创建了几列,因为我不确定最好的方法:

> dates <- read.csv("http://pastebin.com/raw.php?i=sDzXKFxJ", sep=",", header=T)
> head(dates)
       YM       Date Year Month
1 2008-Apr 2008-04-01 2008     4
2 2009-Apr 2009-04-01 2009     4
3 2009-Apr 2009-04-01 2009     4
4 2009-Apr 2009-04-01 2009     4
5 2009-Apr 2009-04-01 2009     4
6 2009-Apr 2009-04-01 2009     4

这是我尝试过的:

library(ggplot2)
library(scales)
dates$converted <- as.Date(dates$Date, format="%Y-%m-%d")

ggplot(dates, aes(x=converted)) + geom_histogram()
+      opts(axis.text.x = theme_text(angle=90))

这产生了这张图。不过,我想要%Y-%b格式化,所以我四处寻找并尝试了以下内容,基于这个 SO

ggplot(dates, aes(x=converted)) + geom_histogram()
+    scale_x_date(labels=date_format("%Y-%b"),
+    breaks = "1 month")
+    opts(axis.text.x = theme_text(angle=90))

stat_bin: binwidth defaulted to range/30. Use 'binwidth = x' to adjust this.

这给了我这张图

  • 正确的 x 轴标签格式
  • 频率分布已改变形状(binwidth 问题?)
  • 刻度线未在条形下方居中显示
  • xlims 也发生了变化

我在该部分的ggplot2 文档中完成了示例,当我将它与相同的 x 轴数据一起使用时,它似乎正确地中断、标记和居中刻度。我不明白为什么直方图不同。scale_x_dategeom_line()


根据 edgester 和 gauden 的回答进行更新

我最初认为gauden的回答帮助我解决了我的问题,但现在仔细观察后感到困惑。请注意代码后两个答案的结果图之间的差异。

假设两者:

library(ggplot2)
library(scales)
dates <- read.csv("http://pastebin.com/raw.php?i=sDzXKFxJ", sep=",", header=T)

根据@edgester 下面的回答,我能够做到以下几点:

freqs <- aggregate(dates$Date, by=list(dates$Date), FUN=length)
freqs$names <- as.Date(freqs$Group.1, format="%Y-%m-%d")

ggplot(freqs, aes(x=names, y=x)) + geom_bar(stat="identity") +
       scale_x_date(breaks="1 month", labels=date_format("%Y-%b"),
                    limits=c(as.Date("2008-04-30"),as.Date("2012-04-01"))) +
       ylab("Frequency") + xlab("Year and Month") +
       theme_bw() + opts(axis.text.x = theme_text(angle=90))

这是我根据高登回答的尝试:

dates$Date <- as.Date(dates$Date)
ggplot(dates, aes(x=Date)) + geom_histogram(binwidth=30, colour="white") +
       scale_x_date(labels = date_format("%Y-%b"),
                    breaks = seq(min(dates$Date)-5, max(dates$Date)+5, 30),
                    limits = c(as.Date("2008-05-01"), as.Date("2012-04-01"))) +
       ylab("Frequency") + xlab("Year and Month") +
       theme_bw() + opts(axis.text.x = theme_text(angle=90))

基于 edgester 方法的绘图:

边缘图

基于高登方法的绘图:

高登图

请注意以下事项:

  • 2009 年 12 月和 2010 年 3 月 gauden 图中的差距;显示数据table(dates$Date)中有 19 个实例2009-12-01和 26 个实例2010-03-01
  • edgester 的情节从 2008 年 4 月开始,到 2012 年 5 月结束。根据 2008 年 4 月 1 日数据中的最小值和 2012 年 5 月 1 日的最大值,这是正确的。由于某种原因,高登的情节从 2008 年 3 月开始,但不知何故仍设法在 2012 年 5 月结束。在计算垃圾箱并阅读月份标签之后,对于我的生活,我无法弄清楚哪个情节有额外的或缺少直方图的垃圾箱!

对这里的差异有什么想法吗?edgester 的创建单独计数的方法


相关参考

顺便说一句,这里还有其他位置,其中包含有关日期和 ggplot2 的信息,供路人寻求帮助:

  • 从learnr.wordpress 开始,这是一个流行的 R 博客它说我需要将我的数据转换为 POSIXct 格式,我现在认为这是错误的并且浪费了我的时间。
  • 另一个学习者帖子在 ggplot2 中重新创建了一个时间序列,但并不真正适用于我的情况。
  • r-bloggers 对此有一个帖子,但它似乎已经过时了。简单的format=选项对我不起作用。
  • 这个 SO 问题正在使用中断和标签。我尝试将我的Date向量视为连续的,但认为它效果不佳。看起来它一遍又一遍地覆盖相同的标签文本,所以这些字母看起来有点奇怪。分布是正确的,但有一些奇怪的休息。我基于接受的答案的尝试是这样的(结果here)。
4

4 回答 4

36

更新

版本 2:使用 Date 类

我更新了示例以演示对齐标签并在图上设置限制。我还证明了as.Date在一致使用时确实有效(实际上它可能比我之前的示例更适合您的数据)。

目标情节 v2

基于日期的直方图

守则 v2

这是(有点过分)注释代码:

library("ggplot2")
library("scales")

dates <- read.csv("http://pastebin.com/raw.php?i=sDzXKFxJ", sep=",", header=T)
dates$Date <- as.Date(dates$Date)

# convert the Date to its numeric equivalent
# Note that Dates are stored as number of days internally,
# hence it is easy to convert back and forth mentally
dates$num <- as.numeric(dates$Date)

bin <- 60 # used for aggregating the data and aligning the labels

p <- ggplot(dates, aes(num, ..count..))
p <- p + geom_histogram(binwidth = bin, colour="white")

# The numeric data is treated as a date,
# breaks are set to an interval equal to the binwidth,
# and a set of labels is generated and adjusted in order to align with bars
p <- p + scale_x_date(breaks = seq(min(dates$num)-20, # change -20 term to taste
                                   max(dates$num), 
                                   bin),
                      labels = date_format("%Y-%b"),
                      limits = c(as.Date("2009-01-01"), 
                                 as.Date("2011-12-01")))

# from here, format at ease
p <- p + theme_bw() + xlab(NULL) + opts(axis.text.x  = theme_text(angle=45,
                                                                  hjust = 1,
                                                                  vjust = 1))
p

版本 1:使用 POSIXct

我尝试了一个解决方案,它可以在ggplot22009 年初和 2011 年底之间完成所有操作,绘制而不进行聚合,并在 x 轴上设置限制。

目标情节 v1

在 ggplot2 中设置限制的绘图

代码 v1

library("ggplot2")
library("scales")

dates <- read.csv("http://pastebin.com/raw.php?i=sDzXKFxJ", sep=",", header=T)
dates$Date <- as.POSIXct(dates$Date)

p <- ggplot(dates, aes(Date, ..count..)) + 
    geom_histogram() +
    theme_bw() + xlab(NULL) +
    scale_x_datetime(breaks = date_breaks("3 months"),
                     labels = date_format("%Y-%b"),
                     limits = c(as.POSIXct("2009-01-01"), 
                                as.POSIXct("2011-12-01")) )

p

当然,它可以通过使用轴上的标签选项来完成,但这是为了在绘图包中使用一个干净的简短例程来完成绘图。

于 2012-05-27T20:18:50.937 回答
6

breaks=我知道这是一个老问题,但对于任何在 2021 年(或以后)提出这个问题的人来说,使用参数 forgeom_histogram()并创建一个小快捷函数来制作所需的序列,这可以更容易地完成。

dates <- read.csv("http://pastebin.com/raw.php?i=sDzXKFxJ", sep=",", header=T)

dates$Date <- lubridate::ymd(dates$Date)

by_month <- function(x,n=1){
  seq(min(x,na.rm=T),max(x,na.rm=T),by=paste0(n," months"))
}

ggplot(dates,aes(Date)) +
  geom_histogram(breaks = by_month(dates$Date)) +
  scale_x_date(labels = scales::date_format("%Y-%b"),
               breaks = by_month(dates$Date,2)) + 
  theme(axis.text.x = element_text(angle=90))

直方图

于 2021-03-11T13:00:52.073 回答
5

我认为关键是您需要在 ggplot 之外进行频率计算。将 aggregate() 与 geom_bar(stat="identity") 一起使用以获得没有重新排序因子的直方图。这是一些示例代码:

require(ggplot2)

# scales goes with ggplot and adds the needed scale* functions
require(scales)

# need the month() function for the extra plot
require(lubridate)

# original data
#df<-read.csv("http://pastebin.com/download.php?i=sDzXKFxJ", header=TRUE)

# simulated data
years=sample(seq(2008,2012),681,replace=TRUE,prob=c(0.0176211453744493,0.302496328928047,0.323054331864905,0.237885462555066,0.118942731277533))
months=sample(seq(1,12),681,replace=TRUE)
my.dates=as.Date(paste(years,months,01,sep="-"))
df=data.frame(YM=strftime(my.dates, format="%Y-%b"),Date=my.dates,Year=years,Month=months)
# end simulated data creation

# sort the list just to make it pretty. It makes no difference in the final results
df=df[do.call(order, df[c("Date")]), ]

# add a dummy column for clarity in processing
df$Count=1

# compute the frequencies ourselves
freqs=aggregate(Count ~ Year + Month, data=df, FUN=length)

# rebuild the Date column so that ggplot works
freqs$Date=as.Date(paste(freqs$Year,freqs$Month,"01",sep="-"))

# I set the breaks for 2 months to reduce clutter
g<-ggplot(data=freqs,aes(x=Date,y=Count))+ geom_bar(stat="identity") + scale_x_date(labels=date_format("%Y-%b"),breaks="2 months") + theme_bw() + opts(axis.text.x = theme_text(angle=90))
print(g)

# don't overwrite the previous graph
dev.new()

# just for grins, here is a faceted view by year
# Add the Month.name factor to have things work. month() keeps the factor levels in order
freqs$Month.name=month(freqs$Date,label=TRUE, abbr=TRUE)
g2<-ggplot(data=freqs,aes(x=Month.name,y=Count))+ geom_bar(stat="identity") + facet_grid(Year~.) + theme_bw()
print(g2)
于 2012-05-27T15:57:14.270 回答
0

标题“基于 Gauden 方法的绘图”下的错误图是由于 binwidth 参数: ... + Geom_histogram (binwidth = 30, color = "white") + ... 如果我们将 30 的值更改为值小于 20,例如 10,您将获得所有频率。

在统计数据中,数值比显示更重要,平淡的图形对非常漂亮的图片但有错误。

于 2015-08-18T22:38:17.277 回答