################ ## Download zip files here: ## http://research.stlouisfed.org/fred2/categories/10 ## http://research.stlouisfed.org/fred2/categories/106 ################ ## Load Data ## Gross Domestic Product, 1 Decimal; Bil. of $; Q; SAAR; 2010-03-26 df.GDP <- read.csv("c:/GDP.csv") ## All Employees: Total Private Industries; Thous.; M; SA; 2010-04-02 df.all.Emp <- read.csv("c:/USPRIV.csv") ## Total Nonfarm Payrolls: All Employees; Thous.; M; SA; 2010-04-02 df.all.Wages <- read.csv("c:/PAYEMS.csv") ## Use merge to make Wages and Emp Quarterly df.Emp <- merge(df.all.Emp ,df.GDP,by.x = "DATE",by.y = "DATE",all.x = FALSE,all.y = FALSE ) df.Wages <- merge(df.all.Wages,df.GDP,by.x = "DATE",by.y = "DATE",all.x = FALSE,all.y = FALSE ) ## Get first and last date df.Emp[1,] df.Wages[1,] df.GDP[1,] df.Emp[NROW(df.Emp),] df.Wages[NROW(df.Wages),] df.GDP[NROW(df.GDP),] ##Get row count Emp.obs <- NROW(df.Emp) Wages.obs <- NROW(df.Wages) GDP.obs <- NROW(df.GDP) ## Take sub samples df.0.Emp <-df.Emp[(Emp.obs -Wages.obs +1):(Emp.obs-2),] df.1.Wages <-df.Wages[ 2 :(Wages.obs -1),] df.1.GDP <-df.GDP[(GDP.obs -Wages.obs +2):(GDP.obs-1),] ## Print first observation df.0.Emp[1,] df.1.Wages[1,] df.1.GDP[1,] ## Build final Data.frame df.Mdl <- data.frame(df.0.Emp$DATE, df.0.Emp$VALUE.x, df.1.GDP$VALUE, df.1.Wages$VALUE.x) names(df.Mdl)<- c('DATE','EMP','L1_GDP','L1_WAGES') ## Get sample SmpEnd<-NROW(df.Mdl)-50 OutSmpStart<-SmpEnd+1 OutSmpEnd <-NROW(df.Mdl) ## Do data transformations df.Mdl.2<-df.Mdl[2:OutSmpEnd,] df.Mdl.2$EMP_log_1D <- log(df.Mdl[2:OutSmpEnd,]$EMP ) - log(df.Mdl[1:(OutSmpEnd-1),]$EMP) df.Mdl.2$L1_GDP_log_1D <- log(df.Mdl[2:OutSmpEnd,]$L1_GDP) - log(df.Mdl[1:(OutSmpEnd-1),]$L1_GDP) df.Mdl.2$L1_WAGES_log_1D<- log(df.Mdl[2:OutSmpEnd,]$L1_WAGES) - log(df.Mdl[1:(OutSmpEnd-1),]$L1_WAGES) ################# ##END OF FILE #################