R na média móvel

A largura da janela em movimento deve ser um número inteiro entre 1 e n uma opção para escolher algoritmos diferentes C - uma versão está escrita em C. Ela pode lidar com números não finitos como NaNs e Infs (como a média (x, na. rm TRUE)). . 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 endrulagem diferente da média. R - código muito mais lento escrito em R. Útil para depuração e como documentação. Exato - o mesmo que C., exceto que todas as adições são realizadas usando o algoritmo que rastreia e corrige a sequência de caracteres de arredondamento de adição, indicando como os valores no início e no final dos dados devem ser tratados. Apenas os primeiros e últimos valores de k2 em ambas as extremidades são afetados, onde k2 é a largura de banda k2 k 2. significa - aplica a função subjacente a seções menores e menores da matriz. Equivalente a: para (i em 1: k2) outi significa (x1: (ik2)). Esta opção é implementada em C se o algC. Caso contrário, é feito em R. trim - trim the ends O comprimento da matriz de saída é igual ao comprimento (x) -2k2 (out out (k21) :( n-k2)). Esta opção imita a saída de aplicar (incorporar (x, k), 1, significar) e outras funções relacionadas. Manter - preencher as extremidades com números de x vetor (out1: k2 x1: k2) constante - preencher as extremidades com o primeiro e último valor calculado na matriz de saída (out1: k2 outk21) NA - preencher as extremidades com NAs (out1: k2 NA ) Func - o mesmo que a média, mas implantado em R. Esta opção pode ser muito lenta e está incluída principalmente para testes semelhantes a endrule na função runmed que possui as seguintes opções: ldquo c (mediana, keep, constant) rdquo. specifica se o resultado Deve ser centrado (padrão), alinhado à esquerda ou alinhado à direita. Se endrule significa que o ajuste alinhado à esquerda ou à direita retornará a implementação mais lenta equivalente a endrule func. 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 esse conjunto de funções foi a lentidão relativa da maioria das funções de janela em movimento disponíveis em R e seus pacotes. Com a exceção de runmed. Uma função de medição de janela de execução, todas as funções listadas na seção veja também são mais lentas do que as ldquo muito ineficientes aplicam (incorporar (x, k), 1, FUN) abordagem rdquo. A velocidade relativa da função runmean é O (n). Function EndRule aplica um dos cinco métodos (veja o argumento endrule) para processar os pontos finais da matriz de entrada x. Na versão atual do código, a opção default endrulemean é calculada dentro do código C. Isso é feito para melhorar a velocidade no caso de grandes janelas em movimento. No caso da função runmean (.gegexact), um algoritmo especial é usado (ver seção de referências) para garantir que os erros de arredondamento não se acumulam. Como resultado, o modo de execução é mais preciso do que as funções filtro (x, rep (1k, k)) e runmean (. AlgC). Retorna um vetor ou matriz numérica do mesmo tamanho que x. Somente em caso de endruletrim, os vetores de saída serão mais curtos e as matrizes de saída terão menos linhas. Função runmean (.algexact) é baseado em código por Vadim Ogranovich, que é baseado no código Python (veja a última referência), apontada por Gabor Grothendieck. Referências Sobre a correção de erro de arredondamento usada no runmean. Shewchuk, Jonathan Adaptive Precision Floating-Point Arithmetic and Fast Robust Geometric Predicates. Www-2.cs. cmu. eduafscsprojectquakepublicpapersrobust-arithmetic. ps Mais informações sobre a correção de erros de arredondamento podem ser encontradas em: aspn. activestateASPNCookbookPythonRecipe393090 Links relacionados a: significa-mudança - significa. Kernapply. filtro. decompor. Stl. Rollmean da biblioteca do zoológico, subsums da biblioteca mágica, Outras funções da janela móvel deste pacote: runmin. Runmax. Runquantile. Runmad e rund funções genéricas de janela de execução: aplique (incorporar (x, k), 1, FUN) (mais rápido), executando o pacote de gtools (extremamente lento para esse propósito), os subsums da biblioteca mágica podem executar operações de janela em execução com dados com Quaisquer dimensões. Pacote caTools versão 1.12 IndexMoving Médias em R No meu melhor conhecimento, R não possui uma função integrada para calcular as 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 (dados, 11) se quisermos especificar um número diferente de pontos de dados Do que o traçado padrão 5 funciona como esperado: plot (mav (data)). Além do número de pontos de dados sobre os quais a média, também podemos alterar o argumento lateral das funções de filtro: sides2 usa ambos os lados, sides1 usa apenas valores passados. Compartilhe isto: Publicar navegação Comentar navegação Comentar navigationgt mav (c (4,5,4,6), 3) Série temporal: Iniciar 1 Final 4 Frequência 1 1 NA 4.333333 5.000000 NA Aqui estava tentando fazer uma média móvel que levou a Conta os últimos 3 números, então eu esperava obter apenas dois números de volta 8211 4.333333 e 5 8211 e se houvesse valores de NA, pensei que eles deveriam estar no início da seqüência. Na verdade, isso é o que o parâmetro 8216sides8217 controla: lados apenas para filtros de convolução. Se os lados 1 os coeficientes de filtro são para valores passados ​​apenas se os lados 2 estiverem centrados em torno de lag 0. Neste caso, o comprimento do filtro deve ser estranho, mas se for uniforme, mais do filtro está em frente do que para trás. Portanto, na nossa função 8216mav8217, a média de rolamento parece ter ambos os lados do valor atual em vez de apenas em valores passados. Nós podemos ajustar isso para obter o comportamento que queremos: gt library (zoo) gt rollmean (c (4,5,4,6), 3) 1 4.333333 5.000000 Eu também percebi que posso listar todas as funções em um pacote com o 8216ls8217 Função, então eu vou estar procurando digitalizar a lista de funções da zoo8217s na próxima vez que eu precisar fazer algumas séries temporais relacionadas 8211 there8217ll provavelmente já é uma função para ela gt ls (quotpackage: zooquot) 1 quotas. Datequot quotas. Date. numericquot quotas. Date. tsquot 4 Quotas. Date. yearmonquot quotas. date. yearqtrquot quotas. yearmonquot 7 quotas. yearmon. defaultquot quotas. yearqtrquot quotas. yearqtr. defaultquot 10 quotas. zooquot quotas. zoo. defaultquot quotas. zooregquot 13 quotas. zooreg. defaultquot quotautoplot. zooquot quotcbind. Zooquot 16 quotcoredataquot quatcoredata. defaultquot quotcoredatalt-quot 19 quotfacetfreequot quotformat. yearqtrquot quotfortify. zooot 22 quotfrequencylt-quet quotifelse. zooot quotindexquot 25 quotindexlt-quotindex2charquot quotis. regularquot 28 quotis. zooquot quotmake. par. listquot q UotMATCHquot 31 quatMATCH. defaultquot quatMATCH. timesquot quimedial. zoootot 34 quotmerge. zooquot quotna. aggregatequot quotna. aggregate. defaultquot 37 quotna. approxquot quotna. approx. defaultquot quotna. fillquot 40 quotna. fill. defaultquot quotna. locfquot quotna. locf. defaultquot 43 Quotna. splinequot quotna. spline. defaultquot quotna. StructTSquot 46 quotna. trimquot quotna. trim. defaultquot quotna. trim. tsquot 49 quotORDERquot quotORDER. defaultquot quotpanel. lines. itsquot 52 quotpanel. lines. tisquot quotpanel. lines. tsquot quotpanel. lines. Zooquot 55 quotpanel. plot. customquot quotpanel. plot. defaultquot quotpanel. points. itsquot 58 quotpanel. points. tisquot quotpanel. points. tsquot quotpanel. points. zooot 61 quotpanel. polygon. itsquot quotpanel. polygon. tisquot quotpanel. polygon. tsquot 64 Quotpanel. polygon. zooquot quotpanel. rect. itsquot quotpanel. rect. tisquot 67 quotpanel. rect. tsquot quotpanel. rect. zooot quotpanel. segments. itsquot 70 quotpanel. segments. tisquot quotpanel. segments. tsquot quotpanel. se Gments. zooot 73 quotpanel. text. itsquot quotpanel. text. tisquot quotpanel. text. tsquot 76 quotpanel. text. zooquot quotplot. zooquot quotquantile. zoootot 79 quotrbind. zooquot quotread. zooquot quotrev. zooot 82 quotrollapplyquot quotrollapplyrquot quotrollmaxquot 85 quotrollmax. defaultquot quotrollmaxrquot quotrollmeanquot 88 quotrollmean. defaultquot quotrollmeanrquot quotrollmedianquot 91 quotrollmedian. defaultquot quotrollmedianrquot quotrollsumquot 94 quotrollsum. defaultquot quotrollsumrquot quotscalexyearmonquot 97 quotscalexyearqtrquot quotscaleyyearmonquot quotscaleyyearqtrquot 100 quotSys. yearmonquot quotSys. yearqtrquot quottimelt-quot 103 quotwrite. zooquot quotxblocksquot quotxblocks. defaultquot 106 quotxtfrm. zooquot quotyearmonquot quotyearmontransquot 109 quotyearqtrquot quotyearqtrtransquot Quotzooot 112 quotzooregquot Be Sociable, Share