元々、ws(風速)とwd(風速)のデータが抜けていて、
timestamp no pm hum
1 1404227061 49 2.5 54
2 1404227121 35 4.5 52
3 1404227181 27 6.6 49
4 1404227241 8 7.1 48
5 1404227301 51 8.4 46
6 1404227361 20 9.0 44
そのギャップを毎日の要約データで埋めなければなりません。
dat <- transform(dat, wd = 230, ws = 2.160000 )
timestamp no pm hum wd ws
1 1404227061 49 2.5 54 230 2.16
2 1404227121 35 4.5 52 230 2.16
3 1404227181 27 6.6 49 230 2.16
4 1404227241 8 7.1 48 230 2.16
5 1404227301 51 8.4 46 230 2.16
6 1404227361 20 9.0 44 230 2.16
次に、この関数を openair パッケージで実行します。
polarPlot(dat, pollutant = "pm", na.rm = TRUE)
結果、
Error in cut.default(mydata[[x]], breaks = seq(0, max.ws, length = 31), :
'breaks' are not unique
どういう意味ですか?これを修正するにはどうすればよいですか? それとも、私のデータに何か問題がありますか?
編集:
ws の値を変更し、
dat[4, "ws"] = 5.6
結果、
Warning message:
Not enough data to fit surface.
Try reducing the value of the smoothing parameter, k to less than 100.