version 8.2 clear set memory 10m set more off capture log close set scheme sj prog drop _all *One-way tables *Example: indicator mode sjlog using examples/example01, replace use drugs describe inco1-inco7 sjlog close, replace sjlog using examples/example02, replace tabstat inco1-inco7, s(sum mean) c(s) sjlog close, replace sjlog using examples/example03, replace mrtab inco1-inco7, sort title(Sources of income) sjlog close, replace *Tech. note: Denominator sjlog using examples/example04, replace mrtab inco1-inco7, include sort title(Sources of income) sjlog close, replace *Tech. note: response() sjlog using examples/example05, replace codebook crime1 sjlog close, replace sjlog using examples/example06, replace mrtab crime1-crime5, include response(1 3) title(Crime (as committer)) /// nonames mrtab crime1-crime5, include response(2 3) title(Crime (as victim)) nonames sjlog close, replace *Example: polytomous mode sjlog using examples/example07, replace codebook pinco1 sjlog close, replace sjlog using examples/example08, replace mrtab pinco1-pinco6, poly response(1/7) include sort abbrev sjlog close, replace sjlog using examples/example09, replace codebook sinco1 mrtab sinco1-sinco6, include sort abbrev sjlog close, replace sjlog using examples/example10, replace mrtab pinco1-pinco6, poly response(1/7) nofreq include generate(_inco) cochran _inco* drop _inco* sjlog close, replace sjlog using examples/example11, replace preserve mrtab pinco1-pinco6, poly response(1/7) nofreq include generate(_inco) qui reshape long _inco, i(id) j(R) xi: logit _inco i.R, cluster(id) nolog restore sjlog close, replace *Two-way tables sjlog using examples/example12, replace mrtab inco1-inco7, include sort title(Sources of income) nonames /// width(28) by(city) column sjlog close, replace sjlog using examples/example13, replace mrtab inco1-inco7, include by(city) nofreq chi2 lrchi2 generate(_inco) sjlog close, replace sjlog using examples/example14, replace list _incorp in 1/6, clean drop _inco* sjlog close, replace sjlog using examples/example15, replace mrtab inco1-inco7, include sort title(Sources of income) nonames /// width(19) by(city) column mtest(bonferroni) sjlog close, replace *Tech. note: significance tests sjlog using examples/example16, replace preserve mrtab inco1-inco7, include nofreq generate(_inco) qui reshape long _inco, i(id) j(R) qui xi: logit _inco i.R*i.city, cluster(id) unab I : _Icity* _IRX* test `I' restore sjlog close, replace preserve mrtab inco1-inco7, include nofreq generate(_inco) reshape long _inco, i(id) j(R) sjlog using examples/example17, replace quietly xi: xtlogit _inco i.R*i.city, re i(id) estimates store A quietly xi: xtlogit _inco i.R, re i(id) lrtest A sjlog close, replace restore *Graphs *Example: one-way graph sjlog using examples/example18, replace mrgraph bar crime1-crime5, include response(2 3) sort width(15) /// title(Criminal experiences (as a victim)) ylabel(,angle(0)) sjlog close, replace graph export examples/example18.eps, replace *Example: two-way graphs sjlog using examples/example19, replace mrgraph bar crime1-crime5, include response(2 3) sort width(15) by(sex) /// stat(column) title(Criminal experiences (as a victim)) /// ylabel(,angle(0)) legend(bmargin(t+1)) sjlog close, replace graph export examples/example19.eps, replace sjlog using examples/example20, replace mrgraph hbar crime1-crime5, include response(2 3) sort by(sex, outboard) /// stat(column) title(Criminal experiences (as a victim)) ylabel(,angle(0)) sjlog close, replace graph export examples/example20.eps, replace *Tech. note: separate panels sjlog using examples/example21, replace mrgraph hbar crime1-crime5, include response(2 3) sort(1) width(16) /// stat(column) by(sex, separate title(Criminal experiences (as a victim))) sjlog close, replace graph export examples/example21.eps, replace *Example: tabplot discard sjlog using examples/example22, replace mrgraph tab crime1-crime5, include response(2 3) sort width(16) /// stat(column) by(sex) rtotal title(Criminal experiences (as a victim)) sjlog close, replace graph export examples/example22.eps, replace *Tech. note: mrsvmat sjlog using examples/example23, replace mrtab crime1-crime5, include response(2 3) by(sex) nofreq preserve _mrsvmat, stat(column) ctotal rtotal clear list, clean string(20) generate x=_n scatter C1 C2 x in 1/5, xtitle(Type of crime) ytitle(Percent) /// title(Criminal experiences (as a victim)) sjlog close, replace restore graph export examples/example23.eps, replace exit