如何计算向量中的符号变化（正数或负数）

Question

我目前正在使用一个大型数据集，在过去20年中，每天有40k km2的每日相对积雪值（每km2），并希望获得积雪在特定时期内增加的日子。 每天进行一次测量，因此，相对的积雪可能会在几天内增加/减少。

但是，我很难计算/指出我的每日值与前几天相比减少/增加的次数。 我创建了一些假的（但可能的）数据，并计算了元素n和元素n + 1之间的差异：

set.seed(134)
c <- c(sort(runif(6, 0, 1), decreasing = FALSE),
  sort(runif(10, 0, 1), decreasing = TRUE),
  sort(runif(4, 0, 1), decreasing = FALSE),
  0.9, 0.9, 0.9, sort(runif(3, 0, 1), decreasing = TRUE),
  sort(runif(8, 0, 1), decreasing = FALSE)) # create fake snow cover fraction data
plot(c)
cdiff <- c [1:length(c) - 1] - c [2:length(c)]

然后我想到了用符号表示代数符号

> sign(cdiff)
 [1] 1  1  1  1  1  1  0  1 -1  1  1  1  1  1  1  1  1  1  1 -1 -1 -1  1  1  1 -1 -1 -1 -1 -1 -1 -1 -1

我现在有两个问题：

1.）如何获取代数符号改变的新向量的位置

2.）甚至可以在（i）两个正数，（ii）两个负数或（iii）混合数之间找到值。 在我的示例中，我想将相同的值视为这些值周围的子向量的一部分（在本例中为正号，但积雪分数减小了）。

我想要的解决方案是接收一个向量，该向量的天数与前一天相比又变为正数

解决方案不必包含sign() ，这只是我认为可以使我更接近解决问题的方法。

提前致谢！

Answer 1

您可以使用diff和which 找到的位置，其中的符号改变或者其变正 。

#Location where sign changes
which(diff(sign(diff(c))) != 0)+1
#[1]  7 16 20 21 23 26

#Turn positive
which(diff(sign(diff(c))) > 0)+1
#[1] 16 21 26

或更健壮的找到所有山谷 ：

peakPosition <- function(x, inclBorders=TRUE) {
  if(inclBorders) {y <- c(min(x), x, min(x))
  } else {y <- c(x[1], x)}
  y <- data.frame(x=sign(diff(y)), i=1:(length(y)-1))
  y <- y[y$x!=0,]
  idx <- diff(y$x)<0
  (y$i[c(idx,F)] + y$i[c(F,idx)] - 1)/2
}

peakPosition(-c) #Get valleys
[1]  1 16 26