*Syntaxis of the article "Can the legal framework for prostitution influence the acceptability of buying sex?"
*
*
*Program STATA 14.2
*
*
*
*
*Import data (indicate direction)
import excel "C:\...\DATA_PROSTITUTION.xls", sheet("DATA") firstrow clear

*Create interactions
gen INT_DECRIM_NOTRAFFICKING = DECRIMINALIZED*NOTRAFFICKING
gen INT_DECRIM_NOPIMPS = DECRIMINALIZED*NOPIMPS
gen INT_DECRIM_PROSVIOLENCE = DECRIMINALIZED*PROSTITUTION_VIOLENCE


*Table 2 ---------------------------------------------------------------------	
*Model 1 (M1)
mixed ACCEPTABILITY DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS EMPOWER || ID_RESPONDENT: DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS, reml stddeviations cov(unstructured)
*standardized effects
estadd beta
mat beta2 = e(beta)'
matrix list beta2
estat icc
estat ic

*Model 2 (M2)
mixed ACCEPTABILITY DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS EMPOWER  WOMAN PROSTITUTION_VIOLENCE INCOMEMORE2000 ECONOMIC_LIBERALISM HEGEMONIC_MASCULINITY ANTIEGALITARIANISM || ID_RESPONDENT: DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS, reml stddeviations cov(unstructured)
*standardized effects
estadd beta
mat beta2 = e(beta)'
matrix list beta2
estat icc
estat ic

*Model 3 (M3)
mixed ACCEPTABILITY DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS EMPOWER  WOMAN PROSTITUTION_VIOLENCE INCOMEMORE2000 ECONOMIC_LIBERALISM HEGEMONIC_MASCULINITY ANTIEGALITARIANISM INT_DECRIM_NOTRAFFICKING INT_DECRIM_NOPIMPS || ID_RESPONDENT: DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS, reml stddeviations cov(unstructured)
*standardized effects
estadd beta
mat beta2 = e(beta)'
matrix list beta2
estat icc
estat ic

*Model 4 (M4)
mixed ACCEPTABILITY DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS EMPOWER  WOMAN PROSTITUTION_VIOLENCE INCOMEMORE2000 ECONOMIC_LIBERALISM HEGEMONIC_MASCULINITY ANTIEGALITARIANISM  INT_DECRIM_NOTRAFFICKING INT_DECRIM_NOPIMPS INT_DECRIM_PROSVIOLENCE || ID_RESPONDENT: DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS, reml stddeviations  cov(unstructured)
*standardized effects
estadd beta
mat beta2 = e(beta)'
matrix list beta2
estat icc
estat ic


*Figure 2 ---------------------------------------------------------------------	
*Model 1, Figure 2.1
mixed ACCEPTABILITY i.DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS EMPOWER || ID_RESPONDENT: DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS, reml stddeviations cov(unstructured)
margins i.DECRIMINALIZED
marginsplot, ytitle("Acceptability of buying sex") xtitle("Decriminalized") name(figure2_1, replace)

*Model 3, Figure 2.2
mixed ACCEPTABILITY DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS EMPOWER  WOMAN PROSTITUTION_VIOLENCE INCOMEMORE2000 ECONOMIC_LIBERALISM HEGEMONIC_MASCULINITY ANTIEGALITARIANISM DECRIMINALIZED##NOTRAFFICKING DECRIMINALIZED##NOPIMPS || ID_RESPONDENT: DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS, reml stddeviations cov(unstructured)
margins DECRIMINALIZED#NOTRAFFICKING
marginsplot, ytitle("Acceptability of buying sex") name(figure2_2, replace)

*Model 3, Figure 2.3
mixed ACCEPTABILITY DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS EMPOWER  WOMAN PROSTITUTION_VIOLENCE INCOMEMORE2000 ECONOMIC_LIBERALISM HEGEMONIC_MASCULINITY ANTIEGALITARIANISM DECRIMINALIZED##NOTRAFFICKING DECRIMINALIZED##NOPIMPS || ID_RESPONDENT: DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS, reml stddeviations cov(unstructured)
margins DECRIMINALIZED#NOPIMPS
marginsplot, ytitle("Acceptability of buying sex") name(figure2_3, replace)

*Model 4, Figure 2.4
mixed ACCEPTABILITY DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS EMPOWER  WOMAN PROSTITUTION_VIOLENCE INCOMEMORE2000 ECONOMIC_LIBERALISM HEGEMONIC_MASCULINITY ANTIEGALITARIANISM DECRIMINALIZED##NOTRAFFICKING DECRIMINALIZED##NOPIMPS  DECRIMINALIZED##PROSTITUTION_VIOLENCE|| ID_RESPONDENT: DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS, reml stddeviations cov(unstructured)
margins DECRIMINALIZED#PROSTITUTION_VIOLENCE
marginsplot, ytitle("Acceptability of buying sex") name(figure2_4, replace)

