0

我正在尝试在 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 数据框一起运行的想法将不胜感激。

谢谢,克雷格。

4

1 回答 1

0

原来功能 calendarPlot 必须将“日期”列标记为日期才能使其工作。我将列名从“日期时间”更改为“日期”,效果很好:

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

日历图的输出

于 2016-03-01T06:24:31.353 回答