# ==============================================================================
# Wooldridge - Capitulo 9 - Exemplo 9.3
# Update: 13/09/2022
# ==============================================================================

rm(list=ls())

# Pacotes
library(wooldridge)
library(tidyverse)
library(stargazer)

# Base de dados
data('crime2')
head(crime2)

# ------------------------------------------------------------------------------

# Definir data frame
df <- crime2[crime2$d87 == 1, ]
# df <- crime2 %>%
#   filter(d87 ==1)

# Regressao 1
reg1<-lm(lcrmrte ~ unem + llawexpc, data=df)
summary(reg1)

# Regressao 2 (DP var como proxy)
reg2<-lm(lcrmrte ~ unem + llawexpc + lcrmrt_1, data=df)
summary(reg2)

# Compara resultados
stargazer(reg1, reg2, type = 'text', keep.stat = c('n', 'rsq'))