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: RE: Plot estimates and ci from two regressions in one graph
From
Nick Cox <[email protected]>
To
[email protected]
Subject
Re: st: RE: Plot estimates and ci from two regressions in one graph
Date
Thu, 21 Jun 2012 09:18:03 +0100
I see here { { } } which looks fine to me.
But the code won't run unless various user-written programs have been
installed previously.
Nick
On Thu, Jun 21, 2012 at 9:13 AM, Muhammad Anees <[email protected]> wrote:
> Thanks for the nice effort,
>
> but the code will work smoothly and only if "}" a used a once.
On Thu, Jun 21, 2012 at 1:06 PM, Nikolaos Kanellopoulos
>> Tirthankar thank you very much for the suggestion to use marginsplot (it never crossed my mind).
>>
>> Since I had something more specific in my mind I tried to utilize the flexibility of Stata and its tw command and I think I found a solution to my problem.
>> If anybody has anything to add in order to improve the following code is more than welcome.
>>
>> Nikos
>>
>>
>> *-----------------------------------------------------------------------------------------------------------------
>>
>> sysuse auto, clear
>> levelsof for, local(a)
>> tempfile file
>>
>> foreach i of local a {
>> sysuse auto, clear
>> count if for==`i'
>>
>> if r(N)>=1 {
>> regress mpg c.(rep78 headroom trunk weight) if for==`i'
>> parmest, norestore
>> egen axis= axis(estimate), label(parm)
>> gen for=`i'
>> capture append using `file'
>> save `file', replace
>> local kk=`i1'+1
>> }
>> }
>>
>> gsort axis -for
>> gen order=_N-_n
>>
>> label define jj 1 "_cons" 3 "rep78" 5 "headroom" ///
>> 7 "weight" 9 "trunk"
>> label value order jj
>>
>> twoway ///
>> (scatter order est if for==1, msymbol(circle)) ///
>> (scatter order est if for==0, msymbol(T)) ///
>> (rcap min95 max95 order, horizontal ), ///
>> ylabel(1(2)9, valuelabel angle(0)) ytitle("") ///
>> legend(order (1 2) label(1 "for==1") label(2 "for==0")) ///
>> xtitle("Coefficients and confidence interval")
>> *-----------------------------------------------------------------------------------------------------------------
From: Tirthankar Chakravarty <[email protected]>
>> -marginsplot- is very powerful for these kinds of problems.
>>
>> *------------------------------------------------------------------------
>> sysuse auto, clear
>> reg mpg i.foreign##c.(rep78 headroom trunk weight)
>> margins, dydx(*) over(foreign)
>> marginsplot, xdim(_deriv)
>> *------------------------------------------------------------------------
>>
>> Note however that without standardizing the variables, you might not
>> be able to make out the variations in coefficients across groups.
>>
>> T
>>
*
* 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/