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/