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]

st: Evaluating the importance of interaction effects in logistic regression


From   Thomas Speidel <[email protected]>
To   <[email protected]>
Subject   st: Evaluating the importance of interaction effects in logistic regression
Date   Tue, 30 Aug 2011 15:31:33 -0600

A recent topic on interaction (http://www.stata.com/statalist/archive/2011-08/msg01455.html), made me wonder what is the correct way to identify important interactions in a logistic model. Maarten in his Stata tip 87

Maarten L. Buis (2010) "Stata tip 87: Interpretation of interactions in non-linear models", The Stata Journal, 10(2), pp. 305-308.)

writes in the context of his first example that the interaction between the two dischotomies black and collgrad "is not significant" (the reported p-value is in fact 0.161). However, another reference Maarten cited:

Edward Norton, Hua Wang, and Chunrong Ai (2004) "Computing interaction effects and standard errors in logit and probit models" The Stata Journal, 4(2): 154-167.

says that "The statistical significance cannot be determined from the z-statistic reported in the regression output" (p.1). I am now confused on the appropriate way of identifying significant interactions. This sentence has confused me. Regardless of interpretation, how does one assess the importance of an interaction?

For example:
********************************************************************************
sysuse nlsw88, clear
gen byte high_occ = occupation < 3 if occupation < .
drop if race==3
logistic high_occ race##collgrad , nolog

Logistic regression Number of obs = 2211 LR chi2(3) = 127.07 Prob > chi2 = 0.0000 Log likelihood = -1199.4399 Pseudo R2 = 0.0503

-------------------------------------------------------------------------------
high_occ | Odds Ratio Std. Err. z P>|z| [95% Conf. Interval]
--------------+----------------------------------------------------------------
2.race | .4194072 .0655069 -5.56 0.000 .3088072 .5696188 1.collgrad | 2.465411 .293568 7.58 0.000 1.952238 3.113478
              |
race#collgrad |
2 1 | 1.479715 .4132536 1.40 0.161 .8559637 2.558003
              |
_cons | .3220524 .0215596 -16.93 0.000 .2824512 .3672059
-------------------------------------------------------------------------------
********************************************************************************

is it correct to say that the interaction race * collgrad is not important because its p-value is 1.161?
What if, for example, we had 3 levels to race:

********************************************************************************
sort idcode
replace race = 3 in 1/300
logistic high_occ race##collgrad, nolog

Logistic regression Number of obs = 2211 LR chi2(5) = 133.72 Prob > chi2 = 0.0000 Log likelihood = -1196.114 Pseudo R2 = 0.0529

-------------------------------------------------------------------------------
high_occ | Odds Ratio Std. Err. z P>|z| [95% Conf. Interval]
--------------+----------------------------------------------------------------
         race |
2 | .3943154 .0652863 -5.62 0.000 .2850437 .5454766 3 | .9770963 .1709298 -0.13 0.895 .6934755 1.376713
              |
1.collgrad | 2.325301 .3051036 6.43 0.000 1.798013 3.007223
              |
race#collgrad |
2 1 | 1.485761 .4431258 1.33 0.184 .828094 2.665744 3 1 | 1.466932 .4419418 1.27 0.203 .812772 2.647591
              |
_cons | .3225389 .0235064 -15.53 0.000 .2796063 .3720635
-------------------------------------------------------------------------------
********************************************************************************

and we want to evaluate the overall importance of the interaction between race and collgrad (i.e. jointly)? Is it approriate to use the likelihood ratio test to compare the model without interaction to the model with interaction, and determine the importance of the interaction effect according based on the results of LR test?

--
Thomas Speidel


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