
. do "C:\Users\csuna\AppData\Local\Temp\STD8d48_000010.tmp"

. clear all

.     cd "D:\Data\AntiSLAPP"
D:\Data\AntiSLAPP

.         
.         set matsize 800
set matsize ignored.
    Matrix sizes are no longer limited by c(matsize) in modern Statas.  Matrix sizes are now limited by edition of
    Stata.  See limits for more details.

.         set more off

.         use finalsample_emission.dta

.         destring, replace 
TRIFD: contains nonnumeric characters; no replace
year already numeric; no replace
P_CO_DB_NUM: contains nonnumeric characters; no replace
P_CO_NAME: contains nonnumeric characters; no replace
state: contains nonnumeric characters; no replace
COUNTY: contains nonnumeric characters; no replace
federal: contains nonnumeric characters; no replace
avgprodratio already numeric; no replace
avgcumprod already numeric; no replace
avgcumprod_n already numeric; no replace
avgtrprod already numeric; no replace
avgtrprod_n already numeric; no replace
nchemicals already numeric; no replace
tr already numeric; no replace
onsite_tr already numeric; no replace
offsite_tr already numeric; no replace
transfer already numeric; no replace
prodwste already numeric; no replace
treat already numeric; no replace
recycle already numeric; no replace
recover already numeric; no replace
air already numeric; no replace
land already numeric; no replace
water already numeric; no replace
other already numeric; no replace
naics4: all characters numeric; replaced as int
PRIMARY_naics: all characters numeric; replaced as long
PRIMARY_SIC: contains nonnumeric characters; no replace
DunsNumber already numeric; no replace
hqduns already numeric; no replace
HQCompany1: contains nonnumeric characters; no replace
City_his: contains nonnumeric characters; no replace
State_his: contains nonnumeric characters; no replace
HQDunsChange: contains nonnumeric characters; no replace
fips already numeric; no replace
emp already numeric; no replace
sales already numeric; no replace
sic2_nets already numeric; no replace
emp_lag1 already numeric; no replace
sales_lag1 already numeric; no replace
hqdunslink already numeric; no replace
logempnum already numeric; no replace
logtremp already numeric; no replace
logtremp1 already numeric; no replace
fips1: all characters numeric; replaced as long
trbp_on already numeric; no replace
a263 already numeric; no replace
a246 already numeric; no replace
trbp already numeric; no replace
population_lag1 already numeric; no replace
population already numeric; no replace
ipc_lag1 already numeric; no replace
IPC already numeric; no replace
Unemployment_Rate_lag1 already numeric; no replace
Unemployment_Rate already numeric; no replace
nccs already numeric; no replace
nccs_lag1 already numeric; no replace
negviewratio already numeric; no replace
posviewratio already numeric; no replace
gvkey: all characters numeric; replaced as long
DATADATE already numeric; no replace
FYEAR already numeric; no replace
INCORP: contains nonnumeric characters; no replace
LOC: contains nonnumeric characters; no replace
SIC: all characters numeric; replaced as int
(20954 missing values generated)
pstate: contains nonnumeric characters; no replace
TIC: contains nonnumeric characters; no replace
CUSIP: contains nonnumeric characters; no replace
CIK: all characters numeric; replaced as long
(22762 missing values generated)
AT already numeric; no replace
CEQ already numeric; no replace
DLC already numeric; no replace
DLTT already numeric; no replace
INVT already numeric; no replace
SEQ already numeric; no replace
COGS already numeric; no replace
DP already numeric; no replace
DVC already numeric; no replace
DVP already numeric; no replace
IB already numeric; no replace
OIBDP already numeric; no replace
XRD already numeric; no replace
OANCF already numeric; no replace
PRCC_F already numeric; no replace
SICH already numeric; no replace
at_lag1 already numeric; no replace
mv already numeric; no replace
size already numeric; no replace
logat already numeric; no replace
lev already numeric; no replace
tang already numeric; no replace
roa already numeric; no replace
assetturn already numeric; no replace
xrdat already numeric; no replace
oi already numeric; no replace
margin already numeric; no replace
sg already numeric; no replace
xrddum already numeric; no replace
headqstate_yb: contains nonnumeric characters; no replace
numlocalnews_lag1 already numeric; no replace
numlocalnews already numeric; no replace
pop already numeric; no replace
logipc already numeric; no replace
employmentrate already numeric; no replace
dpop already numeric; no replace
dipc already numeric; no replace
demp already numeric; no replace
n_county already numeric; no replace
eventdate already numeric; no replace
eventyear already numeric; no replace
eventstate already numeric; no replace
AntiSLAPP already numeric; no replace
avgotherAntiSLAPP already numeric; no replace
avgotherAntiSLAPP_nonHQ already numeric; no replace
wavgotherAntiSLAPP already numeric; no replace
wavgotherAntiSLAPP_nonHQ already numeric; no replace
fips2 already numeric; no replace
RegEPA already numeric; no replace
FedEPA already numeric; no replace
EPAregcode already numeric; no replace
EPAregfips: all characters numeric; replaced as long
(4324 missing values generated)
numsrilocal1 already numeric; no replace
nccs_pop_premean already numeric; no replace
nccs_pop_pre already numeric; no replace
numlocalnews_pre already numeric; no replace
numsrilocal1_pre already numeric; no replace
statepol_pre1 already numeric; no replace
statepro_pre1 already numeric; no replace
statepol_pre2 already numeric; no replace
statepro_pre2 already numeric; no replace
statepol already numeric; no replace
statepro already numeric; no replace
statepol_post1 already numeric; no replace
statepro_post1 already numeric; no replace
statepol_post2 already numeric; no replace
statepro_post2 already numeric; no replace
abpm already numeric; no replace
abgop already numeric; no replace
abprodm already numeric; no replace
abmaterial already numeric; no replace
abother already numeric; no replace
emission_intensity already numeric; no replace

.         drop if year<1990
(0 observations deleted)

.         drop if year>2020
(0 observations deleted)

.         
.         drop if mi(at_lag1)
(21,444 observations deleted)

.         drop if mi(fips)
(0 observations deleted)

.         drop if mi(gvkey)
(0 observations deleted)

.         drop if mi(TRIFD)
(0 observations deleted)

.         drop if mi(year)
(0 observations deleted)

. 
.         isid year fips TRIFD gvkey

.         egen statecode=group(state)

.         gen naics2=int(naics4/100)

.         gen naics3=int(naics4/10)       

.         egen faclld=group(TRIFD)

.         egen gvkeyyear=group(gvkey year)

.         egen fipsyear=group(fips year)

.         drop if SIC>=6000 & SIC<=6999
(1,225 observations deleted)

.         drop if SIC>=4900 & SIC<=4999
(7,007 observations deleted)

.         
.         sort gvkeyyear

.         egen countstate=nvals(AntiSLAPP),by(gvkeyyear) 

.     drop if  countstate<=1 & AntiSLAPP==1       
(7,193 observations deleted)

.         
.     label variable logempnum  "Natural logarithm of the number of employees at plant level "

.     label variable pop  "Natural logarithm of the total population at county level"     

.         label variable dpop  "The growth rate of population at county level"    

.     label variable logipc  "Natural logarithm of the per capita income at county level" 

.         label variable dipc  "The growth rate of per capita income at county level"     

.     label variable employmentrate  "The ratio of employed labors at county level"       

.         label variable demp  "Change of Employment Rate at county level"

.         label variable size  "Natural logarithm of the market value of equity of parent firm"   

.         label variable lev  "Leverage ratio of parent firm"     

.         label variable xrdat  "R&D expense of parent firm, scaled by total asset"       

.         label variable tang  "Net property, plant, and equipment of parent firm, scaled by total asse"  

.         
.         drop if mi(lev)
(0 observations deleted)

.         drop if mi(xrdat)
(0 observations deleted)

.         drop if mi(tang)
(23 observations deleted)

.         drop if mi(size)
(2,513 observations deleted)

.         drop if mi(pop)
(2,873 observations deleted)

.         drop if mi(logipc)
(0 observations deleted)

.         drop if mi(employmentrate)
(35 observations deleted)

.         drop if mi(dpop) 
(0 observations deleted)

.         drop if mi(dipc)
(0 observations deleted)

.         drop if mi(demp)
(16 observations deleted)

.         drop if mi(logempnum)
(0 observations deleted)

.         
.         * Construct variables for Parallel trend and emission ********
. 
.     gen AntiSLAPP_post1=1 if year==(eventyear+1)
(111,482 missing values generated)

.         replace AntiSLAPP_post1=0 if mi(AntiSLAPP_post1)
(111,482 real changes made)

.         gen AntiSLAPP_post2=1 if year==(eventyear+2)
(111,576 missing values generated)

.         replace AntiSLAPP_post2=0 if mi(AntiSLAPP_post2)
(111,576 real changes made)

.         gen AntiSLAPP_post3=1 if year==(eventyear+3)
(111,684 missing values generated)

.         replace AntiSLAPP_post3=0 if mi(AntiSLAPP_post3)
(111,684 real changes made)

.         gen AntiSLAPP_post4=1 if year==(eventyear+4)
(111,739 missing values generated)

.         replace AntiSLAPP_post4=0 if mi(AntiSLAPP_post4)        
(111,739 real changes made)

.         gen AntiSLAPP_post5=1 if year>=(eventyear+5)
(85,090 missing values generated)

