3

我最近找到了该zoo功能na.approx，但想知道是否有一种方法可以在不改变时区显示方式的情况下使用它，如运行我的可重现示例所示...

set.seed(111)
x <- xts(matrix(cumprod(rnorm(100,0,0.001)+1)*100,
                ncol=1,dimnames=list(rep("",100),c("JJ"))),
         Sys.time()-c(100:1), tzone="America/Chicago")

x[30:50,] <- NA

tzone(x)
[1] "America/Chicago"
.index(head(x,1))
[1] 1377270598
attr(,"tzone")
[1] "America/Chicago"
attr(,"tclass")
[1] "POSIXct" "POSIXt" 
head(x)
                           JJ
2013-08-23 10:09:57 100.02352
2013-08-23 10:09:58  99.99044
2013-08-23 10:09:59  99.95928
2013-08-23 10:10:00  99.72914
2013-08-23 10:10:01  99.71210
2013-08-23 10:10:02  99.72609
Warning message:
timezone of object (America/Chicago) is different than current timezone (). 



y <- na.approx(x)

tzone(y)
[1] "America/Chicago"
.index(head(y,1))
[1] 1377270598
head(y)
                           JJ
2013-08-23 16:09:57 100.02352
2013-08-23 16:09:58  99.99044
2013-08-23 16:09:59  99.95928
2013-08-23 16:10:00  99.72914
2013-08-23 16:10:01  99.71210
2013-08-23 16:10:02  99.72609
Warning message:
timezone of object (America/Chicago) is different than current timezone ().

注意我不在芝加哥....而且我Sys.getenv('TZ')没有设置为"America/Chicago"

4

1 回答 1

2

编辑：我可以用 xts 0.9-3 重现这个问题。它已在 xts 0.9-5 中修复。您需要更新 xts。

我相信这可以被认为是一个错误。

head(x)
#                           JJ
#2013-08-23 14:12:19 100.02352
#2013-08-23 14:12:20  99.99044
#2013-08-23 14:12:21  99.95928
#2013-08-23 14:12:22  99.72914
#2013-08-23 14:12:23  99.71210
#2013-08-23 14:12:24  99.72609

methods(na.approx)
#[1] na.approx.default na.approx.ts*     na.approx.zoo*    na.approx.zooreg*

所以，没有 xts 方法。自从

class(x)
#[1] "xts" "zoo"

na.approx.zoo被派遣。在该函数内部，基本上会发生这种情况：

index(x) <- index(x)
head(x)
#                           JJ
#2013-08-23 21:12:19 100.02352
#2013-08-23 21:12:20  99.99044
#2013-08-23 21:12:21  99.95928
#2013-08-23 21:12:22  99.72914
#2013-08-23 21:12:23  99.71210
#2013-08-23 21:12:24  99.72609

请注意，index<-.xts存在并在此处使用。不幸的是，与 相反index<-.zoo，它不保留（索引的）时区。

因此，您可以使用y <- na.approx(as.zoo(x))来避免该问题。

于 2013-08-23T19:28:36.443 回答