时间序列图用x轴表示“年” - “月”表示R

Question

我有一个“月度数据”我想绘制数据，以便在x轴上以%Y-%m" (2001-01)格式获得日期。说我的数据是pcp <- rnorm(24)。

我试过了：

PCP <- ts(pcp, frequency = 12, start = 2001)
plot(PCP) 

但是该图在x轴上只有几年。

如何在x轴上获得我想要的日期格式的图？

Answer 1

这是玩具数据的一个想法，因为这个问题不可重复。希望它有所帮助

R> foo = ts(rnorm(36), frequency = 12, start = 2001)
R> plot(foo, xaxt = "n")
R> tsp = attributes(foo)$tsp
R> dates = seq(as.Date("2001-01-01"), by = "month", along = foo)
R> axis(1, at = seq(tsp[1], tsp[2], along = foo), labels = format(dates, "%Y-%m"))

ggplot版本，其数据与您的数据相似

R> df = data.frame(date = seq(as.POSIXct("2001-01-01"), by = "month", length.out = 36), pcp = rnorm(36))
R> library(ggplot2)
R> library(scales)
R> p = ggplot(data = df, aes(x = date, y = pcp)) + geom_line()
R> p + scale_x_datetime(labels = date_format("%Y-%m"), breaks = date_breaks("months")) + theme(axis.text.x = element_text(angle = 45))

Answer 2

我发现优秀的xts包是存储数据的最佳方式。如果您没有，可以使用install.packages('xts')下载。

让我们从基础开始 - 包括制作pcp，因为你没有提供它。

require(xts)
pcp <- rnorm(24)
PCP <- ts(pcp, frequency = 12, start = 2001)
plot(as.xts(PCP), major.format = "%Y-%m")

这为您提供了如下图表。您可以通过更改传递给major.format的字符串来调整日期。例如，"%b-%y"会以2001年1月的Jan-01格式生成日期。

Answer 3

最好的方法是使用axis.POSIXct {graphics} 这个例子来自这个函数的帮助：

with(beaver1, {
time <- strptime(paste(1990, day, time %/% 100, time %% 100),
                 "%Y %j %H %M")
plot(time, temp, type = "l") # axis at 4-hour intervals.
# now label every hour on the time axis
plot(time, temp, type = "l", xaxt = "n")
r <- as.POSIXct(round(range(time), "hours"))
axis.POSIXct(1, at = seq(r[1], r[2], by = "hour"), format = "%H")
})

在您的情况下，将格式更改为format="%Y-%m"

此致