Question

我的数据看起来像这样：

有10,000行，每行代表一个城市以及自1998-01至2013-9以来的所有月份：

RegionName| State|  Metro|         CountyName|  1998-01|      1998-02|  1998-03

New York|   NY| New York|   Queens|         1.3414|   1.344|             1.3514

Los Angeles|    CA| Los Angeles|    Los Angeles|    12.8841|     12.5466|   12.2737

Philadelphia|   PA| Philadelphia|   Philadelphia|   1.626|    0.5639|   0.2414

Phoenix|            AZ| Phoenix|            Maricopa|    2.7046|       2.5525|  2.3472

我希望自1998年以来可以为任何一个城市或一个以上的城市制作一个地图。

我尝试了这个，但是我收到了一个错误。我不确定我是否尝试过这个权利。任何帮助将不胜感激。谢谢。

forecl <- ts(forecl, start=c(1998, 1), end=c(2013, 9), frequency=12)

plot(forecl)

Error in plots(x = x, y = y, plot.type = plot.type, xy.labels = xy.labels,  : 
  cannot plot more than 10 series as "multiple"

Answer 1

您可以尝试

require(reshape)
require(ggplot2)
forecl <- melt(forecl, id.vars = c("region","state","city"), variable_name = "month")
forecl$month <- as.Date(forecl$month)
ggplot(forecl, aes(x = month, y = value, color = city)) + geom_line()

Answer 2

要添加到@ JLLagrange的答案，如果城市太多且颜色难以区分，您可能希望将city传递给facet_grid()。

ggplot(forecl, aes(x = month, y = value, color = city, group = city)) +
  geom_line() +
  facet_grid( ~ city)

Answer 3

您能否举例说明您的数据，例如： {/ 1}}，之前转换为时间序列对象？问题也可能出在dput(head(forecl))对象上。

无论如何，我认为有两个问题。

首先，数据是宽格式的。我不确定你的专栏名称，因为它们应该以字母开头，但无论如何，一般的想法都是这样的：

ts

其次，我认为你试图同时策划太多的城市。尝试：

test <- structure(list(
  city = structure(1:2, .Label = c("New York", "Philly"), 
  class = "factor"), state = structure(1:2, .Label = c("NY", 
  "PA"), class = "factor"), a2005.1 = c(1, 1), a2005.2 = c(2, 5
  )), .Names = c("city", "state", "a2005.1", "a2005.2"), row.names = c(NA, 
  -2L), class = "data.frame")

test.long <- reshape(test, varying=c(3:4), direction="long")

或

plot(forecl[, 1])

R中的时间序列图

3 个答案: