I hereby claim:
- I am rob-p on github.
- I am rob_p (https://keybase.io/rob_p) on keybase.
- I have a public key whose fingerprint is C7B0 CCB0 FCB8 E9C3 DFB8 A9BF 2EAB 9D25 7707 B487
To claim this, I am signing this object:
from math import * | |
# function [alpha, xmin, L]=plfit(x, varargin) | |
# PLFIT fits a power-law distributional model to data. | |
# Source: http://www.santafe.edu/~aaronc/powerlaws/ | |
# | |
# PLFIT(x) estimates x_min and alpha according to the goodness-of-fit | |
# based method described in Clauset, Shalizi, Newman (2007). x is a | |
# vector of observations of some quantity to which we wish to fit the | |
# power-law distribution p(x) ~ x^-alpha for x >= xmin. |
import requests | |
import argparse | |
import pprint | |
import os | |
def do_upload(args): | |
ACCESS_TOKEN = args.token | |
depid = args.deposition_id | |
url = 'https://zenodo.org/api/deposit/depositions/{}'.format(depid) | |
r = requests.get(url, params={'access_token': ACCESS_TOKEN}) |
import argparse | |
import pandas as pd | |
import numpy as np | |
import sys | |
def main(args): | |
gtable = pd.read_table(args.ginput).set_index('Name') | |
ttable = pd.read_table(args.tinput).set_index('Name') | |
tgmap = pd.read_table(args.tgmap, names=['t', 'g']).set_index('t') | |
gene_lengths = {} |
#!/bin/bash | |
set -e | |
# Activate Holy Build Box environment. | |
source /hbb_exe/activate | |
set -x | |
# Install things we need | |
yum install -y --quiet wget |
def effLen(L, pmf, meanLen): | |
if L < meanLen: | |
return L | |
minVal = 0 | |
maxVal = len(pmf)-1 | |
mval = maxVal | |
clen = minVal | |
maxLen = min(L, mval) | |
effectiveLength = 0 | |
while (clen <= maxLen): |
#!/bin/bash | |
cmd="$@" | |
interleaved_file=`echo $cmd | sed -n 's/.*--interleaved\s\+\(\S\+\)\s\+.*/\1/p'` | |
if [ -z "$interleaved_file" ] | |
then | |
#Run normally in this branch | |
${@} | |
else |
I hereby claim:
To claim this, I am signing this object:
{ | |
"cells": [ | |
{ | |
"cell_type": "code", | |
"execution_count": 21, | |
"metadata": { | |
"collapsed": true | |
}, | |
"outputs": [], | |
"source": [ |
seqtk mergepe ${leftIn} ${rightIn} | pv -r -b | paste - - - - - - - - | shuf | tr "\t" "\n" | tee >(seqtk seq -1 - > ${outBase}_1.fastq) | seqtk seq -2 - > ${outBase}_2.fastq |