使用read.zoo读取时间序列表

时间:2011-01-17 02:59:12

标签: r time-series zoo

我到处都看了,但我以前找不到这个问题的位置。

将这些数据导入适当的动物园系列有什么干净的方法?此版本是复制/粘贴,以使此帖子更容易,但它将始终以下表格式(来自文本文件)。我的read.zoo()语句将Year作为索引读取,但是季度(Qtr1,Qtr2等)被读作列名。我一直试图找出一种非垃圾方式来读取列作为索引的“四分之一”部分,但它很草率(太松散了,不能发布)。我猜这个问题已经解决了,但我找不到了。

> texinp <- "   
+ Year   Qtr1  Qtr2  Qtr3  Qtr4   
+ 1992    566   443   329   341   
+ 1993    344   212   133   112   
+ 1994    252   252   199   207"   
> z <- read.zoo(textConnection(texinp), header=TRUE)
> z  

从as.yearqtr()文档中,目标看起来像:

1992 Q1 1992 Q2 1992 Q3 1992 Q4 1993 Q1 1993 Q2 1993 Q3 1993 Q4    
    566     443     329     341     344     212     133     112      

1994 Q1 1994 Q2 1994 Q3 1994 Q4     
    252     252     199     207    

2 个答案:

答案 0 :(得分:3)

read.zoo假设您的数据最多只有一个时间索引列,因此您必须自己处理。首先使用read.table

阅读
zt <- read.table( textConnection( texinp ), header = TRUE)

然后使用melt包中的reshape函数将其转换为“长表”:

require(reshape)
zt.m <- melt( zt, id = 'Year', variable_name = 'Qtr')

> zt.m
   Year  Qtr value
1  1992 Qtr1   566
2  1993 Qtr1   344
3  1994 Qtr1   252
4  1992 Qtr2   443
5  1993 Qtr2   212
6  1994 Qtr2   252
7  1992 Qtr3   329
8  1993 Qtr3   133
9  1994 Qtr3   199
10 1992 Qtr4   341
11 1993 Qtr4   112
12 1994 Qtr4   207

最后创建所需的zoo对象:

z <- with( zt.m, zoo( value, as.yearqtr(paste(Year, Qtr), format = '%Y Qtr%q')))

> z
1992 Q1 1992 Q2 1992 Q3 1992 Q4 1993 Q1 1993 Q2 1993 Q3 1993 Q4 1994 Q1 1994 Q2 
    566     443     329     341     344     212     133     112     252     252 
1994 Q3 1994 Q4 
    199     207 

答案 1 :(得分:3)

使用read.zoo读取数据,然后将其转换为zooreg时间索引为yearqtr的对象:

texinp <- "Year   Qtr1  Qtr2  Qtr3  Qtr4   
1992    566   443   329   341   
1993    344   212   133   112   
1994    252   252   199   207"

library(zoo)

z <- read.zoo(text = texinp, header=TRUE)
zz <- zooreg(c(t(z)), start = yearqtr(start(z)), freq = 4)

结果如下:

> zz

1992 Q1 1992 Q2 1992 Q3 1992 Q4 1993 Q1 1993 Q2 1993 Q3 1993 Q4 1994 Q1 1994 Q2 1994 Q3 1994 Q4 
    566     443     329     341     344     212     133     112     252     252     199     207