Skip to content

Instantly share code, notes, and snippets.

View squirrelo's full-sized avatar

Joshua Shorenstein squirrelo

View GitHub Profile

Keybase proof

I hereby claim:

  • I am squirrelo on github.
  • I am squirrelo (https://keybase.io/squirrelo) on keybase.
  • I have a public key whose fingerprint is 7E05 CAC7 FFF0 E4C1 9CAC 5611 CB2E 4A92 49B4 62F9

To claim this, I am signing this object:

@squirrelo
squirrelo / scrub.py
Created October 12, 2015 06:41
Scrubber for AG database pulldowns
import psycopg2
from passlib.hash import bcrypt
from random import choice
from sys import argv
def _new_kit_id():
seen = set()
alpha = "abcdefghijklmnopqrstuvwxyz"
alpha += alpha.upper()
while True:
#!/usr/bin/env python
from datetime import datetime
import re
import binascii
import os
from copy import deepcopy
from amgut.lib.data_access.sql_connection import SQLConnectionHandler
from amgut.connections import ag_data
@squirrelo
squirrelo / compare_all_categories
Last active April 6, 2016 18:52
This runs permanova or anosim on all catgories in your qiime mapping file. Simmilar to the compare_category.py script in QIIME, but runs over all categories at once.
#!/usr/bin/env python
import click
from pandas.io.parsers import read_table
from pandas import DataFrame
from skbio.stats.distance import DistanceMatrix, permanova, anosim
from statsmodels.sandbox.stats.multicomp import multipletests
DIST_TESTS = {
@squirrelo
squirrelo / bloom_filter
Created April 16, 2015 17:52
bloom_filter
#!/usr/bin/env python
from os.path import join, basename, splitext
from subprocess import Popen, PIPE
from shlex import split as comm_split
import click
from qiime_default_reference import get_reference_taxonomy
@click.command()
@click.option('--seqs_fp', '-i', type=click.Path(
@squirrelo
squirrelo / create_test_skeletion.py
Last active August 31, 2015 22:03
Test skeleton creator
# with help from http://stackoverflow.com/q/17171878
import ast
from argparse import ArgumentParser
from os.path import basename, splitext
parser = ArgumentParser(description='Creates test code skeleton for a given file')
parser.add_argument('-i', type=str,
help='File to generate test framework for')
parser.add_argument('-o', type=str,
help='File to write test framework to.')
@squirrelo
squirrelo / average_movie.py
Created April 8, 2014 05:41
Simple image parser for averaging frames in a movie
from argparse import ArgumentParser
#requires numpy, moviepy, and PIL libraries
from moviepy.editor import VideoFileClip
from numpy import zeros, array, uint8
from PIL import Image
def frange(x, y, inc):
while x < y:
yield x