Multiple time-series data with missing values in a single plot using R

Question

I have prepared a data frame ordered according to date and time as.POSIXlt:

             Date_Time A1    B2 C1    E2
46 24/06/2012 12:20:00 NA 5.515 20    NA
47 24/06/2012 13:20:00 41    NA NA 3.519
48 25/06/2012 14:00:00 NA    NA NA    NA
49 25/06/2012 14:20:00 30    NA 30    NA
50 27/06/2012 15:20:00 NA    71 NA    NA
51 28/06/2012 18:00:00 11    NA 55    11
... ...

As you can see it contains a lot of 'NA' values. Is it possible to plot all data in a single plot with the 'x-axis' showing month-year and different colors assigned for different data (e.g. green for A1, blue for C1)? I tried 'ggplot2' package which removes all missing values. Trying 'zoo' package doesn't allow x-axis to be plotted as month-year. Any suggestion?

thelatemail · Accepted Answer · 2013-08-23 04:09:39Z

2

This should not be problematic for any of R's plotting packages. I'm not a ggplot user, but here's a base R example using a slightly modified version of your data. Note that using as.POSIXct is preferable to using as.POSIXlt within data.frame objects:

dat <- read.csv(text="Date_Time,A1,B2,C1,E2
24/06/2012 12:20:00,NA,5.515,20,NA
24/06/2012 13:20:00,41,NA,NA,3.519
25/06/2012 14:00:00,NA,NA,NA,NA
25/07/2012 14:20:00,30,NA,30,NA
27/08/2012 15:20:00,NA,71,NA,NA
28/09/2012 18:00:00,11,NA,55,11")

dat$Date_Time <- as.POSIXct(dat$Date_Time,format="%d/%m/%Y %H:%M:%S")

matplot(
  dat[,1],
  dat[2:5],
  type="o",
  pch=19,
  col=c("green","red","blue","black"),
  xaxt="n",
  xlab="Year-Month",
  ylab="Value"
)

Add a POSIXct represented axis with a chosen format:

axis.POSIXct(side=1,x=dat$Date_Time,format="%Y-%m")

enter image description here

answered Aug 23, 2013 at 4:09

thelatemail

94.3k12 gold badges140 silver badges197 bronze badges

Sign up to request clarification or add additional context in comments.

3 Comments

ToNoY Over a year ago

thanks, but say you have time-series data for more than 50 variables, how would you give the legend?

thelatemail Over a year ago

@ToNoY - well, if you were doing this with 50 variables, that would be totally unreadable and a legend wouldn't help at all. Adding a legend is pretty straight-forward otherwise though - see ?legend and the examples at the end of the help page.

ToNoY Over a year ago

There's another problem, due to NA values, I can't plot it as a line which is very important to see the trend when you have at least a number of different variables. Even if I use 'matlines' after plot with the same data, it doesn't seem to work! Any solution??

Collectives™ on Stack Overflow

Multiple time-series data with missing values in a single plot using R

1 Answer 1

3 Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

3 Comments

Your Answer

Sign up or log in

Post as a guest

Related