get file = 'C:\Documents and Settings\Amanda Geller\My Documents\R Files (DP)\replicating mocan\mocan data.sav'.

compute murderrate_m = cmur*1000/pop.
compute yearid = year-1977.
compute timesquared = yearid*yearid.
execute.

*varying effect of time, interacting with executions and custody, varying intercept.
MIXED
  murderrate_m  BY fips year  WITH yearid  exprJLEl amurtlag sentratl2 cust3lag deathrtl perblk repub unemp realinc mortal urban legal18 legal19 legal20 rat2034 rat3544 rat4554 rat55 
  /CRITERIA = CIN(95) MXITER(100) MXSTEP(5) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, ABSOLUTE) LCONVERGE(0, ABSOLUTE)
  PCONVERGE(0.000001, ABSOLUTE)
  /FIXED = yearid  yearid*cust3lag yearid*exprJLEl  exprJLEl amurtlag sentratl2 cust3lag deathrtl perblk repub unemp realinc mortal urban legal18 legal19 legal20 rat2034 rat3544 rat4554 rat55 | SSTYPE(3)
  /METHOD = REML
  /PRINT = solution descriptives
  /RANDOM intercept yearid | SUBJECT(fips) COVTYPE(AR1)
  /REGWGT = pop.

*interacting custody and executions with time and time squared.
MIXED
  murderrate_m  BY fips year  WITH yearid timesquared exprJLEl amurtlag sentratl2 cust3lag deathrtl perblk repub unemp realinc mortal urban legal18 legal19 legal20 rat2034 rat3544 rat4554 rat55 
  /CRITERIA = CIN(95) MXITER(100) MXSTEP(5) SCORING(1) SINGULAR(0.000000000001) HCONVERGE(0, ABSOLUTE) LCONVERGE(0, ABSOLUTE)
  PCONVERGE(0.000001, ABSOLUTE)
  /FIXED = yearid timesquared yearid*cust3lag timesquared*cust3lag yearid*exprJLEl timesquared*exprJLEl  exprJLEl amurtlag sentratl2 cust3lag deathrtl perblk repub unemp realinc mortal urban legal18 legal19 legal20 rat2034 rat3544 rat4554 rat55 | SSTYPE(3)
  /METHOD = REML
  /PRINT = solution descriptive
  /RANDOM intercept yearid timesquared| SUBJECT(fips) COVTYPE(AR1)
  /REGWGT = pop.