1

Tenho o seguinte vetor:

a <- c(1,2,3,1,1,2,3,4,2,5,6,7,8,4,6,4,3,2,1)

Gostaria de achar a maior subsequência crescente e decrescente do mesmo. A saída seria os índices da sequencia dos vetores:

ss_cres <- c(9,10,11,12,13) # ou seja 2,5,6,7 e 8

ss_decres <- c(15,16,17,18,19) # ou seja 6,4,3,2 e 1

Essa é uma operação comum utilizada em análise de séries temporais. Existe algum pacote no R que faça isso? Se não, alguma dica para uma ideia geral para o algoritmo baseado em experiências?

Até o momento tenho feito isso para a maior série decrescente:

long_seq_decres<- function(seq){

seq_inds <- split(1:length(seq), cumsum(c(0, diff(seq)) > 0))  
  
ind_list <- lengths(seq_inds)

long_seq_inds <- seq_inds[which.max(ind_list)]

return(long_seq_inds)

}

Funcionou bem, entretanto achei um pouco lenta a função. Alguma ideia para otimizar esta função?

1 Resposta 1

2

Pode usar rle para poupar a etapa de criar a lista dividida:

long_seq_dec <- function(seq) {
  seq_inds <- split(1:length(seq), cumsum(c(0, diff(seq)) > 0))
  ind_list <- lengths(seq_inds)
  seq_inds[[which.max(ind_list)]]
}

long_seq <- function(seq, type = c("increasing", "decreasing")) {
  type <- match.arg(type)
  dif <- c(0, diff(seq))
  csd <- if (type == "increasing") cumsum(dif <= 0) else cumsum(dif > 0)
  cmps <- rle(csd)
  seq_len(length(seq))[csd == with(cmps, values[lengths == max(lengths)])]
}

a <- c(1,2,3,1,1,2,3,4,2,5,6,7,8,4,6,4,3,2,1)

long_seq_dec(a)
#> [1] 15 16 17 18 19

long_seq(a, "d")
#> [1] 15 16 17 18 19

long_seq(a)
#> [1]  9 10 11 12 13


library(microbenchmark)

x <- sample(1:10, 1e4, TRUE)

microbenchmark(long_seq_dec(x), long_seq(x, "d"))
#> Unit: milliseconds
#>             expr      min       lq     mean   median       uq       max neval
#>  long_seq_dec(x) 5.701912 5.916570 6.247751 5.976668 6.200649 13.209920   100
#> long_seq(x, "d") 1.627581 1.684397 1.919555 1.718479 1.869197  8.719684   100
1
  • Muito bom! Obrigado pela dica! Commented 25/03/2021 às 16:39

Você deve fazer log-in para responder a esta pergunta.

Esta não é a resposta que você está procurando? Pesquise outras perguntas com a tag .