<div dir="ltr">A dica esta' na descricao do proprio Matlab a respeito de como ele faz as coisas (que me faz lembrar o quao importante e' conhecer os detalhes das ferramentas que utilizamos):<div><br></div><div>"If x and y are not the same length, the shorter vector is zero-padded to the length of the longer vector."</div><div><br></div><div>b</div></div><div class="gmail_extra"><br><div class="gmail_quote">Em 23 de setembro de 2014 19:50, Daniel Marcelino <span dir="ltr"><<a href="mailto:dmsilva.br@gmail.com" target="_blank">dmsilva.br@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Paulo,<br>
Não sei se meu comentário vai te ajudar, mas partindo da pretensão de<br>
que as duas funções são equivalentes, valeria a pena colocar uma do<br>
lado da outra para ver no que elas diferem. Aparentemente a função ccf<br>
do R desconsidera o que "sobra".<br>
<br>
n=0:50;<br>
N=sample(length(n))<br>
N2=sample(length(n)/2)<br>
(ccf(N,N))<br>
(ccf(N,N2))<br>
(ccf(N2,N, type="covariance"))<br>
<br>
Daniel<br>
<br>
<br>
2014-09-23 18:11 GMT-03:00 Paulo Nogueira Starzynski <<a href="mailto:paulons@gmail.com">paulons@gmail.com</a>>:<br>
<div><div class="h5">> Olá colegas de lista.<br>
><br>
> Tenho em mãos um programa escrito para Matlab e estou tentando fazer uma<br>
> versão para R. O projeto faz parte de uma migração de longo prazo que talvez<br>
> façamos no longo prazo.<br>
><br>
> Como muitos conhecem "de vista" os registros de eletrocardiograma vou<br>
> contextualizar meu problema, que consiste em detectar o instante em que<br>
> ocorrem os batimentos do coração.<br>
><br>
> Tenho em mãos um script em Matlab que detecta os "picos do ECG" (que na<br>
> verdade se chama complexo QRS). Ele faz isso baseado na correlação cruzada<br>
> entre o vetor contendo o eletrocardiograma completo e um vetor menor<br>
> contendo apenas um complexo QRS.<br>
><br>
> A função xcorr do Matlab faz "correr" o vetor menor através do maior,<br>
> apontando os picos de correlação no instante da coincidência (ocorrência do<br>
> complexo QRS).<br>
> Estou tentando fazer o mesmo no R com através da stats:: ccf, mas sem<br>
> sucesso, porque ela exige séries de tamanhos iguais.<br>
><br>
> Alguém já resolveu esse problema? Procurei bastante e não achei alternativa.<br>
> Aceito sugestões e se for preciso posso enviar um CRM.<br>
><br>
> Obrigado!<br>
><br>
> ==================================<br>
> Paulo Nogueira Starzynski, B.Sc.<br>
> Estudante de Mestrado em Fisiologia Geral<br>
> Laboratório de Fisiologia Teórica (LFT)<br>
> Instituto de Biociências<br>
> Universidade de São Paulo<br>
><br>
><br>
</div></div>> _______________________________________________<br>
> R-br mailing list<br>
> <a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código<br>
> mínimo reproduzível.<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br>
</blockquote></div><br></div>