我最近发现了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 ().
N.B。我不在芝加哥......我的Sys.getenv('TZ')
未设置为"America/Chicago"
答案 0 :(得分: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))
来避免此问题。