data(d.pollZH16)
dp <- d.pollZH16
names(dp)
dp$date <- gendateaxis(date=dp$date, hour=dp$hour)
plyx(O3+NO2~date, data=dp, subset= month=="May", type="l")
dp$summer <- dp$month %in% c("Jun","Jul","Aug")
dp$daylight <- dp$hour>8 & dp$hour<17
plmatrix(O3~temp+logst(rad)+logst(prec), data=dp,
subset = summer & daylight)
Run the code above in your browser using DataLab