3

I'm a newer R user and I need help with a time series plot. I created a time series plot, and cannot figure out how to change my x-axis values to correspond to my sample dates. My data is as follows:

 Year  Month  Level
 2009   8     350
 2009   9     210
 2009   10    173
 2009   11    166
 2009   12    153
 2010   1     141
 2010   2     129
 2010   3     124
 2010   4     103
 2010   5     69
 2010   6     51
 2010   7     49
 2010   8     51
 2010   9     51

Let's say this data is saved as the name "data.csv"

data = read.table("data.csv", sep = ",", header = T)
data.ts = ts(data, frequency = 1)
plot(dat.mission.ts[, 3], ylab = "level", main = "main", axes = T)

I've also tried inputing the start = c(2009, 8) into the ts function but I still get wrong values

When I plot this my x axis does not correlate to August 2009 through Sept. 2010. It will either increase by year or just by decimal. I've looked up many examples online and also through the ? help on R, but cannot find a way to relabel my axis values. Any help would be appreciated.

Spacedman
  • 92,590
  • 12
  • 140
  • 224
Karl Campbell
  • 33
  • 1
  • 1
  • 3

3 Answers3

3

Using base coding, you can accomplish this in a few steps. As described in this SO answer, you can identify your "Month" and "Year" data as a date if you use as.Date and paste functions together and incorporate a day (i.e., first day of the month; "1"). For the purposes of this answer, I will simply refer to the data you provided as df:

df$date<-with(df,as.Date(paste(Year,Month,'1',sep='-'),format='%Y-%m-%d'))
df
   Year Month Level       date
1  2009     8   350 2009-08-01
2  2009     9   210 2009-09-01
3  2009    10   173 2009-10-01
4  2009    11   166 2009-11-01
5  2009    12   153 2009-12-01
6  2010     1   141 2010-01-01
7  2010     2   129 2010-02-01
8  2010     3   124 2010-03-01
9  2010     4   103 2010-04-01
10 2010     5    69 2010-05-01
11 2010     6    51 2010-06-01
12 2010     7    49 2010-07-01
13 2010     8    51 2010-08-01
14 2010     9    51 2010-09-01

Then you can use your basic plot, axis, and mtext functions to control how you want to visualize the data and your axes. For instance:

xmin<-min(df$date,na.rm=T);xmax<-max(df$date,na.rm=T) #ESTABLISH X-VALUES (MIN & MAX)
ymin<-min(df$Level,na.rm=T);ymax<-max(df$Level,na.rm=T) #ESTABLISH Y-VALUES (MIN & MAX)

xseq<-seq.Date(xmin,xmax,by='1 month') #CREATE DATE SEQUENCE THAT INCREASES BY MONTH FROM DATE MINIMUM TO MAXIMUM
yseq<-round(seq(0,ymax,by=50),0) # CREATE SEQUENCE FROM 0-350 BY 50

par(mar=c(1,1,0,0),oma=c(6,5,3,2)) #CONTROLS YOUR IMAGE MARGINS
plot(Level~date,data=df,type='b',ylim=c(0,ymax),axes=F,xlab='',ylab='');box() #PLOT LEVEL AS A FUNCTION OF DATE, REMOVE AXES FOR FUTURE CUSTOMIZATION 
axis.Date(side=1,at=xseq,format='%Y-%m',labels=T,las=3) #ADD X-AXIS LABELS WITH "YEAR-MONTH" FORMAT
axis(side=2,at=yseq,las=2) #ADD Y-AXIS LABELS
mtext('Date (Year-Month)',side=1,line=5) #X-AXIS LABEL
mtext('Level',side=2,line=4) #Y-AXIS LABEL

enter image description here

Community
  • 1
  • 1
ccapizzano
  • 1,556
  • 13
  • 20
1
library(data.table)
library(ggplot2)
library(scales)

data<-data.table(datetime=seq(as.POSIXct("2009/08/01",format="%Y/%m/%d"),
                              as.POSIXct("2010/09/01",format="%Y/%m/%d"),by="1 month"),
                 Level=c(350,210,173,166,153,141,129,124,103,69,51,49,51,51))


ggplot(data)+
 geom_point(aes(x=datetime,y=Level),col="brown1",size=1)+
 scale_x_datetime(labels = date_format("%Y/%m"),breaks = "1 month")+
 theme(axis.text.x = element_text(angle = 90, hjust = 1,vjust=0.3))
Soheil
  • 954
  • 7
  • 20
0

Example using xts package:

library(xts)
ts1 <- xts(data$Level, as.POSIXct(sprintf("%d-%d-01", data$Year, data$Month)))

# or ts1 <- xts(data$Level, as.yearmon(data$Year + (data$Month-1)/12))

plot(ts1)

If you are using ggplot2:

library(ggplot2)
autoplot(ts1)
bergant
  • 7,122
  • 1
  • 20
  • 24