neoclassical.R

plotNeoclassicalGrowth = function()
{
	if ( exists("GDP") == FALSE ) { getSymbols("GDP", src="FRED") }
	if ( exists("GDPC96") == FALSE ) { getSymbols("GDPC96", src="FRED") }
	if ( exists("CLF16OV") == FALSE ) { getSymbols("CLF16OV", src="FRED") }
	if ( exists("GPDI") == FALSE ) { getSymbols("GPDI", src="FRED") }
	if ( exists("COFC") == FALSE ) { getSymbols("COFC", src="FRED") }

	a = na.omit(merge(Delt(CLF16OV, k=12), (GPDI-COFC)/GDP, Delt(GDPC96, k=4) - Delt(CLF16OV, k=12) - ((GPDI-COFC)/CPIAUCSL)/GDP)) * 100
	z = a
	names(a) = c("L", "K", "A")
	a = na.locf(cbind(USRECD, z))
	a = a[ .index(a) %in% .index(z) ]
	par(oma=c(0,0,0,0))
	par(xaxt="n", yaxt="n")
	barplot(a$USRECD, border="grey", col="grey")
	par(new=T, xaxt="s", yaxt="s")
	barplot(z, col=c("blue", "green", "red"))
	legend("topright", c("L(abour)", "K(apital)", "Productivity (A)ugmentation"), col=c("blue","green", "red"), pch=15, lty=0)
}