这是一个我在网上找不到任何帮助的问题。
以下是原始数据的链接(太大而无法粘贴): https://drive.google.com/open?id=0B5KaRU02aEnAbnM5ODNPMkYxN0E
这是它的结构
> str(dfsite)
'data.frame': 11712 obs. of 8 variables:
$ Site_ID : chr "02AB08-001" "02AB08-001" "02AB08-001" "02AB08-001" ...
$ Deployment_ID: num 15352 15352 15352 15352 15352 ...
$ RecDate : chr "01-07-2007" "01-07-2007" "01-07-2007" "01-07-2007" ...
$ RecTime : chr "00:00:00" "00:30:00" "01:00:00" "01:30:00" ...
$ Temperature : num 19.5 19.5 19.6 19.6 19.6 ...
$ TimeStep : num 30 30 30 30 30 30 30 30 30 30 ...
$ RecYr : Factor w/ 19 levels "1996","1997",..: 12 12 12 12 12 12 12 12 12 12 ...
$ RecDT : POSIXct, format: "2007-07-01 00:00:00" "2007-07-01 00:30:00" "2007-07-01 01:00:00" "2007-07-01 01:30:00" ...
RecDT栏包含7月份整个月00:00:00至23:30:00的每日数据: 1)2007-07-01 00:00:00-23:30:00至2007-07-30 00:00:00-23:30:00 2)2009-07-01 00:00:00-23:30:00至2009-07-30 00:00:00-23:30:00 ...
我想为每年创建一个时间序列图,显示特定site_id的7月1日到30日。
我当前的代码创建了一个时间序列图 - 所有的月份,但它应该只有七月 - 这些年来,但它应该创造不同年份的方面
# original data
df <- excel file
# extract year only
df["RecYr"] <- substr(df$RecDate,7,10)
df$RecYr <- as.factor(df$RecYr)
#extract date and time
df["RecDT"] <- paste(df$RecDate,df$RecTime)
df$RecDT <- as.character(df$RecDT)
df$RecDT <- as.POSIXct(df$RecDT,"%d-%m-%Y %H:%M:%S", tz="UTC")
# select specific site
dfsite <- df[df$Site_ID == "02AB08-001",]
# import ggplot2 and create plot
library(ggplot2)
p <- ggplot(dfsite, aes(x=dfsite$RecDT, y=dfsite$Temperature)) +
geom_point() +
print(p)
非常感谢你。
答案 0 :(得分:0)
一个简单的解决方案可能是:
x <- cumsum(table(dfsite$RecYr))
x <- cbind(c(1,x[-length(x)]),x)
xseq <- c(apply(x,1,function(kk) round(seq(kk[1],kk[2],length.out=3))))
p <- ggplot(dfsite, aes(x=1:nrow(dfsite), y=Temperature)) +
geom_point() + facet_grid(.~RecYr, scale="free_x") +
scale_x_continuous("Time", breaks=xseq,
labels = dfsite$RecDT[xseq])+
theme(axis.text.x=element_text(angle=45,hjust=1))
print(p)