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: contrasts


From   Jeffrey Pitblado <[email protected]>
To   [email protected]
Subject   Re: st: contrasts
Date   Mon, 01 Oct 2012 21:32:24 -0500

Nikolaos Pandis <[email protected]> originally asked about
plotting margin contrasts, but I jumped to pairwise comparisons of
margins when I noticed the third comparison in the original email.
Nick's follow-up comment on my reply was:

> Thank you for the suggestion.

marginsplot after :

margins r.wiretype,contrast

looks neater compared to:

margins wiretype, pwcompare

Nick can get -margins- to plot any number of contrasts of margins,
the trick is employing user-defined contrasts.  Here is a rework
of Nick's original example that also includes the 3 vs 2 contrast:

	margins {wiretype -1  1 0} /// 2 vs 1
		{wiretype -1  0 1} /// 3 vs 1
		{wiretype  0 -1 1}  // 3 vs 2

While user-defined contrasts require more from the user in order to
specify them, user-defined contrasts also give the user the most amount
of control.

--Jeff Pitblado
  [email protected]

PS:

Here the example code I was playing with while thinking about
Nick's query:

* BEGIN
sysuse auto

gen fac = (rep<4) + 2*(rep==4) + 3*(rep==5)

keep if fac

regress mpg i.fac

margins r.fac
marginsplot, name(contrast)

margins fac, pwcompare
marginsplot, name(pwcompare)

margins	{fac -1  1 0}	///
	{fac -1  0 1}	///
	{fac  0 -1 1}
marginsplot, name(user)
* END


*
*   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/


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