schochastics

# Recreation of https://twitter.com/stevenstrogatz/status/1295915404574040065/photo/1
library(tidyverse)
library(Polychrome)
fun <- function(x){
  n <- c()
  i <- 2
  r <- x
  while(prod(n)!=x){
    if(!r%%i) {n=c(n,i);r=r/i;i=1}
    i <- i+1

schochastics

library(tidyverse)
library(ggforce)

d <- tibble(x = c(1, 6, 10, 23,18), xpred = c(4, 8, 11, 2,22), y = 1:5, id = factor(1:5))
n <- nrow(d)+1
time_angle <- tibble(hour=1:24,angle=0.25*60*hour)

d$start <- time_angle$angle[d$x]*pi/180
d$end <- time_angle$angle[d$xpred]*pi/180
d$r <- d$y

schochastics

lissajou <- function(a,b,delta){
  t <- seq(0,2*pi,0.01)
  x <- sin(a*t+delta)
  y <- sin(b*t)
  tibble(x,y)
}

# setup parameters ----
n <- 7
params <- tibble(a=sample(1:10,n,replace = FALSE),

schochastics

pyramid <- function(ymax=5,mcol="#CD3333",xoffset=0){
  f <- function(a,b,y){
    (y-b)/a
  }
  
  a1 <- ymax/0.65
  b1 <- 0
  c1 <- -ymax/(1-0.65)
  d1 <- -c1
  a2 <- -ymax/(1-0.65)

schochastics

library(rvest)
library(tidyverse)
library(gganimate)
library(extrafont)
loadfonts()

get_season <- function(x){
 url <- paste0("https://www.basketball-reference.com/leagues/NBA_",x,"_totals.html")
 url %>% 
  read_html() %>% 

schochastics

library(tidyverse)
library(gganimate)
ascii <- c(";", ":", "!", "?", ".", "'", "\"", "(", ")", "[", "]", "{", 
           "}", "@", "*", "/", "\\", "&", "#", "%", "`", "^", "+", "<", 
           "=", ">", "|", "~", "$", "0", "1", "1", "10", "2", "2", "3", 
           "3", "4", "4", "5", "5", "6", "6", "7", "7", "8", "8", "9", "9", 
           "a", "a", "A", "b", "b", "B", "c", "c", "C", "d", "d", "D", "e", 
           "e", "E", "f", "f", "F", "g", "g", "G", "h", "h", "H", "i", "i", 
           "I", "j", "j", "J", "k", "k", "K", "l", "l", "L", "m", "m", "M", 
           "n", "n", "N", "o", "o", "O", "p", "p", "P", "q", "q", "Q", "r", 

schochastics

schochastics

library(igraph)
library(smglr)
library(tidyverse)
library(ggraph)
g <- sample_islands(9,40,0.2,10)
g <- simplify(g)
V(g)$grp <- rep(1:9,each=40)
el <- get.edgelist(g)
E(g)$col <- case_when(
V(g)$grp[el[,1]]==1 & V(g)$grp[el[,2]]  == 1 ~ 1,

schochastics

library(rvest)
library(tidyverse)
# helper functions 

get_theo_refs_cache <- function(url){
  doc <- read_html(url)  
  caps <- doc %>% html_nodes("caption") %>% html_text()
  skippy <- (any(caps=="Hypothesis" | caps=="Hypotheses"))+0
  refs <- doc %>% 
    html_table() %>% 

schochastics

library(rvest)
library(tidyverse)
 
leagues <- "http://www.footballsquads.co.uk/archive.htm" %>% 
   read_html() %>% 
   html_nodes("a") %>% 
   html_attr("href")
 
leagues_tbl <- as_tibble(str_split(leagues,"/",simplify = T)) %>% 
   mutate(link=leagues) %>%