1. Regression with Panel Data(SW Chapter 10)

2. summary statistics

3. sum year

4. sum state

5. create year dummies • gen y82=(year==1982)

6. gen y83=(year==1983) • gen y84=(year==1984) • gen y85=(year==1985) • gen y86=(year==1986) • gen y87=(year==1987) • gen y88=(year==1988)

7. minimum legal drinking age • gen da18=(mlda<19) • gen da19=(mlda>=19)*(mlda<20) • gen da20=(mlda>=20)*(mlda<21) • gen da21=(mlda>=21)

8. fatality rate per 10,000 in the population • gen vfrall=10000*mrall

9. gen incperc=perinc/1000 • gen lincperc = ln(incperc) • gen vmilespd = vmiles/1000 • gen frmall = mrall/(vmiles/100000) • gen jailcom = ((jaild+comserd)>0) • gen mjailcom = (jaild==.)*(comserd==.)

10. mvdecode-Change numeric values to missing values • mvdecode mjailcom,mv(1)

12. replace jailcom = jailcom + mjailcom;

14. reg vfrall beertax if (year==1982), r

18. display

19. display "Adjusted Rsquared = " _result(8)

20. reg vfrall beertax if (year==1988), r • dis "Adjusted Rsquared = " _result(8)

21. preserve • gen dvfrall = vfrall-vfrall[_n-6] • gen dbtax = beertax-beertax[_n-6]

22. keep

24. regdvfralldbtax, r • dis "Adjusted Rsquared = " _result(8) • restore

25. egen-Extensions to generate

27. by state, sort : egen float beertaxtest = total(beertax*(year==1982)) • browse state year beertax beertax1982

28. egen beertax1988 = sum(beertax*(year==1988)) ,by(state) • egen vfrall1982 = sum(vfrall*(year==1982)) ,by(state) • egen vfrall1988 = sum(vfrall*(year==1988)) ,by(state) • gen beertax_diff=beertax1988-beertax1982 • gen vfrall_diff=vfrall1988-vfrall1982

30. areg-Linear regression with a large dummy-variable set

31. areg vfrall beertax, absorb(state) r

33. areg vfrall beertax y82 y83 y84 y85 y86 y87, absorb(state) r

34. reg vfrall beertax, r • dis "Adjusted Rsquared = " _result(8)

35. areg vfrall beertax, absorb(state) r

36. areg vfrall beertax y82 y83 y84 y85 y86 y87, absorb(state) r

37. test y82 y83 y84 y85 y86 y87

39. aregvfrallbeertax da18 da19 da20 jaildcomserdvmilespdunratelincperc y82 y83 y84 y85 y86 y87, absorb(state) r

40. test y82 y83 y84 y85 y86 y87 • test da18 da19 da20

41. test jaild comserd • test unrate lincperc

42. areg vfrall beertax da18 da19 da20 jaild comserd vmilespd y82 y83 y84 y85 y86 y87, absorb(state) r

43. test y82 y83 y84 y85 y86 y87 • test da18 da19 da20 • test jaild comserd

44. areg vfrall beertax mlda jailcom vmilespd unrate lincperc y82 y83 y84 y85 y86 y87, absorb(state) r

45. test y82 y83 y84 y85 y86 y87 • test unrate lincperc

46. aregvfrallbeertax da18 da19 da20 jaildcomserdvmilespdunratelincperc y82 y83 y84 y85 y86 y87, absorb(state) r cluster(state)

47. test y82 y83 y84 y85 y86 y87 • test da18 da19 da20

48. test jaild comserd • test unrate lincperc

49. Example: Traffic deaths and beer taxes in STATA