*Model 7, Figure 2.5
mixed ACCEPTABILITY DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS EMPOWER  WOMAN PROSTITUTION_VIOLENCE INCOMEMORE2000 c.ECONOMIC_LIBERALISM c.HEGEMONIC_MASCULINITY c.ANTIEGALITARIANISM DECRIMINALIZED##NOTRAFFICKING DECRIMINALIZED##NOPIMPS DECRIMINALIZED##PROSTITUTION_VIOLENCE WOMAN##DECRIMINALIZED WOMAN##NOTRAFFICKING WOMAN##INDOOR WOMAN##NOPIMPS WOMAN##EMPOWER   WOMAN##PROSTITUTION_VIOLENCE WOMAN##INCOMEMORE2000 WOMAN##c.ECONOMIC_LIBERALISM WOMAN##c.HEGEMONIC_MASCULINITY WOMAN##c.ANTIEGALITARIANISM WOMAN##DECRIMINALIZED##NOTRAFFICKING WOMAN##DECRIMINALIZED##NOPIMPS WOMAN##DECRIMINALIZED##PROSTITUTION_VIOLENCE || ID_RESPONDENT: DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS, reml stddeviations cov(unstructured)
margins DECRIMINALIZED#WOMAN
marginsplot, ytitle("Acceptability of buying sex") name(figure2_5, replace)



*Table 3 ---------------------------------------------------------------------	
*Model 5 (M5): MEN
mixed ACCEPTABILITY DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS EMPOWER PROSTITUTION_VIOLENCE INCOMEMORE2000 ECONOMIC_LIBERALISM HEGEMONIC_MASCULINITY ANTIEGALITARIANISM  INT_DECRIM_NOTRAFFICKING INT_DECRIM_NOPIMPS INT_DECRIM_PROSVIOLENCE || ID_RESPONDENT: DECRIMINALIZED NOTRAFFICKING INDOOR EMPOWER if WOMAN==0, reml stddeviations  cov(unstructured)
*standardized effects
estadd beta
mat beta2 = e(beta)'
matrix list beta2
estat icc
estat ic

*Model 6 (M6): WOMEN
mixed ACCEPTABILITY DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS EMPOWER PROSTITUTION_VIOLENCE INCOMEMORE2000 ECONOMIC_LIBERALISM HEGEMONIC_MASCULINITY ANTIEGALITARIANISM  INT_DECRIM_NOTRAFFICKING INT_DECRIM_NOPIMPS INT_DECRIM_PROSVIOLENCE || ID_RESPONDENT: DECRIMINALIZED NOTRAFFICKING INDOOR EMPOWER if WOMAN==1, reml stddeviations  cov(unstructured)
*standardized effects
estadd beta
mat beta2 = e(beta)'
matrix list beta2
estat icc
estat ic


*Figure 3 ---------------------------------------------------------------------	
*Standardized variables
foreach var of varlist ACCEPTABILITY DECRIMINALIZED NOTRAFFICKING INDOOR NOPIMPS EMPOWER PROSTITUTION_VIOLENCE INCOMEMORE2000 ECONOMIC_LIBERALISM HEGEMONIC_MASCULINITY ANTIEGALITARIANISM  INT_DECRIM_NOTRAFFICKING INT_DECRIM_NOPIMPS INT_DECRIM_PROSVIOLENCE {
	gen std_`var' = .
	tempvar aux
	levelsof WOMAN, local(VALUES)
	foreach i of local VALUES {
		egen `aux' =  std(`var') if WOMAN == `i'
		replace std_`var' = `aux' if WOMAN == `i'
		drop `aux'
	}
}

*Model 5 (standardized effects): MEN
mixed std_ACCEPTABILITY std_DECRIMINALIZED std_NOTRAFFICKING std_INDOOR std_NOPIMPS std_EMPOWER std_PROSTITUTION_VIOLENCE std_INCOMEMORE2000 std_ECONOMIC_LIBERALISM std_HEGEMONIC_MASCULINITY std_ANTIEGALITARIANISM  std_INT_DECRIM_NOTRAFFICKING std_INT_DECRIM_NOPIMPS std_INT_DECRIM_PROSVIOLENCE || ID_RESPONDENT: std_DECRIMINALIZED std_NOTRAFFICKING std_INDOOR std_EMPOWER if WOMAN==0, reml stddeviations  cov(unstructured)
estimates store men

*Model 6 (standardized effects): WOMEN
mixed std_ACCEPTABILITY std_DECRIMINALIZED std_NOTRAFFICKING std_INDOOR std_NOPIMPS std_EMPOWER std_PROSTITUTION_VIOLENCE std_INCOMEMORE2000 std_ECONOMIC_LIBERALISM std_HEGEMONIC_MASCULINITY std_ANTIEGALITARIANISM  std_INT_DECRIM_NOTRAFFICKING std_INT_DECRIM_NOPIMPS std_INT_DECRIM_PROSVIOLENCE || ID_RESPONDENT: std_DECRIMINALIZED std_NOTRAFFICKING std_INDOOR std_EMPOWER if WOMAN==1, reml stddeviations  cov(unstructured)
estimates store women

*Models M5 and M6, Figure 3
coefplot (men, label(Man) pstyle(p4)) (women, label(Woman)), drop(_cons) xline(0)