.         replace AntiSLAPP_post5=0 if mi(AntiSLAPP_post5)
(85,090 real changes made)

.     gen AntiSLAPP_pre1=1 if year==(eventyear-1)
(111,059 missing values generated)

.         replace AntiSLAPP_pre1=0 if mi(AntiSLAPP_pre1)  
(111,059 real changes made)

.     gen AntiSLAPP_pre2=1 if year==(eventyear-2)
(110,979 missing values generated)

.         replace AntiSLAPP_pre2=0 if mi(AntiSLAPP_pre2)
(110,979 real changes made)

.         gen AntiSLAPP_pre3=1 if year==(eventyear-3)
(111,378 missing values generated)

.         replace AntiSLAPP_pre3=0 if mi(AntiSLAPP_pre3)
(111,378 real changes made)

.         gen AntiSLAPP_pre4=1 if year==(eventyear-4)
(111,615 missing values generated)

.         replace AntiSLAPP_pre4=0 if mi(AntiSLAPP_pre4)
(111,615 real changes made)

.         gen AntiSLAPP_pre5=1 if year<=(eventyear-5)
(51,396 missing values generated)

.         replace AntiSLAPP_pre5=0 if mi(AntiSLAPP_pre5)
(51,396 real changes made)

. 
.     gen emission=ln(trbp)
(9 missing values generated)

.         rename avgcumprod_n production

.         
.         *identify and remove singletons under facility and gvkey#year
.     bysort faclld: gen faclld_count = _N

.         drop if faclld_count == 1
(1,660 observations deleted)

.         drop faclld_count

. 
.         bysort gvkey year: gen gvkey_year_count = _N

.         drop if gvkey_year_count == 1
(4,493 observations deleted)

.         drop gvkey_year_count

. 
.     *construct AntiSlAPP strength  based on current rating in 2021
.         gen Rating_Overall="A+" if inlist(state,"CA")
(101,210 missing values generated)

.         replace Rating_Overall="A" if inlist(state,"GA","NV","OK","TN","VT","NY")
(12,958 real changes made)

.         replace Rating_Overall="A-" if inlist(state,"KS","LA","OR","TX","WA")
(14,076 real changes made)

.         replace Rating_Overall="B+" if inlist(state,"CT","IN")
(6,420 real changes made)

.         replace Rating_Overall="B" if inlist(state,"CO","DC","RI")
(1,197 real changes made)

.         replace Rating_Overall="C+" if inlist(state,"VA")
(1,122 real changes made)

.         replace Rating_Overall="C" if inlist(state,"AR")
(2,504 real changes made)

.         replace Rating_Overall="C-" if inlist(state,"FL")
(2,252 real changes made)

.         replace Rating_Overall="D+" if inlist(state,"IL","MA","MN")
(8,646 real changes made)

.         replace Rating_Overall="D" if inlist(state,"HI","ME","MD","NM")
(1,682 real changes made)

.         replace Rating_Overall="D-" if inlist(state,"AZ","DE","MO","NE","PA","UT")
(11,371 real changes made)

. 
.         gen FedY=0

.         replace FedY=1 if inlist(state,"CA","NV","OR","KS","OK","VT","ME","FL","HI")
(14,489 real changes made)

. 
.         gen FedU=0

.         replace FedU=1 if inlist(state,"WA","UT","AZ","CO","NE","MO","AR","LA","RI")
(12,968 real changes made)

.         replace FedU=1 if inlist(state,"IL","MI","IN","WV","PA","VA","MA","CT","MD")
(26,183 real changes made)

.         replace FedU=1 if inlist(state,"DE")
(359 real changes made)

.         replace FedU=1 if inlist(state,"TN","NY")
(7,071 real changes made)

. 
.         gen FedN=0

.         replace FedN=1 if inlist(state,"NM","TX","GA","MN","DC")
(13,314 real changes made)

.         
.         gen Sall=0

.         gen x=substr(Rating_Overall,1,1)
(38,982 missing values generated)

.         replace Sall=1 if inlist(x,"A","B","C")+inlist(Rating_Overall,"C+","C" ,"C-")
(47,082 real changes made)

. 
.         *Adjust some states with rating changes using their historical rating
.         * Washington
.         replace Sall=0 if state=="WA" & (year<2010 | year>=2015)
(973 real changes made)

.         * Nevada (2013)
.         replace Sall=0 if state=="NV" & year<2013
(205 real changes made)

.         * NewYork (2020)
.         replace Sall=0 if state=="NY" & year<2020
(3,214 real changes made)

.         *Virginia (2017)
.         replace Sall=0 if state=="VA" & year<2017
(1,041 real changes made)

.         * Tennessee (2019)
.     replace Sall=0 if state=="TN" & year<2019
(3,630 real changes made)

. 
.         gen Wall=1-Sall

.         gen FedN_or_Wall=(FedN+Wall>0)

.         gen FedN_and_Wall=FedN*Wall

.         gen FedN_and_Sall=FedN*Sall

.         gen FedY_and_Sall=(FedY+FedU)*Sall

.         gen FedY_and_Wall=(FedY+FedU)*Wall

.         gen FedY1_and_Sall=FedY*Sall

.         gen allother=(1-FedY_and_Sall)

.         replace allother=0 if mi(eventyear)
(38,982 real changes made)

. 
.         winsor2 emission size lev logat  xrddum xrdat tang logipc employmentrate dpop dipc demp logempnum  tr  avgcum
> prod, replace cuts(1 99) by (year)

.         
.         egen fac11d=group(TRIFD fips)

.         
.         save auto1, replace
file auto1.dta saved

.         
.         *summary statistics
.         *Note summary statistics for Local Pollution News, Local Protest News and Air Quality index (AQI) are reporte
> d when factiva dataset and AQI dataset are imported for regression
.         sysuse auto1, clear

.     use auto1.dta 

.         drop if mi(emission)
(8 observations deleted)

.         
.     winsor2 emission production emission_intensity logempnum , replace cuts(1 99) by (year)

.         tabstat AntiSLAPP emission production emission_intensity logempnum   , stat(n mean sd  p25 p50 p75 ) c(s)

    Variable |         N      Mean        SD       p25       p50       p75
-------------+------------------------------------------------------------
   AntiSLAPP |    107755  .3484386   .476478         0         0         1
    emission |    107755  8.794282   3.67427  6.694959  9.548311  11.34451
  production |    107755   1.29562  3.759739  .8966667         1  1.081633
emission_i~y |    107755  8.415711  3.645635  6.282049  9.111735  10.96813
   logempnum |    107755  5.077388   1.48658  4.248495  5.220356  6.070738
--------------------------------------------------------------------------

.         
.         sort fipsyear

.     quietly by fipsyear : gen dup = cond(_N==1,0,_n) 

.         drop if dup>1
(70,403 observations deleted)

.         winsor2 pop logipc employmentrate dpop dipc  demp , replace cuts(1 99) by (year)

.    
.         tabstat pop logipc employmentrate dpop dipc  demp , stat(n mean sd  p25 p50 p75 ) c(s)  

    Variable |         N      Mean        SD       p25       p50       p75
-------------+------------------------------------------------------------
         pop |     37352  11.20977  1.243894  10.30571  11.04777  11.99369
      logipc |     37352  10.21843  .3736137  9.933241  10.21852   10.4919
employment~e |     37352  .9403189  .0238025      .927      .945      .958
        dpop |     37352  .0071378  .0121519 -.0010493  .0052462  .0132303
        dipc |     37352  .0381415  .0302074  .0219407  .0382074  .0551291
        demp |     37352 -.0003521  .0126957     -.004      .002      .007
--------------------------------------------------------------------------

.         
.         sysuse auto1, clear

.     use auto1.dta 

.         sort gvkeyyear

.     quietly by gvkeyyear : gen dup = cond(_N==1,0,_n) 

.         drop if dup>1
(94,852 observations deleted)

.         drop if mi(roa)
(3 observations deleted)

.         drop if mi(margin)
(7 observations deleted)

.         drop if mi(assetturn)
(0 observations deleted)

.         winsor2  roa margin assetturn sg  , replace cuts(1 99) by (year)

.         tabstat size lev xrdat tang roa margin assetturn sg  , stat(n mean sd  p25 p50 p75 ) c(s)

    Variable |         N      Mean        SD       p25       p50       p75
-------------+------------------------------------------------------------
        size |     12901  7.146692  2.018559  5.805227  7.150344  8.429659
         lev |     12901  .3029051  .2045776  .1704656  .2750762  .3977298
       xrdat |     12901  .0186205  .0248615         0    .01003  .0262574
        tang |     12901  .3426799  .1842355  .2036484  .3070043  .4422016
         roa |     12901  .0515499  .0700991  .0195262  .0539209  .0887609
      margin |     12901  .0931547  .0689235   .052095  .0881497  .1273153
   assetturn |     12901  1.279866   .620082  .8686474  1.172698  1.559073
          sg |     12895  .0747739  .2011204 -.0186939  .0556987  .1390802
--------------------------------------------------------------------------

. 
.         *Regression of Table 2
.         
.         sysuse auto1, clear

.     use auto1.dta 

. 
.     est clear

.         local control " logempnum  pop logipc employmentrate dpop dipc  demp " 

.         local control1 "logempnum  pop logipc employmentrate dpop dipc demp size lev xrdat tang "

.         local title " Emission  Emission Emission "

