Openair: calendarPlot没有读取数据

Question

我试图在R中的openair包中使用calendarPlot函数，如果加载包附带的mydata dataframe并运行以下代码，它工作得很好：

calendarPlot(mydata, pollutant = "co", year = 2003)

mydata dataframe如下所示：

'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 ...

我有自己的dataframe，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

我注意到我的数据中的日期格式是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数据文件一起运行的想法都将不胜感激。

谢谢克雷格。

Answer 1

事实证明，函数calendarPlot必须将“date”列标记为date，才能使其正常工作。我把专栏名从“DateTime”改为“date”，效果非常好：

names(R1_Temp_cut)[1] <- "date" 
calendarPlot(R1_Temp_cut, pollutant = "Air.Temp", year =2014)

来自calendarPlot的输出