[英]time series plot with x axis in “year”-“month” in R
我有一个“月度数据”我想绘制数据,以便在x轴上以格式%Y-%m" (2001-01)
得到一个日期。假设我的数据是pcp <- rnorm(24)
。
我试过了:
PCP <- ts(pcp, frequency = 12, start = 2001)
plot(PCP)
但是该图在x轴上只有几年。
如何在x轴上获得我想要的日期格式的图?
这是玩具数据的一个想法,因为这个问题不可重复。 希望它有所帮助
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))
我发现优秀的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
格式的日期。
最好的方法是使用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"
问候,
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.