. 
.         eststo:quietly reghdfe emission AntiSLAPP , absorb(fac11d  gvkey#year ) tol(0.0001) vce(cluster statecode naics2
>  )
(est1 stored)

.     est store Column1

.         
.         eststo:quietly reghdfe emission  AntiSLAPP  `control'   , absorb( fac11d   gvkey#year ) tol(0.0001) vce(cluster 
> statecode  naics2 )
(est2 stored)

.         est store Column2

. 
.         eststo:quietly reghdfe emission  AntiSLAPP  `control1'   , absorb( fac11d   naics2#year ) tol(0.0001) vce(cluste
> r statecode naics2 )
(est3 stored)

.         est store Column3

. 
. *  Display the combined table in Stata's results window
.     esttab Column1 Column2 Column3 , ///
>     b(%9.3f) t(%9.2f)                 ///
>     star(* 0.1 ** 0.05 *** 0.01)     ///
>     order(  AntiSLAPP   ) ///
>     varwidth(20)                        ///
>     stats(N_full r2_a, fmt(%9.0f %9.3f) labels("Observations" "Adj. R-squared")) ///
>     title("Table 2: The Main Effect of the Anti-SLAPP Statutes on Toxic Emissions ") ///
>     mtitles("Model 1" "Model 2" "Model 3")      

Table 2: The Main Effect of the Anti-SLAPP Statutes on Toxic Emissions 
--------------------------------------------------------------------
                              (1)             (2)             (3)   
                          Model 1         Model 2         Model 3   
--------------------------------------------------------------------
AntiSLAPP                  -0.086          -0.085***       -0.099** 
                          (-4.42)         (-3.50)         (-2.32)   

logempnum                                   0.042***        0.053** 
                                           (2.94)          (2.32)   

pop                                         0.108          -0.017   
                                           (0.41)         (-0.31)   

logipc                                      0.702           0.340   
                                           (1.58)          (1.45)   

employmentrate                              1.306**         1.720   
                                           (2.84)          (1.63)   

dpop                                        2.066           1.204   
                                           (1.31)          (0.54)   

dipc                                        0.143           0.029   
                                           (0.38)          (0.08)   

demp                                       -0.027           0.424   
                                          (-0.06)          (1.23)   

size                                                        0.030   
                                                           (1.64)   

lev                                                         0.221** 
                                                           (2.65)   

xrddum                                                      0.018   
                                                           (0.21)   

tang                                                       -0.473***
                                                          (-3.61)   

_cons                       8.826          -1.165           3.485   
                        (1311.49)         (-0.26)          (1.24)   
--------------------------------------------------------------------
Observations               107755          107755          107755   
Adj. R-squared              0.840           0.840           0.816   
--------------------------------------------------------------------
t statistics in parentheses
* p<0.1, ** p<0.05, *** p<0.01

.         
.         *Regression of Table 3
.         
.         sysuse auto1, clear

.     use auto1.dta 

