Skip to content

Instantly share code, notes, and snippets.

@sientifiko
sientifiko / multreg.R
Created September 3, 2024 02:17
Codigo para mostrar interpretación de betas en el contexto de regresión multiple
library(foreign)
dat <- read.dta("http://fmwww.bc.edu/ec-p/data/wooldridge/gpa1.dta")
coef(lm(colGPA ~ ACT + hsGPA, data = dat))
coef(lm(colGPA ~ hsGPA, data = dat))
@sientifiko
sientifiko / stargazertableout.R
Created August 29, 2024 03:02
Función que genera un output similar al de stargazer, pero en formato tabulado para que puedan exportarlo
generate_regression_table <- function(models, se_type = "HC3") {
require(sandwich)
require(lmtest)
# Initialize lists to store results
results_list <- list()
# Loop over each model
for (i in seq_along(models)) {
model <- models[[i]]
@sientifiko
sientifiko / guia_desigualdad1.R
Last active December 19, 2023 08:53
Script para la nota 1 sobre guía de desigualdad en medium
library(tidyverse)
theme_set(theme_bw(base_size = 21))
options(scipen = 999)
dat <- read.csv("dataregionesChile.csv")
# filtrar el año 2020
y2020 <- dat %>% filter(year == 2020)
# generar histograma
@sientifiko
sientifiko / script maddison.R
Created December 11, 2023 20:06
Generar gráficas sobre PIB argentina
library(tidyverse)
library(directlabels)
theme_set(theme_bw(base_size = 21))
options(scipen = 999)
dat <- readxl::read_excel("mpd2020.xlsx", sheet = 3) %>%
filter(country %in% c("Argentina", "Chile", "Brazil","Sweden", "United States",
"Norway", "Uruguay", "United Kingdom", "Panama"),
year >= 1870)
library(tidyverse)
library(countrycode)
library(gganimate)
library(gifski)
theme_set(theme_classic())
tiempo <- read.csv("annual-working-hours-per-worker.csv")
colnames(tiempo)[4] <- "tiempo"
@sientifiko
sientifiko / Text stat prog gob 2021.R
Created November 5, 2021 01:30
Código para replicar la estadística de texto de los programas de gobierno
library(XML)
library(RCurl)
library(tm)
library(wordcloud2)
library(stm)
library(pdftools)
library(tidyverse)
library(patchwork)
library(ggwordcloud)
@sientifiko
sientifiko / scrapservel1.R
Last active November 26, 2021 20:34
Codigo para scrapear el servel sacando votación por comuna
library(tidyverse)
library(rvest)
library(RSelenium)
# funciones pa limpiar
limpieza <- function(objtabla){
objtabla <- objtabla[-which(is.na(objtabla[,1])), ]
objtabla <- objtabla[c(1:(nrow(objtabla)-5)),]
@sientifiko
sientifiko / ineq inevitable.R
Last active June 21, 2021 21:17
Script permite modelar desigualdad "natural"
library(tidyverse)
library(ineq)
options(scipen = 999)
theme_set(theme_classic())
base <- rep(1000, 1000)
lista <- list()
library(tidyverse)
library(patchwork)
theme_set(theme_classic())
paises <- c("Argentina", "Bolivia", "Brasil", "Chile", "Colombia",
"Costa Rica", "Cuba", "Ecuador", "El Salvador", "Guatemala", "Honduras",
"México", "Nicaragua", "Panamá", "Paraguay", "Perú", "Rep Dominicana",
"Uruguay", "Venezuela")