1

我有以下矩阵:

B = as.POSIXct("2015-06-21") + seq(1,18)
P = c(21.675,21.775,21.675,21.725,21.575,21.525,21.625,21.475,21.375,21.325,21.32,21.475,21.375,21.325,21.425,21.525,21.275,21.325)
A = c(-1,1,-1,1,-1,0,1,-1,0,0,0,1,-1,0,1,0,-1,1)
C = c(rep(0,15),1,0,0)
D = c(0,0,0,0,0,-1,0,0,-1,-1,-1,0,0,-1,0,0,0,0)
E = seq(1,18)  

dat = as.data.frame(cbind(P,A,C,D))
rownames(dat) = B

输出矩阵如下所示:

                         P  A C  D  E
2015-06-21 00:00:01 21.675 -1 0  0  1
2015-06-21 00:00:02 21.775  1 0  0  2
2015-06-21 00:00:03 21.675 -1 0  0  3
2015-06-21 00:00:04 21.725  1 0  0  4
2015-06-21 00:00:05 21.575 -1 0  0  5
2015-06-21 00:00:06 21.525  0 0 -1  6
2015-06-21 00:00:07 21.625  1 0  0  7
2015-06-21 00:00:08 21.475 -1 0  0  8
2015-06-21 00:00:09 21.375  0 0 -1  9
2015-06-21 00:00:10 21.325  0 0 -1 10
2015-06-21 00:00:11 21.325  0 0 -1 11
2015-06-21 00:00:12 21.475  1 0  0 12
2015-06-21 00:00:13 21.375 -1 0  0 13
2015-06-21 00:00:14 21.325  0 0 -1 14
2015-06-21 00:00:15 21.425  1 0  0 15
2015-06-21 00:00:16 21.525  0 1  0 16
2015-06-21 00:00:17 21.275 -1 0  0 17
2015-06-21 00:00:18 21.325  1 0  0 18

A 列有交替的 -1 和 1。目标是查看有零的 A 列,看看我是否可以使用 C 和 D 替换它们。

E 列表示行号,P 列包含价格信息。一列表示局部最小和最大点(-1 和 1)。

所以在第 4 行和第 7 行之间,第 6 行有一个零。4 和 7 表示局部最大值点,分别为 21.725 和 21.625。

我想检查第 5 行是否确实是第 4 行和第 7 行之间的最小数字。结果,第 6 行的值较低,21.525 与第 5 行的 21.575 相比。

因此,我想将 A 列的第 5 个条目从 -1 替换为 0,将第 6 个条目从 0 替换为 -1。

我想用上述逻辑遍历矩阵,最后得到以下结果:

                         P  A C  D  E
2015-06-21 00:00:01 21.675 -1 0  0  1
2015-06-21 00:00:02 21.775  1 0  0  2
2015-06-21 00:00:03 21.675 -1 0  0  3
2015-06-21 00:00:04 21.725  1 0  0  4
2015-06-21 00:00:05 21.575  0 0  0  5
2015-06-21 00:00:06 21.525 -1 0 -1  6
2015-06-21 00:00:07 21.625  1 0  0  7
2015-06-21 00:00:08 21.475  0 0  0  8
2015-06-21 00:00:09 21.375  0 0 -1  9
2015-06-21 00:00:10 21.325 -1 0 -1 10
2015-06-21 00:00:11 21.325  0 0 -1 11
2015-06-21 00:00:12 21.475  1 0  0 12
2015-06-21 00:00:13 21.375  0 0  0 13
2015-06-21 00:00:14 21.325 -1 0 -1 14
2015-06-21 00:00:15 21.425  0 0  0 15
2015-06-21 00:00:16 21.525  1 1  0 16
2015-06-21 00:00:17 21.275 -1 0  0 17
2015-06-21 00:00:18 21.325  1 0  0 18

我在编写经过此逻辑的高效代码时遇到了麻烦。任何帮助将不胜感激!

4

0 回答 0