Question

我有一个时间序列如下。我只是想找到t和t+1处的观察结果具有不同符号的指数。

显然，我可以通过将索引放在t*(t-1)<0来实现这一点。但是我不确定实现这种逻辑的最佳功能。

在基础R 中实现此目标的最快方法是什么？

b <- structure(c(NA, NA, NA, NA, NA, NA, NA, NA, NA, 0.00117300000000009, 
            0.00165400000000004, 0.00262999999999991, 0.00236499999999995, 
            0.00571200000000016, 0.00364100000000001, 0.0013479999999999, 
            -0.0013510000000001, -0.00270499999999996, -0.00338699999999981, 
            -0.000166000000000111, -0.00105099999999991, -0.002467, -0.00367300000000004, 
            -0.00569500000000001, -0.00446499999999994, -0.002861, -0.00129699999999988, 
            0.000721000000000194, 0.00351799999999969, 0.0040929999999999, 
            0.00418700000000016, 0.004189, 0.00348199999999999, 0.0015729999999996, 
            5.89999999998092e-05, -0.000866000000000255, -0.00130799999999986, 
            0.00158199999999997, 0.00198000000000009, 0.00180899999999995, 
            0.000526999999999944, -0.000598000000000098, -0.00189899999999987, 
            -0.00285700000000011, -0.00296200000000013, -0.000966000000000244, 
            0.00388200000000016, 0.00311499999999998, 0.00300600000000006
), .Tsp = c(1, 49, 1), class = "ts")

Answer 1

您也可以使用ts个对象执行此操作，但我使用zoo或xts更多。所以这是动物园的简单解决方案：

R> suppressMessages(library(zoo))
R> bz <- as.zoo(b)
R> bz[which(c(NA,diff(sign(bz))) != 0)]
       17        28        36        38        42        47 
-0.001351  0.000721 -0.000866  0.001582 -0.000598  0.003882 
R>

编辑：同样，为了获得索引，您可以使用zoo提供的便利访问器：

R> index(bz[which(c(NA,diff(sign(bz))) != 0)])
[1] 17 28 36 38 42 47
R>

编辑2 ：如果确实想要坚持ts个对象和基础R，请通过tsp()重建索引系列并使用它：

R> tsp(b)[1] + which(diff(sign(b)) != 0) * tsp(b)[3]
[1] 17 28 36 38 42 47
R>

R时间序列发现符号变化的指数

1 个答案: