21

我有一个包含多个时间序列返回的数据框，存储在列中。

第一列包含日期，随后的列是独立的时间序列，每个都有一个名称。列标题是变量名。

## I have a data frame like this
t <- seq(as.Date('2009-01-01'),by='days',length=10)
X <- rnorm(10,0,1)
Y <- rnorm(10,0,2)
Z <- rnorm(10,0,4)

dat <- data.frame(t,X,Y,Z)

## which appears as
           t          X          Y         Z
1 2009-01-01 -1.8763317 -0.1885183 -6.655663
2 2009-01-02 -1.3566227 -2.1851226 -3.863576
3 2009-01-03 -1.3447188  2.4180249 -1.543931

我想将每个时间序列绘制为一个单独的图上的一条线，在一个格子中，每个图都由变量名称标记。

要使用 lattice 绘制此图，数据必须采用 tall 格式，如下所示： 

           t symbol       price
1 2009-01-01      X -1.8763317
2 2009-01-02      Y -0.1885183
2 2009-01-02      Z -6.655663

什么是一个好的函数调用来做到这一点？

4

5 回答 5

18

您还可以使用“reshape”库中的 melt() （我认为它比 reshape() 本身更容易使用） - 这将节省您必须重新添加时间列的额外步骤......

> library(reshape)
> m <- melt(dat,id="t",variable_name="symbol")
> names(m) <- sub("value","price",names(m))
> head(m)
           t symbol       price
1 2009-01-01      X -1.14945096
2 2009-01-02      X -0.07619870
3 2009-01-03      X  0.01547395
4 2009-01-04      X -0.31493143
5 2009-01-05      X  1.26985167
6 2009-01-06      X  1.31492397
> class(m$t)
[1] "Date"
> library(lattice)                                                              
> xyplot( price ~ t | symbol, data=m ,type ="l", layout = c(1,3) )

然而，对于这个特定的任务，我会考虑使用“动物园”库，它不需要你重塑数据框：

> library(zoo)                                                                  
> zobj <- zoo(dat[,-1],dat[,1])                                                 
> plot(zobj,col=rainbow(ncol(zobj)))

R 开发人员/贡献者（在本例中是 Gabor 和 Hadley）为我们提供了许多很棒的选择。（并且不能忘记格子包的 Deepayan）

于 2009-07-28T06:29:20.270 回答
17

从tidyr 收集帮助页面：

例子

library(tidyr)
library(dplyr)
# From http://stackoverflow.com/questions/1181060
stocks <- data.frame(
  time = as.Date('2009-01-01') + 0:9,
  X = rnorm(10, 0, 1),
  Y = rnorm(10, 0, 2),
  Z = rnorm(10, 0, 4)
)

gather(stocks, stock, price, -time)
stocks %>% gather(stock, price, -time)
于 2014-10-08T10:13:36.100 回答
11

如果是多变量时间序列，可以考虑使用同名包将其存储为zoo对象。这使得索引、合并、子集化变得更加容易——参见动物园小插曲。

但是正如你所问的格子图 - 这也可以做到。在本例中，我们构造了一个简单的“long”data.frame，其中包含一个日期列、一个值列“val”和一个变量 id 列“var”：

> set.seed(42)
> D <- data.frame(date=rep(seq(as.Date("2009-01-01"),Sys.Date(),by="week"),2),\
                  val=c(cumsum(rnorm(30)), cumsum(rnorm(30))), \
                  var=c(rep("x1",30), rep("x2",30)))

鉴于该数据集，根据您的描述进行绘图是由 lattice 包中的 xyplot 完成的，方法是要求绘制“按变量分组的值给定数据”的图，我们在每个面板中打开线：

> library(lattice)
> xyplot(val ~ date | var, data=D, panel=panel.lines)
于 2009-07-25T03:07:18.407 回答
5

对于第一列中的日期和其他每一列中的值的数据框“temp”：

> par(mfrow=c(3,4)) # 3x4 grid of plots
> mapply(plot,temp[,-1],main=names(temp)[-1],MoreArgs=list(x=temp[,1],xlab="Date",type="l",ylab="Value") )
于 2009-07-25T03:18:43.100 回答
3

非常感谢大家的回答 - 德克的回答是正确的。

原来缺少的步骤是使用“stack()”函数将数据帧从宽格式转换为长格式。我知道使用 reshape() 函数可能有一种更简单的方法，如果有人想发布它，很高兴看到一个例子。

所以这就是我最终做的事情，使用问题中提到的“dat”数据框：

## use stack() to reshape the data frame to a long format
## <time> <stock> <price>
stackdat <- stack(dat,select=-t) 
names(stackdat) <- c('price','symbol')

## create a column of date & bind to the new data frame
nsymbol <- length(levels(stackdat$symbol))  
date <- rep(dat$t, nsymbol)  
newdat <- cbind(date,stackdat)

## plot it with lattice
library(lattice)
xyplot(price ~ date | symbol,  ## model conditions on 'symbol' to lattice
       data=newdat,            ## data source
       type='l',               ## line
       layout=c(nsymbol,1))    ## put it on a single line

## or plot it with ggplot2
library(ggplot2)
qplot(date, price, data = newdat, geom="line") + facet_grid(. ~ symbol)
于 2009-07-25T22:51:54.827 回答