Largura da janela em movimento deve ser um número inteiro entre 1 e n uma opção para escolher diferentes algoritmos C - uma versão é escrita em C. Pode tratar números não finitos como NaNs e Infs (como mean (x, na. rm TRUE)) . Ele funciona o mais rápido para endrulemean. Rápido - segundo, ainda mais rápido, versão C. Este algoritmo não funciona com números não finitos. Ele também funciona o mais rápido para endrule que não seja médio. R - muito mais lento código escrito em R. Útil para depuração e como documentação. Exato - o mesmo que C. exceto que todas as adições são executadas usando o algoritmo que segue e corrige a adição arredondar-fora da corda de caráter dos erros que indica como os valores no começo e no fim, dos dados, devem ser tratados. Somente os primeiros e últimos valores de k2 em ambas as extremidades são afetados, onde k2 é a metade da largura de banda k2 k / 2. média - aplica a função subjacente a seções menores e menores da matriz. Equivalente a: for (i in 1: k2) outi mean (x1: (ik2)). Esta opção é implementada em C se algC. Caso contrário, é feito em R. trim - trim o comprimento da matriz de saída de termina é igual ao comprimento (x) -2k2 (out out (k21) :( n-k2)). Esta opção imita a saída de apply (embed (x, k), 1, mean) e outras funções relacionadas. - preencher as extremidades com os números do vetor x (out1: k2 x1: k2) constante - preencher as extremidades com o primeiro eo último valor calculado na matriz de saída (out1: k2 outk21) NA - preencher as extremidades com NAs (out1: k2 NA ) Func - igual à média mas implimented em R. Esta opção poderia ser muito lenta, e é incluída principalmente para testar Similar ao endrule na função runmed que tem as seguintes opções: ldquo c (mediana, manter, constante) rdquo. specifies se o resultado Deve ser centralizado (padrão), alinhado à esquerda ou alinhado à direita. Se endrule significa, em seguida, definir alinhamento para a esquerda ou direita vai cair para trás em execução mais lenta equivalente ao func endrule. Além dos valores finais, o resultado de y runmean (x, k) é o mesmo que ldquo para (j (1k2) :( n-k2)) yjmean (x (j-k2) :( jk2)) rdquo. O principal incentivo para escrever este conjunto de funções foi a relativa lentidão da maioria das funções de janelas móveis disponíveis em R e seus pacotes. Com exceção de runmed. Uma função mediana de janela de execução, todas as funções listadas na seção ver também são mais lentas do que muito ineficientes ldquo apply (embed (x, k), 1, FUN) abordagem rdquo. A velocidade relativa da função runmean é O (n). Função EndRule aplica um dos cinco métodos (veja o argumento endrule) para processar pontos finais da matriz de entrada x. Na versão atual do código, a opção endrulemean padrão é calculada dentro do código C. Isso é feito para melhorar a velocidade no caso de grandes janelas móveis. No caso da função runmean (.algexact), um algoritmo especial é usado (veja a seção de referências) para garantir que os erros de arredondamento não se acumulam. Como resultado, runmean é mais preciso do que as funções filter (x, rep (1 / k, k)) e runmean (.AlcC). Retorna um vetor numérico ou matriz do mesmo tamanho de x. Somente no caso de endruletrim os vetores de saída serão mais curtos e as matrizes de saída terão menos linhas. A função runmean (.algexact) é baseada no código de Vadim Ogranovich, que é baseado no código Python (ver última referência), apontado por Gabor Grothendieck. Referências Sobre a correção de erro de arredondamento usada no runmean. Shewchuk, Jonathan Precisão Adaptativa Flutuante-Ponto Aritmética e Fast robusto Geometric Predicates. Www-2.cs. cmu. edu/afs/cs/project/quake/public/papers/robust-arithmetic. ps Mais informações sobre a correção de erros de arredondamento podem ser encontradas em: aspn. activestate / ASPN / Cookbook / Python / Recipe / 393090 Links relacionados a: média móvel - média. Kernapply. filtro. decompor. Stl. Rollmean da biblioteca do jardim zoológico, subsums da biblioteca mágica, outras funções móveis da janela deste pacote: runmin. Runmax. Runquantile Runmax e rund runmed funções genéricas de execução de janela: aplicar (incorporar (x, k), 1, FUN) (mais rápido), executando a partir do pacote gtools (extremamente lento para esta finalidade), subsums da biblioteca mágica pode executar operações de janela em execução em dados com Qualquer dimensão. Pacote caTools versão 1.12 IndexMoving Médias em R Para o melhor de meu conhecimento, R não tem uma função interna para calcular médias móveis. Usando a função de filtro, no entanto, podemos escrever uma função curta para médias móveis: Podemos então usar a função em qualquer dado: mav (dados) ou mav (data, 11) se quisermos especificar um número diferente de pontos de dados Do que o padrão 5 plotando obras como esperado: plot (mav (dados)). Além do número de pontos de dados sobre os quais a média, também podemos alterar o argumento de lados das funções de filtro: sides2 usa ambos os lados, sides1 usa apenas valores passados. Compartilhe: Navegação de posts Navegação de comentários Navegação de comentários Funções incorporadas Extraia ou substitui substrings em um vetor de caracteres. X lt - quotabcdefquot substr (x, 2, 4) é quotbcdquot substr (x, 2, 4) lt - quot22222quot is quota222efquot Busca por padrão em x. Se FALSE fixo, o padrão é uma expressão regular. Se fixedTRUE, então pattern é uma string de texto. Retorna índices correspondentes. Grep (quotAquot, c (quotbquot, quotAquot, quotcquot), fixedTRUE) retorna 2 Encontre padrão em x e substitua por texto de substituição. Se fixedFALSE então pattern é uma expressão regular. Se T fixo, o padrão é uma string de texto. Sub (quotsquot, quot. quot, quotHello Therequot) retorna quotHello. Therequot Dividir os elementos do vetor de caracteres x na divisão. Strsplit (quotabcquot, quotquot) retorna 3 vetor de elemento quotaquot, quotbquot, quotcquot Concatenate strings depois de usar a seqüência de caracteres sep para separá-los. Colar (quotxquot, 1: 3, sepquotquot) retorna c (quotx1quot, quotx2quot quotx3quot) paste (quotxquot, 1: 3, sepquotMquot) retorna c (quotxM1quot, quotxM2quotxM3quot) paste (quotToday isquot, date ()
Comments
Post a Comment