Question

我有一个有三列的动物园对象。我想计算每列的成对差异，以获得随时间的相对变化。我为此使用了for循环。但是，我得到的东西不再是一个动物园对象，而是一个普通的data.frame。

当我想在多个面板中绘制它时，我会松开我的x轴时间标签。当我尝试将数据帧转换回动物园对象时，我收到错误的＆＃34;错误的条目＆＃34;

head(zoo.2017)

          zoo1  zoo2     zoo3
2017-01-01 104.6066 106.6240 112.7336
2017-01-02 103.1034 104.4088 111.8473
2017-01-03 103.4550 107.0606 111.4292
2017-01-04 101.6916 108.0550 111.8183
2017-01-05 101.8274 105.1528 114.3505
2017-01-06 102.9502 107.1381 113.8313

delta<-zoo()

for(n in 2:length(zoo.2017[,1])){

    diff<-as.data.frame(zoo.2017[n,])-as.data.frame(zoo.2017[n-1,])
    delta<-rbind(delta, diff)
}

head(delta)
               Maui6   Olowalu        Pali
 diff       -1.503264 -2.215139  -0.8862847
 2017-01-03 0.3516667  2.651771  -0.4180556
 2017-01-04 -1.763472 0.9944444   0.3890625
 2017-01-05 0.1358681 -2.902257    2.532153
 2017-01-06  1.122743  1.985313  -0.5191319
 2017-01-07  1.430937 0.5648958 -0.01836806

我知道错误可能来自for循环中的as.data.frame（），但我不知道如何计算它。当我在数据帧上尝试read.zoo（）之后，它给了我错误的条目错误...

我仍然很难与动物园合作

Answer 1

如果z是末尾注释中可重复定义的Zoo对象，则要在各列中获得连续的差异，请使用diff.zoo

zd <- diff(z)

给予该动物园对象：

> zd
              zoo1    zoo2    zoo3
2017-01-02 -1.5032 -2.2152 -0.8863
2017-01-03  0.3516  2.6518 -0.4181
2017-01-04 -1.7634  0.9944  0.3891
2017-01-05  0.1358 -2.9022  2.5322
2017-01-06  1.1228  1.9853 -0.5192

关于绘图，我们现在可以使用以下任意一种进行绘图：

plot(zd)

library(ggplot2)
autoplot(zd)

library(lattice)
xyplot(zd)

zoo带有多个小插图（pdf文档）和一本参考手册，您可以阅读以获取更多信息。

注意

Lines <- "date,zoo1,zoo2,zoo3
2017-01-01,104.6066,106.6240,112.7336
2017-01-02,103.1034,104.4088,111.8473
2017-01-03,103.4550,107.0606,111.4292
2017-01-04,101.6916,108.0550,111.8183
2017-01-05,101.8274,105.1528,114.3505
2017-01-06,102.9502,107.1381,113.8313"
library(zoo)
z <- read.csv.zoo(text = Lines)

数学和rbind与动物园对象

1 个答案:

注意