Thursday 28 September 2017

Moving Average Mean Stata


Diese Datenstruktur ist unzweckmäßig. Angenommen, eine ID-ID müssen Sie umformen. z. B. Dann ist ein gleitender Durchschnitt einfach. Verwenden Sie tssmooth oder einfach nur generieren. z. B. Mehr darüber, warum Ihre Datenstruktur ist ziemlich unfit: Nicht nur würde die Berechnung eines gleitenden Durchschnitt benötigen eine Schleife (nicht unbedingt mit egen), aber Sie würden mehrere neue zusätzliche Variablen erstellen. Mit denen in einer späteren Analyse wäre irgendwo zwischen ungeschickt und unmöglich. EDIT Ill geben eine Probe-Schleife, während sich nicht aus meiner Haltung, dass es schlechte Technik. Ich sehe nicht einen Grund für Ihre Namenskonvention, wobei P1947 ist ein Mittel für 1943-1945 Ich nehme an, das ist nur ein Tippfehler. Nehmen wir an, dass wir Daten für 1913-2012 haben. Für Mittel von 3 Jahren verlieren wir ein Jahr an jedem Ende. Das könnte prägnanter geschrieben werden, auf Kosten einer Flut von Makros innerhalb von Makros. Mit ungleichen Gewichten ist einfach, wie oben. Der einzige Grund, egen verwenden ist, dass es nicht aufgeben, wenn es Versäumnisse, die die oben genannten tun wird. Der Vollständigkeit halber sei bemerkt, dass es leicht ist, Fehler zu behandeln, ohne auf egen zurückzugreifen. Und der Nenner Wenn alle Werte fehlen, reduziert sich dies auf 0/0 oder fehlt. Andernfalls, wenn ein Wert fehlt, fügen wir 0 auf den Zähler und 0 auf den Nenner, die die gleiche wie Ignorieren ist. Natürlich ist der Code erträglich wie oben für Mittelwerte von 3 Jahren, aber entweder für diesen Fall oder für die Mittelung über mehr Jahre, würden wir die Zeilen oben durch eine Schleife, die egen tut. Stata: Datenanalyse und statistische Software Nicholas J () Und ihre Einschränkungen Statarsquos offensichtlichste Befehl zur Berechnung der gleitenden Mittelwerte ist die ma () - Funktion von egen. Bei einem Ausdruck wird ein gleitender Durchschnitt für diesen Ausdruck erstellt. Standardmäßig wird als 3. genommen, muss ungerade sein. Allerdings kann, wie der manuelle Eintrag angibt, egen, ma () nicht mit varlist kombiniert werden:. Und aus diesem Grund ist es nicht auf Paneldaten anwendbar. In jedem Fall steht er außerhalb des Satzes von Befehlen, die speziell für Zeitreihen geschrieben werden, siehe Zeitreihen für Details. Alternative Ansätze Zur Berechnung von Bewegungsdurchschnitten für Paneldaten gibt es mindestens zwei Möglichkeiten. Beide hängen davon ab, dass der Dataset zuvor tsset wurde. Das ist sehr viel wert: nicht nur können Sie sich immer wieder spezifizieren Panel variabel und Zeit variabel, aber Stata verhält sich intelligent jede Lücken in den Daten. 1. Schreiben Sie Ihre eigene Definition unter Verwendung von Zeitreihenoperatoren wie L. und F. Geben Sie die Definition des gleitenden Durchschnitts als Argument für eine generierte Anweisung an. Wenn Sie dies tun, sind Sie natürlich nicht auf die gleich gewichteten (ungewichteten) zentrierten Bewegungsdurchschnitte beschränkt, die von egen, ma () berechnet wurden. Zum Beispiel würden gleich gewichtete Dreiphasenbewegungsdurchschnitte gegeben und einige Gewichte können leicht angegeben werden: Sie können natürlich einen Ausdruck wie log (myvar) anstelle eines Variablennamens wie myvar angeben. Ein großer Vorteil dieses Ansatzes ist, dass Stata automatisch das Richtige für Paneldaten macht: führende und nacheilende Werte werden in Panels ausgearbeitet, genauso wie Logik diktiert. Der bemerkenswerteste Nachteil ist, dass die Befehlszeile ziemlich lang werden kann, wenn der gleitende Durchschnitt mehrere Begriffe beinhaltet. Ein anderes Beispiel ist ein einseitiger gleitender Durchschnitt, der nur auf vorherigen Werten basiert. Dies könnte nützlich sein für die Erzeugung einer adaptiven Erwartung dessen, was eine Variable nur auf Informationen basieren wird: was könnte jemand prognostizieren für den aktuellen Zeitraum auf der Grundlage der letzten vier Werte, mit einem festen Gewichtungsschema (A 4-Periode Verzögerung sein könnte Besonders gebräuchlich mit vierteljährlichen Zeitreihen.) 2. Verwenden Sie egen, filter () von SSC Verwenden Sie den benutzerdefinierten egen function filter () aus dem egenmore package auf SSC. In Stata 7 (aktualisiert nach dem 14. November 2001) können Sie dieses Paket installieren, nachdem egenmore auf die Details zu filter () hingewiesen hat. Die beiden obigen Beispiele würden gerendert (In diesem Vergleich ist der generierte Ansatz vielleicht transparenter, aber wir sehen ein Beispiel des Gegenteils in einem Moment.) Die Lags sind eine Numliste. Führt zu negativen Verzögerungen: In diesem Fall erweitert sich -1/1 auf -1 0 1 oder Blei 1, lag 0, lag 1. Die Koeffizienten, eine weitere Numliste, multiplizieren die entsprechenden nacheilenden oder führenden Elemente: In diesem Fall sind diese Elemente F1.myvar Myvar und L1.myvar. Der Effekt der Normalisierungsoption besteht darin, jeden Koeffizienten durch die Summe der Koeffizienten zu skalieren, so daß coef (1 1 1) normalisiert ist, zu Koeffizienten von 1/3 1/3 1/3 äquivalent ist und coef (1 2 1) normalisiert Auf Koeffizienten von 1/4 1/2 1/4. Sie müssen nicht nur die Verzögerungen, sondern auch die Koeffizienten angeben. Da egen, ma () den gleich gewichteten Fall liefert, ist der Hauptgrund für egen, filter (), den ungleich gewichteten Fall zu unterstützen, für den Sie Koeffizienten angeben müssen. Es könnte auch gesagt werden, dass die verpflichtenden Benutzer, um Koeffizienten angeben ist ein wenig mehr Druck auf sie zu denken, welche Koeffizienten sie wollen. Die wichtigste Rechtfertigung für gleiche Gewichte ist, wir schätzen, Einfachheit, aber gleiche Gewichte haben miese Frequenzbereich Eigenschaften, um nur eine Erwägung zu erwähnen. Das dritte Beispiel oben könnte entweder von denen ist nur so kompliziert wie die Generierung Ansatz. Es gibt Fälle, in denen egen, filter () eine einfachere Formulierung ergibt als erzeugen. Wenn Sie einen neun-term-Binomialfilter suchen, der von den Klimatologen als nützlich empfunden wird, dann sieht es vielleicht weniger schrecklich aus und ist leichter zurecht zu kommen. Genau wie beim generierten Ansatz funktioniert egen, filter () ordnungsgemäß mit Panel-Daten. Tatsächlich hängt es, wie oben erwähnt, davon ab, daß der Dataset vorher tsset wurde. Eine grafische Spitze Nach der Berechnung Ihrer gleitenden Durchschnitte werden Sie wahrscheinlich einen Graphen betrachten wollen. Der benutzerdefinierte Befehl tsgraph ist schlau um Tsset-Datasets. Installieren Sie es in einem up-to-date Stata 7 von ssc inst tsgraph. Was ist mit der Teilmenge mit if Keine der obigen Beispiele verwenden, wenn Einschränkungen. In der Tat egen, ma () wird nicht zulassen, wenn angegeben werden. Gelegentlich Menschen wollen verwenden, wenn bei der Berechnung der gleitenden Durchschnitte, aber seine Verwendung ist ein wenig komplizierter als es normalerweise ist. Was würden Sie von einem gleitenden Durchschnitt erwarten? Lassen Sie uns zwei Möglichkeiten identifizieren: Schwache Interpretation: Ich möchte keine Ergebnisse für die ausgeschlossenen Beobachtungen sehen. Starke Interpretation: Ich möchte nicht, dass Sie die Werte für die ausgeschlossenen Beobachtungen verwenden. Hier ist ein konkretes Beispiel. Angenommen, infolge einer Bedingung sind die Beobachtungen 1-42 enthalten, aber nicht die Beobachtungen 43 an. Aber der gleitende Durchschnitt für 42 wird unter anderem von dem Wert für die Beobachtung 43 abhängen, wenn der Mittelwert sich nach hinten und vorne erstreckt und eine Länge von mindestens 3 hat, und er wird in einigen Fällen von einigen der Beobachtungen 44 abhängen. Unsere Vermutung ist, dass die meisten Menschen für die schwache Interpretation gehen würde, aber ob das korrekt ist, egen, filter () nicht unterstützt, wenn entweder. Sie können immer ignorieren, was Sie donrsquot wollen oder sogar unerwünschte Werte auf fehlende danach mit replace setzen. Eine Notiz über fehlende Ergebnisse an den Enden der Serie Da gleitende Mittelwerte Funktionen von Lags und Leads sind, erzeugt eMe () fehlende Stellen, wo die Lags und Leads nicht existieren, am Anfang und Ende der Reihe. Eine Option nomiss zwingt die Berechnung der kürzeren, nicht beanspruchten gleitenden Mittelwerte für die Schwänze. Im Gegensatz dazu weder erzeugen noch egen, filter () macht oder erlaubt, etwas Besonderes, um fehlende Ergebnisse zu vermeiden. Wenn einer der für die Berechnung benötigten Werte fehlt, fehlt dieses Ergebnis. Es ist bis zu den Benutzern zu entscheiden, ob und welche Korrekturchirurgie für solche Beobachtungen erforderlich ist, vermutlich nach dem Betrachten des Datensatzes und unter Berücksichtigung jeder zugrunde liegenden Wissenschaft, die gebracht werden kann, zu tragen. Announcement Es ist mein erster Beitrag und ich werde versuchen, so klar zu sein wie möglich. Der Link für die Hauptdatenbank befindet sich am Ende des Beitrags. Kontext Ich benutze Stata / SE 12.0 unter Windows 10. Ich habe mit Stata erst vor ein paar Wochen begonnen und ich versuche, selbstständig für einen Auftrag zu lernen, der in ein paar Tagen fällig wird (weil jeder Tisch oder jede Figur mir Tage und Stunden brachte) Tage): Replikation der Papier-QuoteDoes obligatorische Schulbesuch Affect Schoolig und Earningsquot: web. stanford. edu/pista / angrist. pdf, dass die Menschen geboren in den letzten Quartalen der Jahre haben mehr Bildung im Durchschnitt als die in der ersten geboren Aufgrund der Pflichtschulgesetze. In den ersten Zahlen wird für alle Personen, die in einem bestimmten Jahr (Variable YOB für Geburtsjahr) in einem bestimmten Quartal (QOB) geboren wurden, ein Schaubild der durchschnittlichen Erziehungsdauer (Variable EDUC) erstellt. Es gibt eine allmähliche Tendenz, und die Daten zu schädigen, verwenden sie einen gleitenden Durchschnitt (Abbildung IV), wo ich seit den letzten 5 Tagen blockiert bin. Problem In der Datenbank gibt es 27 Variablen, darunter v4 umbenannt EDUC, v27 umbenannt YOB (Geburtsjahr), und v18 umbenannt QOB (Viertel der Geburt). Was für den gleitenden Durchschnitt gebraucht wird, ist für alle Personen, die im Jahr c und Viertel j geboren sind, die durchschnittliche Anzahl der Jahre der Bildung nicht für dieses Jahr und Quartal, sondern für das Quartal kurz vor, 2 Quartale vorher, ein Viertel berechnet Später und 2 Quartale später (s. 985 der Zeitung). Wenn ich zum Beispiel die von 1930 bis 1939 geborenen Männer in dieser Figur betrachte (Abbildung IV des Artikels: onedrive. live/redirresi. ntphoto2cpng), muss ich mit der im Jahr 1930, dem 3. Quartal geborenen Kohorte beginnen und die (Jahrgang 1930), 2. Quartal (geb. ein Viertel vor der angegebenen Kohorte), gleiches Geburtsjahr 1930, 1. Quartal (geb. 2 Quartale vor der angegebenen Kohorte), gleiches Geburtsjahr 1930, 4. Quartal (ein Viertel nach der angegebenen Kohorte) und dieselben für die Jahrgänge 1931, 1. Quartal (2 Quartale nach der angegebenen Kohorte). Dann wird der gleitende Durchschnitt durch Addition dieser 4 Werte und Division durch 4 erhalten. Dieser ganze Prozess sollte für jede Kohorte zwischen 1930, 3. Quartal und 1939, 2. Quartal wiederholt werden. Do-Datei für den Do-Datei (onedrive. live/redirresid6919D329B3BF1EF23227ampauthkeyAO2cxEN AGpZMgsMampithintfile2cdo), begann ich mit dem Modell der anderen Figuren und versuchte foreach-Schleife und viele andere Dinge zu verwenden (Sie die Fehlermeldungen nicht erinnern, / wusste nicht, ich wollte Um hier zu posten), aber immer noch nicht herauszufinden, wie zu erzählen Stata: quotfor jeder YOBQn. (EDUC) von YOBQ n-1, YOBQ n-2, YOBQ n1, YOBQ n2quot. Um die Summe und Division durch 4 nach, dass es einfacher sein sollte. Ich habe einen außergewöhnlichen Hinweis auf die Lehr-Assistent: quot versuchen Sie die tssmooth Befehl. Sie werden zunächst eine Zeitvariable zu schaffen haben, für die die egen Gruppe Befehl sehr useful. quot sein wird, aber nach meinen Recherchen über quotegenquot und quottssetquot in den Datenhandbüchern und im Buch Cameron amp Trivedi, quotEconometrics mit Stataquot (letzter Link): Stata / manuals14 / degen. Tfolder2cdta stata / handbücher14 / gsw11.pdf stata / manuals14 / u11.p. SprachenSyntax stata / manuals14 / u13.p. Itsubscripting onedrive. live/redirresi. intfile2cpdf sollte ich die Daten vor tssmooth tsset aber ich habe nicht über dieses Stadium erhalten, da es scheint, ist die Notation n nicht mit quotegenquot erlaubt (Fehler r (101) quotweights nicht allowedquot) und ich bin immer noch sehr verwirrt mit wie egen zu kombinieren, tsset Und tssmooth. Es wäre großartig, wenn jemand könnte mir helfen, mit, wie zu lösen, die quotweights not allowedquot Fehler und wie die Befehle quotegenquot, quottssetquot und quottssmoothquot kombinieren. Danke so viel Postscript: hier ist die Datenbank übrigens onedrive. live/redirresi. Tfolder2cdta Hinweis: Ich habe die do-Datei für die wichtigsten anderen Figuren und Tabellen des Artikels, außer Tabelle I, aber diese Datei ist wahrscheinlich nicht notwendig / nur für info: onedrive. live/redirresi. Hintfile2cdo Zuletzt bearbeitet von Amarylis Durand 25 Mär 2016, 01:55. 25. März 2016, 22.50 Uhr Hier ist eine kürzere Version meiner Frage: Wie kann der Fehler unter (R451 fett gedruckt) zu vermeiden und wie Stata zu sagen, dass der gleitende Durchschnitt von quotmedstay1quot sollte für jeden Wert von tps berechnet werden / von YOB QOB Für alle, die im selben Jahr und im gleichen Quartal geboren wurden, nach Jahren und Vierteln sortieren und die durchschnittliche Anzahl der Jahre der Ausbildung / von YOB QOB berechnen. sort: egen medstay1 Mittelwert (EDUC) / eine neue Variable YOBNew erzeugen, da der Befehl yq das erste Argument erfordert zwischen 1000 und 9999 und unsere Daten für YOB in der Volkszählung von 1980 zu sein, ist zwischen 30 und 49 anstelle von 1930 und 1949 / gen YOBNewYOB YOBNew YOB1900 ersetzen, wenn CENSUS80 / eine Zeitvariable erzeugen, die das Format in Hilfe tsset / gen tpsyq (YOBNew, QOB) Format TQ tps / die folgende Anweisung liefert R451 benötigt hat. Wiederholt Zeitwerte im Panel, wahrscheinlich, weil es Tausende von Menschen geboren im selben Jahr und im gleichen Quartal, offensichtlich mit dem gleichen durchschnittlichen Anzahl von Jahren der Bildung. Wie zu vermeiden, diese Fehler / tsset medstay1 tps / Anweisung haben die gleitende durchschnittliche MA, aber ich möchte die Verschiebung des Mittels von medstay1 für jeden Wert des timevariable quottpsquot berechnet werden. Was Befehl würde mir erlauben, dies zu tun oder ist es automatisch getan / tssmooth ma MA medstay1. Fenster (2 0 2) Ich hoffe, dass jemand helfen kann.

No comments:

Post a Comment