Um método simples e geral para preencher os dados em falta, se tiver execuções de dados completos, é usar regressão linear. Digamos que você tem 1000 execuções de 5 em uma linha com nenhuma falta. Configurar o vetor 1000 x 1 y e 1000 x 4 matriz X: regressão lhe dará 4 números a b c d que dão uma melhor correspondência para suas 1000 linhas de dados mdash dados diferentes, diferentes a b c d. Então você usa estas a b c d para estimar (prever, interpolar) wt0 ausente. (Para os pesos humanos, eu esperava que abcd fosse em torno de 14.) (Há centenas de livros e artigos sobre regressão, em todos os níveis. Para a conexão com a interpolação, porém, eu não sei de uma boa introdução ninguém) Arnaud Dowkiw escreveu: Eu tenho um arquivo de dados (x, y, z) onde xey são coordenadas espaciais e variável za. Eu tenho alguns valores em falta na coluna z e eu gostaria de substituí-los com uma estimativa otimizada a partir das células vizinhas. Eu não poderia encontrar qualquer função em R para fazer isso. Alguém sabe de tal função Talvez alguém saiba como implementar o método Papadakis em R. Obrigado pela sua ajuda, não tenho certeza sobre o método Papadakis, eu não acho que ele está lá, mas poderia ser implementado usando ferramentas padrão. Ripley (1981) relaciona o método aos modelos SAR e CAR, mas prever com eles está longe de ser fácil (há artigos de Martin a partir da década de 1980 sobre isso). São os métodos geoestatísticos uma opção para você Que é, o que é o apoio de suas células - - - - - - - - - - - - - - - - - - - - - Arnaud DOWKIW Unidade INRA. Amlioration Gntique et Physiologie Forestires 2163 avenue de la Pomme de Pin BP 20619 ARDON 45166 OLIVET CEDEX FRANÇA Tel. 33 2 38 41 78 00 Fax. 33 2 38 41 48 09 Roger Bivand Secção de Geografia Económica, Departamento de Economia, Escola Norueguesa de Economia e Administração de Empresas, Helleveien 30, N-5045 Bergen, Noruega. Voz: 47 55 95 93 55 fax 47 55 95 95 43 e-mail: Roger. Bivand em nhh. noCalculando Médias Móveis com Valores Faltos set obs 10 obs foi 0, agora 10 tsset t variável de tempo: t, 1 a 10 delta: 1 tssmooth ma y3 y, w (1 1 1) O mais suave aplicado foi (13) x (t-1) 1x (t) x (t1) x (t) y tssmooth ma y5 y, w (2 1 2) O mais liso aplicado foi (15) x (t-2) x (t-1) 1x (t) x (t1) x (t2) x (t) y-lisa mapeia o melhor que pode nas extremidades da série (O que não é o mesmo que tratar falhas como 0s, caso contrário os resultados acima não seriam obtidos). Jeff pode substituir os resultados de - ssmooth - se ele não quer isso. Substitua y3. Se ausente (L1.y, F1.y) substituir y5. Se ausente (L1.y, L2.y, F1.y, F2.y) Por exemplo, L1.y está ausente no início de uma série y e F1.y no final. Isso também funciona com dados de painel. Uma alternativa escrita por usuário para - tssmooth - é a função - egen - filter () - de - egenmore-. Em alternativa, especifique qualquer forma de lidar com os fins que preferir. Quando uso a função tssmooth para criar médias móveis de 2 ou 5 anos, a Stata calcula estas médias móveis até ao final da série temporal em vez de parar 1 ano (no caso da média móvel de 2 anos) ou 4 anos Caso da média móvel de 5 anos) antes do final da série temporal. (C (4,5,4,6), 3) Séries Temporais: Início 1 Fim 4 Freqüência 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 assim que eu esperava obter apenas dois números de volta 8211 4,333333 e 5 8211 e se houvesse ser NA valores eu pensei que 8217d estar no início da seqüência. Na verdade, verifica-se isto é 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 os valores passados apenas se os lados 2 eles são centrados em torno de atraso 0. Neste caso, o comprimento do filtro deve ser ímpar, mas se é mesmo, mais do filtro é para a frente no tempo do que para trás. Assim, na nossa função 8216mav8217, a média de rolamento considera ambos os lados do valor atual em vez de apenas 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 para I8217ll ser varredura zoo8217s lista de funções próxima vez que eu preciso fazer algo série de tempo relacionados 8211 there8217ll provavelmente já ser uma função para ele 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 quotcoredata. defaultquot quotcoredatalt-quot 19 quotfacetfreequot quotformat. yearqtrquot quotfortify. zooquot 22 quotfrequencylt-quot quotifelse. zooquot quotindexquot 25 quotindexlt-quotindex2charquot quotis. regularquot 28 quotis. zooquot quotmake. par. listquot q UotMATCHquot 31 quotMATCH. defaultquot quotMATCH. timesquotquest. medicam. zooquot 34 quotmerge. zooquot quotna. aggregatequot quotna. aggregate. defaultquot 37 quotna. approxquot quotna. approx. defaultquot quotna. fillquot 40 quotna. fill. defaultquotquotna. locfquot quotna. locf. defaultquot 43 Quotna. splinequot quotna. spline. defaultquot 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. zooquot 61 quotpanel. polygon. itsquot quotpanel. polygon. tisquot quotpanel. polygon. tsquot 64 Quotpanel. polygon. zooquot quotpanel. rect. itsquot quotpanel. rect. tisquot 67quotpanel. rect. tsquotquotpanel. rect. zooquotquotpanel. segments. itsquot 70quotpanel. segments. tisquotquopanel. segments. tsquotquotpanel. se O. Quotquot. Quotquot. Quotquot. Quotquot. Quotquot. Quotquot. Quotquot. Quotquot. Quotquot. Quotquot. Quotquot. 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 Quotzooquot 112 quotzooregquot Be Sociable, Share
No comments:
Post a Comment