sientifiko sientifiko
sientifiko
/ visual stm.R
Created
October 18, 2024 23:09
Script para generar una visualización de las distribuciones de efectos condicionales a posterior a partir de un objeto de tipo estimateEffect (básicamente el que creas cuando haces la regresión)
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| library(stm) | |
| library(tidytext) | |
| library(ggridges) | |
| n_topicos = 20 | |
| n_nuestras = 100 | |
| # iteración sobre tópicos | |
| lapply(1:n_topicos, function(i){ |
sientifiko
/ pib_imacec.R
Created
November 5, 2024 23:57
Script para obtener la proyección de crecimiento tendencial (asumiendo linealidad) a partir del PIB y el IMACEC
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| # "pib trim.xlsx" es el PIB trimestral a precios encadenados de este enlace: | |
| #https://si3.bcentral.cl/Siete/ES/Siete/Cuadro/CAP_CCNN/MN_CCNN76/CCNN2018_P0_V2/637801082315858005?cbFechaInicio=2000&cbFechaTermino=2024&cbFrecuencia=QUARTERLY&cbCalculo=NONE&cbFechaBase= | |
| pibtrim <- readxl::read_excel("pib trim.xlsx",skip = 1) | |
| colnames(pibtrim) <- c("time", "pib") | |
| pibtrim <- pibtrim %>% | |
| mutate(time2 = ymd(str_sub(time, 1, 10))) %>% |
OlderNewer