replace mjtimes1 = 0 if mjnumb==0 egen totalnumb=sum(mjnumb) if mjnumb<996 summarize totalnumb egen total1x=sum(mjtimes1) if mjtimes1<996 summarize total1x centile mjnumb if mjnumb<996, centile(50 80 90) egen totalnumb90=sum(mjnumb) if mjnumb>=15 & mjnumb<996 summarize totalnumb90 generate mjprop1 = mjtimes1/mjnumb if mjtimes1<996 & mjnumb<996 summarize mjprop1 if mjnumb>0 & mjnumb<=5 & mjtimes1<996 summarize mjprop1 if mjnumb>5 & mjnumb<=10 & mjtimes1<996 summarize mjprop1 if mjnumb>10 & mjnumb<=20 & mjtimes1<996 summarize mjprop1 if mjnumb>20 & mjnumb<=30 & mjtimes1<996 summarize mjprop1 if mjnumb>30 & mjnumb<=40 & mjtimes1<996 summarize mjprop1 if mjnumb>40 & mjtimes1<996 & mjnumb<996 generate hetmale = 2 replace hetmale = 1 if gender==1 & sxident==1 replace hetmale = 1 if gender==1 & sxident==5 generate nonhetmale = 2 replace nonhetmale = 1 if gender==1 & sxident==2 replace nonhetmale = 1 if gender==1 & sxident==3 replace nonhetmale = 1 if gender==1 & sxident==4 generate hetfemale = 2 replace hetfemale = 1 if gender==2 & sxident==1 replace hetfemale = 1 if gender==2 & sxident==5 generate nonhetfemale = 2 replace nonhetfemale = 1 if gender==2 & sxident==2 replace nonhetfemale = 1 if gender==2 & sxident==3 replace nonhetfemale = 1 if gender==2 & sxident==4 generate mjminus1 = mjnumb-mjtimes1 if mjnumb<996 & mjtimes1<996 summarize mjnumb if mjnumb<996 summarize mjnumb if mjnumb<996 & hetmale==1 summarize mjnumb if mjnumb<996 & nonhetmale==1 summarize mjnumb if mjnumb<996 & hetfemale==1 summarize mjnumb if mjnumb<996 & nonhetfemale==1 summarize mjminus1 if mjnumb<996 & mjtimes1<996 summarize mjminus1 if mjnumb<996 & mjtimes1<996 & hetmale==1 summarize mjminus1 if mjnumb<996 & mjtimes1<996 & nonhetmale==1 summarize mjminus1 if mjnumb<996 & mjtimes1<996 & hetfemale==1 summarize mjminus1 if mjnumb<996 & mjtimes1<996 & nonhetfemale==1 generate mjminus100 = mjnumb-mjtimes1 if mjnumb<996 & mjtimes1<996 summarize mjminus100 generate mjminus50 = mjnumb-(mjtimes1*.5) if mjnumb<996 & mjtimes1<996 summarize mjminus50 generate mjminus20 = mjnumb-(mjtimes1*.2) if mjnumb<996 & mjtimes1<996 summarize mjminus20 generate mjminus10 = mjnumb-(mjtimes1*.1) if mjnumb<996 & mjtimes1<996 summarize mjminus10 replace v18frqs1=0 if v18num==0 replace v18frqs2=0 if v18num==0 replace v18frqs3=0 if v18num==0 drop if v18num~=(v18frqs1+v18frqs2+v18frqs3) egen total18num=sum(v18num) if v18num<98 summarize total18num egen total181x=sum(v18frqs1) if v18num<98 & v18frqs1<97 summarize total181x summarize v18num if v18num<98 generate v18minus1 = v18num-v18frqs1 if v18num<98 & v18frqs1<96 summarize v18minus1