Bookmark and Share

Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.


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

RE: st: Graphic displays or results from margins


From   "Martin Weiss" <[email protected]>
To   <[email protected]>
Subject   RE: st: Graphic displays or results from margins
Date   Mon, 20 Sep 2010 19:50:21 +0200

<>

Not sure whether I am being helpful, but Roger`s -ssc d parmest-is also
available after -mergins, post-:


*************
cap which parmest
if _rc ssc inst parmest

webuse nhanes2f, clear
logit diabetes i.black i.female age i.female#c.age, nolog
margins female#black, at(age = (20(10)70)) post

parmest, saving(myfile, replace)

u myfile, clear
l
*************


HTH
Martin


-----Original Message-----
From: [email protected]
[mailto:[email protected]] On Behalf Of Phil Schumm
Sent: Montag, 20. September 2010 17:31
To: [email protected]
Subject: Re: st: Graphic displays or results from margins

On Sep 20, 2010, at 10:32 AM, Richard Williams wrote:
> I am warming up to the new margins command.  However, unlike the  
> older adjust command, margins does not include the ability to  
> generate new variables. I am therefore finding it hard to create  
> graphic displays of the results.
>
> Consider the following:
>
> webuse nhanes2f, clear
> logit diabetes i.black i.female age i.female#c.age, nolog
> margins female#black, at(age = (20(10)70))
>
> I would love to create a graphic where the Y axis was the  
> probability of diabetes, the X axis was age, and there were 4  
> separate lines, for black males, black females, nonblack males,  
> nonblack females.


I believe that the only option currently is to go manually through  
e(b).  For example, if you add the -post- option to your -margins-  
command above, you could do something like this:


     mat A = e(b)'
     svmat A
     egen _age = fill(10 10 10 10 20 20 20 20)
     egen _grp = fill(1 2 3 4 1 2 3 4)
     separate A1, by(_grp)
     line A11-A14 _age if !mi(A1)


It wouldn't be difficult to write a little command to do this in a  
more general way (indeed, someone may have done it already).


-- Phil

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/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/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/


© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index