r - R中时间序列数据的箱线图

Question

这是我的示例数据，我一直到 2018 年都有月度数据，不同的 state_ID 与不同的 city_ID 相关联。

structure(list(state_ID = c("00001", "00001", "00001", "00001", 
"00001", "00001", "00001", "00001", "00001", "00001", "00001", 
"00001", "00001", "00001", "00001", "00001", "00001", "00001", 
"00001", "00001", "00001", "00001", "00001", "00001", "00001", 
"00001", "00001", "00001", "00001", "00001", "00001", "00001", 
"00001", "00001", "00001", "00001", "00001", "00001", "00001", 
"00001", "00001", "00001", "00001", "00001", "00001", "00001", 
"00002", "00002", "00002", "00002", "00002", "00002", "00002", 
"00002", "00002", "00002", "00002", "00002", "00002", "00002", 
"00002", "00002", "00002", "00002", "00002", "00002", "00002", 
"00002", "00002", "00002", "00002", "00002", "00002", "00002", 
"00002", "00002", "00002", "00002", "00002", "00002", "00002", 
"00002", "00002", "00002", "00002", "00002", "00002", "00002", 
"00002", "00002", "00002", "00002"), city_ID = c(436L, 436L, 
436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 
436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 
436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 
436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 
436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 
436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 
436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 
436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 436L, 
436L, 436L), DATE_ = structure(c(1.404e+09, 1406419200, 1409443200, 
1411862400, 1414281600, 1417305600, 1419724800, 1422144000, 1424563200, 
1427587200, 1430006400, 1433030400, 1435449600, 1437868800, 1440892800, 
1443312000, 1445731200, 1448755200, 1451174400, 1454198400, 1456617600, 
1459036800, 1461456000, 1464480000, 1466899200, 1469923200, 1472342400, 
1474761600, 1477785600, 1480204800, 1482624000, 1485648000, 1488067200, 
1490486400, 1493510400, 1495929600, 1498348800, 1501372800, 1503792000, 
1506211200, 1509235200, 1511654400, 1514678400, 1517097600, 1519516800, 
1521331200, 1.404e+09, 1406419200, 1409443200, 1411862400, 1414281600, 
1417305600, 1419724800, 1422144000, 1424563200, 1427587200, 1430006400, 
1433030400, 1435449600, 1437868800, 1440892800, 1443312000, 1445731200, 
1448755200, 1451174400, 1454198400, 1456617600, 1459036800, 1461456000, 
1464480000, 1466899200, 1469923200, 1472342400, 1474761600, 1477785600, 
1480204800, 1482624000, 1485648000, 1488067200, 1490486400, 1493510400, 
1495929600, 1498348800, 1501372800, 1503792000, 1506211200, 1509235200, 
1511654400, 1514678400, 1517097600, 1519516800, 1521331200), class = c("POSIXct", 
"POSIXt"), tzone = "UTC"), value = c(175.81640448, 175.63913856, 
176.406913536, 189.90377088, 203.23923648, 232.85606087743, 257.99215104, 
283.66071936, 271.38357504, 267.250732032, 268.928352, 275.858909184, 
296.18742528, 296.67188352, 321.60720384, 323.3726592, 304.4304192, 
317.916162048, 339.82880256, 353.526595584, 334.11075456, 374.4671808, 
384.64440576, 374.4813312, 386.65659264, 309.52533504, 338.709529860565, 
335.95442304, 291.20721344, 189.42136016, 345.55611264, 284.379576832, 
364.26294144, 348.0561024, 356.95953408, 381.45053184, 399.6471744, 
411.383721984, 341.01576384, 169.67886144, 166.71692544, 229.002507780565, 
210.709393009355, 129.83239632, 186.77904096, 199.03727138742, 
217.75355904, 202.74088512, 214.35128832, 206.65861632, 207.05547072, 
230.90021831743, 260.32053504, 280.698912, 278.942976, 292.585402368, 
295.38779904, 309.370555392, 309.17826432, 322.06284672, 317.232003072, 
348.15669888, 328.85735424, 440.925229056, 452.2647936, 505.882254336, 
512.65638528, 533.06049024, 563.41078272, 563.734646784, 555.21101184, 
459.154031616, 511.35196690356, 518.53857792, 237.594366592, 
265.52476272, 26.70911584, 570.316075008, 560.95298688, 531.25721472, 
540.957133824, 534.50769024, 579.9181248, 591.968553984, 318.88770064, 
244.57011072, 229.676840448, 301.25467474356, 314.062592568545, 
152.38488576, 114.50917472, 170.496469876475)), row.names = c(NA, 
-92L), class = c("grouped_df", "tbl_df", "tbl", "data.frame"), vars = c("state_ID", 
"city_ID"), drop = TRUE, indices = list(0:45, 46:91), group_sizes = c(46L, 
46L), biggest_group_size = 46L, labels = structure(list(state_ID = c("00001", 
"00002"), city_ID = c(436L, 436L)), row.names = c(NA, -2L), class = "data.frame", vars = c("state_ID", 
"city_ID"), drop = TRUE))

想获取每个组的箱线图（按 state_ID 分组）并检测异常值
有人可以帮我在 R 中做吗

标签： r