2017-10-18 86 views
2

这是我无法在网上找到任何帮助的问题。在R:如何在使用日期时间列的posixct时删除日期 - 时间序列图中的空白?

这里是链接到原始数据(太大,贴在这里): 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列由每天的数据从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 to 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) 

here is my plot

创建面

非常感谢您提前。

+0

任何机会,你可以使用'dput()'共享数据?如果您提供一个[完整的最小可重现示例](http://stackoverflow.com/help/mcve)来解决您的问题,我们就更有可能帮助您。我们可以从中学习并使用它来向您展示如何回答您的问题。另外,欢迎! –

+0

@EricFail请参阅数据集的谷歌驱动器链接。谢谢。 –

+0

如果有任何更改,您可以创建一个较小的(即_minimal_)数据集,以便可以复制该问题? –

回答

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) 

enter image description here

相关问题