Grafiken aus mehreren Margen in einem marginsplot in Stata

Ich würde gerne zeichnen Margen, die aus der margins Befehl in eine marginsplot, aber aus verschiedenen margins Schätzungen. Wichtige Einschränkung: Diese Koeffizienten sind innerhalb der gleichen min-und max-und daher vergleichbar. Wie mache ich das?

Hier ist ein code Beispiel:

webuse nhanes2, clear

tnbreg psu weight hdresult iron, iterate(5) //I am using this regression type so I stick with it here

Weiß ich, dass ich alle Ränder Reaktion Graphen in einem plot

margins, dydx(*)
marginsplot, horizontal xline(0) yscale(reverse) recast(scatter)

Aber in der Tat, ich bin mit drei margins Befehl für jeden der Regressoren getrennt, weil ich möchte, zu vergleichen, die Auswirkungen, wenn das regressor würde variieren. Der code ist daher

foreach var in weight hdresult iron {
  * Procedure to get the numbers for margins right
  quietly summarize `var '
  local max = r(max)
  local step = round(r(max)/6)

  quietly margins, at(`cvar'=(1(`step')`max'))
  marginsplot, title("") ytitle("")
}

Dies gibt mir drei separaten Dateien. Aber ich möchte, dass alle Zeilen in einer einzigen Figur, in verschiedenen Farben, natürlich.

Irgendwelche Vorschläge, wie das zu tun?

InformationsquelleAutor MERose | 2014-09-24

  1. 5

    Mit combomarginsplot (und von der Hilfe-Datei):

    sysuse auto, clear

oprobit rep78 i.foreign mpg price weight
margins foreign, at(mpg=(10(5)50)) predict(outcome(3)) saving(file1, replace)

oprobit rep78 i.foreign mpg
margins foreign, at(mpg=(10(5)50)) predict(outcome(3)) saving(file2, replace)

oprobit rep78 i.foreign mpg gear
margins foreign, at(mpg=(10(5)50)) predict(outcome(3)) saving(file3, replace)

combomarginsplot file1 file2 file3, ///
    labels("Full model" "Restricted model" "Gear Model") noci

    combomarginsplot ist ein user-written command von Nicholas Winter. Man kann es installieren, läuft 

    ssc install combomarginsplot

    InformationsquelleAutor Roberto Ferrer

  2. 1

    Basierend auf @RobertoFerrer Vorschlag zur Verwendung combomarginsplot ich bin jetzt trickst, Paket (danke an Nicholas Winter):

    webuse nhanes2, clear

* Run regressions
foreach var in weight hdresult iron {
  * Trick: always regress on the same variable
  gen testvar = `var'

  * Any regression where testvar enters first - the identical variable will be omitted
  tnbreg psu ///
     testvar weight hdresult iron, iterate(5)

  * Procedure to get the numbers for margins right
  quietly summarize testvar
  local max = r(max)
  local step = round(r(max)/6)

  * Margins post estimation
  quietly margins, at(testvar=(1(`step')`max')) saving(margins_`var', replace)

  * Drop testvar so that it can be reassigned within the loop
  drop testvar
}

* Combine the margins graph information
combomarginsplot margins_weight margins_hdresult margins_iron, labels("Weight" "HDrestul" "Iron")

    Natürlich macht es nur Sinn zu vergleichen, die Koeffizienten der Variablen, die alle im gleichen Bereich. Diese Einschränkung war nicht Teil meiner ursprünglichen Antwort - sorry dafür.

    InformationsquelleAutor MERose

Schreibe einen Kommentar