# load data for Germany
data("gap")
data("indicator")
country <- "Germany"
tsList <- amecoData2input(gap[[country]], alpha = 0.65)
# compute cubs indicator
namesCubs <- c("indu", "serv", "buil")
namesVACubs <- paste0("va", namesCubs)
tscubs <- cubs(
tsCU = gap[[country]][, namesCubs],
tsVA = gap[[country]][, namesVACubs]
)
tsList <- c(tsList, tscubs)
# define tfp model
model <- TFPmodel(tsl = tsList, cycle = "RAR2", cubsErrorARMA = c(1,0))
Run the code above in your browser using DataLab