**charts **Use the second sheet in "classprobabilities.xslx" **Figure 3.2 label define varcode 1 "Member Charitable Society" 2 "Laws in Accordance with Islam" 3 "Apply Sharia Strictly" 4 "Less Rights for non-Muslims" 5 "Religious Practice Private" 6 "Sharia in Marriage" 7 "Sharia Word of God" 8 "Religious People in Public Office" 9 "Differences among Scholars" 10 "Religious Minority Practice" 11 "Prefer Religious Parties" 12 "No Influence Religious Institutions" label values varcode varcode label define prob 1 "Religious Individualist" 2 "Post-Islamist" 3 "Religious Communitarian" 4 "Status Quoist" label values prob prob graph bar (mean) var5, over(varval) over(varcode) stack asyvars by(prob) bar(1, bcolor(gs14)) bar(2, bcolor(gs10)) bar(3, bcolor(gs6)) bar(4, bcolor(gs2)) scheme(s1mono) legend(order(1 "Weak Preference" 5 "Strong Preference")) ***class vs religiosity **Figure 3.3 **use ABIII after running “abIIIcode.do” and " LCAandmodels.do" bysort religindex: egen Group1 = mean(V1classpost1) bysort religindex: egen Group2 = mean(V1classpost2) bysort religindex: egen Group3 = mean(V1classpost3) bysort religindex: egen Group4 = mean(V1classpost4) keep religindex Group1 Group2 Group3 Group4 filter gen case=_n reshape long Group, i(case) j(group) label define group 1 "Religious Individualist" 2 "Post-Islamist" 3 "Religious Communitarian" 4 "Status Quoist" label values group group graph twoway (histogram religindex if Group!=0 & filter==1, frac) (line Group religindex, lcolor(black) ), scheme(s1mono) ytitle(Posterior Probabilities) legend (label(1 "Fraction of Religious") label(2 "Class Probabilities") ) by(group) ***Figure 4 **use ABIII after running Figure 3.3“abIIIcode.do” and " LCAandmodels.do" reg diff i.class2 i.class1 i.class3 i.trust c.poltrust c.polinterest c.egalindex c.pecon c.secon c.age i.gender2 c.education c.income i.country if filter==1 margins class1 class2 class3 reg diff i.class2 i.class1 i.class3 religindex i.trust c.poltrust c.polinterest c.egalindex c.pecon c.secon c.age i.gender2 c.education c.income i.country if filter==1 margins class1 class2 class3 reg diff i.v2class1 i.v2class2 i.v2class3 i.trust c.poltrust c.polinterest c.egalindex c.pecon c.secon c.age i.gender2 c.education c.income i.country margins v2class1 v2class2 v2class3 reg diff i.v2class1 i.v2class2 i.v2class3 religindex i.trust c.poltrust c.polinterest c.egalindex c.pecon c.secon c.age i.gender2 c.education c.income i.country margins v2class1 v2class2 v2class3 twoway (connected margin classval if class==1, mcolor(black) msize(small) msymbol(circle) lcolor(black) lpattern(solid) connect(direct)) (rcap lower upper classval if class==1, lcolor(black)) (connected margin classval if class==2, mcolor(black) msize(small) msymbol(circle) lcolor(black) lwidth(medthick) lpattern(dash_dot) connect(direct)) (rcap lower upper classval if class==2, lcolor(black) lwidth(thin) lpattern(solid)) (connected margin classval if class==3, mcolor(black) msize(small) msymbol(circle) lcolor(black) lwidth(thick) lpattern(tight_dot) connect(direct)) (rcap lower upper classval if class==3 , lcolor(black)) if withrelig==0, ytitle(Marginal Effects) ytitle(, size(large)) ylabel(#3, labsize(medium)) xtitle(Class) xtitle(, size(medium)) xlabel(#2, labsize(medium)) legend(label(1 "Religious Individualist") label(3 "Post-Islamist") label(5 "Religious Communitarian") ) scheme(s1mono) graphregion(fcolor(white)) plotregion(fcolor(white)) by(sample)