.         
.         qui foreach x in FedY_and_Sall  allother FedY_and_Wall FedN_and_Sall FedN_and_Wall {

. 
.     est clear

.         local control " logempnum pop logipc employmentrate dpop dipc demp "

.         local title " (1) (2)  "

. 
.         eststo:quietly reghdfe emission FedY_and_Sall allother   `control', absorb( fac11d gvkey#year) tol(0.0001) vc
> e(cluster statecode naics2)
(est1 stored)

.         est store Column1

.         test FedY_and_Sall=allother     

 ( 1)  FedY_and_Sall - allother = 0

       F(  1,    18) =    4.41
            Prob > F =    0.0502

.         
.         eststo:quietly reghdfe emission FedY_and_Sall  FedY_and_Wall FedN_and_Sall FedN_and_Wall   `control', absorb(
>  fac11d gvkey#year)  tol(0.0001) vce(cluster statecode naics2)
(est2 stored)

.         est store Column2

.         test FedY_and_Sall=FedY_and_Wall

 ( 1)  FedY_and_Sall - FedY_and_Wall = 0

       F(  1,    18) =    0.57
            Prob > F =    0.4597

.     test FedY_and_Sall=FedN_and_Sall

 ( 1)  FedY_and_Sall - FedN_and_Sall = 0

       F(  1,    18) =   23.42
            Prob > F =    0.0001

.     test FedY_and_Sall=FedN_and_Wall

 ( 1)  FedY_and_Sall - FedN_and_Wall = 0

       F(  1,    18) =    0.63
            Prob > F =    0.4393

. 
. *  Display the combined table in Stata's results window
.     esttab Column1 Column2 , ///
>     b(%9.3f) t(%9.2f)                 ///
>     star(* 0.10 ** 0.05 *** 0.01)     ///
>     order(FedY_and_Sall allother FedY_and_Wall FedN_and_Sall FedN_and_Wall  ) ///
>     varwidth(20)                        ///
>     stats(N_full r2_a, fmt(%9.0f %9.3f) labels("Observations" "Adj. R-squared")) ///
>     title("Table 3: The Role of Strength of Anti-SLAPP Statutes and Applicability in Federal Courts ") ///
>     mtitles("Model 1" "Model 2" )

Table 3: The Role of Strength of Anti-SLAPP Statutes and Applicability in Federal Courts 
----------------------------------------------------
                              (1)             (2)   
                          Model 1         Model 2   
----------------------------------------------------
FedY_and_Sall              -0.161***       -0.162***
                         (-11.16)         (-7.90)   

allother                   -0.045                   
                          (-1.04)                   

FedY_and_Wall                              -0.078   
                                          (-0.84)   

FedN_and_Sall                               0.051   
                                           (0.85)   

FedN_and_Wall                               0.029   
                                           (0.12)   

logempnum                   0.043***        0.041** 
                           (2.91)          (2.37)   

pop                         0.096          -0.081   
                           (0.37)         (-1.71)   

logipc                      0.686           0.671   
                           (1.57)          (1.44)   

employmentrate              1.267**         0.746   
                           (2.72)          (1.19)   

dpop                        2.033           1.594   
                           (1.31)          (0.90)   

dipc                        0.138           0.162   
                           (0.37)          (0.43)   

demp                       -0.023          -0.193   
                          (-0.05)         (-0.45)   

_cons                      -0.815           1.990   
                          (-0.18)          (0.41)   
----------------------------------------------------
Observations               107755          107755   
Adj. R-squared              0.840           0.839   
----------------------------------------------------
t statistics in parentheses
* p<0.10, ** p<0.05, *** p<0.01

.                 
.   *Regression of Table 4
.         
.     *Column  1 of Table 4***********
. 
.         sysuse auto1, clear

.     use auto1.dta 

.         
.         gen dum=1 if posviewratio>0.56 & posviewratio~=.
(67,718 missing values generated)

.     replace dum=0 if posviewratio<=0.56 & posviewratio~=.
(67,718 real changes made)

.         gen dum1=1 if negviewratio<0.38 & negviewratio~=.
(70,122 missing values generated)

.     replace dum1=0 if negviewratio>=0.38 & negviewratio~=.
(70,122 real changes made)

. 
.         gen AntiSLAPP_HighE= AntiSLAPP *dum1

.         
.     est clear

.         local control " logempnum  pop logipc employmentrate dpop dipc  demp " 

. 
.         eststo:quietly reghdfe emission   AntiSLAPP  AntiSLAPP_HighE   `control', absorb( fac11d  gvkey#year) tol(0.0
> 001) vce(cluster statecode naics2)
(est1 stored)

.     est store Column1

. 
. *  Display in Stata's results window
.     esttab Column1  , ///
>     b(%9.3f) t(%9.2f)                 ///
>     star(* 0.10 ** 0.05 *** 0.01)     ///
>     order(  AntiSLAPP  AntiSLAPP_HighE ) ///
>     varwidth(20)                        ///
>     stats(N_full r2_a, fmt(%9.0f %9.3f) labels("Observations" "Adj. R-squared")) ///
>     title("Table 4-Column1: The Role of the Abilities and Incentives of Local Stakeholders to Pressure Polluters ") /
> //
>     mtitles("Emission" )        

Table 4-Column1: The Role of the Abilities and Incentives of Local Stakeholders to Pressure Polluters 
------------------------------------
                              (1)   
                         Emission   
------------------------------------
AntiSLAPP                  -0.045   
                          (-1.07)   

AntiSLAPP_HighE            -0.117** 
                          (-2.16)   

logempnum                   0.042***
                           (2.91)   

pop                         0.084   
                           (0.32)   

logipc                      0.695   
                           (1.57)   

employmentrate              1.238** 
                           (2.18)   

dpop                        1.948   
                           (1.27)   

dipc                        0.153   
                           (0.41)   

demp                       -0.036   
                          (-0.08)   

_cons                      -0.740   
                          (-0.17)   
------------------------------------
Observations               107755   
Adj. R-squared              0.840   
------------------------------------
t statistics in parentheses
* p<0.10, ** p<0.05, *** p<0.01

.         
.         **Column 2 of Table 4***********        
.         
.         sysuse auto1, clear

.     use auto1.dta 

.     *median value of the unique county level pre-NGOs scaled by county population is 20.32616
.         gen dum=1 if nccs_pop_premean> 20.32616 & nccs_pop_premean!=.
(64,956 missing values generated)

.         replace dum=0 if mi(dum)
(64,956 real changes made)

.         gen AntiSLAPP_HighNGOs= AntiSLAPP*dum
.         drop if state=="WA" // drop it because no pre-NGO data for WA adoption           

.                 
.     est clear

.         local control " logempnum  pop logipc employmentrate dpop dipc  demp " 

. 
.         eststo:quietly  reghdfe emission  AntiSLAPP   AntiSLAPP_HighNGOs   `control', absorb( fac11d  gvkey#year ) to
> l(0.0001) vce(cluster  statecode naics2)
(est1 stored)

.     est store Column1

.         
. *  Display in Stata's results window
.     esttab Column1  , ///
>     b(%9.3f) t(%9.2f)                 ///
>     star(* 0.10 ** 0.05 *** 0.01)     ///
>     order(  AntiSLAPP  AntiSLAPP_HighNGOs ) ///
>     varwidth(20)                        ///
>     stats(N_full r2_a, fmt(%9.0f %9.3f) labels("Observations" "Adj. R-squared")) ///
>     title("Table 4-Column2: The Role of the Abilities and Incentives of Local Stakeholders to Pressure Polluters ") /
> //
>     mtitles("Emission" )        

Table 4-Column2: The Role of the Abilities and Incentives of Local Stakeholders to Pressure Polluters 
------------------------------------
                              (1)   
                         Emission   
------------------------------------
AntiSLAPP                   0.073   
                           (1.33)   

AntiSLAPP_HighNGOs         -0.195** 
                          (-2.28)   

logempnum                   0.039** 
                           (2.63)   

pop                        -0.079   
                          (-1.44)   

logipc                      0.677   
                           (1.37)   

employmentrate              0.761   
                           (0.95)   

dpop                        1.437   
                           (0.81)   

dipc                        0.091   
                           (0.23)   

demp                       -0.167   
                          (-0.31)   

_cons                       1.874   
                           (0.37)   
------------------------------------
Observations               106616  
Adj. R-squared              0.838   
------------------------------------
t statistics in parentheses
* p<0.10, ** p<0.05, *** p<0.01

.         
.    *Column 3 of Table 4
.         
.         sysuse auto1, clear

.     use auto1.dta 

. 
.     drop if year>2015
(12,805 observations deleted)

. 
.     gen dum=1 if numlocalnews_pre>1 & numlocalnews_pre!=.
(77,270 missing values generated)

.         replace dum=0 if  mi(dum)
(77,270 real changes made)

.         gen AntiSLAPP_LocalNewspaper = AntiSLAPP *dum

. 
.     est clear

.         local control " logempnum  pop logipc employmentrate dpop dipc  demp " 

. 
.         eststo:quietly  reghdfe emission  AntiSLAPP   AntiSLAPP_LocalNewspaper    `control', absorb( fac11d gvkey#yea
> r) tol(0.00001) vce(cluster statecode naics2)
(est1 stored)

.     est store Column1

.         
. *  Display in Stata's results window
.     esttab Column1  , ///
>     b(%9.3f) t(%9.2f)                 ///
>     star(* 0.10 ** 0.05 *** 0.01)     ///
>     order( AntiSLAPP  AntiSLAPP_LocalNewspaper ) ///
>     varwidth(20)                        ///
>     stats(N_full r2_a, fmt(%9.0f %9.3f) labels("Observations" "Adj. R-squared")) ///
>     title("Table 4-Column 3: The Role of the Abilities and Incentives of Local Stakeholders to Pressure Polluters ") 
> ///
>     mtitles("Emission" )        

Table 4-Column 3: The Role of the Abilities and Incentives of Local Stakeholders to Pressure Polluters 
------------------------------------
                              (1)   
                         Emission   
------------------------------------
AntiSLAPP                  -0.061   
                          (-1.59)   

AntiSLAPP_LocalNew~r       -0.118***
                          (-3.73)   

logempnum                   0.039** 
                           (2.48)   

pop                         0.008   
                           (0.04)   

logipc                      0.810   
                           (1.57)   

employmentrate              1.031** 
                           (2.28)   

dpop                        1.212   
                           (0.76)   

dipc                        0.254   
                           (0.63)   

demp                       -0.543   
                          (-1.18)   

_cons                      -0.620   
                          (-0.13)   
------------------------------------
Observations                94956   
Adj. R-squared              0.841   
------------------------------------
t statistics in parentheses
* p<0.10, ** p<0.05, *** p<0.01

.         
.  *Column 4 of Table 4
. 
.         sysuse auto1, clear

.     use auto1.dta 

. 
.         drop if FYEAR>2016
(7,452 observations deleted)

.         drop if FYEAR<1992
(12,576 observations deleted)

. ***1992-2016 because the SRI fund construction is based KLD rating which is available with large coverage from 1992-2
> 016 ****************
.         gen dum=1 if numsrilocal1_pre>=1 & numsrilocal1_pre!=.  
(63,635 missing values generated)

.         replace dum=0 if  mi(dum)
(63,635 real changes made)

.         gen AntiSLAPP_LocalSRMF = AntiSLAPP *dum

. 
.         est clear

.         local control " logempnum  pop logipc employmentrate dpop dipc  demp " 

. 
.         eststo:quietly  reghdfe emission  AntiSLAPP   AntiSLAPP_LocalSRMF    `control', absorb( fac11d gvkey#year) to
> l(0.00001) vce(cluster statecode naics2)
(est1 stored)

.     est store Column1

.         
. *  Display in Stata's results window
.     esttab Column1  , ///
>     b(%9.3f) t(%9.2f)                 ///
>     star(* 0.10 ** 0.05 *** 0.01)     ///
>     order(  AntiSLAPP  AntiSLAPP_LocalSRMF ) ///
>     varwidth(20)                        ///
>     stats(N_full r2, fmt(%9.0f %9.3f) labels("Observations" "Adj. R-squared")) ///
>     title("Table 4-Column 4: The Role of the Abilities and Incentives of Local Stakeholders to Pressure Polluters ") 
> ///
>     mtitles("Emission" )        

Table 4-Column 4: The Role of the Abilities and Incentives of Local Stakeholders to Pressure Polluters 
------------------------------------
                              (1)   
                         Emission   
------------------------------------
AntiSLAPP                   0.012   
                           (0.19)   

AntiSLAPP_LocalSRMF        -0.155** 
                          (-2.86)   

logempnum                   0.021   
                           (1.59)   

pop                        -0.070   
                          (-0.35)   

logipc                      0.680*  
                           (1.83)   

employmentrate              1.840*  
                           (1.83)   

dpop                        0.419   
                           (0.28)   

dipc                        0.254   
                           (0.63)   

demp                        0.104   
                           (0.18)   

_cons                       0.715   
                           (0.23)   
------------------------------------
Observations                87734   
Adj. R-squared              0.887   
------------------------------------
t statistics in parentheses
* p<0.10, ** p<0.05, *** p<0.01

.         
. * Regression of Table 5
. 
.         clear all

.         cd "D:\Data\AntiSLAPP"
D:\Data\AntiSLAPP

.         set matsize 800
set matsize ignored.
    Matrix sizes are no longer limited by c(matsize) in modern Statas.  Matrix sizes are now limited by edition of
    Stata.  See limits for more details.

.         set more off

.         use finalsample_factivanewstest.dta     

.         destring, replace 
TRIFD: contains nonnumeric characters; no replace
year already numeric; no replace
P_CO_DB_NUM: contains nonnumeric characters; no replace
P_CO_NAME: contains nonnumeric characters; no replace
state: contains nonnumeric characters; no replace
COUNTY: contains nonnumeric characters; no replace
federal: contains nonnumeric characters; no replace
avgprodratio already numeric; no replace
avgcumprod already numeric; no replace
avgcumprod_n already numeric; no replace
avgtrprod already numeric; no replace
avgtrprod_n already numeric; no replace
nchemicals already numeric; no replace
tr already numeric; no replace
onsite_tr already numeric; no replace
offsite_tr already numeric; no replace
transfer already numeric; no replace
prodwste already numeric; no replace
treat already numeric; no replace
recycle already numeric; no replace
recover already numeric; no replace
air already numeric; no replace
land already numeric; no replace
water already numeric; no replace
other already numeric; no replace
naics4: all characters numeric; replaced as int
PRIMARY_naics: all characters numeric; replaced as long
PRIMARY_SIC: contains nonnumeric characters; no replace
DunsNumber already numeric; no replace
hqduns already numeric; no replace
HQCompany1: contains nonnumeric characters; no replace
City_his: contains nonnumeric characters; no replace
State_his: contains nonnumeric characters; no replace
HQDunsChange: contains nonnumeric characters; no replace
fips already numeric; no replace
emp already numeric; no replace
sales already numeric; no replace
sic2_nets already numeric; no replace
emp_lag1 already numeric; no replace
sales_lag1 already numeric; no replace
hqdunslink already numeric; no replace
logempnum already numeric; no replace
logtremp already numeric; no replace
logtremp1 already numeric; no replace
fips1: all characters numeric; replaced as long
trbp_on already numeric; no replace
a263 already numeric; no replace
a246 already numeric; no replace
trbp already numeric; no replace
nccs already numeric; no replace
nccs_lag1 already numeric; no replace
negviewratio already numeric; no replace
posviewratio already numeric; no replace
gvkey: all characters numeric; replaced as long
DATADATE already numeric; no replace
FYEAR already numeric; no replace
INCORP: contains nonnumeric characters; no replace
LOC: contains nonnumeric characters; no replace
SIC: all characters numeric; replaced as int
(17461 missing values generated)
pstate: contains nonnumeric characters; no replace
TIC: contains nonnumeric characters; no replace
CUSIP: contains nonnumeric characters; no replace
CIK: all characters numeric; replaced as long
(18780 missing values generated)
AT already numeric; no replace
CEQ already numeric; no replace
DLC already numeric; no replace
DLTT already numeric; no replace
INVT already numeric; no replace
SEQ already numeric; no replace
COGS already numeric; no replace
DP already numeric; no replace
DVC already numeric; no replace
DVP already numeric; no replace
IB already numeric; no replace
OIBDP already numeric; no replace
XRD already numeric; no replace
OANCF already numeric; no replace
PRCC_F already numeric; no replace
SICH already numeric; no replace
at_lag1 already numeric; no replace
mv already numeric; no replace
size already numeric; no replace
logat already numeric; no replace
lev already numeric; no replace
tang already numeric; no replace
xrdat already numeric; no replace
oi already numeric; no replace
xrddum already numeric; no replace
headqstate_yb: contains nonnumeric characters; no replace
numlocalnews_lag1 already numeric; no replace
numlocalnews already numeric; no replace
n_county already numeric; no replace
eventdate already numeric; no replace
eventyear already numeric; no replace
eventstate already numeric; no replace
AntiSLAPP already numeric; no replace
avgotherAntiSLAPP already numeric; no replace
avgotherAntiSLAPP_nonHQ already numeric; no replace
wavgotherAntiSLAPP already numeric; no replace
wavgotherAntiSLAPP_nonHQ already numeric; no replace
fips2 already numeric; no replace
RegEPA already numeric; no replace
FedEPA already numeric; no replace
EPAregcode already numeric; no replace
EPAregfips: all characters numeric; replaced as long
(3919 missing values generated)
numsrilocal1 already numeric; no replace
nccs_pop_premean already numeric; no replace
nccs_pop_pre already numeric; no replace
numlocalnews_pre already numeric; no replace
numsrilocal1_pre already numeric; no replace
statepol_pre1 already numeric; no replace
statepro_pre1 already numeric; no replace
statepol_pre2 already numeric; no replace
statepro_pre2 already numeric; no replace
statepol already numeric; no replace
statepro already numeric; no replace
statepol_post1 already numeric; no replace
statepro_post1 already numeric; no replace
statepol_post2 already numeric; no replace
statepro_post2 already numeric; no replace
count already numeric; no replace
Population already numeric; no replace
ipc already numeric; no replace
employmentrate already numeric; no replace
population_lag1 already numeric; no replace
ipc_lag1 already numeric; no replace
employmentrate_lag1 already numeric; no replace
numplant already numeric; no replace
sumtrbp already numeric; no replace
sumemp already numeric; no replace
sumemp_lag1 already numeric; no replace
sumtrbp_lag1 already numeric; no replace
firmemp already numeric; no replace
firmemp_lag1 already numeric; no replace
polnews1 already numeric; no replace
polnews already numeric; no replace
pronews already numeric; no replace
dpop already numeric; no replace
dipc already numeric; no replace
demp already numeric; no replace
pop already numeric; no replace
logipc already numeric; no replace
logsumemp already numeric; no replace
factivafirm already numeric; no replace
Initialyear already numeric; no replace
rank_sumtrbp_lag1 already numeric; no replace
margin already numeric; no replace
assetturn already numeric; no replace
roa already numeric; no replace
sg already numeric; no replace
pollutiondum1 already numeric; no replace
protestdum1 already numeric; no replace

.         isid year state gvkey

.          
.         drop if year<1990
(0 observations deleted)

.         drop if year>2020
(0 observations deleted)

. 
.     egen statecode=group(state )

.         gen naics2=int(naics4/100)

.     drop if mi(logat)
(17,461 observations deleted)

.         drop if mi(at_lag1)
(477 observations deleted)

.         drop if mi(size)
(2,277 observations deleted)

.         drop if mi(lev)
(0 observations deleted)

.         drop if mi(xrddum)
(0 observations deleted)

.         drop if mi(tang)
(115 observations deleted)

.         drop if mi(pop)
(0 observations deleted)

.         drop if mi(logipc)
(0 observations deleted)

.         drop if mi(employmentrate)
(0 observations deleted)

.         drop if mi(dpop)
(0 observations deleted)

.         drop if mi(dipc)
(0 observations deleted)

.         drop if mi(demp)
(0 observations deleted)

.         drop if mi(logsumemp)
(0 observations deleted)

.         
.     save auto, replace
file auto.dta saved

.                 
.         sysuse auto, clear
(1978 automobile data)

.     use auto.dta 

.     drop if mi(factivafirm)
(21,006 observations deleted)

.     drop if year<Initialyear
(849 observations deleted)

.     label variable logsumemp "Natural logarithm of the total number of employees at firm-state level "

.     label variable pop  "Natural logarithm of the total population at state level"      

.         label variable dpop  "The growth rate of population at state level"     

.     label variable logipc  "Natural logarithm of the per capita income at state level"  

.         label variable dipc  "The growth rate of per capita income at state level"      

.     label variable employmentrate  "The ratio of employed labors at state level"        

.         label variable demp  "Change of Employment Rate at state level" 

.         
.     winsor2  logsumemp  logipc employmentrate dpop dipc  demp , replace cuts(1 99) 

. 
.         gen HighEmissionDum=rank_sumtrbp_lag1-1
(37,635 missing values generated)

.         replace HighEmissionDum=0 if mi(rank_sumtrbp_lag1)
(37,635 real changes made)

. 
.         gen AntiSLAPP_HighEmission= AntiSLAPP*HighEmissionDum

. 
.         gen pollutiondum=1 if polnews>=1
(74,492 missing values generated)

.         replace pollutiondum=0 if mi(pollutiondum)      
(74,492 real changes made)

.         
.         gen protestdum=1 if pronews>=1
(75,345 missing values generated)

.         replace protestdum=0 if mi(protestdum)
(75,345 real changes made)

. 
.     tabstat pollutiondum1 protestdum1 , stat(n mean sd p25 p50 p75 ) c(s)

    Variable |         N      Mean        SD       p25       p50       p75
-------------+------------------------------------------------------------
pollutiond~1 |     75458  .0638033  .2444039         0         0         0
 protestdum1 |     75458  .0239279  .1528256         0         0         0
--------------------------------------------------------------------------

. 
.     est clear

.     local control " logsumemp pop logipc employmentrate dpop dipc  demp   "  

.         local fix "   statecode#gvkey year"

. 
.         eststo:quietly reghdfe pollutiondum  AntiSLAPP      `control', absorb( `fix') tol(0.0001)   vce(cluster  stat
> ecode)
(est1 stored)

