/*
Replication Program:
	Green, Strolovitch, and Wong (in press)
	"Defended Neighborhoods, Integration, and Racially Motivated Crime"
	American Journal of Sociology

Poisson & Negative binomial regression analyses
*/

@-----------------------------------------------@
@ Read in GAUSS Count Module and Gauss dataset  @
@-----------------------------------------------@
library count;
#include count.ext;
countset;

open f1 = "temp";
                     x = readr(f1,rowsf(f1));
                     makevars(x,0,getname("temp"));
                     f1 = close(f1);

print "
@-----------------------------------------------@
@ Replicate Table 5: Anti-Asian Incidents       @
@   with white or unknown offenders             @
@-----------------------------------------------@";

 dep1 = { newahc };
 ind1 = { pw80,chngasn,xasn };

_cnstart=2;
@ poisson regression provides start values for negative binomial regression @
screen off;
{ b,vc,llik } = poisson("temp",dep1,ind1);
call countprt(b,vc,llik);
_cnsval=b|-2;
@ negative binomial regression @
{ b,vc,llik } = negbin("temp",dep1,ind1,0);
screen on;
call countprt(b,vc,llik);

print "
    Gauss likelihood minus sumc(lnfact(y)) " (llik - sumc(lnfact(newahc)));

print "
@-----------------------------------------------@
@ Replicate Table 6: Anti-Latino Incidents      @
@   with white or unknown offenders             @
@-----------------------------------------------@";

 dep1 = { newlhc };
 ind1 = { pw80,chnghsp,xhsp };

_cnstart=2;
@ poisson regression provides start values for negative binomial regression @
screen off;
{ b,vc,llik } = poisson("temp",dep1,ind1);
call countprt(b,vc,llik);
_cnsval=b|-2;
@ negative binomial regression @
{ b,vc,llik } = negbin("temp",dep1,ind1,0);
screen on;
call countprt(b,vc,llik);

print "
    Gauss likelihood minus sumc(lnfact(y)) " (llik - sumc(lnfact(newlhc)));

print "
@-----------------------------------------------@
@ Replicate Table 7: Anti-Black Incidents       @
@   with white or unknown offenders             @
@-----------------------------------------------@";

 dep1 = { newbhc };
 ind1 = { pw80,chngblk,xblk };

_cnstart=2;
@ poisson regression provides start values for negative binomial regression @
screen off;
{ b,vc,llik } = poisson("temp",dep1,ind1);
call countprt(b,vc,llik);
_cnsval=b|-2;
@ negative binomial regression @
{ b,vc,llik } = negbin("temp",dep1,ind1,0);
screen on;
call countprt(b,vc,llik);

print "
    Gauss likelihood minus sumc(lnfact(y)) " (llik - sumc(lnfact(newbhc)));

print "
@-----------------------------------------------@
@ Replicate Table 5: Anti-Asian Incidents       @
@   with white or unknown offenders             @
@ INCLUDES LOG OF 1990 POPULATION AS CONTROL    @
@-----------------------------------------------@";

 dep1 = { newahc };
 ind1 = { pw80,chngasn,xasn,lntot90 };

_cnstart=2;
@ poisson regression provides start values for negative binomial regression @
screen off;
{ b,vc,llik } = poisson("temp",dep1,ind1);
call countprt(b,vc,llik);
_cnsval=b|-2;
@ negative binomial regression @
{ b,vc,llik } = negbin("temp",dep1,ind1,0);
screen on;
call countprt(b,vc,llik);

print "
@-----------------------------------------------@
@ Replicate Table 6: Anti-Latino Incidents      @
@   with white or unknown offenders             @
@ INCLUDES LOG OF 1990 POPULATION AS CONTROL    @
@-----------------------------------------------@";

 dep1 = { newlhc };
 ind1 = { pw80,chnghsp,xhsp,lntot90 };

_cnstart=2;
@ poisson regression provides start values for negative binomial regression @
screen off;
{ b,vc,llik } = poisson("temp",dep1,ind1);
call countprt(b,vc,llik);
_cnsval=b|-2;
@ negative binomial regression @
{ b,vc,llik } = negbin("temp",dep1,ind1,0);
screen on;
call countprt(b,vc,llik);

print "
@-----------------------------------------------@
@ Replicate Table 7: Anti-Black Incidents       @
@   with white or unknown offenders             @
@ INCLUDES LOG OF 1990 POPULATION AS CONTROL    @
@-----------------------------------------------@";

 dep1 = { newbhc };
 ind1 = { pw80,chngblk,xblk,lntot90 };

_cnstart=2;
@ poisson regression provides start values for negative binomial regression @
screen off;
{ b,vc,llik } = poisson("temp",dep1,ind1);
call countprt(b,vc,llik);
_cnsval=b|-2;
@ negative binomial regression @
{ b,vc,llik } = negbin("temp",dep1,ind1,0);
screen on;
call countprt(b,vc,llik);