Data Downloaded from http://www.eia.doe.gov/cneaf/coal/ctrdb/database.html ------------------------------------------------------------------------------- log: z:/hedonic_coal/industrial.log log type: text opened on: 16 Jun 2007, 15:06:32 ****************txt to stata formatted data; . insheet using "z:/CTRDB79_92.txt"; (92 vars, 12132 obs) . quietly compress; . save z:/hedonic_coal/CTRDB79_92.dta, replace; file z:/hedonic_coal/CTRDB79_92.dta saved . clear; . insheet using "z:/CTRDB93_01.txt"; (92 vars, 8830 obs) . quietly compress; . save z:/hedonic_coal/CTRDB93_01.dta, replace; file z:/hedonic_coal/CTRDB93_01.dta saved . clear; . *****************merging data sets; . use z:/hedonic_coal/CTRDB79_92.dta; . append using z:/hedonic_coal/CTRDB93_01.dta; btulow was int now long btusshipped was int now long (note: distanceofmode1 is float in using data but will be str6 now) . save z:/hedonic_coal/indeconsulfprem.dta, replace; file z:/hedonic_coal/indeconsulfprem.dta saved . sort plantcode; . clear; . insheet using z:/hedonic_coal/region.txt; (4 vars, 769 obs) . quietly compress; . sort plantcode; . save z:/hedonic_coal/region.dta, replace; file z:/hedonic_coal/region.dta saved . clear; . use z:/hedonic_coal/indeconsulfprem.dta; . sort plantcode; . merge plantcode using z:/hedonic_coal/region.dta; variable plantcode does not uniquely identify observations in the master data variable plantcode does not uniquely identify observations in z:/hedonic_coal/region.dta . save z:/hedonic_coal/indeconsulfprem.dta, replace; file z:/hedonic_coal/indeconsulfprem.dta saved . *****Price Deflator, PPI Crude Materials (table b-66 economic report of the p > resident); . gen pricefactor=0; . replace pricefactor=54.9 if year==1979; (930 real changes made) . replace pricefactor=73.1 if year==1980; (886 real changes made) . replace pricefactor=97.7 if year==1981; (870 real changes made) . replace pricefactor=100.0 if year==1982; (769 real changes made) . replace pricefactor=98.7 if year==1983; (736 real changes made) . replace pricefactor=98.0 if year==1984; (793 real changes made) . replace pricefactor=93.3 if year==1985; (791 real changes made) . replace pricefactor=71.8 if year==1986; (826 real changes made) . replace pricefactor=75.0 if year==1987; (816 real changes made) . replace pricefactor=67.7 if year==1988; (871 real changes made) . replace pricefactor=75.9 if year==1989; (883 real changes made) . replace pricefactor=85.9 if year==1990; (984 real changes made) . replace pricefactor=80.4 if year==1991; (968 real changes made) . replace pricefactor=78.8 if year==1992; (1009 real changes made) . replace pricefactor=76.7 if year==1993; (992 real changes made) . replace pricefactor=72.1 if year==1994; (1286 real changes made) . replace pricefactor=69.4 if year==1995; (1196 real changes made) . replace pricefactor=85.0 if year==1996; (947 real changes made) . replace pricefactor=87.3 if year==1997; (956 real changes made) . replace pricefactor=68.6 if year==1998; (1040 real changes made) . replace pricefactor=78.5 if year==1999; (1026 real changes made) . *****************deriving scrubber dummies; . gen scrubber_dummy = 0; . replace scrubber_dummy =1 if > plantcode == 51| plantcode == 56| plantcode == 87| plantcode == 108| plantcod > e == 113| plantcode == 127| plantcode == 130 > | plantcode == 136| plantcode == 160| plantcode == 165| plantcode == 207| pla > ntcode == 298| plantcode == 465| plantcode == 469 > | plantcode == 477| plantcode == 525| plantcode == 564| plantcode == 645| pla > ntcode == 676| plantcode == 728| plantcode == 961 > | plantcode == 963| plantcode == 976| plantcode == 994| plantcode == 995| pla > ntcode == 1012| plantcode == 1040| plantcode == 1082; (922 real changes made) . replace scrubber_dummy =1 if > plantcode == 1167| plantcode == 1241| plantcode == 1250| plantcode == 1356| p > lantcode == 1357| plantcode == 1363| plantcode == 1364 > | plantcode == 1374| plantcode == 1378| plantcode == 1382| plantcode == 1825| > plantcode == 1843| plantcode == 1866| plantcode == 1891 > | plantcode == 1893| plantcode == 1926| plantcode == 1927| plantcode == 1934| > plantcode == 2039| plantcode == 2079| plantcode == 2324 > | plantcode == 2378| plantcode == 2442| plantcode == 2451| plantcode == 2535| > plantcode == 2823| plantcode == 2824| plantcode == 2840; (1107 real changes made) . replace scrubber_dummy =1 if > plantcode == 2861| plantcode == 2917| plantcode == 3098| plantcode == 3099| > plantcode == 3118| plantcode == 3122| plantcode == 3159 > | plantcode == 3161| plantcode == 3181| plantcode == 3399| plantcode == 3470| > plantcode == 3845| plantcode == 3944| plantcode == 3954 > | plantcode == 4005| plantcode == 4040| plantcode == 4158| plantcode == 4162| > plantcode == 4259| plantcode == 4941| plantcode == 6004; (1055 real changes made) . replace scrubber_dummy =1 if > plantcode == 6016| plantcode == 6018| plantcode == 6019| plantcode == 6021| > plantcode == 6030| plantcode == 6041| plantcode == 6061 > | plantcode == 6068| plantcode == 6071| plantcode == 6076| plantcode == 6082| > plantcode == 6085| plantcode == 6090| plantcode == 6094 > | plantcode == 6101| plantcode == 6113| plantcode == 6136| plantcode == 6137| > plantcode == 6146| plantcode == 6147| plantcode == 6165; (790 real changes made) . replace scrubber_dummy =1 if > plantcode == 6177| plantcode == 6179| plantcode == 6183| plantcode == 6195| > plantcode == 6204| plantcode == 6213| plantcode == 6238 > | plantcode == 6249| plantcode == 6288| plantcode == 6469| plantcode == 6481| > plantcode == 6639| plantcode == 6648| plantcode == 6761 > | plantcode == 6768| plantcode == 6823| plantcode == 7097| plantcode == 7210| > plantcode == 7213| plantcode == 7504| plantcode == 7737; (213 real changes made) . replace scrubber_dummy =1 if > plantcode == 7790| plantcode == 7902| plantcode == 8066| plantcode == 8069| > plantcode == 8102| plantcode == 8222| plantcode == 8223 > | plantcode == 8224; (227 real changes made) . ***********gen voluntary dummies; . gen vol =0; . replace vol= 1 if plantcode == 7| plantcode == 641| plantcode == 642| plantco > de == 645| plantcode == 699| plantcode == 709| plantcode == 727 > | plantcode == 733| plantcode == 863| plantcode == 891| plantcode == 897| pla > ntcode == 898| plantcode == 991| plantcode == 1241 > | plantcode == 1570| plantcode == 1571| plantcode == 1606| plantcode == 1723| > plantcode == 1912| plantcode == 2079| plantcode == 2094 > | plantcode == 2104| plantcode == 2161| plantcode == 2168| plantcode == 2516| > plantcode == 2836| plantcode == 2838| plantcode == 2844; (1334 real changes made) . replace vol =1 if > plantcode == 2850| plantcode == 2857| plantcode == 2858| plantcode == 2864| p > lantcode == 2867| plantcode == 2877| plantcode == 2878 > | plantcode == 3138| plantcode == 3181| plantcode == 3942| plantcode == 3945| > plantcode == 4050| plantcode == 4057| plantcode == 4072 > | plantcode == 4078| plantcode == 4140| plantcode == 4271| plantcode == 6017| > plantcode == 6018| plantcode == 6061| plantcode == 6090 > | plantcode == 6094| plantcode == 6124| plantcode == 6155| plantcode == 6195| > plantcode == 6639; (1064 real changes made) . label variable vol "voluntary plants throughout the whole Phase I period"; . ***********gen voluntary dummies; . gen vol1 =0; . replace vol= 1 if plantcode == 7| plantcode == 642 | plantcode == 699| plantc > ode == 709| plantcode == 727 > | plantcode == 733| plantcode == 863| plantcode == 891| plantcode == 898| pla > ntcode == 1241 > | plantcode == 1570| plantcode == 1606| plantcode == 1723| plantcode == 2079 > | plantcode == 2104 | plantcode == 2836| plantcode == 2838; (0 real changes made) . replace vol1 =1 if > plantcode == 2850| plantcode == 2857| plantcode == 2858 | plantcode == 2867| > plantcode == 2877| plantcode == 2878 > | plantcode == 3181| plantcode == 3945| plantcode == 4057 > | plantcode == 4078| plantcode == 4140| plantcode == 4271| plantcode == 6017| > plantcode == 6018| plantcode == 6061| plantcode == 6090 > | plantcode == 6094| plantcode == 6124| plantcode == 6155| plantcode == 6195| > plantcode == 6639; (833 real changes made) . *************converting prices; . gen converted_mine_rate=.; (21368 missing values generated) . replace converted_mine_rate = fobminerate if fobmineunit==1; (27 real changes made) . replace converted_mine_rate = ((2*fobminerate)*btusshipped)/100000 if fobmine > unit==2; (14108 real changes made) . replace converted_mine_rate = ((2*fobminerate)*btusshipped)/1000 if fobmineun > it==3; (932 real changes made) . gen converted_plant_rate=.; (21368 missing values generated) . replace converted_plant_rate = fobplantrate if fobplantunit==1; (188 real changes made) . replace converted_plant_rate = ((2*fobplantrate)*btusshipped)/100000 if fobpl > antunit==2; (18967 real changes made) . replace converted_plant_rate = ((2*fobplantrate)*btusshipped)/1000 if fobplan > tunit==3; (1021 real changes made) . gen converted_rate_mode_1=.; (21368 missing values generated) . replace converted_rate_mode_1 = rateofmode1 if unitofrate1==1; (3725 real changes made) . replace converted_rate_mode_1= ((2*rateofmode1)*btusshipped)/100000 if unitof > rate1==2; (12193 real changes made) . replace converted_rate_mode_1 = ((2*rateofmode1)*btusshipped)/1000 if unitofr > ate1==3; (982 real changes made) . gen converted_rate_mode_2=.; (21368 missing values generated) . replace converted_rate_mode_2 = rateofmode2 if unitofrate2==1; (857 real changes made) . replace converted_rate_mode_2 = ((2*rateofmode2)*btusshipped)/100000 if unito > frate2==2; (876 real changes made) . replace converted_rate_mode_2 = ((2*rateofmode2)*btusshipped)/1000 if unitofr > ate2==3; (103 real changes made) . gen converted_rate_mode_3=.; (21368 missing values generated) . replace converted_rate_mode_3 = rateofmode3 if unitofrate3==1; (330 real changes made) . replace converted_rate_mode_3 = ((2*rateofmode3)*btusshipped)/100000 if unito > frate3==2; (225 real changes made) . replace converted_rate_mode_3 = ((2*rateofmode3)*btusshipped)/1000 if unitofr > ate3==3; (3 real changes made) . gen converted_rate_mode_4=.; (21368 missing values generated) . replace converted_rate_mode_4 = rateofmode4 if unitofrate4==1; (313 real changes made) . replace converted_rate_mode_4 = ((2*rateofmode4)*btusshipped)/100000 if unito > frate4==2; (71 real changes made) . replace converted_rate_mode_4 = ((2*rateofmode4)*btusshipped)/1000 if unitofr > ate4==3; (0 real changes made) . gen converted_transfer_fee_1=.; (21368 missing values generated) . replace converted_transfer_fee_1 = transferfeemode1 if transferfeeunit1==1; (419 real changes made) . replace converted_transfer_fee_1 = ((2*transferfeemode1)*btusshipped)/100000 > if transferfeeunit1==2; (450 real changes made) . replace converted_transfer_fee_1 = ((2*transferfeemode1)*btusshipped)/1000 if > transferfeeunit1==3; (2 real changes made) . gen converted_transfer_fee_2=.; (21368 missing values generated) . replace converted_transfer_fee_2 = transferfeemode2 if transferfeeunit2==1; (58 real changes made) . replace converted_transfer_fee_2 = ((2*transferfeemode2)*btusshipped)/100000 > if transferfeeunit2==2; (17 real changes made) . replace converted_transfer_fee_2 = ((2*transferfeemode2)*btusshipped)/1000 if > transferfeeunit2==3; (0 real changes made) . gen converted_transfer_fee_3=.; (21368 missing values generated) . replace converted_transfer_fee_3 = transferfeemode3 if transferfeeunit3==1; (51 real changes made) . replace converted_transfer_fee_3= ((2*transferfeemode3)*btusshipped)/100000 i > f transferfeeunit3==2; (5 real changes made) . replace converted_transfer_fee_3 = ((2*transferfeemode3)*btusshipped)/1000 if > transferfeeunit3==3; (0 real changes made) . gen converted_transfer_fee_4=.; (21368 missing values generated) . replace converted_transfer_fee_4 = transferfeemode4 if transferfeeunit4==1; (0 real changes made) . replace converted_transfer_fee_4= ((2*transferfeemode4)*btusshipped)/100000 i > f transferfeeunit4==2; (0 real changes made) . replace converted_transfer_fee_4 = ((2*transferfeemode4)*btusshipped)/1000 if > transferfeeunit4==3; (0 real changes made) . ****************calculating mine rate; . gen predminerate = .; (21368 missing values generated) . replace predminerate = converted_mine_rate if converted_mine_rate>0 & convert > ed_mine_rate~=.; (7348 real changes made) . *replace predminerate = calcmine if predminerate==. & calcmine~=.; . gen realpredminerate = (predminerate/pricefactor)*100; (14081 missing values generated) . gen realplantrate = (converted_plant_rate/pricefactor)*100; (2120 missing values generated) . replace realplantrate=. if realplantrate >=100 | realplantrate<=5; (104 real changes made, 104 to missing) . label variable realpredminerate "Real price at mine"; . label variable realplantrate "Real price at plant"; . *******************seperating time periods; . gen pre =0; . replace pre=1 if(yrsigned >0 & yrsigned<=1988 & (year ar>=yrmod) & (yrmod>0 & yrmod<=1988)); (10465 real changes made) . gen mid=0; . replace mid=1 if (yrsigned>=1989 & yrsigned<=1990 & (year ((yrmod>=1989 & yrmod<=1990) & (year>=yrmod)); (1680 real changes made) . gen tran=0; . replace tran=1 if (yrsigned>=1991 & yrsigned<=1994 & (year | ((yrmod>=1991 & yrmod<=1994) & (year>=yrmod)); (3543 real changes made) . gen post=0; . replace post=1 if (yrsigned>=1995 & yrsigned<=1999 & (year ((yrmod>=1995 & yrmod<=1999) & (year>=yrmod)); (3064 real changes made) . gen ptwo=0; . replace ptwo=1 if (yrsigned>=2000 & yrsigned<=2001 & (year ((yrmod>=2000 & yrmod<=2001) & (year>=yrmod)); (335 real changes made) . gen noauction=0; . replace noauction=1 if (yrsigned>=1990 & yrsigned<=1993 & (year =0)) | ((yrmod>=1990 & yrmod<=1993) & (year>=yrmod)); (3442 real changes made) . gen noauction2=0; . replace noauction2=1 if (yrsigned>=1990 & yrsigned<=1992 & (year ==0)) | ((yrmod>=1990 & yrmod<=1992) & (year>=yrmod)); (2446 real changes made) . label variable pre "contracts signed or renegotiated before 1989"; . label variable mid "contracts signed or renegotiated between 1989 and 1990"; . label variable tran "contracts signed or renegotiated between 1991 and 1994"; . label variable post "contracts signed or renegotiated after 1994"; . label variable noauction "contracts signed or renegotiated between 1990 and 1 > 993"; . gen before=.; (21368 missing values generated) . replace before=1 if (pre==1|mid==1); (12145 real changes made) . replace before=0 if (post==1|tran==1); (6607 real changes made) . gen after=.; (21368 missing values generated) . replace after=1 if (post==1|tran==1); (6607 real changes made) . replace after=0 if (pre==1|mid==1); (12145 real changes made) . **************generating basin dummies (based on EIA classifications); . gen app =.; (21368 missing values generated) . replace app = 1 if (bomdistrict >=1 & bomdistrict <=8) | (bomdistrict == 13) > | (bomdistrict == 24); (13071 real changes made) . replace app = 0 if (bomdistrict >=9 & bomdistrict<=12) | (bomdistrict >=14 & > bomdistrict <=23); (7732 real changes made) . gen interior =.; (21368 missing values generated) . replace interior = 1 if (bomdistrict >=9 & bomdistrict <=12) | (bomdistrict = > = 14) | (bomdistrict == 15); (3054 real changes made) . replace interior = 0 if (bomdistrict >=1 & bomdistrict<=8) | (bomdistrict >=1 > 6 & bomdistrict <=23) | (bomdistrict == 13); (17747 real changes made) . gen west =.; (21368 missing values generated) . replace west = 1 if (bomdistrict >=16 & bomdistrict <=23); (4678 real changes made) . replace west = 0 if (bomdistrict >=1 & bomdistrict<=15) | (bomdistrict == 24) > ; (16125 real changes made) . label variable app "coal mined from appalachian mine"; . label variable interior "coal mined from interior mine"; . label variable west "coal from western mine"; . replace bomdistrict=. if bomdistrict ==0 | bomdistrict>24; (159 real changes made, 159 to missing) . ****************generating plant region dummies (based on bureau of census cl > assification); . gen eastp=.; (21368 missing values generated) . replace eastp = 1 if (region == 1) | (region == 2); (2590 real changes made) . replace eastp = 0 if (region>=3); (18778 real changes made) . gen midp=.; (21368 missing values generated) . replace midp = 1 if (region == 3) | (region == 4); (8164 real changes made) . replace midp =0 if (region<=2) |(region>=5); (13204 real changes made) . gen southp=.; (21368 missing values generated) . replace southp=1 if (region>=5) & (region<=7); (8842 real changes made) . replace southp=0 if (region<=4)|(region>=8); (12526 real changes made) . gen westp=.; (21368 missing values generated) . replace westp = 1 if (region == 8)| (region >= 9); (1772 real changes made) . replace westp=0 if (region<=7); (19596 real changes made) . label variable eastp"coal to eastern plant"; . label variable midp "coal to midwest plant"; . label variable southp "coal to southern plant"; . label variable westp "coal to western plant"; . *****************calculating quadratic terms; . gen EF=.; (21368 missing values generated) . replace EF = 38 if btusshipped>=9750; (17296 real changes made) . replace EF = 35 if btusshipped<9750 & btusshipped>=7500; (3473 real changes made) . replace EF = 30 if btusshipped<7500; (599 real changes made) . label variable EF "emissions factor"; . gen sulfoxmm =.; (21368 missing values generated) . replace sulfoxmm = (EF*sulfurshipped*500)/btusshipped; (20530 real changes made) . replace sulfoxmm=. if sulfoxmm==0; (3 real changes made, 3 to missing) . replace sulfoxmm=. if sulfoxmm>7.5; (118 real changes made, 118 to missing) . label variable sulfoxmm "so2 content in pounds per million btu"; . gen sulfoxpt=.; (21368 missing values generated) . replace sulfoxpt = (EF* sulfurshipped)/2000; (20962 real changes made) . label variable sulfoxpt "so2 content per ton"; . gen dirty =.; (21368 missing values generated) . replace dirty=1 if sulfoxmm >1.2; (14361 real changes made) . replace dirty=0 if sulfoxmm<=1.2; (7007 real changes made) . label variable dirty "so2 content above 1.2 lbs per million btu"; . gen dirtysox =.; (21368 missing values generated) . replace dirtysox = sulfoxpt * dirty; (20962 real changes made) . gen dirtyphase1 =.; (21368 missing values generated) . replace dirtyphase1=1 if sulfoxmm >2.5; (6899 real changes made) . replace dirtyphase1=0 if sulfoxmm<=2.5; (14469 real changes made) . label variable dirty "so2 content above 2.5 lbs per million btu"; . gen dirtysoxphase1 =.; (21368 missing values generated) . replace dirtysoxphase1 = sulfoxpt * dirtyphase1; (20962 real changes made) . gen sulfoxpt1=.; (21368 missing values generated) . replace sulfoxpt1 = ((sulfurshipped*2)/100); (20962 real changes made) . gen sulfoxpt2=.; (21368 missing values generated) . replace sulfoxpt2 = ((sulfurshipped*1.75)/100); (20962 real changes made) . gen ashpt = (ashshipped/100); (406 missing values generated) . gen moispt = (moistureshipped/100); (406 missing values generated) . gen btupt = (btusshipped*2000)/1000000; (406 missing values generated) . gen caaa = 0; . replace caaa=1 if caaboiler>0; (3461 real changes made) . label variable ashpt "ash content per ton"; . label variable moispt "moisture content per ton"; . label variable btupt "btu content per ton"; . label variable caaa "table a boiler dummy"; . gen durationm =0; . replace durationm = yrexp-yrsigned; (19196 real changes made, 406 to missing) . replace durationm = . if durationm<0; (120 real changes made, 120 to missing) . replace durationm = . if durationm>100; (8 real changes made, 8 to missing) . gen mmpl=.; (21368 missing values generated) . replace mmpl=1 if mmpflag1==1; (305 real changes made) . replace mmpl=0 if mmpflag1==0; (20657 real changes made) . gen ct1 =0; . replace ct1 = 1 if contracttype=="1"; (12607 real changes made) . gen ct2= 0; . replace ct2=1 if contracttype =="2"; (1285 real changes made) . gen ct3=0; . replace ct3=1 if contracttype=="3"; (258 real changes made) . gen ct4=0; . replace ct4=1 if contracttype=="4"; (462 real changes made) . gen ct5=0; . replace ct5=1 if contracttype=="5"; (205 real changes made) . gen ct6=0; . replace ct6=1 if contracttype=="6"; (1114 real changes made) . gen ct7=0; . replace ct7=1 if contracttype=="7"; (2217 real changes made) . gen ct12=0; . replace ct12=1 if contracttype=="1,2"; (222 real changes made) . gen ct13=0; . replace ct13=1 if contracttype=="1,3"; (543 real changes made) . gen ct16=0; . replace ct16=1 if contracttype=="1,6"; (4 real changes made) . gen ct23=0; . replace ct23=1 if contracttype=="2,3"; (48 real changes made) . gen ct37=0; . replace ct37=1 if contracttype=="3,7"; (24 real changes made) . gen ct71=0; . replace ct71=1 if contracttype=="7,1"; (18 real changes made) . ***************year dummies; . tab year, gen (dyear); Year | Freq. Percent Cum. ------------+----------------------------------- 1979 | 930 4.44 4.44 1980 | 886 4.23 8.66 1981 | 870 4.15 12.81 1982 | 769 3.67 16.48 1983 | 736 3.51 19.99 1984 | 793 3.78 23.78 1985 | 791 3.77 27.55 1986 | 826 3.94 31.49 1987 | 816 3.89 35.38 1988 | 871 4.16 39.54 1989 | 883 4.21 43.75 1990 | 984 4.69 48.44 1991 | 968 4.62 53.06 1992 | 1,009 4.81 57.88 1993 | 992 4.73 62.61 1994 | 1,286 6.13 68.74 1995 | 1,196 5.71 74.45 1996 | 947 4.52 78.97 1997 | 956 4.56 83.53 1998 | 1,040 4.96 88.49 1999 | 1,026 4.89 93.38 2000 | 667 3.18 96.57 2001 | 720 3.43 100.00 ------------+----------------------------------- Total | 20,962 100.00 . tab bomdistrict, gen (dbomdistrict); BOM | District | Freq. Percent Cum. ------------+----------------------------------- 1 | 1,455 6.99 6.99 2 | 935 4.49 11.49 3 | 1,336 6.42 17.91 4 | 833 4.00 21.92 6 | 87 0.42 22.33 7 | 88 0.42 22.76 8 | 7,753 37.27 60.02 9 | 596 2.86 62.89 10 | 1,145 5.50 68.39 11 | 1,035 4.98 73.37 12 | 29 0.14 73.51 13 | 582 2.80 76.31 14 | 6 0.03 76.34 15 | 243 1.17 77.50 16 | 3 0.01 77.52 17 | 450 2.16 79.68 18 | 190 0.91 80.59 19 | 2,834 13.62 94.22 20 | 303 1.46 95.67 21 | 105 0.50 96.18 22 | 770 3.70 99.88 23 | 23 0.11 99.99 24 | 2 0.01 100.00 ------------+----------------------------------- Total | 20,803 100.00 . gen phase1 =0; . replace phase1 = 1 if (vol==1 & caaa==0) | (vol==0 & caaa==1)| (vol==1 & caaa > ==1); (5473 real changes made) . replace phase1 =. if vol==. |caaa==.; (0 real changes made) . gen phase11 =0; . replace phase11 = 1 if (vol1==1 & caaa==0) | (vol1==0 & caaa==1)| (vol1==1 & > caaa==1); (4291 real changes made) . replace phase11 =. if vol1==. |caaa==.; (0 real changes made) . gen sulfcaaa = sulfoxpt*caaa; (406 missing values generated) . gen sulfph1 = sulfoxpt*phase1; (406 missing values generated) . gen sulfph11 = sulfoxpt*phase11; (406 missing values generated) . gen sulfcaaaafter = sulfoxpt*caaa* after; (2616 missing values generated) . gen sulfph1after = sulfoxpt*phase1* after; (2616 missing values generated) . gen sulfph11after = sulfoxpt*phase11* after; (2616 missing values generated) . gen sulfcaaanoau = sulfoxpt*caaa*noauction; (406 missing values generated) . gen sulfnoau = sulfoxpt*noauction; (406 missing values generated) . gen sulfafter = sulfoxpt*after; (2616 missing values generated) . gen sulfph1post = sulfoxpt*phase1*post; (406 missing values generated) . gen sulfpost = sulfoxpt*post; (406 missing values generated) . gen sulfcaaapost = sulfoxpt*caaa*post; (406 missing values generated) . gen dirtyafter = dirty*after; (2616 missing values generated) . gen dirtypost = dirty*post; . gen sulfapp = sulfoxpt*app; (565 missing values generated) . gen sulfint = sulfoxpt*int; (567 missing values generated) . gen sulfwest = sulfoxpt*west; (565 missing values generated) . gen sulfeastp = sulfoxpt*eastp; (406 missing values generated) . gen sulfmidp = sulfoxpt*midp; (406 missing values generated) . gen sulfwestp = sulfoxpt*westp; (406 missing values generated) . gen sulfsouthp = sulfoxpt*southp; (406 missing values generated) . gen dissq=.; (21368 missing values generated) . replace dissq = totaldistance*totaldistance; (20962 real changes made) . gen dissqafter=.; (21368 missing values generated) . replace dissqafter = dissq*after; (18752 real changes made) . gen trend=.; (21368 missing values generated) . replace trend = year- 1978; (20962 real changes made) . replace trend = . if trend<0 | trend>25; (0 real changes made) . gen multipletrip =0; . replace multipletrip =1 if (mode1== "B" |mode1=="T") |(mode2=="B"|mode2=="T") > |(mode3=="B" |mode3=="T")|(mode4=="B"|mode4=="T"); (6422 real changes made) . label variable multipletrip "more than one company transporting coal"; . save z:/hedonic_coal/indeconsulfprem.dta, replace; file z:/hedonic_coal/indeconsulfprem.dta saved . log close; log: z:/hedonic_coal/industrial.log log type: text closed on: 16 Jun 2007, 15:07:54 -------------------------------------------------------------------------------