log using c:\abenv\program.log, replace
set mem 32m
set matsize 200

use c:\abenv\dataset
describe
generate yr96=year==1996
generate yr97=year==1997
generate yr98=year==1998
generate yr99=year==1999
generate yr00=year==2000
generate yr01=year==2001
generate yr02=year==2002
generate qtr1=month==1 | month==2 | month==3
generate qtr2=month==4 | month==5 | month==6
generate qtr3=month==7 | month==8 | month==9
generate qtr4=month==10 | month==11 | month==12

generate pubcompreports=pubcomplaints*number
generate chempc=pubcomplaints*chemical
generate foodpc=pubcomplaints*foodandbeverage
generate fertpc=pubcomplaints*fertilizer
generate pulppc=pubcomplaints*pulppaper
generate petpc=pubcomplaints*oilandgas
generate incpc=averageincome*pubcomplaints
generate denpc=density*pubcomplaints
generate penfpc=pastenforcement*pubcomplaints
generate pincpc=pastincidents*pubcomplaints
generate imppc=impact*pubcomplaints
generate legpc=legislation*pubcomplaints
generate relpc=release*pubcomplaints
generate q1pc=qtr1*pubcomplaints
generate q2pc=qtr2*pubcomplaints
generate q3pc=qtr3*pubcomplaints
generate q4pc=qtr4*pubcomplaints
generate pastinvpc=pastinvestigations*pubcomplaints
generate linvestpc=lastinvestigated*pubcomplaints
generate lenfpc=lastenforced*pubcomplaints

probit resolution number pastincidents pastinvestigations pastenforcement pincpc pastinvpc penfpc lastinvestigated lastenforced linvestpc lenfpc density  denpc averageincome incpc albertaenvironment agency pubcomplaints chemical foodandbeverage fertilizer pulppaper chempc foodpc fertpc pulppc impact imppc legislation legpc pubcompreports yr97 yr98 yr99 yr00 yr01 yr02 qtr2 qtr3 qtr4

test pincpc pastinvpc penfpc linvestpc lenfpc pubcomplaints chempc foodpc fertpc pulppc imppc legpc pubcompreports denpc incpc

test pubcomplaints+pubcompreports+pincpc+pastinvpc+penfpc+linvestpc+lenfpc+denpc+incpc=0
test pubcomplaints+pubcompreports+imppc+pincpc+pastinvpc+penfpc+linvestpc+lenfpc+denpc+incpc=0
test pubcomplaints+pubcompreports+legpc+pincpc+pastinvpc+penfpc+linvestpc+lenfpc+denpc+incpc=0
test pubcomplaints+pubcompreports+chempc+pincpc+pastinvpc+penfpc+linvestpc+lenfpc+denpc+incpc=0
test pubcomplaints+pubcompreports+foodpc+pincpc+pastinvpc+penfpc+linvestpc+lenfpc+denpc+incpc=0
test pubcomplaints+pubcompreports+fertpc+pincpc+pastinvpc+penfpc+linvestpc+lenfpc+denpc+incpc=0
test pubcomplaints+pubcompreports+pincpc+pastinvpc+penfpc+linvestpc+lenfpc+denpc+incpc+pulppc=0

*tests within public complaints
test pastincidents+pincpc=0
test pastenforcement+penfpc=0
test pastinvestigations+pastinvpc=0
test lastinvestigated+linvestpc=0
test lastenforced+lenfpc=0
test number+pubcompreports=0

test impact+imppc=0
test legislation+legpc=0
test chemical+chempc=0
test foodandbeverage+foodpc=0
test fertilizer+fertpc=0
test pulppaper+pulppc=0

test impact-legislation=0
test chemical-foodandbeverage=0
test chemical-fertilizer=0
test chemical-pulppaper=0
test foodandbeverage-fertilizer=0
test foodandbeverage-pulppaper=0
test fertilizer-pulppaper=0

test impact+imppc-legislation-legpc=0
test chemical+chempc-foodandbeverage-foodpc=0
test chemical+chempc-fertilizer-fertpc=0
test chemical+chempc-pulppaper-pulppc=0
test foodandbeverage+foodpc-fertilizer-fertpc=0
test foodandbeverage+foodpc-pulppaper-pulppc=0
test fertilizer+fertpc-pulppaper-pulppc=0

*enforcement
probit enforcement pastincidents pastinvestigations pastenforcement lastinvestigated lastenforced density averageincome albertaenvironment agency pubcomplaints chemical foodandbeverage fertilizer pulppaper impact legislation number pubcompreports yr97 yr98 yr99 yr00 yr01 yr02 qtr2 qtr3 qtr4 if resolution==1
test pubcomplaint 
test albertaenvironment-pubcomplaints=0
test agency-pubcomplaints=0
test agency-albertaenvironment=0
test number+pubcompreports=0
test pubcomplaints+pubcompreports=0

drop _all 
log close
exit, clear