Pessoal, consegui resolver aqui com uma pesquisa mais aprofundada. Segue o código de leitura dos picos e vales para a amostra utilizada: maximo = data.frame() for(i in 1:length(amostra)){ if( max(amostra[i+1], amostra[i-1]) < amostra[i]) maximo = rbind(maximo, data.frame(print(amostra[i]))) } #LEITURA DOS MINIMOS minimo = data.frame() for(i in 1:length(amostra)){ if( min(amostra[i+1], amostra[i-1]) > amostra[i]) minimo = rbind(minimo, data.frame(print(amostra[i]))) } Atenciosamente, José Lucas. Message: 1
Date: Sat, 19 Dec 2015 19:04:52 -0200 From: José Lucas Safanelli <zecojls@gmail.com> To: r-br@listas.c3sl.ufpr.br Subject: [R-br] Loop para leitura de picos Message-ID: < CADni8G5B-N9M4LqpZTVkAOgoSkTO56aK0icdg-PjVf0Us8PcFQ@mail.gmail.com> Content-Type: text/plain; charset="utf-8"
Boa tarde, estou com o seguinte caso:
DADO: amostra = c(49.936157, 48.133488, 45.775238, 43.229397, 40.021607, 35.927558, 31.365912, 26.849848, 23.233177, 20.770278, 18.572067, 15.590341, 11.936978, 8.478520, 5.572115, 3.272254, 1.642041, 0.999997, 1.473886, 2.616919, 4.169691, 6.464019, 9.790810, 13.820101, 17.978423, 21.674281, 24.065947, 24.875069, 24.590811, 23.620608, 22.240339, 20.676521, 19.088229, 17.689105, 16.720912, 16.358739, 16.792281, 18.321545, 21.053896, 24.777756, 29.121638, 33.551339, 37.628649, 41.272700, 44.405581, 46.899820, 48.602420, 49.253714, 48.816692, 47.380737, 44.920718) wavelength = seq(400, 450, by=1) plot(wavelenght, amostra, xlab = "Comp. Onda (nm)", ylab = 'Valor', type = "l") axis(side=1, at=seq(400, 450, by=5))
PROBLEMA: Estou tentando ler os valores mínimos e máximos de cada pico para ter o valor da amplitude. Consigo fazer isso delimitando o intervalo de comp. de ondas da posição do pico, mas é muito laborioso EX: max(amostra[,1:20]) e min(amostraamostra[,1:20]). Amplitude seria a diferença deles.
Para várias amostras tenho utilizado a função apply para dentro da linha, PORÉM decorrente da heterogeneidade das amostras, as posições variam causando erro de leitura e modelagem.
Outro fato é que existem inúmeros picos dentro de uma faixa de interesse, e delimitando visualmente torna-se inviável. Para isso, estou tentando criar uma rotina que leia todos os picos dentro da faixa de interesse, para que possa utiliza-los diferentemente nas modelagens.
Tentei o código abaixo, mas me retorna o erro que o argumento é de comprimento zero!
for(i in 1:length(amostra)){ if((amostra[i+1] > amostra[i])&(amostra[i-1]>amostra[i])) print(amostra[i]) }
(Error in if ((amostra[i + 1] > amostra[i]) & (amostra[i - 1] > amostra[i])) print(amostra[i]) : argument is of length zero)
Um resultado que observei foi de do seguinte código, somente com uma condição colocada nele:
for(i in 1:length(amostra)){ if(amostra[i+1] > amostra[i]) print(amostra[i]) }
Nesse código, consegue-se ler os valores crescentes dos picos, desconsiderando os valores do objeto quando os valores são decrescentes. A partir dele que tentei colocar duas condições depois do IF para que ele lê-se somente os valores mínimos, mas sem sucesso.
Portanto, peço ajuda para confecção de um código que faça a leitura dos min e max, ou amplitude, de todos os picos da amostra, retornando em um objeto/dataframe para modelagem.
Grato pela atenção! José Lucas. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: < http://listas.inf.ufpr.br/pipermail/r-br/attachments/20151219/add2cc3b/attac...