Eu lutava pela busca de uma função simples para mover médias que tinham alguma flexibilidade para fazer o que eu precisava. Eu finalmente escrevi algumas funções estendendo a baseada na função de filtro que o rinni dá acima no comentário (mas que não funcionará porque incluirá a observação atual na média de 3 períodos). Função de média móvel que inclui a observação atual Função de média móvel que não inclui a observação atual Função média móvel de visualização para trás, não incluindo obss atuais, com base em leituras h2 iniciando os períodos h1 retornados resposta 24 de agosto 16 às 2:25 Sua resposta 2017 Stack Exchange , A média móvel IncR e a biblioteca de valores de NA (zoológico) x lt-zoo (1:10) x5 lt - NA rollapply (x, 3, mean, na. rm TRUE) 2 3 4 5 6 7 8 9 2.0 3.0 3.5 5.0 6.5 7.0 8.0 9.0 xm lt - rollapply (x, 3, mean, na. rm TRUE) xm 2 3 4 5 6 7 8 9 2.0 3.0 3.5 5.0 6.5 7.0 8.0 9.0 coredata (xm) unzoo it 1 2.0 3.0 3.5 5.0 6.5 7.0 8.0 9.0 Veja as duas vinhetas do zoológico para obter mais informações. Em 10 de dezembro de 2007 6:45 da manhã, Cornelis de Gier escreveu: A função S-plus moving. ave (data, span 2) calcula a média móvel, mas não possui um argumento para dizer como lidar com os valores de NA , Então retornará NA para todas as médias, como mostrado abaixo. Existe uma função de média móvel R ou S que é capaz de omitir alguns valores de NA no conjunto de dados. Na amostra simples mostrada abaixo, seria possível simplesmente remover as linhas com valores de NA. O conjunto de dados no qual eu quero usar a função de média móvel com um intervalo de 270 é um conjunto de dados de séries temporais, apenas remover linhas corromperam este conjunto de dados e torná-lo impróprio para plotar. T lt - (1:10) moving. ave (t, 2) 36aves: 1 1.0 1.5 2.5 3.5 4.5 5.5 6.5 7.5 8.5 9.5 36sizes: 1 1 2 2 2 2 2 2 2 2 2 t5 lt - NA moving. ave ( T, 2) 36aves: 1 NA NA NA NA NA NA NA NA NA NA 36sizes: 1 1 2 2 2 2 2 2 2 2 2 R-help na lista de correspondência r-project. org stat. ethz. chmailmanlistinfor-help POR FAVOR faça Leia o guia de publicação R-project. orgposting-guide. html e forneça o código comentável, mínimo, autônomo e reprodutível. gt mav (c (4,5,4,6), 3) Série de tempo: Iniciar 1 Freqüência de fim de 4 1 1 NA 4.333333 5.000000 NA Aqui eu estava tentando fazer uma média móvel que levou em 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 NA eu pensei que eles deveriam ser 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
No comments:
Post a Comment