Notice: On March 31, it was **announced** that Statalist is moving from an email list to a **forum**. The old list will shut down on April 23, and its replacement, **statalist.org** is already up and running.

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

From |
<michael.beck@edu.unibe.ch> |

To |
<statalist@hsphsun2.harvard.edu> |

Subject |
AW: st: Conditional Effect Plots after mi estimate |

Date |
Sun, 17 Feb 2013 13:00:17 +0000 |

Thank you very much for the answer. Actually the margins and marginsplot commands do not work after MI. But after a little search I found this website: http://www.ats.ucla.edu/stat/stata/faq/ologit_mi_marginsplot.htm It contains a manual for a workaround, which actually after a few modifications worked out fine for me. Here is my version of the code: use mi1502.dta program emargins, eclass properties(mi) version 12 args outcome logit uebertritt zh sex i.egp4 migerf i.mig schnitt d_cost inc, vce(cluster klasse) margins egp4, at (schnitt=(4(.1)6)) atmeans asbalanced post end program myret, rclass return add return matrix b = b return matrix V= V end mi estimate: logit uebertritt zh sex i.egp4 migerf i.mig schnitt d_cost inc, vce(cluster klasse) mi estimate, cmdok: emargins uebertritt // emargins is defined above mat b= e(b_mi) // save mi point estimates mat V = e(V_mi) // save mi vce * run ologit and margins on the _mi_m==0 data quietly logit uebertritt zh sex i.egp4 migerf i.mig schnitt d_cost inc, vce(cluster klasse) quietly margins egp4, at (schnitt=(4(.1)6)) atmeans asbalanced post myret // myret is defined above mata: st_global("e(cmd)", "margins") // set previous cmd to margins marginsplot, x(schnitt) recast(line) name(logit, replace) xlabel(4 4.5 5 5.5 6) ylabel(0 .25 .50 .75 1) ________________________________________ Von: owner-statalist@hsphsun2.harvard.edu [owner-statalist@hsphsun2.harvard.edu]" im Auftrag von "Ulrich Kohler [ukohler@uni-potsdam.de] Gesendet: Sonntag, 17. Februar 2013 11:40 An: statalist@hsphsun2.harvard.edu Betreff: Re: st: Conditional Effect Plots after mi estimate I haven't tried it, but my guess would be that the official Stata commands -margins- and -marginsplot- work after -mi-: . sysuse auto . logit for weight length . margins, at(weight =(1760(100)4840) length=(150(25)225)) . marginsplot, by(length) Am Samstag, den 16.02.2013, 11:22 +0000 schrieb michael.beck@edu.unibe.ch: > Dear Stata List Community > > After perfoming a Logistic Regression on Educational Transitions like the following: > > logit y x1 x2 x3 ... xi, vce(cluster school) > > where y is the choice between two different Educational Tracks and xi are the predictor variables (like social class, GPA and migration background), i want to produce conditional effect plots, that are showing the change in transition probabilities for indivuals given certain characteristics for different values of GPA. Thanks to the command 'prgen' this ist pretty easy, my command for that looks like this: > > prgen gpa, from (4) to (6) gen (pr1) x (x1 =0 x2=1) rest (mean) > prgen gpa, from (4) to (6) gen (pr2) x (x1 =1 x2=1) rest (mean) > > twoway (scatter pr1p1 pr2p1 pr2px) > > Now if I want to do the same with imputed data, the command prgen doesn't work with results from regressions like: > > mi estimate: logit y x1 x2 x3 ... xi, vce(cluster school) > > There seem to be solutions for similar problems, but not the same problem (e.g. for svy data) and the more I'm searching, the more confused I actually get. Is there a "simple" solution for plotting Conditional Probabilities after mi estimate: logit? > > Thanks a lot > Michael > * > * For searches and help try: > * http://www.stata.com/help.cgi?search > * http://www.stata.com/support/faqs/resources/statalist-faq/ > * http://www.ats.ucla.edu/stat/stata/ * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/faqs/resources/statalist-faq/ * http://www.ats.ucla.edu/stat/stata/ * * For searches and help try: * http://www.stata.com/help.cgi?search * http://www.stata.com/support/faqs/resources/statalist-faq/ * http://www.ats.ucla.edu/stat/stata/

**References**:**st: Conditional Effect Plots after mi estimate***From:*<michael.beck@edu.unibe.ch>

**Re: st: Conditional Effect Plots after mi estimate***From:*Ulrich Kohler <ukohler@uni-potsdam.de>

- Prev by Date:
**Re: st: Nonlinear constraints in -ml- programming** - Next by Date:
**Re: st: Need help with regular expression** - Previous by thread:
**Re: st: Conditional Effect Plots after mi estimate** - Next by thread:
**st: Nonlinear constraints in -ml- programming** - Index(es):