# Não se esqueça de instalar o pacote antes de executar esse roteiro,
# o que pode ser feito pelo seguinte comando:
## install.packages("ElastH")

# Para maiores informações sobre o funcionamento deste pacote leia o
# tutorial disponibilizado no sítio da SPE, e que também pode ser acessado
# utilizando o seguinte comando:
## vignette("ElastH")

library(ElastH) # Para chamar a biblioteca de funções de estimação das elasticidades

# Para executar os comandos a seguir é necessário definir as variáveis de input,
# a saber:
#    arquivo
#    Receitas
#    Hpib
#    Hpet
#    tce
#    comeco
#    fim

# arquivo deve conter o caminho e nome do arquivo onde será salvo os resultados.
# Caso nenhum caminho seja fornecido o arquivo será salvo na pasta de trabalho
# do R.
arquivo <- "elasticidades.csv"

# Receitas deve conter todos os 11 grupos de receitas, identificados
# com os nomes corretos: "TRT", "TFP", "TRC", "TI", "TM", "TGC", 
# "ROY", "PE", "TRAN", "ICMS", "ISS".
# A classe dessa variável deve ser c("mts", "ts", "matrix")
Receitas  <- Exemplo$receitas

# Hpib deve ser conter a série temporal de ln(Y/Y*), em que Y e Y* são
# respectivamente o PIB e o PIB tendencial, ou por aproximação o hiato
# do produto. 
# A classe dessa variável deve ser "ts"
Hpib      <- Exemplo$Hpib

# Hpet deve ser conter a série temporal de ln(Pet/Pet*), em que Pet e Pet* são
# respectivamente o preço do petróleo e o preço do petróleo tendencial, ou por
# aproximação o hiato do produto. 
# A classe dessa variável deve ser "ts"
Hpet      <- Exemplo$Hpet

# tce deve conter a série temporal de taxa de câmbio efetiva, que será utilizada
# como variável de controle na estimação da elasticidade das receitas de 
# importação. Esta variável será ignorada se for definida como NULL.
tce       <- NULL

# fim deve conter o período final de estimação, que ordinariamente é o último
# dado disponível para todos os grupos de receita. Por padrão, esta variável é
# definida como o último trimestre do último ano em qua há dados disponíveis,
# no caso:
fim       <- c(2015,4)

# comeco deve conter o período inicial de estimação, e podem variar para cada
# grupo de receita. Por padrão, esta variável é definida para os valores
# utilizados pela metodologia de SPE.
comeco <- list(
  trt  = 2000,
  tfp  = 1997,
  trc  = 1997,
  ti   = 2001,
  tm   = 2005,
  tgc  = 2005,
  roy  = 2000,
  pe   = 2000,
  tran = 2002,
  icms = 1997,
  iss  = 2002)


