source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
View(log.likelihood.wrapper)
View(kalman.log.likelihood)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
print(Q)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
print(Q)
View(unpack.parameters.stage1)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
View(P.tp1t)
print(Q)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
print(theta[1])
print(theta[8])
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
Q
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/unpack.parameters.stage1.R')
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/unpack.parameters.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
x.data
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
x.data
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
y.data
plot(y.data[,1])
print(Q)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/unpack.parameters.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/unpack.parameters.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
plot(r.ph)
b.ph
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
b.ph
b2.constraint
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/unpack.parameters.stage1.R')
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
y.data
plot(y.data[,1])
plot(y.data[,2])
y.data[, 2] <- y.data[, 2] - 1:dim(y.data)[1] * (-0.1196)
plot(y.data[,2])
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/unpack.parameters.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
#------------------------------------------------------------------------------#
# World: Read in data, run estimation, and save output
#------------------------------------------------------------------------------#
# Read in the input data
wd.data <- read.table("inputData/rstar.data.wd.csv",
sep = ',', na.strings = ".", header=TRUE, stringsAsFactors=FALSE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
View(y.data)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/unpack.parameters.stage1.R')
View(rstar.stage3)
View(unpack.parameters.stage1)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
g.pot
View(log.output.hp.trend)
print(g.pot)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/unpack.parameters.stage1.R')
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
print(Q)
P.tp1t
View(P.ttm1.array)
View(Q)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
View(rstar.stage3)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
View(format.output)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
? nloptr
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
print(Q)
F
View(P.tm1tm1)
View(P.ttm1)
View(P.tt)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
xr
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
print(t)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
print(Q)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.R', echo=TRUE)
print(Q)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.R', echo=TRUE)
initial.parameters
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
initial.parameters
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
# Lower bound on b_2 parameter (slope of the Phillips curve)
b2.constraint <- 0.0002
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
initial.parameters
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
install.packages(seasonal)
seasonal
install.packages("seasonal")
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
seas(wd.data.gdp.log)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
library("x13binary", lib.loc="/Library/Frameworks/R.framework/Versions/3.4/Resources/library")
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
detach("package:seasonal", unload=TRUE)
library("seasonal", lib.loc="/Library/Frameworks/R.framework/Versions/3.4/Resources/library")
detach("package:x13binary", unload=TRUE)
library("x13binary", lib.loc="/Library/Frameworks/R.framework/Versions/3.4/Resources/library")
outputseasonal <- seas(
x = wd.log.output,
x11 = "",
arima.model = "(0 1 1)"
)
install.packages("seasonal")
library(seasonal)
install.packages("seasonal")
library(seasonal)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
library("x13binary", lib.loc="/Library/Frameworks/R.framework/Versions/3.4/Resources/library")
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
outputseasonal <- seas(
x = wd.log.output,
x11 = "",
arima.model = "(0 1 1)"
)
outputseasonal <- seas(
x = wd.log.output,
x11 = "",
arima.model = "(0 1 1)"
)
checkX13()
vignette(seasonal)
vignette("seas")
install.packages("seasonal")
detach("package:seasonal", unload=TRUE)
detach("package:x13binary", unload=TRUE)
install.packages("seasonal")
outputseasonal <- seas(wd.log.output)
library(seasonal)
outputseasonal <- seas(wd.log.output)
library("x13binary", lib.loc="/Library/Frameworks/R.framework/Versions/3.4/Resources/library")
library(seasonal)
outputseasonal <- seas(wd.log.output)
remove.packages("x13binary")
remove.packages("seasonal")
install.packages("seasonal")
install.packages("seasonal")
library(seasonal)
library(x13binary)
library(seasonal)
outputseasonal <- seas(wd.log.output)
x13path()
checkX13binary()
Sys.setenv(X13_PATH = "/Library/Frameworks/R.framework/Versions/3.4/Resources/library/x13binary/bin")
outputseasonal <- seas(wd.log.output)
timeseriesoutput <- ts(wd.log.output)
outputseasonal <- seas(timeseriesoutput)
timeseriesoutput <- ts(wd.log.output, frequency = 4, start = c(1980, 1))
outputseasonal <- seas(timeseriesoutput)
?seas
timeseriesoutput <- ts(wd.log.output, frequency = 4, start = c(1980, 1))
outputseasonal <- robust.seas(timeseriesoutput)
timeseriesoutput <- ts(wd.log.output, frequency = 4, start = c(1980, 1))
outputseasonal <- seas(timeseriesoutput, x11="")
View(outputseasonal)
timeseriesoutput <- ts(wd.log.output, frequency = 4, start = c(1980, 1))
os <- seas(timeseriesoutput, x11="")
plot(os)
ts.output <- ts(wd.log.output, frequency = 4, start = c(1980, 1))
os <- seas(ts.output, x11="")
fos <- final(os)
print(fos)
wd.log.outputsa <- data.frame(final(os))
View(wd.log.outputsa)
View(wd.log.outputsa)
plot(wd.log.outputsa)
ts.inflation <- ts(wd.log.inflation, frequency = 4, start = c(1980, 1))
inf.s <- seas(ts.inflation, x11="")
plot(inf.s)
ts.inflation <- ts(wd.log.inflation, frequency = 4, start = c(1980, 1))
ts.inflation <- ts(wd.inflation, frequency = 4, start = c(1980, 1))
inf.s <- seas(ts.inflation, x11="")
plot(inf.s)
view(inf.s)
install.packages(("seasonalview"))
view(inf.s)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
plot(wd.log.output)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
View(run.hlw.estimation)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.estimation.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
out.stage1.theta
out.stage1$theta
plot(out.stage1$potential.smoothed)
plot(out.stage1$potential.gap.smoothed)
out.stage1$potential.gap.smoothed
plot(out.stage1$output.gap.smoothed)
plot(out.stage1$output.gap.smoothed, type='o')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
out.stage1.theta
out.stage1$theta
plot(out.stage1$output.gap.smoothed, type='o')
plot(out.stage1$potential.smoothed, type='o')
View(rstar.stage2)
View(unpack.parameters.stage2)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.estimation.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
theta
View(unpack.parameters.stage2)
View(rstar.stage1)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
theta
? nloptr
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/rstar.stage1.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
theta
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
theta
theta
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
theta
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
theta
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
theta
theta
theta
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R')
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
b.oil <- solve(t(wd.oil) %*% wd.oil, t(wd.oil) %*% wd.inflation)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R')
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R')
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
wd.data <- read.table("inputData/rstar.data.wd.csv",
sep = ',', na.strings = ".", header=TRUE, stringsAsFactors=FALSE)
wd.data <- read.table("inputData/rstar.data.wd.csv",
sep = ',', na.strings = ".", header=TRUE, stringsAsFactors=FALSE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R')
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
inf.no.oil <- x.oil*b.oil
inf.no.oil <- x.oil %*% b.oil
plot(inf.no.oil)
inf.no.oil <- y.oil - x.oil %*% b.oil
plot(inf.no.oil)
inf.fitted = x.oil %*% b.oil
plot(inf.fitted)
plot(inf.no.oil)
plot(wd.oil)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
View(median.unbiased.estimator.stage2)
View(unpack.parameters.stage2)
View(rstar.stage2)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.estimation.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.estimation.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.wd.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
View(rstar.stage3)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
theta
theta
theta
Q
rm(list=ls())
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
plot(inf.s)
plot(out.s)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/unpack.parameters.stage1.R')
plot(y.data[,1])
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/unpack.parameters.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/unpack.parameters.stage1.R', echo=TRUE)
plot(y.data[,1])
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
plot(y.data[,1])
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/unpack.parameters.stage1.R')
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/unpack.parameters.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
plot(y.data[,1])
plot(x.data[,1])
plot(x.data[,2])
plot(ae.inflation)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/unpack.parameters.stage1.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
theta
theta
plot(inf.s)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
theta
plot(inf.s)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
theta
theta
plot(out.s)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
theta
pllot(inf.s)
plot(inf.s)
plot(out.s)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
theta
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
theta
theta
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
theta
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/calculate.covariance.R')
debugSource('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
View(format.output)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
write.table(output.ae, 'output/output.ae.king.csv', col.names = output.col.names, quote=FALSE, row.names=FALSE, sep = ',', na = '')
#------------------------------------------------------------------------------#
# File:        run.hlw.ae.R
#
# Description: This the main file for HLW, which rund the HLW procedure with Adavanced Economies data
#              (1) Prepares data to be used in estimation
#              (2) Runs the three-stage HLW estimation for each economy
#              (3) Saves output.
#------------------------------------------------------------------------------#
rm(list=ls())
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.kingshort.R', echo=TRUE)
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.kingshort.R', echo=TRUE)
clear all
clear
ae.log.output
source('~/Documents/Research/Research ideas/Laubach Williams for the World/HLW_Code/run.hlw.ae.kingshort.R', echo=TRUE)
source('~/Documents/Research/Research_ideas/Laubach Williams for the World/BPEA/BPEA Final Draft/Prog/HLW_Code/run.hlw.ae.R', echo=TRUE)
cd
pwd
'pwd'
disp pwd
getwd
getwd()
source('~/Documents/Research/Research_ideas/Laubach Williams for the World/BPEA/BPEA Final Draft/Prog/HLW/run.hlw.ae.R', echo=TRUE)
source('~/Documents/Research/Research_ideas/Laubach Williams for the World/BPEA/BPEA Final Draft/Prog/HLW/run.hlw.ae.R', echo=TRUE)
