Question

使用strptime将日期/时间字符串转换为POSIXlt后，我留下以下内容（数据在此处被截断）：

         DateTime  North  South   West   East  Seast  System
1 2008-09-12 01:00:00 1919.9 3721.4 2085.9 2565.5 2571.1 12863.8

2 2008-09-12 02:00:00 1827.0 3518.1 1965.3 2396.9 2410.7 12118.0

3 2008-09-12 03:00:00 1755.4 3388.4 1866.8 2338.7 2335.2 11684.5

4 2008-09-12 04:00:00 1733.5 3327.1 1810.0 2295.6 2290.2 11456.4

5 2008-09-12 05:00:00 1742.7 3327.3 1831.4 2314.2 2302.3 11517.9

6 2008-09-12 06:00:00 1912.2 3504.4 1986.7 2515.0 2502.6 12420.9

然后，我使用以下代码片段将数据（看似正确）汇总为年度平均值： North_Monthly_Avg <- aggregate(North, list(Date=format(DateTime, "%Y-%m")),mean) 产生以下结果：

 Date        x
1 2008-09 2192.066

2 2008-10 1885.074

3 2008-11 1675.373

4 2008-12 1637.231

5 2009-01 1752.693

6 2009-02 1743.393

我可以绘制'x'值，但不能在x轴上正确标记年份 - 因为它只绘制索引。不知道我错过了什么......我玩过axis.POSIXct，但没有运气。

Answer 1

尝试zoo和lattice：

library(zoo)
library(lattice)

dat <- 'Date  Time  North  South   West   East  Seast  System
2008-09-12 01:00:00 1919.9 3721.4 2085.9 2565.5 2571.1 12863.8
2008-09-12 02:00:00 1827.0 3518.1 1965.3 2396.9 2410.7 12118.0
2008-09-12 03:00:00 1755.4 3388.4 1866.8 2338.7 2335.2 11684.5
2008-09-12 04:00:00 1733.5 3327.1 1810.0 2295.6 2290.2 11456.4
2008-09-12 05:00:00 1742.7 3327.3 1831.4 2314.2 2302.3 11517.9
2008-09-12 06:00:00 1912.2 3504.4 1986.7 2515.0 2502.6 12420.9'

z <- read.zoo(text = dat,  header = TRUE,  index.column = 1:2, tz = "")
xyplot(z)

zAgg <- aggregate(z$North,  by = as.yearmon,  FUN = mean)

dat2 <- 'Date        x
2008-09 2192.066
2008-10 1885.074
2008-11 1675.373
2008-12 1637.231
2009-01 1752.693
2009-02 1743.393'

zAgg <- read.zoo(text = dat2, header = TRUE,  FUN = as.yearmon)

plot(zAgg, xaxt = "n")

tt <- time(zAgg)
m <- format(tt, "%m")
axis(side = 1, at = tt, labels = ifelse(m == "01", trunc(tt), m), cex.axis = .7)

Answer 2

尝试在日期

上使用 as.integer（）

North_Monthly_Avg <- aggregate(North, list(Date=as.integer(format(DateTime, "%Y-%m"))),mean)

Answer 3

@ user1062431，要将刻度名称编辑为首选格式，请在Oscar答案中编辑m＆lt; - format（tt，“％m”）行。

要获得格式12 - 2008，您需要修改：

m＆lt; - format（tt，“％m”）到m＆lt; - format（tt，“％m - ％Y”）

要获得2008年12月的格式，您需要修改：

m＆lt; - format（tt，“％m”）到m＆lt; - format（tt，“％b％Y”）

Answer 4

我认为问题在于没有约会。您必须在当月的第1天或每月的15日结算并将其应用于汇总表。

我想出了这个：

North_Monthly_Avg=aggregate(North,by=list(format(DateTime,'%Y-%m')),mean)
names(North_Monthly_Avg)=c('Month','North')
North_Monthly_Avg$day=15
North_Monthly_Avg$Date=paste(North_Monthly_Avg$Month,North_Monthly_Avg$day,sep='-')
North_Monthly_Avg$Date=strptime(North_Monthly_Avg$Date,'%Y-%m-%d')
plot(m$Date,m$North,xaxt='n') # the xaxt='n' removes any ticks on the x axis
axis(1,as.numeric(m$Date),labels=format(m$Date,'%Y-%m')) # formats the x axis to your liking

我对R来说相当新，所以这可能不是最优雅的解决方案，但它会起作用。

如果您更喜欢本月的第一天，请将{15}替换为$day行中的1，并将sep中的paste更改为'-0'。

Answer 5

您可以尝试使用该套件openair并使用它的功能timeAverage

每小时到每月

library(openair)  
mydata$date <- as.POSIXct(strptime(mydata$date, format = "%d/%m/%Y %H:%M", tz = "GMT"))
hourly<-timeAverage(mydata, average.time = "day")

R：计算每小时数据的月平均值，然后绘制

5 个答案: