R moving average missing values

Prostą i ogólną metodą uzupełniania brakujących danych, jeśli masz kompletne dane, jest użycie regresji liniowej. Załóżmy, że masz 1000 przebiegów po 5 z rzędu, których nie brakuje. Skonfiguruj macierz 1000 x 1 yi matrycę 1000 x 4 X: Regresja da ci 4 cyfry a b c d, które najlepiej pasują do 1000 wierszy danych, mdasz inne dane, inne a b c d. Następnie używa się tych b c d do oszacowania (przewidywania, interpolacji) brakującego wt0. (W przypadku ludzkich ciężarów, Spodziewam się, że abcd będzie w okolicy 14.) (Na wszystkich poziomach znajduje się miliard książek i artykułów na temat regresji, jednak w odniesieniu do związku z interpolacją nie wiem o dobrym wprowadzeniu) R Zamień brak wartości w projektowaniu przestrzennym za pomocą średniej ruchomej W środę, 9 listopada 2005 r., Arnaud Dowkiw napisał: Mam plik danych (x, y, z), w którym xiy są współrzędnymi przestrzennymi i zmienną za. Mam kilka brakujących wartości w kolumnie z i chciałbym zastąpić je zoptymalizowaną estymacją z sąsiednich komórek. Nie mogłem znaleźć żadnej funkcji w R, aby to zrobić. Czy ktoś jest świadomy takiej funkcji? Może ktoś wie, jak wdrożyć metodę Papadakis w R. Dzięki za twoją pomoc, nie jestem pewien co do metody Papadakisa, nie sądzę, że ona istnieje, choć mogłaby być zaimplementowana przy użyciu standardowych narzędzi. Ripley (1981) opisuje tę metodę w modelach SAR i CAR, ale przewidywanie z nimi jest daleka od łatwości (są to prace Martina z lat 40. XX wieku na ten temat). Czy metody geostatystyczne są dla ciebie opcją? To jest, co jest wsparciem twoich komórek - - - - - - - - - - - - - - - - - - - - Jednostka Arnaud DOWKIW INRA. Amlioration Gntique et Physiologie Forestires 2163 avenue de la Pomme de Pin BP 20619 ARDON 45166 OLIVET CEDEX FRANCE Tel. 33 2 38 41 78 00 Faks. 33 2 38 41 48 09 - - - - - - - - - - - - - - - - - - - - Roger Bivand Wydział Geografii Ekonomicznej, Wydział Ekonomii, Norweska Szkoła Ekonomii i Administracji Biznesowej, Helleveien 30, N-5045 Bergen, Norwegia. głos: 47 55 95 93 55 faks 47 55 95 95 43 e-mail: Roger. Bivand w nhh. no Obliczanie średnich kroczących z brakującymi wartościami zestaw obs 10 obs było 0, teraz 10 tsset t zmienna czasowa: t, 1 do 10 delta: 1 jednostka tssmooth ma y3 y, w (1 1 1) Zastosowano gładszą (13) x (t-1) 1x (t) x (t1) x (t) y tssmooth ma y5 y, w (2 1 2) Zastosowano gładszą (15) x (t-2) x (t-1) 1 x (t) x (t1) x (t2) x (t) y - gładkie wypełnienia najlepiej jak to możliwe na końcach serii (co nie jest tym samym, co traktowanie błędów jako 0, ponieważ w przeciwnym razie powyższe wyniki nie zostałyby uzyskane). Jeff może zignorować wyniki - pociechy - jeśli tego nie chce. zastąp y3. jeśli brakuje (L1.y, F1.y), zastąp y5. w przypadku braku (L1.y, L2.y, F1.y, F2.y) Na przykład brakuje L1.y na początku serii yi F1.y na końcu. Działa to również z danymi panelu. Alternatywą dla - tssmooth napisaną przez użytkownika jest - egen-function - filter () - from - egenmore-. Alternatywnie, proszę przeliterować jakikolwiek sposób radzenia sobie z końcami, które preferujesz. Kiedy używam funkcji tssmooth do tworzenia średnich kroczących z 2 lub 5 lat, Stata oblicza te ruchome średnie aż do końca szeregu czasowego zamiast zatrzymywania 1 roku (w przypadku 2-letniej średniej ruchomej) lub 4 lat (w przypadek 5-letniej średniej kroczącej) przed końcem szeregów czasowych. Tym samym zlicza do przodu brakujące wartości jako zero. Czy znasz możliwą strategię naprawienia tego problemu mav (c (4,5,4,6), 3) Szeregi czasowe: Start 1 Koniec 4 Częstotliwość 1 1 NA 4.333333 5.000000 NA Tutaj próbowałem zrobić średnią kroczącą, która zajęła pod uwagę ostatnie 3 cyfry, więc spodziewałem się uzyskać tylko dwie liczby z powrotem 8211 4.333333 i 5 8211 i gdyby istniały wartości NA, pomyślałem, że są na początku sekwencji. W rzeczywistości okazuje się, że to właśnie kontroluje parametr 8216sides8217: boki tylko dla filtrów splotowych. Jeśli boki 1 współczynniki filtra odnoszą się do wartości z przeszłości tylko wtedy, gdy boki 2 są wyśrodkowane wokół opóźnienia 0. W tym przypadku długość filtra powinna być nieparzysta, ale jeśli jest równa, większa część filtra jest przesunięta do przodu w czasie niż do tyłu. Tak więc w naszej funkcji 8216mav8217 średnia krocząca wyszukuje obie strony bieżącej wartości, a nie tylko poprzednie wartości. Możemy to zmienić, aby uzyskać pożądane zachowanie: gt biblioteka (zoo) gt rollmean (c (4,5,4,6), 3) 1 4,333333 5.000000 Zrozumiałem też, że mogę wymienić wszystkie funkcje w pakiecie z 8216ls8217 funkcja, więc I8217ll będzie skanowaniem zoo8217s lista funkcji następnym razem muszę zrobić coś z szeregów czasowych związanych 8211 there8217ll prawdopodobnie już będzie funkcją dla niego gt ls (quotpackage: zooquot) 1 przydziały. Dataquot quas. Date. numericquot quotas. Date. tsquot 4 quotas. Date. yearmonquot quotas. Date. yearqtrquot quotas. yearmonquot 7 quotas. yearmon. defaultquot. yearqtrquot quotas. yearqtr. defaultquot 10 kwot. zooquot quotas. zoo. defaultquot quotas. zooregquot 13 quotas. zooreg. defaultquot quotautoplot. zooquot quotcbind. zooquot 16 quotcoredataquot quotcoredata. defaultquot quotoredatt-quot 19 quotacacetfreequot quotformat. yearqtrquot quotfortify. zooquot 22 quotfrequencylt-quot quotifelse. zooquot quotindexquot 25 quotindexlt-quot quotindex2charquot quotis. regularquot 28 quotis. zooquot quotmake. par. listquot q uotMATCHquot 31 quotMATCH. defaultquotMATCH. timesquot quotmedian. zooquot 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 quotERER. 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 67 quotpanel. rect. tsquot quotpanel. rect. zooquot quotpanel. segments. itsquot 70 quotpanel. segments. tisquot quotpanel. segments. tsquot quotpanel. se gments. zooquot 73 quotpanel. text. itsquot quotpanel. textquit. tablet. text. tsquot 76 quotpanel. text. zooquot quotplot. zooquot quotquantile. zooquot 79 quotrbind. zooquot quotread. zooquot quotrev. zooquot 82 quotrollapplyquot quotrollapplyrquot quotrollmaxquot 82 quotrollmax. defaultquotmaxrquot 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, Udostępnij