Wie schreibt man den code für den median einer variable in stata

Frage ich mich, wie man die codes für die Suche nach dem median für eine variable in stata, ohne mit sort, egen, summarize. Dies ist, was ich so weit gekommen:

capture program drop find_median
program find_median
    local n = _N
    gen ord=0
    forvalues i= 0/`n' {
        replace ord = `i' if [`1']> [`1'][_n-1] & [`1']> [`1'][_n+1] 
    }
end

find_median (the variables...)
  • Beachten Sie, dass Sie Ihre eigenen Schleifen über Beobachtungen zum scheitern verurteilt ist, zu viel, viel langsamer als Stata-internen code.
InformationsquelleAutor kkl | 2013-09-11
  1. 2

    Den centile Befehl haben Sie den median direkt.

    Wenn Sie darauf bestehen, auf die Wiederherstellung des Rades, dann können Sie sort in Ihrem Programm, aber lassen Sie die sortorder der Daten unverändert, nachdem das Programm beendet, indem die sortpreserve option, um die Zeile program find_median finden Sie unter: http://www.stata.com/help.cgi?program. Das sollte das Programm viel einfacher und damit leichter zu Debuggen.

    • Insbesondere müssen Sie "Perzentile myvar", dann haben Sie Zugriff zum "r(c_1)", ist der median, oder nahe genug, um es. In meinem Daten-set, das eigentliche Zentrum Wert ist genau 7.3733, und die Perzentile-basiert median wird berichtet, wie 7.373301. Seltsam, aber nahe genug.
    InformationsquelleAutor Maarten Buis
Schreibe einen Kommentar