.     est store Column1

.         
.         eststo:quietly reghdfe pollutiondum  AntiSLAPP  AntiSLAPP_HighEmission    `control'               , absorb( `
> fix') tol(0.0001)  vce(cluster statecode)
(est2 stored)

.         est store Column2

. 
.         eststo:quietly reghdfe protestdum   AntiSLAPP      `control', absorb(`fix') tol(0.0001) vce(cluster statecode
> )  
(est3 stored)

.         est store Column3

.         
.         eststo:quietly reghdfe protestdum   AntiSLAPP  AntiSLAPP_HighEmission     `control'               , absorb(`f
> ix') tol(0.0001)  vce(cluster statecode)
(est4 stored)

.         est store Column4

.         
. *  Display the combined table in Stata's results window
.     esttab Column1 Column2 Column3 Column4, ///
>     b(%9.3f) t(%9.2f)                 ///
>     star(* 0.10 ** 0.05 *** 0.01)     ///
>     order(  AntiSLAPP  AntiSLAPP_HighEmission ) ///
>     varwidth(30)                        ///
>     stats(N_full r2_a,fmt(%9.0f %9.3f) labels("Observations" "Adj. R-squared")) ///
>     title("Table 5: Anti-SLAPP Statutes and Pollution News and Protest News Reported by Local Newspapers ") ///
>     mtitles("Local Pollution News" "Local Pollution News" "Local Protest News" "Local Protest News")    

Table 5: Anti-SLAPP Statutes and Pollution News and Protest News Reported by Local Newspapers 
----------------------------------------------------------------------------------------------
                                        (1)             (2)             (3)             (4)   
                               Local Poll~s    Local Poll~s    Local Prot~s    Local Prot~s   
----------------------------------------------------------------------------------------------
AntiSLAPP                             0.006*         -0.006           0.002          -0.001   
                                     (1.83)         (-1.31)          (0.62)         (-0.52)   

AntiSLAPP_HighEmission                                0.018***                        0.006** 
                                                     (3.44)                          (2.17)   

logsumemp                            -0.006***       -0.006***       -0.004***       -0.004***
                                    (-4.76)         (-4.66)         (-3.21)         (-3.16)   

pop                                   0.066           0.039          -0.005          -0.007   
                                     (1.58)          (1.14)         (-0.22)         (-0.32)   

logipc                               -0.005          -0.042           0.007           0.001   
                                    (-0.09)         (-0.82)          (0.23)          (0.02)   

employmentrate                        0.192           0.258**         0.162**         0.180** 
                                     (1.21)          (2.14)          (2.13)          (2.40)   

dpop                                  0.064          -0.223          -0.128          -0.159   
                                     (0.23)         (-0.83)         (-0.52)         (-0.64)   

