Skip to content

Instantly share code, notes, and snippets.

@andybega
andybega / chart-model.R
Created January 11, 2019 12:36
RCT-A subset machine/human relative performance by question group
library("tidyverse")
arima_qs <- readr::read_csv("/path/to/forecast_rct_a_subset.csv",
col_types = cols(user_id = col_integer())) %>%
mutate(Forecaster = factor(Forecaster))
# Add # of questions to each Data source label
arima_qs <- arima_qs %>%
group_by(Data_source) %>%
mutate(by_Data_source_n_ifps = length(unique(ifp_id))) %>%
@andybega
andybega / plot-yemen.R
Created February 5, 2019 10:47
cshapes with sf
library(cshapes)
library(sf)
states <- cshp(as.Date("2016-01-01"), useGW = TRUE)
states <- st_as_sf(states)
yemen <- states[states$CNTRY_NAME=="Yemen", ]
plot(yemen[, 1])
@andybega
andybega / line-segment-colors.R
Created April 1, 2019 14:58
Line segment colors in base plot
plot(c(1, 2), c(1, 1), col = c("blue", "red"), pch= 19, cex=4)
lines(c(1, 2), c(1, 1), col = c("blue", "red"), lwd=3)
@andybega
andybega / us48-spatial-lag.R
Last active August 7, 2019 17:09
Fraction of neighboring US states with a policy
library("splm")
data("usaww")
?usaww
# Only Alabama has this policy
policy_x <- c(1, rep(0, 47))
names(policy_x) <- rownames(usaww)
head(policy_x)
@andybega
andybega / runme.R
Last active October 2, 2019 19:58
Scrape the Wikipedia terror incidents pages/lists
#
# copy and run this code in R to setup and source `wikipedia-terrorism.R` from this gist
#
packs <- c("dplyr", "rvest", "httr", "XML", "stringr", "purrr", "xml2", "tibble", "devtools")
need <- packs[!packs %in% rownames(installed.packages())]
install.packages(need)
devtools::source_gist("https://gist.github.com/andybega/0f16b6a0de0185c157e2a280b5e05074",
filename = "wikipedia-terrorism.R")
@andybega
andybega / add-countries-to-cshapes
Last active November 24, 2020 18:27
Example of how to add countries (polygons) to cshapes
# Example of how to add missing countries/territories to cshapes using another source for
# the map data/polygon
#
library(cshapes)
library(rgdal) # need this for projection transform (spTransform)
# this returns all polygons, with start and end dates for indepdence
# we'll have to manually subset for dates below
all <- cshp()