2

I have these data:

> dput(dataF)
structure(list(key_date = c("2016-09-01", "2016-09-01", "2016-09-01", 
"2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", 
"2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", 
"2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", 
"2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", 
"2016-09-01"), date = c("2016-09-01 00:00:00", "2016-09-01 01:00:00", 
"2016-09-01 02:00:00", "2016-09-01 03:00:00", "2016-09-01 04:00:00", 
"2016-09-01 05:00:00", "2016-09-01 06:00:00", "2016-09-01 07:00:00", 
"2016-09-01 08:00:00", "2016-09-01 09:00:00", "2016-09-01 10:00:00", 
"2016-09-01 11:00:00", "2016-09-01 12:00:00", "2016-09-01 13:00:00", 
"2016-09-01 14:00:00", "2016-09-01 15:00:00", "2016-09-01 16:00:00", 
"2016-09-01 17:00:00", "2016-09-01 18:00:00", "2016-09-01 19:00:00", 
"2016-09-01 20:00:00", "2016-09-01 21:00:00", "2016-09-01 22:00:00", 
"2016-09-01 23:00:00"), PM = c("10.9", "9.8", "9.9", "14.1", 
"13.6", "16.5", "15", "18.5", "18", "17", "16.6", "12", "12.1", 
"18.1", "15.9", "15.9", "16.9", "21.6", "23.5", "40.7", "16.6", 
"12.7", "12.4", "12.2"), site = c("LW2", "LW2", "LW2", "LW2", 
"LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", 
"LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", 
"LW2", "LW2"), code = c("LW2", "LW2", "LW2", "LW2", "LW2", "LW2", 
"LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", 
"LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2"
), airport = structure(c(NA_integer_, NA_integer_, NA_integer_, 
NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, 
NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, 
NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, 
NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, 
NA_integer_), .Label = character(0), class = "factor"), ws = c(NA_real_, 
NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
NA_real_, NA_real_), wd = c(NA_real_, NA_real_, NA_real_, NA_real_, 
NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), 
    tempi = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_), humidity = c(NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
    NA_real_, NA_real_)), .Names = c("key_date", "date", "PM", 
"site", "code", "airport", "ws", "wd", "tempi", "humidity"), row.names = c(NA, 
-24L), class = "data.frame")

It seems fine but why do I get this error when I try to plot the data with openair?

timePlot(dataF, pollutant = c("PM"), group = TRUE, avg.time = "default")

Error in object[[name, exact = TRUE]] : subscript out of bounds
In addition: Warning message:
In checkPrep(mydata, vars, type, remove.calm = FALSE) :
  Detected data with Daylight Saving Time, converting to UTC/GMT

What is wrong with the data? Any ideas?

zx8754
  • 52,746
  • 12
  • 114
  • 209
Run
  • 54,938
  • 169
  • 450
  • 748

1 Answers1

3

I suppose you want to plot the PM values against the date ones, using timeplot(). All date-time values in date belong to the same day "2016-09-01", which is specified under key_date, instead.

The timeplot() description https://www.rdocumentation.org/packages/openair/versions/0.3-8/topics/timePlot specifies that data have to be a data.frame, with a date column, which is date in the example under consideration.

Therefore I use the fucntion ymd_hms from the package lubridate for handling date-time format in date and I force the values in PM to be numeric with as.numeric arriving at

require(lubridate)
library(lubridate)
dataF <-   data.frame(key_date = c("2016-09-01", "2016-09-01", "2016-09-01", 
                            "2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", 
                            "2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", 
                            "2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", 
                            "2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", "2016-09-01", 
                            "2016-09-01"), 
                    date =  ymd_hms( c("2016-09-01 00:00:00", "2016-09-01 01:00:00",
                            "2016-09-01 02:00:00", "2016-09-01 03:00:00", "2016-09-01 04:00:00",
                            "2016-09-01 05:00:00", "2016-09-01 06:00:00", "2016-09-01 07:00:00",
                            "2016-09-01 08:00:00", "2016-09-01 09:00:00", "2016-09-01 10:00:00",
                            "2016-09-01 11:00:00", "2016-09-01 12:00:00", "2016-09-01 13:00:00",
                            "2016-09-01 14:00:00", "2016-09-01 15:00:00", "2016-09-01 16:00:00",
                            "2016-09-01 17:00:00", "2016-09-01 18:00:00", "2016-09-01 19:00:00",
                            "2016-09-01 20:00:00", "2016-09-01 21:00:00", "2016-09-01 22:00:00",
                            "2016-09-01 23:00:00")), 
               PM = as.numeric(c("10.9", "9.8", "9.9", "14.1", 
                     "13.6", "16.5", "15", "18.5", "18", "17", "16.6", "12", "12.1", 
                     "18.1", "15.9", "15.9", "16.9", "21.6", "23.5", "40.7", "16.6", 
                     "12.7", "12.4", "12.2")), 
               site = c("LW2", "LW2", "LW2", "LW2", 
                       "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", 
                       "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", 
                       "LW2", "LW2"), 
               code = c("LW2", "LW2", "LW2", "LW2", "LW2", "LW2", 
                       "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", 
                       "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2", "LW2"),             
         airport = c(NA_integer_, NA_integer_, NA_integer_, 
                        NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, 
                        NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, 
                        NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, 
                        NA_integer_, NA_integer_, NA_integer_, NA_integer_, NA_integer_, 
                        NA_integer_), 
          ws = c(NA_real_, 
                 NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
                 NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
                 NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
                 NA_real_, NA_real_), 
          wd = c(NA_real_, NA_real_, NA_real_, NA_real_, 
                 NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
                 NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
                 NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_), 
           tempi = c(NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
                         NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
                         NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
                         NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
                         NA_real_), 
          humidity = c(NA_real_, NA_real_, NA_real_, NA_real_, 
                       NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
                       NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
                       NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, NA_real_, 
                       NA_real_, NA_real_)
        )

We check the structure of dataF:

str(dataF)
'data.frame':   24 obs. of  10 variables:
 $ key_date: Factor w/ 1 level "2016-09-01": 1 1 1 1 1 1 1 1 1 1 ...
 $ date    : POSIXct, format: "2016-09-01 00:00:00" "2016-09-01 01:00:00" "2016-09-01 02:00:00" ...
 $ PM      : num  10.9 9.8 9.9 14.1 13.6 16.5 15 18.5 18 17 ...
 $ site    : Factor w/ 1 level "LW2": 1 1 1 1 1 1 1 1 1 1 ...
 $ code    : Factor w/ 1 level "LW2": 1 1 1 1 1 1 1 1 1 1 ...
 $ airport : Factor w/ 0 levels: NA NA NA NA NA NA NA NA NA NA ...
 $ ws      : num  NA NA NA NA NA NA NA NA NA NA ...
 $ wd      : num  NA NA NA NA NA NA NA NA NA NA ...
 $ tempi   : num  NA NA NA NA NA NA NA NA NA NA ...
 $ humidity: num  NA NA NA NA NA NA NA NA NA NA ...

Now we plot

timePlot(dataF, pollutant = "PM", group = TRUE, avg.time = "default", lty = 1)
dev.off()
Avitus
  • 734
  • 2
  • 14
  • 27
  • thanks for answer. I thought i have converted the PM data to numberic with `dataF[, columns] <- lapply(dataF[, columns, drop = FALSE], function(x) as.numeric(as.character(x)))` – Run May 30 '17 at 03:05
  • 1
    @teelou thank you for your comment. Please upvote and flag the answer if you believe it can be useful – Avitus May 30 '17 at 05:15