dipc                                  0.119*          0.031           0.033           0.030   
                                     (1.92)          (0.40)          (0.87)          (0.77)   

demp                                  0.044           0.408*          0.057           0.058   
                                     (0.16)          (1.72)          (0.41)          (0.41)   

_cons                                -1.049          -0.266          -0.051           0.039   
                                    (-1.15)         (-0.41)         (-0.14)          (0.11)   
----------------------------------------------------------------------------------------------
Observations                          75458           75458           75458           75458   
Adj. R-squared                        0.255           0.255           0.256           0.256   
----------------------------------------------------------------------------------------------
t statistics in parentheses
* p<0.10, ** p<0.05, *** p<0.01

.                 
.   *Regression of Table 6
.   
.         sysuse auto1, clear

.     use auto1.dta 

. 
.         gen delta1=statepol_post1+statepol_post2-statepol_pre1-statepol_pre2
(38,982 missing values generated)

.         gen delta2=statepro_post1+statepro_post2-statepro_pre1-statepro_pre2
(38,982 missing values generated)

. 
.         /* median for news change is 5 and 19*******/
.         gen dum1=1 if delta1>5 & delta1!=.
(73,408 missing values generated)

.         gen dum2=1 if delta2>19 & delta2!=.
(75,281 missing values generated)

.         replace dum1=0 if mi(eventyear) | delta1<=5
(73,408 real changes made)

.         replace dum2=0 if mi(eventyear) | delta2<=19 
(75,281 real changes made)

. 
.         gen AntiSLAPP_HighPollutionNewsInc= AntiSLAPP *dum1

.         gen AntiSLAPP_HighProtestNewsInc= AntiSLAPP *dum2

. 
.     est clear

.         local control " logempnum  pop logipc employmentrate dpop dipc  demp " 

. 
.         eststo:quietly reghdfe emission  AntiSLAPP  AntiSLAPP_HighPollutionNewsInc `control' , absorb( fac11d   gvkey
> #year)  vce(cluster  statecode naics2 )
(est1 stored)

.     est store Column1

.         
.         eststo:quietly reghdfe emission  AntiSLAPP  AntiSLAPP_HighProtestNewsInc `control' , absorb( fac11d   gvkey#y
> ear)  vce(cluster  statecode naics2 )
(est2 stored)

.         est store Column2

. 
. *  Display the combined table in Stata's results window
.     esttab Column1 Column2 , ///
>     b(%9.3f) t(%9.2f)                 ///
>     star(* 0.10 ** 0.05 *** 0.01)     ///
>     order(  AntiSLAPP  AntiSLAPP_HighPollutionNewsInc AntiSLAPP_HighProtestNewsInc ) ///
>     varwidth(20)                        ///
>     stats(N_full r2_a, fmt(%9.0f %9.3f) labels("Observations" "Adj. R-squared")) ///
>     title("Table 6: Increase in Local Public Pressure and the Effect of Anti-SLAPP Statutes on Toxic Emissions") ///
>     mtitles("Emission" "Emission" )     

Table 6: Increase in Local Public Pressure and the Effect of Anti-SLAPP Statutes on Toxic Emissions
----------------------------------------------------
                              (1)             (2)   
                         Emission        Emission   
----------------------------------------------------
AntiSLAPP                   0.014           0.011   
                           (0.35)          (0.43)   

AntiSLAPP_HighPoll~c       -0.192***                
                          (-5.70)                   

AntiSLAPP_HighProt~c                       -0.193***
                                          (-4.29)   

logempnum                   0.042***        0.042** 
                           (2.95)          (2.79)   

pop                         0.102           0.094   
                           (0.41)          (0.38)   

logipc                      0.698           0.684   
                           (1.57)          (1.55)   

employmentrate              1.409**         1.197** 
                           (2.60)          (2.45)   

dpop                        2.073           1.983   
                           (1.26)          (1.23)   

dipc                        0.163           0.146   
                           (0.44)          (0.40)   

demp                        0.042          -0.079   
                           (0.10)         (-0.18)   

_cons                      -1.144          -0.702   
                          (-0.26)         (-0.16)   
----------------------------------------------------
Observations               107755          107755   
Adj. R-squared              0.840           0.840   
----------------------------------------------------
t statistics in parentheses
* p<0.10, ** p<0.05, *** p<0.01

.                  
. *regression of Table 7
.         
.         sysuse auto1, clear

.     use auto1.dta 

