libname system04 "c:\sasdat\system04\"; libname revise04 "c:\sasdat\system04\revise"; options linesize=256; options pagesize=200; data system04.cd1(compress=yes); infile "c:\sasdat\system04\rawcd1.csv" delimiter="," missover; length respid $ 10; input year respid $ agency $ ltype prop $ lpurp occup lamt preapp $ actype msa $ state $ county $ tract $ apeth $ caeth $ aprace1 $ aprace2 $ aprace3 $ aprace4 $ aprace5 $ carace1 $ carace2 $ carace3 $ carace4 $ carace5 apsex casex apinc $ ptype $ den1 $ den2 $ den3 $ ratespr $ hoepa $ lien $ edit $ sequence $ pop $ minpct $ hudmfi $ trmsa $ ownocc $ sfunits $ apflag; if (apflag=0 and msa ne "NA"); run; data system04.cd2(compress=yes); infile "c:\sasdat\system04\rawcd2.csv" delimiter="," missover; length respid $ 10; input year respid $ agency $ ltype prop $ lpurp occup lamt preapp $ actype msa $ state $ county $ tract $ apeth $ caeth $ aprace1 $ aprace2 $ aprace3 $ aprace4 $ aprace5 $ carace1 $ carace2 $ carace3 $ carace4 $ carace5 apsex casex apinc $ ptype $ den1 $ den2 $ den3 $ ratespr $ hoepa $ lien $ edit $ sequence $ pop $ minpct $ hudmfi $ trmsa $ ownocc $ sfunits $ apflag; if (apflag=0 and msa ne "NA"); run; data system04.cd3(compress=yes); infile "c:\sasdat\system04\rawcd3.csv" delimiter="," missover; length respid $ 10; input year respid $ agency $ ltype prop $ lpurp occup lamt preapp $ actype msa $ state $ county $ tract $ apeth $ caeth $ aprace1 $ aprace2 $ aprace3 $ aprace4 $ aprace5 $ carace1 $ carace2 $ carace3 $ carace4 $ carace5 apsex casex apinc $ ptype $ den1 $ den2 $ den3 $ ratespr $ hoepa $ lien $ edit $ sequence $ pop $ minpct $ hudmfi $ trmsa $ ownocc $ sfunits $ apflag; if (apflag=0 and msa ne "NA"); run; data revise04.revised(compress=yes); set system04.cd1 system04.cd2 system04.cd3; where (apinc ne "na" and apinc ne "NA" and tract ne "na" and tract ne "NA"); nstate=0; nstate=state; ncounty=0; ncounty=county; ntract=0; ntract=tract; nfips=(nstate*10000000)+(ncounty*10000)+ntract; fips=state||county||tract; rate=.; if (ratespr ne 'na' and ratespr ne "NA") then rate=ratespr; trxmsa=.; if trmsa ne "NA" then trxmsa=trmsa; hicost=0; if rate gt 0 then hicost=1; rvhicost=1; if rate gt 0 then rvhicost=0; msainc=0; if hudmfi ne "NA" then msainc=hudmfi; napinc=(apinc*1000); napinc2=napinc*napinc; l_napinc=log(napinc); nlamt=(lamt*1000); vx1=0; if (prop ne '3') then vx1=1; label vx1="all sf records"; vx2=0; if (prop ne '3') and (ltype=2 or ltype=3 or ltype=4) then vx2=1; label vx2="all sf fha records"; vx3=0; if (prop ne '3') and (ltype=1) and (actype lt 6) then vx3=1; label vx3="sf conventional excl purch preapp"; vx4=0; if (prop ne '3') and (ltype=1) and (actype = 3) then vx4=1; label vx4="sf conventional denials"; vx5=0; if (prop ne '3') and (ltype=1) and (actype lt 6) and (lpurp=1) then vx5=1; label vx5="sf conventional hp apps"; vx6=0; if (prop ne '3') and (ltype=1) and (actype=1) and (lpurp=1) then vx6=1; label vx6="sf conventional hp origs"; vx7=0; if (prop ne '3') and (ltype=1) and (actype=1) and (lpurp=1) and hicost=1 then vx7=1; label vx7="hicost sf conventional hp origs"; vx8=0; if (prop ne '3') and (ltype=1) and (actype lt 6) and (lpurp=2) then vx8=1; label vx8="sf conventional hi apps"; vx9=0; if (prop ne '3') and (ltype=1) and (actype=1) and (lpurp=2) then vx9=1; label vx9="sf conventional hi orig"; vx10=0; if (prop ne '3') and (ltype=1) and (actype=1) and (lpurp=2) and hicost=1 then vx10=1; label vx10="hicost sf conventional hi origs"; vx11=0; if (prop ne '3') and (ltype=1) and (actype lt 6) and (lpurp=3) then vx11=1; label vx11="sf conventional refi apps"; vx12=0; if (prop ne '3') and (ltype=1) and (actype=1) and (lpurp=3) then vx12=1; label vx12="sf conventional refi origs"; vx13=0; if (prop ne '3') and (ltype=1) and (actype=1) and (lpurp=3) and hicost=1 then vx13=1; label vx13="hicost sf conventional refi origs"; vx14=0; if (prop = '2') then vx14=1; label vx14="all mobile home records"; vx15=0; if (prop ne '3') and (apeth='2') and (aprace1='3') then vx15=1; label vx15="nonhisp black records"; vx16=0; if (prop ne '3') and (apeth='1') then vx16=1; label vx16="hispanic records"; vx17=0; if (prop ne '3') and (apeth='2') and (aprace1='5') then vx17=1; label vx17="nonhispanic white records"; vx18=0; if (prop ne '3') and (apeth='2') and (aprace1='2') then vx18=1; label vx18="nonhispanic asian records"; vx19=0; if (prop ne '3') and (apeth='2') and (aprace1='1' or aprace1='4') then vx19=1; label vx19="nonhispanic native records"; vx20=0; if (prop ne '3') and (ltype=1) and (actype=1) and (apeth='2') and (aprace1='5') then vx20=1; label vx20="NH white sf originations"; vx21=0; if (prop ne '3') and (ltype=1) and (actype=1) and (apeth='2') and (aprace1='5') and (hicost=1) then vx21=1; label vx21="NH white hicost originations"; vx22=0; if (prop ne '3') and (ltype=1) and (actype=1) and (apeth='2') and (aprace1='3') then vx22=1; label vx22="NH black sf originations"; vx23=0; if (prop ne '3') and (ltype=1) and (actype=1) and (apeth='2') and (aprace1='3') and (hicost=1) then vx23=1; label vx23="NH black hicost originations"; vx24=0; if (prop ne '3') and (ltype=1) and (actype=1) and (apeth='1') then vx24=1; label vx24="Hispanic sf originations"; vx25=0; if (prop ne '3') and (ltype=1) and (actype=1) and (apeth='1') and (hicost=1) then vx25=1; label vx25="Hispanic hicost originations"; vx26=0; if (prop ne '3') and (agency='1' or agency='4' or agency='5') then vx26=1; label vx26="All SF Applications to Reason-required Lenders"; vx27=0; if (prop ne '3') and actype=3 and (agency='1' or agency='4' or agency='5') then vx27=1; label vx27="All Denials at Reason-Required Lenders"; vx28=0; if (prop ne '3') and actype=3 and (den1='3' or den2='3' or den3='3') and (agency='1' or agency='4' or agency='5') then vx28=1; label vx28="Bad-Credit Denials at Reason-Required Lenders"; ransel=rannor(0); *Denial Reverse Dummy Codes. Note that denial reason reporting is only required for lenders supervised; *by OCC, OTS, and NCUA. So we only define reverse dummies for them, and we only model the credit history; *regression model for applications at these types of lenders; rvcred1=1; if actype=3 and (den1='3' or den2='3' or den3='3') and (agency='1' or agency='4' or agency='5') then rvcred1=0; label rvcred1="rv dummy1 bad credit"; rvcred2=1; if actype=3 and (den1='3') and (agency='1' or agency='4' or agency='5') then rvcred2=0; label rvcred2="rv dummy2 bad credit"; rvemploy=1; if actype=3 and (den1='2' or den2='2' or den3='2') and (agency='1' or agency='4' or agency='5') then rvemploy=0; label rvemploy="rv dummy emply hist"; rvcollat=1; if actype=3 and (den1='4' or den2='4' or den3='4') and (agency='1' or agency='4' or agency='5') then rvcollat=0; label rvcollat="rv dummy collateral"; rvcash=1; if actype=3 and (den1='5' or den2='5' or den3='5') and (agency='1' or agency='4' or agency='5') then rvcash= 0; label rvcash="rv dummy cash"; rvinfo=1; if actype=3 and (den1='6' or den1='7') and (agency='1' or agency='4' or agency='5') then rvinfo=0; label rvinfo="rv dummy info"; *For this one we define it for all lenders, because it is always required reporting; rvaband=1; if actype=5 then rvaband=0; label rvaband="rv dummy file incomplete"; incrat=napinc/msainc; lratio=napinc/nlamt; lratio2=lratio*lratio; l_ratio=log(lratio*100); label l_ratio="log of (loan ratio as percentile)"; *note that preapproval reporting is only required for certain 'covered' preapproval programs; *therefore, preapp can be 1, yes, 2, no, or 3, not applicable, if the lender does not have a 'covered'; *pre-approval procedure; xpre=0; if preapp='1' then xpre=1; *Note, demunk is coded 1 if any application information is unreported -- sex, race, or ethnicity; demunk=0; if (apeth='3' or apeth='4') or (aprace1='6' or aprace1='7') or (apsex=3 or apsex=4) then demunk=1; female=0; if apsex=2 then female=1; hispan=0; if apeth='1' then hispan=1; native=0; if (apeth='2' and aprace1='1') then native=1; asian=0; if (apeth='2' and (aprace1='2' or aprace1='4')) then asian=1; black=0; if (apeth='2' and aprace1='3') then black=1; xoccup=0; if occup=1 then xoccup=1; conv=0; if ltype=1 then conv=1; lien2=0; if lien='2' then lien2=1; lien0=0; if lien='3' then lien0=1; r_occ=0; if agency='1' then r_occ=1; r_fdc=0; if agency='3' then r_fdc=1; r_ots=0; if agency='4' then r_ots=1; r_ncu=0; if agency='5' then r_ncu=1; r_hud=0; if agency='7' then r_hud=1; *reference category for regulator is hud; imprv=0; if lpurp=2 then imprv=1; refin=0; if lpurp=3 then refin=1; s_port=0; if actype=1 and ptype='0' then s_port=1; s_gses=0; if actype=1 and ptype='1' or ptype='2' or ptype='3' or ptype='4' then s_gses=1; s_priv=0; if actype=1 and ptype='5' then s_priv=1; s_bank=0; if actype=1 and ptype='6' then s_bank=1; s_fcom=0; if actype=1 and ptype='7' then s_fcom=1; s_affl=0; if actype=1 and ptype='8' then s_affl=1; s_othr=0; if actype=1 and ptype='9' then s_othr=1; cr_logit= -3.5010 +(napinc*(-6.43E-6)) +(lratio*0.1828) +(lratio2*-0.00007) +(xoccup*0.3193) +(lien2*0.3587) +(lien0*1.2814) +(xpre*-0.8372) +(r_ots*-0.0331) +(r_ncu*-1.4243) +(imprv*1.6149) +(refin*0.8124) +(lratio*imprv*-0.1254) +(lratio*refin*-0.0724) +(lratio*xoccup*-0.0252) +(demunk*0.4055) +(female*0.1141) +(hispan*0.8116) +(native*0.7220) +(asian*0.0131) +(black*1.1351); cr_ins=(exp(cr_logit))/(1+exp(cr_logit)); cr_ins2=cr_ins*cr_ins; l_crins=log(cr_ins*100); label l_crins="natural log of credit instrument as PERCENTILE"; newlogit= 3.1634 +(l_napinc * (-0.7327)) +(l_ratio * 0.3601) +(xoccup * 0.3038) +(lien2 * (-0.0420)) +(lien0 * 0.6230) +(xpre * (-1.5046)) +(r_ots * (-0.2655)) +(r_ncu * (-1.9481)) +(imprv * 1.2751) +(refin * 0.5715) +(demunk * 0.2545) +(female * 0.0275) +(hispan * 0.7325) +(native * 0.9236) +(asian * 0.2988) +(black * 0.8319); ncr_ins=(exp(newlogit))/(1+exp(newlogit)); l_ncr=log(ncr_ins*100); label l_ncr="natural log of NEW credit instrument as PERCENTILE"; **********************************; *Code distance of each state from New Mexico in the seven-dimensional space of predatory lending; *regulation; **********************************; mdsdist=.; if nstate= 1 then mdsdist= 0.8571 ; if nstate= 2 then mdsdist= 0.7820 ; if nstate= 4 then mdsdist= 0.9248 ; if nstate= 5 then mdsdist= 0.4290 ; if nstate= 6 then mdsdist= 0.9248 ; if nstate= 8 then mdsdist= 0.8571 ; if nstate= 9 then mdsdist= 0.7451 ; if nstate= 10 then mdsdist= 0.9718 ; if nstate= 11 then mdsdist= 0.3779 ; if nstate= 12 then mdsdist= 0.8290 ; if nstate= 13 then mdsdist= 0.3515 ; if nstate= 15 then mdsdist= 0.8571 ; if nstate= 16 then mdsdist= 0.8361 ; if nstate= 17 then mdsdist= 0.5432 ; if nstate= 18 then mdsdist= 0.9718 ; if nstate= 19 then mdsdist= 0.7143 ; if nstate= 20 then mdsdist= 0.8831 ; if nstate= 21 then mdsdist= 0.8621 ; if nstate= 22 then mdsdist= 1.0000 ; if nstate= 23 then mdsdist= 0.9248 ; if nstate= 24 then mdsdist= 0.7143 ; if nstate= 25 then mdsdist= 0.0724 ; if nstate= 26 then mdsdist= 0.7143 ; if nstate= 27 then mdsdist= 0.5553 ; if nstate= 28 then mdsdist= 0.8571 ; if nstate= 29 then mdsdist= 0.7685 ; if nstate= 30 then mdsdist= 0.9248 ; if nstate= 31 then mdsdist= 0.8571 ; if nstate= 32 then mdsdist= 0.8571 ; if nstate= 33 then mdsdist= 0.8571 ; if nstate= 34 then mdsdist= 0.2809 ; if nstate= 35 then mdsdist= 0.0000 ; if nstate= 36 then mdsdist= 0.1754 ; if nstate= 37 then mdsdist= 0.0911 ; if nstate= 38 then mdsdist= 1.0000 ; if nstate= 39 then mdsdist= 0.8831 ; if nstate= 40 then mdsdist= 0.9248 ; if nstate= 41 then mdsdist= 0.9248 ; if nstate= 42 then mdsdist= 0.9718 ; if nstate= 44 then mdsdist= 0.8571 ; if nstate= 45 then mdsdist= 0.3322 ; if nstate= 46 then mdsdist= 0.9248 ; if nstate= 47 then mdsdist= 0.8571 ; if nstate= 48 then mdsdist= 0.6391 ; if nstate= 49 then mdsdist= 0.8571 ; if nstate= 50 then mdsdist= 0.8290 ; if nstate= 51 then mdsdist= 0.6438 ; if nstate= 53 then mdsdist= 0.9248 ; if nstate= 54 then mdsdist= 0.0182 ; if nstate= 55 then mdsdist= 1.0000 ; if nstate= 56 then mdsdist= 0.8571 ; mdsdist2=mdsdist*mdsdist; ***********************************; *Code the individual msas that are; *part of the black belt exploitation; *cluster 7, in order of distance; *from cluster centroid; m33860=0; if msa='33860' then m33860=1; m43340=0; if msa='43340' then m43340=1; m33740=0; if msa='33740' then m33740=1; m10780=0; if msa='10780' then m10780=1; m40580=0; if msa='40580' then m40580=1; m12260=0; if msa='12260' then m12260=1; m44940=0; if msa='44940' then m44940=1; m19260=0; if msa='19260' then m19260=1; m22500=0; if msa='22500' then m22500=1; m47580=0; if msa='47580' then m47580=1; m17980=0; if msa='17980' then m17980=1; m17900=0; if msa='17900' then m17900=1; m27140=0; if msa='27140' then m27140=1; m13820=0; if msa='13820' then m13820=1; m46220=0; if msa='46220' then m46220=1; m10500=0; if msa='10500' then m10500=1; m15500=0; if msa='15500' then m15500=1; m12940=0; if msa='12940' then m12940=1; m33660=0; if msa='33660' then m33660=1; m27180=0; if msa='27180' then m27180=1; m20020=0; if msa='20020' then m20020=1; m19460=0; if msa='19460' then m19460=1; m31420=0; if msa='31420' then m31420=1; m30020=0; if msa='30020' then m30020=1; m22180=0; if msa='22180' then m22180=1; m24140=0; if msa='24140' then m24140=1; m38220=0; if msa='38220' then m38220=1; m32820=0; if msa='32820' then m32820=1; m25980=0; if msa='25980' then m25980=1; m19804=0; if msa='19804' then m19804=1; **************************************; *Define the codes for the lenders identified in Scheesselle's 2004 Subprime Lender List; ***************************************************************************************; spcode=0; if agency='7' and respid='72-0506617' then spcode=1; if agency='7' and respid='1403800009' then spcode=1; if agency='2' and respid='0001844818' then spcode=1; if agency='2' and respid='0001844827' then spcode=1; if agency='4' and respid='0000008183' then spcode=1; if agency='1' and respid='52-2113031' then spcode=1; if agency='7' and respid='7884800009' then spcode=1; if agency='2' and respid='0001843606' then spcode=1; if agency='2' and respid='0001844845' then spcode=1; if agency='2' and respid='0001844966' then spcode=1; if agency='2' and respid='0001847305' then spcode=1; if agency='2' and respid='0001847163' then spcode=1; if agency='2' and respid='0002040989' then spcode=1; if agency='2' and respid='0003144576' then spcode=1; if agency='2' and respid='0002040381' then spcode=1; if agency='2' and respid='0002041007' then spcode=1; if agency='7' and respid='1085800002' then spcode=1; if agency='2' and respid='0002817118' then spcode=1; if agency='2' and respid='0001844265' then spcode=1; if agency='2' and respid='0001844331' then spcode=1; if agency='2' and respid='0001844975' then spcode=1; if agency='7' and respid='23-2425397' then spcode=1; if agency='2' and respid='0001847323' then spcode=1; if agency='2' and respid='0003106181' then spcode=1; if agency='2' and respid='0001844359' then spcode=1; if agency='2' and respid='0003236677' then spcode=1; if agency='2' and respid='0001844957' then spcode=1; if agency='2' and respid='0001844984' then spcode=1; if agency='2' and respid='0001426027' then spcode=1; if agency='2' and respid='0001844340' then spcode=1; if agency='2' and respid='0001845011' then spcode=1; if agency='7' and respid='75-2958981' then spcode=1; if agency='2' and respid='0002750532' then spcode=1; if agency='1' and respid='0000014740' then spcode=1; if agency='7' and respid='7062300004' then spcode=1; if agency='2' and respid='0001845039' then spcode=1; if agency='2' and respid='0001844948' then spcode=1; if agency='2' and respid='0002040402' then spcode=1; if agency='2' and respid='0002751762' then spcode=1; if agency='2' and respid='0002752022' then spcode=1; if agency='2' and respid='0001845002' then spcode=1; if agency='2' and respid='0002751922' then spcode=1; if agency='7' and respid='1206600009' then spcode=1; if agency='2' and respid='0001845020' then spcode=1; if agency='7' and respid='41-1810165' then spcode=1; if agency='2' and respid='0002751968' then spcode=1; if agency='7' and respid='2130700003' then spcode=1; if agency='2' and respid='0002134936' then spcode=1; if agency='2' and respid='0002752013' then spcode=1; if agency='2' and respid='0002751847' then spcode=1; if agency='2' and respid='0001847378' then spcode=1; if agency='2' and respid='0002751995' then spcode=1; if agency='2' and respid='0001844836' then spcode=1; if agency='7' and respid='7624900003' then spcode=1; if agency='2' and respid='0002751801' then spcode=1; if agency='2' and respid='0002041034' then spcode=1; if agency='2' and respid='0002751810' then spcode=1; if agency='2' and respid='0001843624' then spcode=1; if agency='2' and respid='0001844920' then spcode=1; if agency='2' and respid='0000860473' then spcode=1; if agency='2' and respid='0002039488' then spcode=1; if agency='7' and respid='22-3887207' then spcode=1; if agency='2' and respid='0001844993' then spcode=1; if agency='2' and respid='0001844939' then spcode=1; if agency='3' and respid='58-2257419' then spcode=1; if agency='7' and respid='1463300003' then spcode=1; if agency='2' and respid='0000902270' then spcode=1; if agency='7' and respid='1596800005' then spcode=1; if agency='2' and respid='0002750242' then spcode=1; if agency='2' and respid='0002751986' then spcode=1; if agency='2' and respid='0002135401' then spcode=1; if agency='7' and respid='0542046140' then spcode=1; if agency='7' and respid='56-2270727' then spcode=1; if agency='2' and respid='0002041098' then spcode=1; if agency='7' and respid='1549600001' then spcode=1; if agency='7' and respid='02-0560839' then spcode=1; if agency='2' and respid='0002752321' then spcode=1; if agency='2' and respid='0001950308' then spcode=1; if agency='3' and respid='0000025653' then spcode=1; if agency='7' and respid='1059700002' then spcode=1; if agency='7' and respid='62-1717999' then spcode=1; if agency='2' and respid='0001847369' then spcode=1; if agency='7' and respid='7775100007' then spcode=1; if agency='7' and respid='7673400003' then spcode=1; if agency='2' and respid='0001868177' then spcode=1; if agency='7' and respid='2064900008' then spcode=1; if agency='2' and respid='0002752684' then spcode=1; if agency='2' and respid='0002040961' then spcode=1; if agency='7' and respid='7069000008' then spcode=1; if agency='2' and respid='0001844283' then spcode=1; if agency='7' and respid='75-2921540' then spcode=1; if agency='7' and respid='1534900004' then spcode=1; if agency='7' and respid='1289800005' then spcode=1; if agency='7' and respid='1394900006' then spcode=1; if agency='7' and respid='13-2607259' then spcode=1; if agency='7' and respid='1438500002' then spcode=1; if agency='2' and respid='0002752077' then spcode=1; if agency='7' and respid='7177000002' then spcode=1; if agency='2' and respid='0002041418' then spcode=1; if agency='7' and respid='0232646780' then spcode=1; if agency='7' and respid='1851400008' then spcode=1; if agency='7' and respid='33-0806023' then spcode=1; if agency='2' and respid='0002861595' then spcode=1; if agency='7' and respid='7652800004' then spcode=1; if agency='7' and respid='22-3191019' then spcode=1; if agency='7' and respid='1958100003' then spcode=1; if agency='7' and respid='2033300001' then spcode=1; if agency='3' and respid='38-3589255' then spcode=1; if agency='2' and respid='0002752527' then spcode=1; if agency='7' and respid='2147800003' then spcode=1; if agency='7' and respid='1483300003' then spcode=1; if agency='7' and respid='1999200003' then spcode=1; if agency='7' and respid='1958200004' then spcode=1; if agency='7' and respid='1770800008' then spcode=1; if agency='7' and respid='7885300000' then spcode=1; if agency='2' and respid='0001942602' then spcode=1; if agency='7' and respid='05-0560686' then spcode=1; if agency='7' and respid='7373400006' then spcode=1; if agency='4' and respid='0000014470' then spcode=1; if agency='4' and respid='0640713034' then spcode=1; if agency='7' and respid='1843200000' then spcode=1; if agency='7' and respid='7756600001' then spcode=1; if agency='7' and respid='13-2868346' then spcode=1; if agency='7' and respid='1512400000' then spcode=1; if agency='7' and respid='7642300001' then spcode=1; if agency='7' and respid='27-0011311' then spcode=1; if agency='7' and respid='2007800007' then spcode=1; if agency='7' and respid='0458600405' then spcode=1; if agency='7' and respid='04-3697858' then spcode=1; if agency='7' and respid='95-2622032' then spcode=1; if agency='7' and respid='1830200003' then spcode=1; if agency='7' and respid='1611300007' then spcode=1; if agency='7' and respid='1281200009' then spcode=1; if agency='7' and respid='1917700009' then spcode=1; if agency='7' and respid='1707500002' then spcode=1; if agency='7' and respid='51-0113718' then spcode=1; if agency='7' and respid='1857500007' then spcode=1; if agency='7' and respid='7854200009' then spcode=1; if agency='2' and respid='0001847332' then spcode=1; if agency='7' and respid='1795300001' then spcode=1; if agency='7' and respid='1718800007' then spcode=1; if agency='7' and respid='7185300006' then spcode=1; if agency='7' and respid='1557900001' then spcode=1; if agency='7' and respid='51-0003820' then spcode=1; if agency='7' and respid='7083400004' then spcode=1; if agency='7' and respid='05-0402708' then spcode=1; if agency='7' and respid='7650700000' then spcode=1; if agency='7' and respid='1809900002' then spcode=1; if agency='2' and respid='0001847275' then spcode=1; if agency='7' and respid='1259600009' then spcode=1; if agency='4' and respid='0000012504' then spcode=1; if agency='7' and respid='7943800003' then spcode=1; if agency='4' and respid='0061104540' then spcode=1; if agency='7' and respid='95-4855230' then spcode=1; if agency='7' and respid='1194900009' then spcode=1; if agency='7' and respid='7900200006' then spcode=1; if agency='7' and respid='7431100008' then spcode=1; if agency='7' and respid='1061700007' then spcode=1; if agency='7' and respid='7479800008' then spcode=1; if agency='7' and respid='7604800006' then spcode=1; if agency='7' and respid='22-3740390' then spcode=1; if agency='1' and respid='0000024095' then spcode=1; if agency='2' and respid='0000612618' then spcode=1; if agency='7' and respid='1003800004' then spcode=1; if agency='7' and respid='36-1239445' then spcode=1; if agency='7' and respid='1750100005' then spcode=1; if agency='7' and respid='1345700002' then spcode=1; if agency='7' and respid='1349000009' then spcode=1; if agency='7' and respid='1786900005' then spcode=1; if agency='7' and respid='6480209999' then spcode=1; if agency='2' and respid='0003103489' then spcode=1; if agency='7' and respid='1991500005' then spcode=1; if agency='7' and respid='7892800004' then spcode=1; if agency='7' and respid='1345300006' then spcode=1; if agency='7' and respid='7718200008' then spcode=1; if agency='7' and respid='7935900009' then spcode=1; if agency='4' and respid='0000014939' then spcode=1; if agency='7' and respid='31-1480758' then spcode=1; if agency='7' and respid='0351887380' then spcode=1; if agency='7' and respid='0330756645' then spcode=1; if agency='7' and respid='7621700002' then spcode=1; if agency='7' and respid='33-0668958' then spcode=1; if agency='7' and respid='7542800002' then spcode=1; if agency='7' and respid='1864800009' then spcode=1; if agency='7' and respid='22-3397885' then spcode=1; if agency='7' and respid='7343000001' then spcode=1; if agency='7' and respid='7177900003' then spcode=1; if agency='7' and respid='1891200009' then spcode=1; if agency='7' and respid='1422400005' then spcode=1; if agency='2' and respid='0003073694' then spcode=1; if agency='7' and respid='1775300005' then spcode=1; if agency='7' and respid='75-1794257' then spcode=1; if agency='7' and respid='7323800008' then spcode=1; if agency='7' and respid='1118100001' then spcode=1; if agency='7' and respid='33-0967630' then spcode=1; if agency='7' and respid='7152500000' then spcode=1; if agency='7' and respid='1665100001' then spcode=1; if agency='7' and respid='1736300009' then spcode=1; if agency='7' and respid='1207500004' then spcode=1; if agency='7' and respid='1390800005' then spcode=1; if agency='7' and respid='1474600000' then spcode=1; if agency='7' and respid='7439600003' then spcode=1; if agency='7' and respid='77-0463885' then spcode=1; if agency='7' and respid='2033000009' then spcode=1; if agency='7' and respid='7786700007' then spcode=1; if agency='7' and respid='7756300009' then spcode=1; if agency='7' and respid='1596200005' then spcode=1; if agency='7' and respid='7770700000' then spcode=1; if agency='7' and respid='05-0527593' then spcode=1; if agency='7' and respid='11-2815564' then spcode=1; if agency='7' and respid='7042100008' then spcode=1; if agency='7' and respid='7619100004' then spcode=1; if agency='7' and respid='0041198897' then spcode=1; if agency='7' and respid='58-2509828' then spcode=1; if agency='7' and respid='1250200008' then spcode=1; if agency='7' and respid='1283400003' then spcode=1; if agency='7' and respid='1809800001' then spcode=1; if agency='7' and respid='0510312284' then spcode=1; if agency='7' and respid='30-0237185' then spcode=1; if agency='7' and respid='1291700001' then spcode=1; if agency='7' and respid='7086400004' then spcode=1; if agency='7' and respid='23-3051009' then spcode=1; if agency='7' and respid='7875200001' then spcode=1; if agency='3' and respid='0000026799' then spcode=1; if agency='4' and respid='0000002630' then spcode=1; if agency='4' and respid='0000008266' then spcode=1; if agency='7' and respid='2120900003' then spcode=1; if agency='3' and respid='0000000178' then spcode=1; seg_type=.; if spcode=1 and hicost=1 then seg_type=0; *High cost originations at subprime lenders; if spcode=0 and actype=3 then seg_type=1; *Rejections at mainstream lenders; if spcode=1 and actype=3 then seg_type=2; *Rejections at subprime lenders; if spcode=0 and hicost=1 then seg_type=3; *High cost originations at mainstream lenders; *Note, all other types - withdrawals, etc. - are not included in these four categories; *For these cases, seg_type will be set to missing; run; *************************************************; *Revised Credit History Instrument Model; *************************************************; proc logistic data=revise04.revised; where (actype < 6) and (ransel gt 3.15) and (conv=1) and (nstate<72) and (nfips ne .) and (edit ne '5') and (edit ne '6') and (edit ne '7') and (agency='1' or agency='4' or agency='5'); model rvcred1=l_napinc l_ratio xoccup lien2 lien0 xpre r_ots r_ncu imprv refin demunk female hispan native asian black / expb rsquare; output out=revise04.pred(compress=yes) xbeta=xb_01 p=pbad; title 'NEW Credit model, only OCC, OTS, NCUA lenders'; units l_napinc=SD l_ratio=SD ; run; data revise04.ctemp(compress=yes); set revise04.pred; dec1=.;if 0.00<=pbad<0.05 then dec1=1; dec2=.;if 0.05<=pbad<0.10 then dec2=1; dec3=.;if 0.10<=pbad<0.15 then dec3=1; dec4=.;if 0.15<=pbad<0.20 then dec4=1; dec5=.;if 0.20<=pbad<0.25 then dec5=1; dec6=.;if 0.25<=pbad<0.30 then dec6=1; dec7=.;if 0.30<=pbad<0.35 then dec7=1; dec8=.;if 0.35<=pbad<0.40 then dec8=1; dec9=.;if 0.40<=pbad<0.45 then dec9=1; dec10=.;if 0.45<=pbad<0.50 then dec10=1; dec11=.;if 0.50<=pbad<=0.55 then dec11=1; dec12=.;if 0.55<=pbad<0.60 then dec12=1; dec13=.;if 0.60<=pbad<0.65 then dec13=1; dec14=.;if 0.65<=pbad<0.70 then dec14=1; dec15=.;if 0.70<=pbad<0.75 then dec15=1; dec16=.;if 0.75<=pbad<0.80 then dec16=1; dec17=.;if 0.80<=pbad<0.85 then dec17=1; dec18=.;if 0.85<=pbad<0.90 then dec18=1; dec19=.;if 0.90<=pbad<0.95 then dec19=1; dec20=.;if 0.95<=pbad<1.00 then dec20=1; p_dec1=pbad*dec1; p_dec2=pbad*dec2; p_dec3=pbad*dec3; p_dec4=pbad*dec4; p_dec5=pbad*dec5; p_dec6=pbad*dec6; p_dec7=pbad*dec7; p_dec8=pbad*dec8; p_dec9=pbad*dec9; p_dec10=pbad*dec10; p_dec11=pbad*dec11; p_dec12=pbad*dec12; p_dec13=pbad*dec13; p_dec14=pbad*dec14; p_dec15=pbad*dec15; p_dec16=pbad*dec16; p_dec17=pbad*dec17; p_dec18=pbad*dec18; p_dec19=pbad*dec19; p_dec20=pbad*dec20; y_dec1=(1-rvcred1)*dec1; y_dec2=(1-rvcred1)*dec2; y_dec3=(1-rvcred1)*dec3; y_dec4=(1-rvcred1)*dec4; y_dec5=(1-rvcred1)*dec5; y_dec6=(1-rvcred1)*dec6; y_dec7=(1-rvcred1)*dec7; y_dec8=(1-rvcred1)*dec8; y_dec9=(1-rvcred1)*dec9; y_dec10=(1-rvcred1)*dec10; y_dec11=(1-rvcred1)*dec11; y_dec12=(1-rvcred1)*dec12; y_dec13=(1-rvcred1)*dec13; y_dec14=(1-rvcred1)*dec14; y_dec15=(1-rvcred1)*dec15; y_dec16=(1-rvcred1)*dec16; y_dec17=(1-rvcred1)*dec17; y_dec18=(1-rvcred1)*dec18; y_dec19=(1-rvcred1)*dec19; y_dec20=(1-rvcred1)*dec20; run; proc means data=revise04.ctemp n nmiss min mean max; var p_dec1 p_dec2 p_dec3 p_dec4 p_dec5 p_dec6 p_dec7 p_dec8 p_dec9 p_dec10 p_dec11 p_dec12 p_dec13 p_dec14 p_dec15 p_dec16 p_dec17 p_dec18 p_dec19 p_dec20 y_dec1 y_dec2 y_dec3 y_dec4 y_dec5 y_dec6 y_dec7 y_dec8 y_dec9 y_dec10 y_dec11 y_dec12 y_dec13 y_dec14 y_dec15 y_dec16 y_dec17 y_dec18 y_dec19 y_dec20; run; proc reg data=revise04.revised; where (actype < 6) and (ransel gt 3.15) and (conv=1) and (nstate<72) and (nfips ne .) and (edit ne '5') and (edit ne '6') and (edit ne '7') and (agency='1' or agency='4' or agency='5'); model rvcred1=l_napinc l_ratio xoccup lien2 lien0 xpre r_ots r_ncu imprv refin demunk female hispan native asian black / tol; title 'Tolerance Tests for New Credit Model'; run; *****************************************************************************; *Standard Denial and Segmentation Models to Test for Robust Credit Instrument; *****************************************************************************; proc logistic data=revise04.revised; where(ransel gt 3.15) and (conv=1) and (nstate<72) and (nfips ne .) and (edit ne '5') and (edit ne '6') and (edit ne '7') and (prop ne '3') and actype=1; model rvhicost=l_napinc l_ratio xoccup lien2 lien0 xpre r_occ r_ots r_fdc r_hud imprv refin demunk female hispan native asian black ncr_ins / expb rsquare; units l_napinc=SD l_ratio=SD ncr_ins=SD ; title 'Hicost Segmentation Model'; run; proc reg data=revise04.revised; where(ransel gt 3.15) and (conv=1) and (nstate<72) and (nfips ne .) and (edit ne '5') and (edit ne '6') and (edit ne '7') and (prop ne '3') and actype=1; model rvhicost=l_napinc l_ratio xoccup lien2 lien0 xpre r_occ r_ots r_fdc r_hud imprv refin demunk female hispan native asian black ncr_ins / tol; title 'Tolerance Tests for Hicost Segmentation Model'; run; proc logistic data=revise04.revised; where(ransel gt 3.15) and (conv=1) and (nstate<72) and (nfips ne .) and (edit ne '5') and (edit ne '6') and (edit ne '7') and (prop ne '3') and (actype=1 or actype=3); model actype=l_napinc l_ratio xoccup lien2 lien0 xpre r_occ r_ots r_fdc r_hud imprv refin demunk female hispan native asian black ncr_ins / expb rsquare; units l_napinc=SD l_ratio=SD ncr_ins=SD ; title 'Approval Model'; run; proc reg data=revise04.revised; where(ransel gt 3.15) and (conv=1) and (nstate<72) and (nfips ne .) and (edit ne '5') and (edit ne '6') and (edit ne '7') and (prop ne '3') and (actype=1 or actype=3); model actype=l_napinc l_ratio xoccup lien2 lien0 xpre r_occ r_ots r_fdc r_hud imprv refin demunk female hispan native asian black ncr_ins / tol; title 'Tolerance Tests for Approval Model'; run; ******************************************************************; *Create summary-level tabulations of HMDA mortgage market measures; *and merge with metropolitan area labels and with SF 3 sample data; *from the 2000 Census. Note that the metromod database will not ; *include data for metropolitan areas where 2000 Census data could ; *not be matched to the new MA codes used in the 2004 HMDA ; ******************************************************************; proc means data=revise04.revised noprint; class msa; var vx1-vx28; output out=revise04.sums(compress=yes) sum=vx1-vx28; run; data revise04.msa(compress=yes); infile "c:\sasdat\system04\msas.csv" delimiter="," missover; length msaname $ 75; input msa $ msaname $; run; proc sort data=revise04.msa; by msa; run; proc sort data=revise04.sums; by msa; run; data revise04.msasum(compress=yes); merge revise04.msa revise04.sums; by msa; if (msa ne '99999' and _type_=1); m_fha=0; m_fha=vx2/vx1; label m_fha="fha share"; m_den=0; m_den=vx4/vx3; label m_den="denial rate, conventional sf"; m_hi=0; m_hi=vx8/vx3; label m_hi="hi share, conventional sf"; m_rf=0; m_rf=vx11/vx3; label m_rf="rf share, conventional sf"; m_mob=0; m_mob=vx14/vx1; label m_mob="mobile home share, all sf"; m_blk=0; m_blk=vx15/vx1; label m_blk="nh black share, all sf"; m_hsp=0; m_hsp=vx16/vx1; label m_hsp="hispanic share, all sf"; m_wht=0; m_wht=vx17/vx1; label m_wht="nh white share, all sf"; m_asn=0; m_asn=vx18/vx1; label m_asn="nh asian share, all sf"; m_ntv=0; m_ntv=vx19/vx1; label m_ntv="nh natve share, all sf"; m_hihp=0; m_hihp=vx7/vx6; label m_hihp="hicost share home purchase"; m_hihi=0; m_hihi=vx10/vx9; label m_hihi="hicost share home imprvemt"; m_hirf=0; m_hirf=vx13/vx12;label m_hirf="hicost share refinance"; m_sseg=0; m_sseg=(vx7+vx10+vx13)/(vx6+vx9+vx12); label m_sseg="overall conventional hicost share"; m_cred=0; m_cred=(vx28/vx27); label m_cred="bad-credit denial rate"; run; data revise04.cens(compress=yes); infile "c:\sasdat\system04\revise\newmsa.csv" delimiter="," missover; input msa $ mx_pop mx_nhw mx_nhb mx_nhn mx_hsp mx_for mx_frc mx_mil mx_une mx_pci mx_pov mx_btm mx_btu mx_htm mx_htu mx_wtm mx_wtu mx_bpq mx_hpq mx_own mx_p50 mx_a95 mx_rnt mx_val mx_mfree mx_2mtg mx_oeld mx_bfree mx_hfree; label mx_pop = "Total metropolitan population"; label mx_nhw = "Non-hispanic white only share"; label mx_nhb = "Non-hispanic black share"; label mx_nhn = "Non-hispanic native share"; label mx_hsp = "Hispanic share"; label mx_for = "Foreign-born share"; label mx_frc = "Share of foreign born recent"; label mx_mil = "Share of labor force in armed forces"; label mx_une = "Civilian unemployment rate"; label mx_pci = "Per capita income 1999"; label mx_pov = "Person poverty rate 1999"; label mx_btm = "Proportion of Black teens in the military"; label mx_btu = "Black teen non school civilian unemployment rate"; label mx_htm = "Proportion of Hispanic teens in the military"; label mx_htu = "Hispanic teen non school civilian unemployment rate"; label mx_wtm = "Proportion of white teens in the military"; label mx_wtu = "White teen non school civilian unemployment rate"; label mx_bpq = "Ratio, White to Black PCI"; label mx_hpq = "Ratio, White to Hispanic PCI"; label mx_own = "Owner occupied as share of all OHU"; label mx_p50 = "Percent of own occ hu built before 1950"; label mx_a95 = "Percent of own occ hu built after 1995"; label mx_rnt = "MGRent as share of household income"; label mx_val = "Ratio, median value to median household income"; label mx_mfree = "Share of spec ohu with no mortgage"; label mx_2mtg = "Share of mtged units with 2ndmtg"; label mx_oeld = "Ratio: Share of elderly living in pre 1950 owned homes / Share of owned homes that are pre 1950"; label mx_bfree = "Ratio: Black-White share mortgage free"; label mx_hfree = "Ratio: Hispanic-White share mortgage free"; lmx_pop=log(mx_pop); label lmx_pop = "Log of Metro Population"; run; proc sort; by msa; run; proc sort data=revise04.msasum; by msa; run; data revise04.metromod(compress=yes); merge revise04.msasum revise04.cens; by msa; if mx_pop ne .; run; *************************************; *Merge the Individual Loan-Level Database; *with the Metropolitan-level measures derived from a) HMDA themselves,; *and b) SF 3 of the 2000 Census of Population and Housing; *EXCLUDES applications in areas where census data could not be matched; *to the new MA codes in HMDA; proc sort data=revise04.metromod; by msa; run; proc datasets library=revise04; modify revised; index create msa; run; data revise04.multimod(compress=yes); merge revise04.revised revise04.metromod; by msa; if mx_pop ne .; run; **************************************; *Metropolitan Level Models; **************************************; proc reg data=revise04.metromod simple; title 'Metro Subprime Models, MaxR Selections'; id msaname; model m_sseg = lmx_pop mx_nhb mx_nhn mx_hsp mx_frc mx_mil mx_une mx_pci mx_bpq mx_hpq mx_own mx_p50 mx_a95 mx_rnt mx_val mx_mfree mx_2mtg mx_oeld mx_bfree mx_hfree m_fha m_hi m_rf m_mob m_cred m_den / selection=maxr; run; proc reg data=revise04.metromod simple; title 'Metro Subprime Model, Parsimonious'; id msaname; model m_sseg = lmx_pop mx_nhb mx_nhn mx_hsp mx_frc mx_mil mx_pci mx_bpq mx_hpq mx_own mx_p50 mx_rnt mx_mfree mx_2mtg m_fha m_hi m_rf m_mob m_cred m_den / tol stb scorr2 pcorr2 p; run; proc reg data=revise04.metromod simple; title 'Metro Denial Models, MaxR Selections'; id msaname; model m_den = lmx_pop mx_nhb mx_nhn mx_hsp mx_frc mx_mil mx_une mx_pci mx_bpq mx_hpq mx_own mx_p50 mx_a95 mx_rnt mx_val mx_mfree mx_2mtg mx_oeld mx_bfree mx_hfree m_fha m_hi m_rf m_mob m_cred m_sseg / selection=maxr; run; proc reg data=revise04.metromod simple; title 'Metro Denial Model, Parsimonious'; id msaname; model m_den = lmx_pop mx_nhb mx_nhn mx_hsp mx_frc mx_mil mx_pci mx_bpq mx_hpq mx_own mx_p50 mx_rnt mx_mfree mx_2mtg m_fha m_hi m_rf m_mob m_cred m_sseg / tol stb scorr2 pcorr2 p; run; **************************************; *New Segmentation Models; **************************************; proc logistic data=revise04.multimod; where(ransel gt 2.8) and (conv=1) and (nstate<72) and (nfips ne .) and (edit ne '5') and (edit ne '6') and (edit ne '7') and (prop ne '3') and (seg_type=0 or seg_type=1); model seg_type=l_napinc l_ratio xoccup lien2 lien0 xpre r_occ r_ots r_fdc r_hud imprv refin demunk female hispan native asian black ncr_ins / expb rsquare; units l_napinc=SD l_ratio=SD ncr_ins=SD ; title '0=hicost at spl 1=deny prime'; run; proc logistic data=revise04.multimod; where(ransel gt 2.8) and (conv=1) and (nstate<72) and (nfips ne .) and (edit ne '5') and (edit ne '6') and (edit ne '7') and (prop ne '3') and (seg_type=0 or seg_type=1); model seg_type=l_napinc l_ratio xoccup lien2 lien0 xpre r_occ r_ots r_fdc r_hud imprv refin demunk female hispan native asian black ncr_ins m_den / expb rsquare; units l_napinc=SD l_ratio=SD ncr_ins=SD m_den=SD ; title '0=hicost at spl 1=deny prime'; run; proc logistic data=revise04.multimod; where(ransel gt 2.8) and (conv=1) and (nstate<72) and (nfips ne .) and (edit ne '5') and (edit ne '6') and (edit ne '7') and (prop ne '3') and (seg_type=1 or seg_type=2); model seg_type=l_napinc l_ratio xoccup lien2 lien0 xpre r_occ r_ots r_fdc r_hud imprv refin demunk female hispan native asian black ncr_ins / expb rsquare; units l_napinc=SD l_ratio=SD ncr_ins=SD ; title '1=deny prime 2=deny subprime'; run; proc logistic data=revise04.multimod; where(ransel gt 2.8) and (conv=1) and (nstate<72) and (nfips ne .) and (edit ne '5') and (edit ne '6') and (edit ne '7') and (prop ne '3') and (seg_type=0 or seg_type=3); model seg_type=l_napinc l_ratio xoccup lien2 lien0 xpre r_occ r_ots r_fdc r_hud imprv refin demunk female hispan native asian black ncr_ins / expb rsquare; units l_napinc=SD l_ratio=SD ncr_ins=SD ; title '0=hicost at spl 3=hicost at prime'; run; ***************************************************************************************; *Old Version of Segmentation Models applied to new database with Metropolitan variables; ***************************************************************************************; proc logistic data=revise04.multimod; where(ransel gt 2.8) and (conv=1) and (nstate<72) and (nfips ne .) and (edit ne '5') and (edit ne '6') and (edit ne '7') and (prop ne '3') and actype=1; model rvhicost=l_napinc l_ratio xoccup lien2 lien0 xpre r_occ r_ots r_fdc r_hud imprv refin demunk female hispan native asian black ncr_ins / expb rsquare; units l_napinc=SD l_ratio=SD ncr_ins=SD ; title 'Hicost Segment & Denial Model 1'; run; proc logistic data=revise04.multimod; where(ransel gt 2.8) and (conv=1) and (nstate<72) and (nfips ne .) and (edit ne '5') and (edit ne '6') and (edit ne '7') and (prop ne '3') and actype=1; model rvhicost=l_napinc l_ratio xoccup lien2 lien0 xpre r_occ r_ots r_fdc r_hud imprv refin demunk female hispan native asian black ncr_ins m_den / expb rsquare; units l_napinc=SD l_ratio=SD ncr_ins=SD m_den=SD ; title 'Hicost Segment & Denial Model 2'; run; proc logistic data=revise04.multimod; where(ransel gt 2.8) and (conv=1) and (nstate<72) and (nfips ne .) and (edit ne '5') and (edit ne '6') and (edit ne '7') and (prop ne '3') and actype=1; model rvhicost=l_napinc l_ratio xoccup lien2 lien0 xpre r_occ r_ots r_fdc r_hud imprv refin demunk female hispan native asian black ncr_ins m_den hispan*m_den native*m_den asian*m_den black*m_den / expb rsquare; units l_napinc=SD l_ratio=SD ncr_ins=SD m_den=SD ; title 'Hicost Segment & Denial Model 3'; run; proc logistic data=revise04.multimod; where(ransel gt 2.8) and (conv=1) and (nstate<72) and (nfips ne .) and (edit ne '5') and (edit ne '6') and (edit ne '7') and (prop ne '3') and actype=1; model rvhicost=l_napinc l_ratio xoccup lien2 lien0 xpre r_occ r_ots r_fdc r_hud imprv refin demunk female hispan native asian black ncr_ins mdsdist mdsdist2 / expb rsquare; units l_napinc=SD l_ratio=SD ncr_ins=SD mdsdist=SD ; title 'Hicost Segment Model with MDSdist'; run; **********; *Segmentation model to test the effects of the Black Belt; *Choose all applications in these metros, and a random sample of all other applications; proc logistic data=revise04.revised; where (ransel gt 3.0 or m33860=1 or m43340=1 or m33740=1 or m10780=1 or m40580=1 or m12260=1 or m44940=1 or m19260=1 or m22500=1 or m47580=1 or m17980=1 or m17900=1 or m27140=1 or m13820=1 or m46220=1 or m10500=1 or m15500=1 or m12940=1 or m33660=1 or m27180=1 or m20020=1 or m19460=1 or m31420=1 or m30020=1 or m22180=1 or m24140=1 or m38220=1 or m32820=1 or m25980=1 or m19804=1) and (conv=1) and (nstate<72) and (nfips ne .) and (edit ne '5') and (edit ne '6') and (edit ne '7') and (prop ne '3') and actype=1; model rvhicost=l_napinc l_ratio xoccup lien2 lien0 xpre r_occ r_ots r_fdc r_hud imprv refin demunk female hispan native asian black ncr_ins m33860 m43340 m33740 m10780 m40580 m12260 m44940 m19260 m22500 m47580 m17980 m17900 m27140 m13820 m46220 m10500 m15500 m12940 m33660 m27180 m20020 m19460 m31420 m30020 m22180 m24140 m38220 m32820 m25980 m19804 / expb rsquare; units l_napinc=SD l_ratio=SD ncr_ins=SD ; title 'Hicost Segment Model with Black Belt Metros'; run;