baseDir = "/Users/andreamanera/Dropbox (MIT)/Research/Brookings20"
# data sources
rawDir = paste0(baseDir,"/DepreciationData/R/raw")
# source for Tax BEA data
taxDir = paste0(rawDir,"/TaxSeries")
# out figures
figureDir = rawDir = paste0(baseDir,"/Docs/Figures")
# out other
outDir = paste0(baseDir,"/DepreciationData/R/out")
# McDaniel dir: where dataset McDanielMine is
McDanielDir = paste0(baseDir,"/McDanielComp/out")
#IMPORT AND MERGE DATASETS----------
tableGDP = setDT(read_csv(paste0(taxDir,"/Table_1_1_5_GDP.csv")))
tableGDI = setDT(read_csv(paste0(taxDir,"/Table_1_10_GDI.csv")))
tableVADD = setDT(read_csv(paste0(taxDir,"/Table_1_14_Value_Add.csv")))
tableGCURR = setDT(read_csv(paste0(taxDir,"/Table_3_1_Govt_Current_BEA.csv")))
# Select interesting variables and rename
tableGDP = tableGDP[Line == 1 | Line == 2 | Line == 8, ]
tableGDP$Item = c("GDP", "C", "I")
tableGDI = tableGDI[Line == 2 | Line == 10 | Line == 20 | Line == 22 | Line == 23 , ]
tableGDI$Item = c("W", "NOSPCU", "NOSGOV", "CFCPCU", "CFCGOV")
tableVADD = tableVADD[Line == 2 | Line == 8, ]
tableVADD$Item = c("CFCCORP", "NOSCORP")
tableGCURR = tableGCURR[(Line >= 3 & Line <= 5 )| Line == 8 | Line == 30 ,]
tableGCURR$Item = c("HHT", "TPI", "CT", "SS", "Sub")
tableBEA = rbind(tableGDP, tableGDI, tableVADD, tableGCURR)
tableBEA = tableBEA[,Line:= NULL]
# reshape to long
tableBEA = melt(tableBEA,
id.vars = c("Item"),
measure.vars = colnames(tableBEA)[colnames(tableBEA) !="Item"],
variable.name = "Year",
value.name = "Value")
tableBEA$Value = as.numeric(tableBEA$Value)
tableBEA$Year = as.numeric(as.character(tableBEA$Year))
tableBEA = dcast(tableBEA, Year~ Item, value.var="Value")
# Import McDaniel series on old vintage of data
tableMcDaniel<- setDT(read_csv(paste0(McDanielDir,"/McDanielMine.csv")))
# Get Business Share of TPI
mu = tableMcDaniel[1,]$mean_mu
# Get Share of net TPI that fall on consumption only
lambda = tableMcDaniel[1,]$mean_lambda
#COMPUTE AVERAGE TAXES USING MCDANIEL FORMULAS ON BEA DATA
tableBEA = tableBEA[, CONS_SHARE:= C/(C+I)]
tableBEA = tableBEA[, T_INC:=  HHT/(GDP - (TPI - Sub))]
# GET GROSS OPERATING SURPLUSES: Net + Consumption of fixed capital
tableBEA = tableBEA[, OSPCU:= NOSPCU + CFCPCU]
tableBEA = tableBEA[, OSCORP:= NOSCORP + CFCCORP]
tableBEA = tableBEA[, OSGOV:= NOSGOV + CFCGOV]
# OSPUE
tableBEA = tableBEA[, OSPUE:= OSPCU - OSCORP]
tableBEA = tableBEA[, LABOR_SHARE:= W/(GDP - (TPI - Sub) - OSPUE)]
#HH direct taxes on labor
tableBEA = tableBEA[, HHT_L:= T_INC * LABOR_SHARE * (GDP - (TPI - Sub))]
# Avreage tax on labor income
tableBEA = tableBEA[, T_H:= (SS + HHT_L) /(LABOR_SHARE * (GDP - (TPI - Sub)))]
# Revenue from taxes that fall exclusively on consumption
tableBEA = tableBEA[, TPI_C:= (lambda + (1 - lambda) * CONS_SHARE) * ((1 - mu) * TPI - Sub)]
# Average consumption tax
tableBEA = tableBEA[, T_C:= TPI_C/(C - TPI_C) ]
# Revenue from taxes on investment expenditures
tableBEA = tableBEA[, TPI_X:= (1 - mu) * TPI - Sub -TPI_C]
# Average investment tax
tableBEA = tableBEA[, TPI_X:= TPI_X/(I - TPI_X)]
#Household direct taxes on capital
tableBEA = tableBEA[, HHT_C_ALT:= T_INC * ((1 - LABOR_SHARE)*
(GDP - (TPI - Sub)) - OSGOV - CT)  ]
tableBEA = tableBEA[, HHT_C:= T_INC * ((1 - LABOR_SHARE)*
OSPUE + OSCORP - CT)  ]
#Average Capital Tax
tableBEA = tableBEA[, T_K:= (HHT_C + CT + mu * TPI)/((1 - LABOR_SHARE)*
OSPUE + OSCORP - CT)   ]
tableBEA = tableBEA[, T_K_ALT:= (HHT_C_ALT + CT + mu * TPI)/((1 - LABOR_SHARE)*
(GDP - (TPI - Sub)) - OSGOV - CT)   ]
tableBEA$Year = as.numeric(tableBEA$Year)
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y=T_K,
colour="K tax, BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=t_k,
colour="K tax, OECD"))+
geom_line(data=tableBEA,aes(x=Year,y=T_K_ALT,
colour="K tax Alt, BEA"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Capital taxes")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y=CT,
colour="BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=ct/1000 ,
colour="OECD"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Oscorp")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y=(1 - LABOR_SHARE)*
(GDP - (TPI - Sub)) - OSGOV,
colour="BEA, OECD"))+
geom_line(data=tableBEA,aes(x=Year,y=(1 - LABOR_SHARE)*
OSPUE + OSCORP ,
colour="BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=  ( (1- labor_share)*ospue +(os - ospue- osgov))/1000 ,
colour="OECD, BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=  ( (1- labor_share)*(gdp -(tpi - sub)) -osgov)/1000 ,
colour="OECD"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Oscorp")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y=(1 - LABOR_SHARE)*
(GDP - (TPI - Sub)) - OSGOV,
colour="BEA, OECD"))+
geom_line(data=tableBEA,aes(x=Year,y=(1 - LABOR_SHARE)*
OSPUE + OSCORP ,
colour="BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=  ( (1- labor_share)*ospue +(os - ospue- osgov))/1000 ,
colour="OECD, BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=  ( (1- labor_share)*(gdp -(tpi - sub)) -osgov)/1000 ,
colour="OECD"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Oscorp")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y= OSGOV,
colour="BEA, "))+
geom_line(data=tableMcDaniel,aes(x=year,y=  osgov/1000 ,
colour="OECD "))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Oscorp")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y=(1 - LABOR_SHARE)*
(GDP - (TPI - Sub)) - OSGOV,
colour="BEA, OECD"))+
geom_line(data=tableBEA,aes(x=Year,y=(1 - LABOR_SHARE)*
OSPUE + OSCORP ,
colour="BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=  ( (1- labor_share)*ospue +(os - ospue- osgov))/1000 ,
colour="OECD, BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=  ( (1- labor_share)*(gdp -(tpi - sub)) -osgov)/1000 ,
colour="OECD"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Oscorp")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y=T_K,
colour="K tax, BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=t_k,
colour="K tax, OECD"))+
geom_line(data=tableBEA,aes(x=Year,y=T_K_ALT,
colour="K tax Alt, BEA"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Capital taxes")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y=(1 - LABOR_SHARE)*
(GDP - (TPI - Sub)) - OSGOV,
colour="BEA, OECD"))+
geom_line(data=tableBEA,aes(x=Year,y=(1 - LABOR_SHARE)*
OSPUE + OSCORP ,
colour="BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=  ( (1- labor_share)*ospue +(os - ospue- osgov))/1000 ,
colour="OECD, BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=  ( (1- labor_share)*(gdp -(tpi - sub)) -osgov)/1000 ,
colour="OECD"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Oscorp")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
mergedata = tableBEA[tableMcDaniel, on=.(Year = year)]
ggplot(mergedata)+
geom_line(aes(x=Year,y=(HHT_C_ALT + CT + mu * TPI)/(hht_c + ct + mu * tpi),
colour="num, BEA/McD"))+
geom_line( aes(x=Year,y=((1 - LABOR_SHARE)*
(GDP - (TPI - Sub)) - OSGOV - CT)/
((1 - labor_share)*
(gdp - (tpi - sub)) - osgov - ct),
colour="denom, BEA/McD"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="K tax comparison")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
ggplot(mergedata)+
geom_line(aes(x=Year,y=(HHT_C_ALT + CT + mu * TPI)/(hht_c + ct + mu * tpi)*1000,
colour="num, BEA/McD"))+
geom_line( aes(x=Year,y=((1 - LABOR_SHARE)*
(GDP - (TPI - Sub)) - OSGOV - CT)/
((1 - labor_share)*
(gdp - (tpi - sub)) - osgov - ct)*1000,
colour="denom, BEA/McD"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="K tax comparison")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y=T_K,
colour="K tax, BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=t_k,
colour="K tax, OECD"))+
geom_line(data=tableBEA,aes(x=Year,y=T_K_ALT,
colour="K tax Alt, BEA"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Capital taxes")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
ggplot(mergedata)+
geom_line(aes(x=Year,y=(HHT_C_ALT + CT + mu * TPI)/(hht_c + ct + mu * tpi)*1000,
colour="num, BEA/McD"))+
geom_line( aes(x=Year,y=((1 - LABOR_SHARE)*
(GDP - (TPI - Sub)) - OSGOV - CT) /
( (1- labor_share)*ospue +(os - ospue- osgov)) *1000,
colour="denom, BEA/McD"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="K tax comparison")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
ggplot(mergedata)+
geom_line(aes(x=Year,y=(HHT_C_ALT + CT + mu * TPI)/(hht_c + ct + mu * tpi)*1000,
colour="num, BEA/McD"))+
geom_line( aes(x=Year,y=((1 - LABOR_SHARE)*
(GDP - (TPI - Sub)) - OSGOV - CT) /
( (1- labor_share)*(gdp -(tpi - sub)) -osgov) *1000,
colour="denom, BEA/McD"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="K tax comparison")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
tableBEA = tableBEA[, T_X:= TPI_X/(I - TPI_X)]
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y=T_C,
colour="C tax, BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=t_c,
colour="C tax, OECD"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Consumption taxes")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
#LABOR INCOME TAXES
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y=T_H,
colour="L tax, BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=t_h,
colour="L tax, OECD"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Labor taxes")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
#INVESTMENT TAXES
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y=T_X,
colour="I tax, BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=t_x,
colour="I tax, OECD"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Investment taxes")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
library(readr)
library(dplyr)
library(reshape2)
library(tibble)
library(data.table)
library(ggplot2)
# Change this to run in another directory with same folder structure
baseDir = "/Users/andreamanera/Dropbox (MIT)/Research/Brookings20"
# data sources
rawDir = paste0(baseDir,"/DepreciationData/R/raw")
# source for Tax BEA data
taxDir = paste0(rawDir,"/TaxSeries")
# out figures
figureDir = rawDir = paste0(baseDir,"/Docs/Figures")
# out other
outDir = paste0(baseDir,"/DepreciationData/R/out")
# McDaniel dir: where dataset McDanielMine is
McDanielDir = paste0(baseDir,"/McDanielComp/out")
#IMPORT AND MERGE DATASETS----------
tableGDP = setDT(read_csv(paste0(taxDir,"/Table_1_1_5_GDP.csv")))
tableGDI = setDT(read_csv(paste0(taxDir,"/Table_1_10_GDI.csv")))
tableVADD = setDT(read_csv(paste0(taxDir,"/Table_1_14_Value_Add.csv")))
tableGCURR = setDT(read_csv(paste0(taxDir,"/Table_3_1_Govt_Current_BEA.csv")))
# Select interesting variables and rename
tableGDP = tableGDP[Line == 1 | Line == 2 | Line == 8, ]
tableGDP$Item = c("GDP", "C", "I")
tableGDI = tableGDI[Line == 2 | Line == 10 | Line == 20 | Line == 22 | Line == 23 , ]
tableGDI$Item = c("W", "NOSPCU", "NOSGOV", "CFCPCU", "CFCGOV")
tableVADD = tableVADD[Line == 2 | Line == 8, ]
tableVADD$Item = c("CFCCORP", "NOSCORP")
tableGCURR = tableGCURR[(Line >= 3 & Line <= 5 )| Line == 8 | Line == 30 ,]
tableGCURR$Item = c("HHT", "TPI", "CT", "SS", "Sub")
tableBEA = rbind(tableGDP, tableGDI, tableVADD, tableGCURR)
tableBEA = tableBEA[,Line:= NULL]
# reshape to long
tableBEA = melt(tableBEA,
id.vars = c("Item"),
measure.vars = colnames(tableBEA)[colnames(tableBEA) !="Item"],
variable.name = "Year",
value.name = "Value")
tableBEA$Value = as.numeric(tableBEA$Value)
tableBEA$Year = as.numeric(as.character(tableBEA$Year))
tableBEA = dcast(tableBEA, Year~ Item, value.var="Value")
# Import McDaniel series on old vintage of data
tableMcDaniel<- setDT(read_csv(paste0(McDanielDir,"/McDanielMine.csv")))
# Get Business Share of TPI
mu = tableMcDaniel[1,]$mean_mu
# Get Share of net TPI that fall on consumption only
lambda = tableMcDaniel[1,]$mean_lambda
#COMPUTE AVERAGE TAXES USING MCDANIEL FORMULAS ON BEA DATA
tableBEA = tableBEA[, CONS_SHARE:= C/(C+I)]
tableBEA = tableBEA[, T_INC:=  HHT/(GDP - (TPI - Sub))]
# GET GROSS OPERATING SURPLUSES: Net + Consumption of fixed capital
tableBEA = tableBEA[, OSPCU:= NOSPCU + CFCPCU]
tableBEA = tableBEA[, OSCORP:= NOSCORP + CFCCORP]
tableBEA = tableBEA[, OSGOV:= NOSGOV + CFCGOV]
# OSPUE
tableBEA = tableBEA[, OSPUE:= OSPCU - OSCORP]
tableBEA = tableBEA[, LABOR_SHARE:= W/(GDP - (TPI - Sub) - OSPUE)]
#HH direct taxes on labor
tableBEA = tableBEA[, HHT_L:= T_INC * LABOR_SHARE * (GDP - (TPI - Sub))]
# Avreage tax on labor income
tableBEA = tableBEA[, T_H:= (SS + HHT_L) /(LABOR_SHARE * (GDP - (TPI - Sub)))]
# Revenue from taxes that fall exclusively on consumption
tableBEA = tableBEA[, TPI_C:= (lambda + (1 - lambda) * CONS_SHARE) * ((1 - mu) * TPI - Sub)]
# Average consumption tax
tableBEA = tableBEA[, T_C:= TPI_C/(C - TPI_C) ]
# Revenue from taxes on investment expenditures
tableBEA = tableBEA[, TPI_X:= (1 - mu) * TPI - Sub -TPI_C]
# Average investment tax
tableBEA = tableBEA[, T_X:= TPI_X/(I - TPI_X)]
#Household direct taxes on capital
tableBEA = tableBEA[, HHT_C:= T_INC * ((1 - LABOR_SHARE)*
OSPUE + OSCORP - CT)  ]
tableBEA = tableBEA[, HHT_C_ALT:= T_INC * ((1 - LABOR_SHARE)*
(GDP - (TPI - Sub)) - OSGOV - CT)  ]
#MAIN GRAPHS-----------
#Average Capital Tax
tableBEA = tableBEA[, T_K:= (HHT_C + CT + mu * TPI)/((1 - LABOR_SHARE)*
OSPUE + OSCORP - CT)   ]
tableBEA = tableBEA[, T_K_ALT:= (HHT_C_ALT + CT + mu * TPI)/((1 - LABOR_SHARE)*
(GDP - (TPI - Sub)) - OSGOV - CT)   ]
tableBEA$Year = as.numeric(tableBEA$Year)
# CAPITAL TAXES
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y=T_K,
colour="K tax, BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=t_k,
colour="K tax, OECD"))+
geom_line(data=tableBEA,aes(x=Year,y=T_K_ALT,
colour="K tax Alt, BEA"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Capital taxes")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
#CONSUMPTION TAXES
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y=T_C,
colour="C tax, BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=t_c,
colour="C tax, OECD"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Consumption taxes")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
#LABOR INCOME TAXES
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y=T_H,
colour="L tax, BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=t_h,
colour="L tax, OECD"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Labor taxes")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
#INVESTMENT TAXES
ggplot()+
geom_line(data=tableBEA,aes(x=Year,y=T_X,
colour="I tax, BEA"))+
geom_line(data=tableMcDaniel,aes(x=year,y=t_x,
colour="I tax, OECD"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Investment taxes")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
tableAllow = merge(read_csv(paste0(taxDir,"/CCAll_Corp.csv")), read_csv(paste0(taxDir,"/CCAll_NonCorp.csv")))
View(tableAllow)
tableAllow = setDT(
merge(read_csv(paste0(taxDir,"/CCAll_Corp.csv")),
read_csv(paste0(taxDir,"/CCAll_NonCorp.csv"))))
tableBEA = tableBEA[tableAllow, on.(Year = DATE)]
tableBEA = tableBEA[tableAllow, on=.(Year = DATE)]
View(tableBEA)
tableBEA = tableBEA[, T_K_taxbase_nopp:= (T_INC * ((1 - LABOR_SHARE)*
NOSPUE + NOSCORP - CT) + CT) /((1 - LABOR_SHARE)*
OSPUE + OSCORP - CT)  ]
tableBEA = tableBEA[, T_K_taxbase_pp:= (T_INC * ((1 - LABOR_SHARE)*
NOSPUE + NOSCORP - CT) + CT + mu * TPI) /((1 - LABOR_SHARE)*
OSPUE + OSCORP - CT)  ]
tableBEA = tableBEA[, T_K_taxbase_nopp:= (T_INC * ((1 - LABOR_SHARE)*
NOSPUE_IRS + NOSCORP_IRS - CT) + CT) /((1 - LABOR_SHARE)*
NOSPUE_IRS + NOSCORP_IRS  - CT)  ]
tableBEA = tableBEA[, T_K_taxbase_pp:= (T_INC * ((1 - LABOR_SHARE)*
NOSPUE_IRS + NOSCORP_IRS - CT) + CT + mu * TPI) /((1 - LABOR_SHARE)*
NOSPUE_IRS + NOSCORP_IRS  - CT)  ]
tableBEA = tableBEA[, NOSPUE_IRS:= OSPUE - NonCorporateAllowance  ]
tableBEA = tableBEA[, NOSCORP_IRS:= OSCORP - CorporateAllowance  ]
tableBEA = tableBEA[, HHT_NC:= T_INC * ((1 - LABOR_SHARE)*
NOSPUE + NOSCORP - CT)   ]
tableBEA = tableBEA[, HHT_NC:= T_INC * ((1 - LABOR_SHARE)*
NOSPUE_IRS + NOSCORP_IRS - CT)   ]
tableBEA = tableBEA[, T_K_taxbase_nopp:= (T_INC * ((1 - LABOR_SHARE)*
NOSPUE_IRS + NOSCORP_IRS - CT) + CT)
/((1 - LABOR_SHARE)*NOSPUE_IRS + NOSCORP_IRS  - CT)  ]
tableBEA = tableBEA[, T_K_taxbase_pp:= (T_INC * ((1 - LABOR_SHARE)* NOSPUE_IRS + NOSCORP_IRS - CT)
+ CT + mu * TPI) /((1 - LABOR_SHARE)*NOSPUE_IRS + NOSCORP_IRS  - CT)]
ggplot(tableBEA)+
geom_line(aes(x=Year,y=T_K,
colour="K tax "))+
geom_line(aes(x=Year,y=T_K_taxbase_nopp,
colour="K tax base, no prop. taxes"))+
geom_line(aes(x=Year,y=T_K_taxbase_pp,
colour="K tax base,  prop. taxes"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Capital taxes")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
tableBEA = tableBEA[, T_K_corp_net:= CT/NOSCORP]
ggplot(tableBEA)+
geom_line(aes(x=Year,y=T_K,
colour="K tax "))+
geom_line(aes(x=Year,y=T_K_taxbase_nopp,
colour="K tax base, no prop. taxes"))+
geom_line(aes(x=Year,y=T_K_taxbase_pp,
colour="K tax base,  prop. taxes"))+
geom_line(aes(x=Year,y=T_K_corp,
colour="K tax base,  corp. only"))+
geom_line(aes(x=Year,y=T_K_corp_net,
colour="K tax base,  corp. only"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Capital taxes")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
tableBEA = tableBEA[, T_K_corp:= CT/OSCORP]
tableBEA = tableBEA[, T_K_corp_net:= CT/NOSCORP]
ggplot(tableBEA)+
geom_line(aes(x=Year,y=T_K,
colour="K tax "))+
geom_line(aes(x=Year,y=T_K_taxbase_nopp,
colour="K tax base, no prop. taxes"))+
geom_line(aes(x=Year,y=T_K_taxbase_pp,
colour="K tax base,  prop. taxes"))+
geom_line(aes(x=Year,y=T_K_corp,
colour="K tax base,  corp. only"))+
geom_line(aes(x=Year,y=T_K_corp_net,
colour="K tax base,  corp. only"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Capital taxes")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
ggplot(tableBEA)+
geom_line(aes(x=Year,y=T_K,
colour="K tax "))+
geom_line(aes(x=Year,y=T_K_taxbase_nopp,
colour="K tax base, no prop. taxes"))+
geom_line(aes(x=Year,y=T_K_taxbase_pp,
colour="K tax base,  prop. taxes"))+
geom_line(aes(x=Year,y=T_K_corp,
colour="K tax base,  corp. only"))+
geom_line(aes(x=Year,y=T_K_corp_net,
colour="K tax base,  corp. only net"))+
scale_colour_discrete("")+
labs(x="Year",y="$bn",title="Capital taxes")+
theme(plot.title = element_text(hjust = 0.5),legend.position = "bottom")
mu
tableInv <- setDT(read_csv(paste0(invDir,"/Table_2_7_INV.csv")))
tableInv_Q <- setDT(read_csv(paste0(invDir,"/Table_2_8_Chain_Inv_q.csv")))
invDir = paste0(rawDir,"/InvPrices")
tableInv <- setDT(read_csv(paste0(invDir,"/Table_2_7_INV.csv")))
tableInv_Q <- setDT(read_csv(paste0(invDir,"/Table_2_8_Chain_Inv_q.csv")))
# data sources
rawDir = paste0(baseDir,"/DepreciationData/R/raw")
invDir = paste0(rawDir,"/InvPrices")
invDir = paste0(rawDir,"/InvPrices")
tableInv <- setDT(read_csv(paste0(invDir,"/Table_2_7_INV.csv")))
tableInv_Q <- setDT(read_csv(paste0(invDir,"/Table_2_8_Chain_Inv_q.csv")))
tableInv = cbind(tableInv[,1:2], tableInv[,3:end]/tableInv_Q )
tableInv = cbind(tableInv[,1:2], tableInv[,3:ncol(tableInv)]/tableInv_Q[,3:ncol(tableInv)] )
View(tableInv)
tableInv = melt(tableInv,
measure.vars = colnames(tableInv)[colnames(tableInv) !="Item" & colnames(tableInv) !="Line"],
variable.name = "Year",
value.name = "Value")
tableInv = cbind(tableInv[,1:2], tableInv[,3:ncol(tableInv)]/tableInv_Q[,3:ncol(tableInv)] )
# Reshape to Long
tableInv = melt(tableInv,
measure.vars = colnames(tableInv)[colnames(tableInv) !="Type" & colnames(tableInv) !="Line"],
variable.name = "Year",
value.name = "Value")
View(tableInv)
tableInv = melt(tableInv,
id.vars = c("Type", "Line"),
measure.vars = colnames(tableInv)[colnames(tableInv) !="Type" & colnames(tableInv) !="Line"],
variable.name = "Year",
value.name = "Value")
tableInv = tableInv[, Line = NULL]
tableInv = tableInv[, .Line = NULL]
tableInv = tableInv[, .Line := NULL]
tableInv = tableInv[, Line:= NULL]
tableInv = melt(tableInv,
id.vars = c("Type"),
measure.vars = colnames(tableInv)[colnames(tableInv) !="Type" & colnames(tableInv) !="Line"],
variable.name = "Year",
value.name = "Value")
tableInv = cbind(tableInv[,1:2], tableInv[,3:ncol(tableInv)]/tableInv_Q[,3:ncol(tableInv)] )
tableInv = tableInv[, Line:= NULL]
# Reshape to Long
tableInv = melt(tableInv,
id.vars = c("Type"),
measure.vars = colnames(tableInv)[colnames(tableInv) !="Type" & colnames(tableInv) !="Line"],
variable.name = "Year",
value.name = "Value")
tableInv <- setDT(read_csv(paste0(invDir,"/Table_2_7_INV.csv")))
tableInv_Q <- setDT(read_csv(paste0(invDir,"/Table_2_8_Chain_Inv_q.csv")))
tableInv = cbind(tableInv[,1:2], tableInv[,3:ncol(tableInv)]/tableInv_Q[,3:ncol(tableInv)] )
tableInv = tableInv[, Line:= NULL]
tableInv = melt(tableInv,
id.vars = c("Type"),
measure.vars = colnames(tableInv)[colnames(tableInv) !="Type" & colnames(tableInv) !="Line"],
variable.name = "Year",
value.name = "Value")
View(tableInv)