resultado <- list(
  trt  = decompor.todos(Receitas[, "TRT"] , X = Hpib, comeco = comeco$trt , fim = fim),
  trt1 = decompor.todos(Receitas[, "TRT"] , X = lag(Hpib,-1), comeco = comeco$trt , fim = fim),
  trt2 = decompor.todos(Receitas[, "TRT"] , X = cbind(Hpib, lag(Hpib, -1)), comeco = comeco$trt , fim = fim),
  
  tfp  = decompor.todos(Receitas[, "TFP"] , X = Hpib, comeco = comeco$tfp , fim = fim),
  tfp1 = decompor.todos(Receitas[, "TFP"] , X = lag(Hpib,-1), comeco = comeco$tfp , fim = fim),
  tfp2 = decompor.todos(Receitas[, "TFP"] , X = cbind(Hpib, lag(Hpib, -1)), comeco = comeco$tfp , fim = fim),

  trc  = decompor.todos(Receitas[, "TRC"] , X = Hpib, comeco = comeco$trc , fim = fim),
  trc1 = decompor.todos(Receitas[, "TRC"] , X = lag(Hpib,-1), comeco = comeco$trc , fim = fim),
  trc2 = decompor.todos(Receitas[, "TRC"] , X = cbind(Hpib, lag(Hpib, -1)), comeco = comeco$trc , fim = fim),

  ti   = decompor.todos(Receitas[, "TI"]  , X = Hpib, comeco = comeco$ti  , fim = fim),
  ti1  = decompor.todos(Receitas[, "TI"]  , X = lag(Hpib,-1), comeco = comeco$ti  , fim = fim),
  ti2  = decompor.todos(Receitas[, "TI"]  , X = cbind(Hpib, lag(Hpib, -1)), comeco = comeco$ti  , fim = fim),

  tm   = decompor.todos(Receitas[, "TM"]  , X = cbind(Hpib, log(tce)), comeco = comeco$tm, fim = fim),
  tm1  = decompor.todos(Receitas[, "TM"]  , X = cbind(lag(Hpib,-1), log(tce)), comeco = comeco$tm, fim = fim),
  tm2  = decompor.todos(Receitas[, "TM"]  , X = cbind(Hpib, lag(Hpib, -1), log(tce)), comeco = comeco$tm, fim = fim),

  tgc  = decompor.todos(Receitas[, "TGC"] , X = Hpib, comeco = comeco$tgc , fim = fim),
  tgc1 = decompor.todos(Receitas[, "TGC"] , X = lag(Hpib, -1), comeco = comeco$tgc , fim = fim),
  tgc2 = decompor.todos(Receitas[, "TGC"] , X = cbind(Hpib, lag(Hpib, -1)), comeco = comeco$tgc , fim = fim),

  roy  = decompor.todos(Receitas[, "ROY"] , X = Hpet, comeco = comeco$roy , fim = fim),
  roy1 = decompor.todos(Receitas[, "ROY"] , X = lag(Hpet, -1), comeco = comeco$roy , fim = fim),
  roy2 = decompor.todos(Receitas[, "ROY"] , X = cbind(Hpet, lag(Hpet, -1)), comeco = comeco$roy , fim = fim),
  
  pe   = decompor.todos(Receitas[, "PE"]  , X = Hpet, comeco = comeco$pe  , fim = fim),
  pe1  = decompor.todos(Receitas[, "PE"]  , X = lag(Hpet, -1), comeco = comeco$pe  , fim = fim),
  pe2  = decompor.todos(Receitas[, "PE"]  , X = cbind(Hpet, lag(Hpet, -1)), comeco = comeco$pe  , fim = fim),

  tran = decompor.todos(Receitas[, "TRAN"], X = Hpib, comeco = comeco$tran, fim = fim),
  tran1= decompor.todos(Receitas[, "TRAN"], X = lag(Hpib, -1), comeco = comeco$tran, fim = fim),
  tran2= decompor.todos(Receitas[, "TRAN"], X = cbind(Hpib, lag(Hpib, -1)), comeco = comeco$tran, fim = fim),

  icms = decompor.todos(Receitas[, "ICMS"], X = Hpib, comeco = comeco$icms, fim = fim),
  icms1= decompor.todos(Receitas[, "ICMS"], X = lag(Hpib,-1), comeco = comeco$icms, fim = fim),
  icms2= decompor.todos(Receitas[, "ICMS"], X = cbind(Hpib, lag(Hpib, -1)), comeco = comeco$icms, fim = fim),

  iss  = decompor.todos(Receitas[, "ISS"] , X = Hpib, comeco = comeco$iss , fim = fim),
  iss1 = decompor.todos(Receitas[, "ISS"] , X = lag(Hpib, -1), comeco = comeco$iss , fim = fim),
  iss2 = decompor.todos(Receitas[, "ISS"] , X = cbind(Hpib, lag(Hpib, -1)), comeco = comeco$iss , fim = fim)
)

ret <- exportar(resultado)
write.csv2(ret, file=arquivo)