12

我在这里有点困惑，一定缺少一些非常基本的东西。我想从多个时间序列对象中提取列。也许可以通过将 ts 对象放入数据框然后提取它们来完成，但是有一种直接的方法可以对时间序列进行子集化，如本问题链接中所给出

要弄清楚如何为 ts 对象分配名称， ?ts 显示

ts(data = NA, start = 1, end = numeric(), frequency = 1,
   deltat = 1, ts.eps = getOption("ts.eps"), class = , names = )

但是在 ?ts 帮助中运行代码时，无论有无名称参数，都会给出 names(z)NULL。

z <- ts(matrix(rnorm(300), 100, 3), start=c(1961, 1), frequency=12) 或者z <- ts(matrix(rnorm(300), 100, 3), start=c(1961, 1), frequency=12, names=c("x1", "x2", "x3"))

主要问题：1.如何为时间序列对象 ts 和 mts 中的列分配名称？2. 随时间索引直接从时间序列对象中提取列的方法有哪些？是否有必要将其转换为 zoo 或 xts 类？还是单独添加时间索引？

为了给出我正在尝试解决的问题的想法：

# using inbuilt ldeaths time series dataset
ldeaths
d <- diff(ldeaths)
percen <- quantile(d, 0.9)
i <- ifelse(d>percen, 1,0)
signal <- cbind(d,i)

现在要提取指标为 1 的数据集，使用 time index，我不知道如何继续。str(signal) 是一个 mts 对象，但打印信号显示没有时间索引。

非常感谢。

4

2 回答 2

14

提取"mts"对象中的系列名称

你想要colnames()：

> colnames(z)
[1] "x1" "x2" "x3"

这是因为z实际上是一个带有额外属性的矩阵，而矩阵没有，names但它们确实有colnames。

分配/更改"mts"对象中系列的名称

colnames要事后分配或更改它们，请使用替换功能'colnames<-'

> colnames(z) <- paste0("a", 1:3)
> colnames(z)
[1] "a1" "a2" "a3"

"mts"从对象中提取特定系列

至于提取列，[适用于"ts"和"mts"对象。例如：

> z[,1]
             Jan         Feb         Mar         Apr         May
1961  0.81800833 -0.30852155  0.05915071  0.14937058  0.67734362
1962  1.12993606 -0.81176485 -0.51903387  1.12527537 -0.34377553
1963  1.30469813  0.32486340  0.01029512 -1.13631688 -1.22013150
1964  0.72449621 -0.88704234  0.78834391 -0.92956537 -0.31584252
1965  0.24610412  0.97980266  0.17136276  2.45216318  0.15846038
1966 -0.48891587 -0.62820331  0.33190472  2.14094813  1.32389152
1967  0.49120472 -0.10149521 -0.39070688 -0.78743955 -1.20563040
1968 -0.70749150  0.52333087 -0.51991721  0.02037504 -0.59848254
1969 -0.80156968 -1.38172513  0.09400527  0.66966443            
....

OP问题的胆量

最后一点，我不确定你希望得到什么。一个"ts"或对象是一个规则的时间序列，并在一个向量而不是一个时间序列中"mts"提取结果的位。signal[, "d"]

> signal[signal[,2] == 1, 1]
[1]  761 1104  810  653  522  956  593

没有任何时间索引，因为这不再是一个对象"ts"。如果你想这样做，动物园包可能是你想要的方式。这是一个示例，我们使用以下方法转换为 zoo 对象as.zoo()

require(zoo)
sz <- as.zoo(signal)

然后我们可以提取我们想要的观察结果（名称是时间索引的有用指标）

> sz[sz[, "i"] == 1, "d"]
1975(12)  1976(2) 1976(12) 1977(12)  1978(1) 1978(12)  1979(1) 
     761     1104      810      653      522      956      593

然后一个类似的子集调用，但使用index()返回整个动物园对象的时间索引，然后选择我们想要的位

> index(sz)[sz[, "i"] == 1]
[1] 1975.917 1976.083 1976.917 1977.917 1978.000 1978.917
[7] 1979.000
于 2012-11-07T22:00:04.580 回答
1

如果要重命名mts（矩阵时间序列）对象的列。只是简单地使用colnames()函数。

例如：您有一个包含三列的时间序列数据，并且您想要更改该列名称：

colnames(data) <- c('newname1','newname2','newname3')

希望这可以帮助。

于 2017-09-16T07:11:56.617 回答