我正在尝试在 R 的 openair 包中使用 calendarPlot 函数。如果我加载包附带的 mydata 数据框并运行以下代码,它可以正常工作:
calendarPlot(mydata, pollutant = "co", year = 2003)
mydata 数据框如下所示:
'data.frame': 65533 obs. of 10 variables:
$ date: POSIXt, format: "1998-01-01 00:00:00" "1998-01-01 01:00:00" "1998- 01-01 02:00:00" "1998-01-01 03:00:00" ...
$ ws : num 0.6 2.16 2.76 2.16 2.4 3 3 3 3.36 3.96 ...
$ wd : int 280 230 190 170 180 190 140 170 170 170 ...
$ nox : int 285 NA NA 493 468 264 171 195 137 113 ...
$ no2 : int 39 NA NA 52 78 42 38 51 42 39 ...
$ o3 : int 1 NA 3 3 2 0 0 0 1 2 ...
$ pm10: int 29 37 34 35 34 16 11 12 12 12 ...
$ so2 : num 4.72 NA 6.83 7.66 8.07 ...
$ co : num 3.37 NA 9.6 10.22 8.91 ...
$ pm25: int NA NA NA NA NA NA NA NA NA NA ...
我有自己的数据框 R1_Temp_cut,如下所示:
'data.frame': 722880 obs. of 3 variables:
$ DateTime: POSIXct, format: "2013-10-01 00:01:00" "2013-10-01 00:02:00" "2013-10-01 00:03:00" "2013-10-01 00:04:00" ...
$ Air.Temp: num 17.5 17.4 17.4 17.7 17.7 17.7 17.7 17.6 17.6 17 ...
$ AirTemp : num 17.5 17.4 17.4 17.7 17.7 17.7 17.7 17.6 17.6 17 ...
当我对我的数据运行相同的代码时,我收到以下错误:
> calendarPlot(R1_Temp_cut, pollutant = "AirTemp", year = 2014)
Error in calendarPlot(R1_Temp_cut, pollutant = "AirTemp", year = 2014) :
No data to plot - check year chosen
In addition: Warning message:
In match(x, table, nomatch = 0L) : NAs introduced by coercion
我注意到 mydata 中的日期格式是 POSIXt,而 R1_Temp_cut 是 POSIXct。我将格式更改为 POSIXct,它仍然可以正常工作。帮助文件说日期格式应该是 POSIXct。
数据框 R1_Temp_cut 总结:
> summary(R1_Temp_cut)
DateTime Air.Temp AirTemp
Min. :2013-10-01 00:01:00 Min. :-4.10 Min. : 0.00
1st Qu.:2014-02-03 11:46:45 1st Qu.:14.40 1st Qu.:14.40
Median :2014-06-09 05:15:30 Median :18.20 Median :18.20
Mean :2014-06-09 07:11:51 Mean :18.57 Mean :18.57
3rd Qu.:2014-10-13 01:41:15 3rd Qu.:22.00 3rd Qu.:22.00
Max. :2015-02-15 23:59:00 Max. :99.00 Max. :99.00
任何关于如何让脚本与 R1_Temp_cut 数据框一起运行的想法将不胜感激。
谢谢,克雷格。