airalcorn2’s gists

airalcorn2 / gdp.csv

Created November 2, 2016 17:29

Comparing fiscal policies and GDP to poverty rates in different countries.

airalcorn2 / d_stat_experiment.py

Last active December 22, 2016 00:05

Kolmogorov–Smirnov statistic experiment.

	# Michael A. Alcorn ([email protected])

	import matplotlib.pyplot as plt
	import numpy as np
	import seaborn as sns

	import numpy as np
	import statsmodels.api as sm # recommended import according to the docs
	import matplotlib.pyplot as plt

airalcorn2 / difference_of_biased_estimators.py

Last active December 22, 2016 00:12

	# Michael A. Alcorn ([email protected])

	import matplotlib.pyplot as plt
	import pandas as pd
	import seaborn as sns
	import random

	fixed_rates = {"p_1": 0.3, "p_2": 0.5}

	find_rates = {"biased": {"fixed": 0.7, "non_fixed": 0.9},

airalcorn2 / pokemon_filters.py

Last active March 21, 2017 00:20

Visualize filters from: https://jgeekstudies.org/2017/03/12/who-is-that-neural-network/.

	kernels = model.layers[0].get_weights()[0].T

	(rows, cols) = (8, 4)
	plt.figure(figsize = (16, 24))
	gs = gridspec.GridSpec(rows, cols, wspace = 0.1, hspace = 0.1)

	for i in range(kernels.shape[0]):
	row = i // cols
	col = i % cols
	ax = plt.subplot(gs[row, col])

airalcorn2 / not_space_filling.py

Created May 9, 2017 21:55

Investigating http://koaning.io/homogenous-autoencoder-ui.html.

	import matplotlib.pyplot as plt
	import numpy as np
	import pandas as pd
	import seaborn as sns

	from keras import backend
	from keras.layers import Input, Dense
	from keras.models import Model

	# Autoencoder.

airalcorn2 / species_plots.py

Last active April 13, 2018 17:40

Python code to generate the plots used to make this GIF --> http://imgur.com/BjgN6UA.

	# Michael A. Alcorn

	import matplotlib
	matplotlib.use("Agg")
	import matplotlib.pyplot as plt
	import multiprocessing
	import numpy as np
	import pandas as pd
	import seaborn as sns
	import time

airalcorn2 / clade_counts.py

Created April 13, 2018 18:54

A script for counting the total species in each clade in a Newick formatted file.

	# Michael A. Alcorn

	from Bio import Phylo

	tree = Phylo.read("species_newick.txt", "newick")
	clade_counts = {}
	for clade in tree.find_clades():
	if not clade.is_terminal():
	clade_counts[str(clade)] = clade.count_terminals()

airalcorn2 / species_count.R

Created April 13, 2018 19:28

Code to generate the plot seen here --> https://www.facebook.com/photo.php?fbid=10103790197221251.

	# Michael A. Alcorn

	library(ggplot2)
	library(plyr)

	data <- read.csv("species.csv")
	data$Date <- as.Date(data$Date, "%Y-%m-%d")

	ggplot(data, aes(Date, Total)) + geom_line() + scale_x_date(date_breaks = "1 month")

airalcorn2 / attention_grok.py

Last active March 21, 2020 14:10

A barebones PyTorch implementation of a seq2seq model with attention.

	# Michael A. Alcorn

	import torch
	import torch.autograd as autograd
	import torch.nn as nn


	def create_lstm(params):
	"""Create a LSTM from a dictionary of parameters.

airalcorn2 / rtw_fatalities.py

Created July 3, 2018 16:36

Plotting fatality rates for states that had a change in right-to-work laws from this paper --> https://oem.bmj.com/content/early/2018/06/13/oemed-2017-104747.

	# Michael A. Alcorn

	import matplotlib.pyplot as plt
	import pandas as pd

	df = pd.read_csv("repdata.tab", sep = "\t", index_col = False)
	no_law = set(df[df["RTW"] == 0]["state"])
	yes_law = set(df[df["RTW"] == 1]["state"])
	states_that_changed = list(no_law & yes_law)
	states_that_changed.sort()

Michael A. Alcorn airalcorn2