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?