#######################################
# Fixed Effects Panel Data Model in R #
#######################################

# clean all variables
rm(list=ls())

# options and libraries
library(readxl)     # required for opening          
library(plm)        # required for panel data manipulation and  modeling

# Read data
setwd("c:/Documents/Classes/Econ4230/R_and_Stata/")
data = read_excel("data/murder.xlsx")

# Organizing the data as a panel
murderpd = pdata.frame(data, index=c("StateNo", "Year"))

# FE Regression
fit = lm(Murders ~ Execs + Unemp + factor(StateNo) -1, data=murderpd)
summary(fit)