. 
.         gen AntiSLAPP_hishead=1 if headqstate_yb=="AZ" & year>2006
(107,499 missing values generated)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="AR" & year>2005
(671 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="CA" & year>1992
(3,575 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="CO" & year>2019
(21 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="CT" & year>2018
(212 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="DC" & year>2011
(147 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="DE" & year>1992
(1,382 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="FL" & year>2000
(438 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="GA" & year>1996
(1,919 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="HI" & year>2002
(0 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="IL" & year>2007
(3,261 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="IN" & year>1998
(1,021 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="KS" & year>2016
(0 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="LA" & year>1999
(164 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="MA" & year>1994
(1,359 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="MD" & year>2004
(358 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="ME" & year>1995
(0 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="MN" & year>1994
(3,042 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="MO" & year>2004
(1,428 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="NE" & year>1994
(2,287 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="NM" & year>2001
(0 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="NV" & year>1993
(0 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="NY" & year>1993
(5,065 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="OK" & year>2014
(59 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="OR" & year>2002
(619 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="PA" & year>2001
(5,025 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="RI" & year>1993
(563 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="TN" & year>1997
(1,205 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="TX" & year>2011
(2,225 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="UT" & year>2001
(179 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="VT" & year>2006
(0 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="VA" & year>2007
(545 real changes made)

.         replace AntiSLAPP_hishead=1 if headqstate_yb=="WA" & year>1989
(1,377 real changes made)

.         replace AntiSLAPP_hishead=0 if mi(AntiSLAPP_hishead)
(69,352 real changes made)

.         
.         gen samestate_flag=1 if state==headqstate_yb
(90,095 missing values generated)

.         replace samestate_flag=0 if mi(samestate_flag)
(90,095 real changes made)

.         gen AntiSLAPP_headdif=AntiSLAPP_hishead*(1-samestate_flag)

. 
.         replace wavgotherAntiSLAPP=0 if mi(wavgotherAntiSLAPP)
(1,418 real changes made)

.         replace wavgotherAntiSLAPP_nonHQ=0 if mi(wavgotherAntiSLAPP_nonHQ)
(1,418 real changes made)

. 
.         rename wavgotherAntiSLAPP AntiSLAPP_Other

.         rename AntiSLAPP_headdif AntiSLAPP_Other_HQ     

.         rename wavgotherAntiSLAPP_nonHQ AntiSLAPP_Other_nonHQ   

.         
.     est clear

.         local control1 " logempnum  pop logipc employmentrate dpop dipc demp size lev xrdat tang"

. 
.         eststo:quietly reghdfe emission AntiSLAPP AntiSLAPP_Other `control1', absorb( fac11d  naics2#year ) tol(0.000
> 1) vce(cluster  statecode naics2 )
(est1 stored)

.     est store Column1

.         
.         eststo:quietly reghdfe emission AntiSLAPP AntiSLAPP_Other_HQ AntiSLAPP_Other_nonHQ `control1', absorb( fac11d
>   naics2#year ) tol(0.0001) vce(cluster  statecode naics2 )
(est2 stored)

.         est store Column2

. 
. *  Display the combined table in Stata's results window
.     esttab Column1 Column2 , ///
>     b(%9.3f) t(%9.2f)                 ///
>     star(* 0.10 ** 0.05 *** 0.01)     ///
>     order( AntiSLAPP AntiSLAPP_Other AntiSLAPP_Other_HQ AntiSLAPP_Other_nonHQ) ///
>     varwidth(20)  ///
>         stats(N_full r2_a, fmt(%9.0f %9.3f) labels("Observations" "Adj. R-squared")) ///
>     title("Table 7: Effects of Anti-SLAPP Statutes in Other States") ///
>     mtitles("Emission" "Emission")      

Table 7: Effects of Anti-SLAPP Statutes in Other States
----------------------------------------------------
                              (1)             (2)   
                         Emission        Emission   
----------------------------------------------------
AntiSLAPP                  -0.099**        -0.096*  
                          (-2.17)         (-2.13)   

AntiSLAPP_Other             0.004                   
                           (0.03)                   

AntiSLAPP_Other_HQ                          0.038   
                                           (1.40)   

AntiSLAPP_Other_no~Q                       -0.024   
                                          (-0.19)   

logempnum                   0.053**         0.053** 
                           (2.36)          (2.33)   

pop                        -0.018          -0.018   
                          (-0.32)         (-0.32)   

logipc                      0.341           0.338   
                           (1.45)          (1.43)   

employmentrate              1.726           1.727   
                           (1.50)          (1.52)   

dpop                        1.189           1.192   
                           (0.54)          (0.54)   

dipc                        0.028           0.028   
                           (0.08)          (0.08)   

demp                        0.428           0.426   
                           (1.45)          (1.31)   

size                        0.030           0.029   
                           (1.71)          (1.64)   

lev                         0.220**         0.221** 
                           (2.49)          (2.45)   

xrdat                      -0.396          -0.381   
                          (-0.25)         (-0.23)   

tang                       -0.465***       -0.463***
                          (-4.15)         (-4.18)   

_cons                       3.489           3.511   
                           (1.24)          (1.23)   
----------------------------------------------------
Observations               107755          107755   
Adj. R-squared              0.816           0.816   
----------------------------------------------------
t statistics in parentheses
* p<0.10, ** p<0.05, *** p<0.01

. 
.         *Regression of Table 8
. 
.         sysuse auto1, clear

.     use auto1.dta 

. 
.     drop if mi(emission)
(8 observations deleted)

.     winsor2 production emission emission_intensity, replace cuts(1 99) by (state)

. 
.     est clear

.         local control "  logempnum pop logipc employmentrate dpop dipc demp " 

. 
.         eststo:quietly reghdfe production AntiSLAPP `control', absorb( fac11d gvkey#year) tol(0.0001)  vce(cluster  s
> tatecode naics2)
(est1 stored)

.     est store Column1

.         
.         eststo:quietly reghdfe emission_intensity AntiSLAPP `control', absorb( fac11d gvkey#year) tol(0.0001) vce(clu
> ster  statecode naics2)
(est2 stored)

.         est store Column2

. 
. *  Display the combined table in Stata's results window
.     esttab Column1 Column2 , ///
>     b(%9.3f) t(%9.2f)                 ///
>     star(* 0.10 ** 0.05 *** 0.01)     ///
>     order(AntiSLAPP) ///
>     varwidth(20)                        ///
>     stats(N_full r2_a, fmt(%9.0f %9.3f) labels("Observations" "Adj. R-squared")) ///
>     title("Table 8: Source of the Emissions Decrease") ///
>     mtitles("Production" "Emission Intensity" )         

Table 8: Source of the Emissions Decrease
----------------------------------------------------
                              (1)             (2)   
                       Production    Emission I~y   
----------------------------------------------------
AntiSLAPP                   0.041          -0.079***
                           (0.91)         (-3.22)   

logempnum                  -0.007           0.047** 
                          (-0.60)          (2.17)   

pop                         0.008           0.107   
                           (0.10)          (0.61)   

logipc                     -0.078           0.521   
                          (-0.56)          (1.08)   

employmentrate             -2.573           1.894***
                          (-1.52)          (6.02)   

dpop                       -1.314           1.720   
                          (-0.87)          (0.91)   

dipc                       -0.138          -0.109   
                          (-0.59)         (-0.39)   

demp                       -0.941          -0.871   
                          (-0.70)         (-0.89)   

_cons                       4.333          -0.238   
                           (1.24)         (-0.04)   
----------------------------------------------------
Observations               107755          107755   
Adj. R-squared              0.287           0.825   
----------------------------------------------------
t statistics in parentheses
* p<0.10, ** p<0.05, *** p<0.01

.         
.    *Regression of Table 9
. 
.     clear all

. 
.         cd "D:\Data\AntiSLAPP"
D:\Data\AntiSLAPP

.         set matsize 800
set matsize ignored.
    Matrix sizes are no longer limited by c(matsize) in modern Statas.  Matrix sizes are now limited by edition of
    Stata.  See limits for more details.

.         set more off

.         use finalsample_factivanews.dta

.         
.         destring, replace 
TRIFD: contains nonnumeric characters; no replace
year already numeric; no replace
P_CO_DB_NUM: contains nonnumeric characters; no replace
P_CO_NAME: contains nonnumeric characters; no replace
state: contains nonnumeric characters; no replace
COUNTY: contains nonnumeric characters; no replace
federal: contains nonnumeric characters; no replace
avgprodratio already numeric; no replace
avgcumprod already numeric; no replace
avgcumprod_n already numeric; no replace
avgtrprod already numeric; no replace
avgtrprod_n already numeric; no replace
nchemicals already numeric; no replace
tr already numeric; no replace
onsite_tr already numeric; no replace
offsite_tr already numeric; no replace
transfer already numeric; no replace
prodwste already numeric; no replace
treat already numeric; no replace
recycle already numeric; no replace
recover already numeric; no replace
air already numeric; no replace
land already numeric; no replace
water already numeric; no replace
other already numeric; no replace
naics4: all characters numeric; replaced as int
PRIMARY_naics: all characters numeric; replaced as long
PRIMARY_SIC: contains nonnumeric characters; no replace
DunsNumber already numeric; no replace
hqduns already numeric; no replace
HQCompany1: contains nonnumeric characters; no replace
City_his: contains nonnumeric characters; no replace
State_his: contains nonnumeric characters; no replace
HQDunsChange: contains nonnumeric characters; no replace
fips already numeric; no replace
emp already numeric; no replace
sales already numeric; no replace
sic2_nets already numeric; no replace
emp_lag1 already numeric; no replace
sales_lag1 already numeric; no replace
hqdunslink already numeric; no replace
logempnum already numeric; no replace
logtremp already numeric; no replace
logtremp1 already numeric; no replace
fips1: all characters numeric; replaced as long
trbp_on already numeric; no replace
a263 already numeric; no replace
a246 already numeric; no replace
trbp already numeric; no replace
nccs already numeric; no replace
nccs_lag1 already numeric; no replace
negviewratio already numeric; no replace
posviewratio already numeric; no replace
gvkey: all characters numeric; replaced as long
DATADATE already numeric; no replace
FYEAR already numeric; no replace
INCORP: contains nonnumeric characters; no replace
LOC: contains nonnumeric characters; no replace
SIC: all characters numeric; replaced as int
(17461 missing values generated)
pstate: contains nonnumeric characters; no replace
TIC: contains nonnumeric characters; no replace
CUSIP: contains nonnumeric characters; no replace
CIK: all characters numeric; replaced as long
(18780 missing values generated)
AT already numeric; no replace
CEQ already numeric; no replace
DLC already numeric; no replace
DLTT already numeric; no replace
INVT already numeric; no replace
SEQ already numeric; no replace
COGS already numeric; no replace
DP already numeric; no replace
DVC already numeric; no replace
DVP already numeric; no replace
IB already numeric; no replace
OIBDP already numeric; no replace
XRD already numeric; no replace
OANCF already numeric; no replace
PRCC_F already numeric; no replace
SICH already numeric; no replace
at_lag1 already numeric; no replace
mv already numeric; no replace
size already numeric; no replace
logat already numeric; no replace
lev already numeric; no replace
tang already numeric; no replace
xrdat already numeric; no replace
oi already numeric; no replace
xrddum already numeric; no replace
headqstate_yb: contains nonnumeric characters; no replace
numlocalnews_lag1 already numeric; no replace
numlocalnews already numeric; no replace
n_county already numeric; no replace
eventdate already numeric; no replace
eventyear already numeric; no replace
eventstate already numeric; no replace
AntiSLAPP already numeric; no replace
avgotherAntiSLAPP already numeric; no replace
avgotherAntiSLAPP_nonHQ already numeric; no replace
wavgotherAntiSLAPP already numeric; no replace
wavgotherAntiSLAPP_nonHQ already numeric; no replace
fips2 already numeric; no replace
RegEPA already numeric; no replace
FedEPA already numeric; no replace
EPAregcode already numeric; no replace
EPAregfips: all characters numeric; replaced as long
(3919 missing values generated)
numsrilocal1 already numeric; no replace
nccs_pop_premean already numeric; no replace
nccs_pop_pre already numeric; no replace
numlocalnews_pre already numeric; no replace
numsrilocal1_pre already numeric; no replace
statepol_pre1 already numeric; no replace
statepro_pre1 already numeric; no replace
statepol_pre2 already numeric; no replace
statepro_pre2 already numeric; no replace
statepol already numeric; no replace
statepro already numeric; no replace
statepol_post1 already numeric; no replace
statepro_post1 already numeric; no replace
statepol_post2 already numeric; no replace
statepro_post2 already numeric; no replace
count already numeric; no replace
Population already numeric; no replace
ipc already numeric; no replace
employmentrate already numeric; no replace
population_lag1 already numeric; no replace
ipc_lag1 already numeric; no replace
employmentrate_lag1 already numeric; no replace
numplant already numeric; no replace
sumtrbp already numeric; no replace
sumemp already numeric; no replace
sumemp_lag1 already numeric; no replace
sumtrbp_lag1 already numeric; no replace
firmemp already numeric; no replace
firmemp_lag1 already numeric; no replace
polnews1 already numeric; no replace
polnews already numeric; no replace
pronews already numeric; no replace
dpop already numeric; no replace
dipc already numeric; no replace
demp already numeric; no replace
pop already numeric; no replace
logipc already numeric; no replace
logsumemp already numeric; no replace
factivafirm already numeric; no replace
Initialyear already numeric; no replace
rank_sumtrbp_lag1 already numeric; no replace
margin already numeric; no replace
assetturn already numeric; no replace
roa already numeric; no replace
sg already numeric; no replace
pollutiondum1 already numeric; no replace
protestdum1 already numeric; no replace

.         isid year state gvkey

.          
.         drop if year<1990
(0 observations deleted)

.         drop if year>2020
(0 observations deleted)

.  
.     egen statecode=group(state )

.         gen naics2=int(naics4/100)

.         gen w=sumemp/firmemp_lag1
(4,437 missing values generated)

. 
.     gen logfirmemp=ln(1+firmemp_lag1)
(4,437 missing values generated)

.     drop if mi(logat)
(17,461 observations deleted)

.         drop if mi(at_lag1)
(477 observations deleted)

.         drop if mi(size)
(2,277 observations deleted)

.         drop if mi(lev)
(0 observations deleted)

.         drop if mi(xrdat)
(0 observations deleted)

.         drop if mi(tang)
(115 observations deleted)

.         drop if mi(pop)
(0 observations deleted)

.         drop if mi(logipc)
(0 observations deleted)

.         drop if mi(employmentrate)
(0 observations deleted)

.         drop if mi(dpop)
(0 observations deleted)

.         drop if mi(dipc)
(0 observations deleted)

.         drop if mi(demp)
(0 observations deleted)

.         drop if mi(logfirmemp)
(3,771 observations deleted)

.         drop if PRCC_F<=2       
(2,361 observations deleted)

.         drop if mi(roa)
(4,753 observations deleted)

.         drop if mi(margin)
(52 observations deleted)

.         drop if mi(assetturn)
(0 observations deleted)

.         
.         label variable logfirmemp  "Natural logarithm of the number of employees at parent firm level "

.     label variable pop  "Natural logarithm of the total population at state level"      

.         label variable dpop  "The growth rate of population at state level"     

.     label variable logipc  "Natural logarithm of the per capita income at state level"  

.         label variable dipc  "The growth rate of per capita income at state level"      

.     label variable employmentrate  "The ratio of employed labors at state level"        

.         label variable demp  "Change of Employment Rate at state level" 

.         
.     winsor2  logfirmemp logat lev  xrdat tang size  pop  logipc employmentrate dpop dipc  demp  , replace cuts(1 99) 

.     winsor2  margin roa assetturn  sg , replace cuts(1 99) 

.                 
.     est clear

.         local control " logfirmemp pop logipc employmentrate dpop dipc demp  size lev xrdat tang  " 

.         local title " roa margin  assetturn  sg"

.         local fix "  gvkey#statecode  year "

. 
.         eststo:quietly reghdfe roa AntiSLAPP `control' [aweight=w], absorb( `fix')  vce(cluster   statecode)
(est1 stored)

.     est store Column1

.         
.         eststo:quietly reghdfe margin AntiSLAPP `control'  [aweight=w], absorb( `fix')  vce(cluster  statecode )
(est2 stored)

.         est store Column2

. 
.         eststo:quietly reghdfe assetturn AntiSLAPP `control' [aweight=w], absorb( `fix') vce(cluster  statecode)     
>    
(est3 stored)

.         est store Column3

.         
.         eststo:quietly reghdfe sg AntiSLAPP  `control' [aweight=w], absorb( `fix')  vce(cluster   statecode )
(est4 stored)

.         est store Column4

.         
. *  Display the combined table in Stata's results window
.     esttab Column1 Column2 Column3 Column4, ///
>     b(%9.3f) t(%9.2f)                 ///
>     star(* 0.10 ** 0.05 *** 0.01)     ///
>     order(AntiSLAPP ) ///
>     varwidth(30)                        ///
>     stats(N_full r2_a, fmt(%9.0f %9.3f) labels("Observations" "Adj. R-squared")) ///
>     title("Table 9: Anti-SLAPP Statutes and Firm Performance ") ///
>     mtitles("ROA" "Profit Margin" "Asset Turnover" "Sales Growth")      

Table 9: Anti-SLAPP Statutes and Firm Performance 
----------------------------------------------------------------------------------------------
                                        (1)             (2)             (3)             (4)   
                                        ROA    Profit Mar~n    Asset Turn~r    Sales Growth   
----------------------------------------------------------------------------------------------
AntiSLAPP                            -0.007**        -0.006**        -0.012          -0.003   
                                    (-2.34)         (-2.11)         (-0.86)         (-0.39)   

logfirmemp                           -0.004***       -0.004***        0.000          -0.013** 
                                    (-2.86)         (-3.14)          (0.02)         (-2.60)   

pop                                   0.028*          0.015           0.057          -0.009   
                                     (1.76)          (0.70)          (0.50)         (-0.18)   

logipc                               -0.020          -0.015          -0.307*         -0.115   
                                    (-0.53)         (-0.32)         (-1.81)         (-1.40)   

employmentrate                        0.114           0.125           0.789           0.007   
                                     (1.21)          (1.00)          (1.51)          (0.03)   

dpop                                 -0.143          -0.319          -1.733          -0.718   
                                    (-0.61)         (-1.16)         (-1.26)         (-1.12)   

dipc                                  0.103**         0.108           0.940***        0.097   
                                     (2.33)          (1.38)          (4.06)          (0.65)   

demp                                 -0.075          -0.106          -1.381**        -0.150   
                                    (-0.47)         (-0.37)         (-2.59)         (-0.44)   

size                                  0.007***        0.015***       -0.089***       -0.034***
                                     (4.42)          (7.45)        (-12.76)         (-7.60)   

lev                                  -0.016**         0.012**        -0.060**        -0.091***
                                    (-2.64)          (2.42)         (-2.43)         (-8.16)   

xrdat                                 0.110*         -0.107*          1.150***        0.091   
                                     (1.78)         (-1.76)          (4.99)          (0.66)   

tang                                 -0.006           0.004           0.015           0.038*  
                                    (-0.75)          (0.30)          (0.35)          (1.81)   

_cons                                -0.287          -0.142           3.517           1.725*  
                                    (-0.73)         (-0.28)          (1.42)          (1.85)   
----------------------------------------------------------------------------------------------
Observations                          87356           87356           87356           87326   
Adj. R-squared                        0.509           0.687           0.868           0.369   
----------------------------------------------------------------------------------------------
t statistics in parentheses
* p<0.10, ** p<0.05, *** p<0.01

.         
. *Regression of Table 10
.    clear all

.         cd "D:\Data\AntiSLAPP"
D:\Data\AntiSLAPP

.         set matsize 800
set matsize ignored.
    Matrix sizes are no longer limited by c(matsize) in modern Statas.  Matrix sizes are now limited by edition of
    Stata.  See limits for more details.

.         set more off

. 
.         use finalsample_AQI.dta

. 
.         destring, replace 
year already numeric; no replace
state: contains nonnumeric characters; no replace
fips already numeric; no replace
pop already numeric; no replace
logipc already numeric; no replace
employmentrate already numeric; no replace
dpop already numeric; no replace
dipc already numeric; no replace
demp already numeric; no replace
AntiSLAPP already numeric; no replace
aqi_ozone already numeric; no replace
aqi already numeric; no replace

.         isid year fips

.         drop if mi(pop)
(667 observations deleted)

.         drop if mi(logipc)
(0 observations deleted)

.         drop if mi(employmentrate)
(12 observations deleted)

.         drop if mi(dpop)
(0 observations deleted)

.         drop if mi(dipc)
(0 observations deleted)

.         drop if mi(demp)
(6 observations deleted)

.     gen logaqi=ln(aqi)
(46 missing values generated)

.         
.         label variable logaqi   "The natural logarithm of the median value of the daily Air Quality Index at county l
> evel "     
note: label truncated to 80 characters

.         
.         winsor2  logipc employmentrate dpop dipc demp , replace cuts(1 99)

.         
.         tabstat logaqi  , stat(n mean sd  p25 p50 p75   ) c(s)

    Variable |         N      Mean        SD       p25       p50       p75
-------------+------------------------------------------------------------
      logaqi |     22841   4.27672  .8050568  3.713572  4.369448   4.85203
--------------------------------------------------------------------------

. 
.     est clear

.         local control " pop logipc employmentrate dpop dipc  demp " 

. 
.         eststo:quietly reghdfe logaqi AntiSLAPP, absorb( fips year ) tol(0.0001) vce(cluster  statecode)
(est1 stored)

.     est store Column1

.         
.         eststo:quietly reghdfe logaqi AntiSLAPP `control', absorb( fips year) tol(0.0001) vce(cluster  statecode)
(est2 stored)

.         est store Column2

. 
. *  Display the combined table in Stata's results window
.     esttab Column1 Column2 , ///
>     b(%9.3f) t(%9.2f)                 ///
>     star(* 0.10 ** 0.05 *** 0.01)     ///
>     order(AntiSLAPP ) ///
>     varwidth(20)                        ///
>     stats(N_full r2_a, fmt(%9.0f %9.3f) labels("Observations" "Adj. R-squared")) ///
>     title("Table 10: Using the Emission data from the EPA Air Quality Monitoring Stations") ///
>     mtitles("Air Quality Index" "Air Quality Index" )                   

Table 10: Using the Emission data from the EPA Air Quality Monitoring Stations
----------------------------------------------------
                              (1)             (2)   
                     Air Qualit~x    Air Qualit~x   
----------------------------------------------------
AntiSLAPP                  -0.053**        -0.062** 
                          (-1.97)         (-2.30)   

pop                                         0.542***
                                           (4.77)   

logipc                                      0.240   
                                           (1.40)   

employmentrate                              1.192*  
                                           (1.81)   

dpop                                        0.360   
                                           (0.48)   

dipc                                        0.019   
                                           (0.11)   

demp                                        0.464   
                                           (0.90)   

_cons                       4.300***       -5.707***
                         (396.02)         (-2.60)   
----------------------------------------------------
Observations                22841           22841   
Adj. R-squared              0.762           0.766   
----------------------------------------------------
t statistics in parentheses
* p<0.10, ** p<0.05, *** p<0.01

.                 
. 
end of do-file

. 
