*/ Replication of results presented Michael C. Horowitz, “Nonstate Actors and the Diffusion of Innovations: The Case of Suicide Terrorism,” International Organization 64:1 (Winter 2010), pp. 33-64.*/

*/ If you find any mistakes or issues, please contact the author at horom (at) sas.upenn.edu */

*/ Table 1 Model 1 */
logit suicide orgage, robust

*/ Table 1 Model 2 */
logit suicide orgage  religious communistsocialist leftist nationalistseparatist other, robust

*/ Table 1 Model 3 */
logit suicide orgage  religious communistsocialist leftist nationalistseparatist other lebanon iraqwar israel aqlink, robust

*/ Table 1 Model 4 */
logit suicide orgage religious religorg communistsocialist leftist nationalistseparatist other lebanon iraqwar israel aqlink, robust

*/ Figure 1 */
postgr3 orgage, by(religious) asis(orgage religious religorg) gen(suicidefigure) table

*/ Figure 2 */
twoway (line suicidefigure0 orgage, sort lcolor(black) lpattern(solid) lwidth(medthick)) (line suicidefigure1 orgage, sort lcolor(black) lpattern(dash) lwidth(medthick)) (lowess suicide orgage if religious==1 & syria==0, sort mean bwidth(.6)), ytitle(Probability Group Adopts Suicide Terrorism,  justification(right) margin(medsmall)) xtitle(Organizational Age, margin(medsmall)) legend(label(1 "Not Religiously-Motivated") label(2 "Religiously-Motivated") label(3 "Lowess Curve if Religious = 1"))graphregion(fcolor(white) lcolor(white) ilcolor(white)) plotregion(fcolor(white) lcolor(white) ilcolor(white))

*/ Table 3 */
table firstsuicideyear firstsuicideage

*/ Table 2 (values for religion and religorg reversed to ease in interpretation of relative risks) */

estsimp logit suicide orgage religious2 religorg2 communistsocialist leftist nationalistseparatist other lebanon iraqwar israel aqlink, robust

setx mean
summarize orgage, meanonly        
local mx1=`r(mean)'          
summarize religious2, meanonly       
local mx2=`r(mean)'           
setx religorg2 `mx1'*`mx2'       
*/ producing mean value */
simqi, pr

*/ org age low, religious low */
setx orgage min religious2 max religorg2 1

simqi, pr

setx mean
summarize orgage, meanonly        
local mx1=`r(mean)'          
summarize religious2, meanonly        
local mx2=`r(mean)'           
setx religorg2 `mx1'*`mx2'         
*/ producing mean value */
simqi, pr

*/ org age high, religious low  */
setx orgage max religious2 max religorg2 51

simqi, pr

setx mean
summarize orgage, meanonly        
local mx1=`r(mean)'           
summarize religious2, meanonly        
local mx2=`r(mean)'           
setx religorg2 `mx1'*`mx2'         
*/ producing mean value */
simqi, pr

*/ org age low, religious high */
setx orgage min religious2 min religorg2 0

simqi, pr

setx mean
summarize orgage, meanonly        
local mx1=`r(mean)'          
summarize religious2, meanonly        
local mx2=`r(mean)'         
setx religorg2 `mx1'*`mx2'         
*/ producing mean value */
simqi, pr

*/ org age high, religious high */
setx orgage max religious2 min religorg2 0

simqi, pr
clear