Skip to content

Instantly share code, notes, and snippets.

import subprocess
import shutil
import argparse
import pathlib
parser = argparse.ArgumentParser(description='Download all SRAs from a GSM.')
parser.add_argument('GSMs', metavar='GSMs', type=str, nargs='+',
help='GSMs to download, can be more than one')
parser.add_argument('-g', '--guess-title',
help='If specified, guess the title of each GSMs and use it to store SRRs',
@golobor
golobor / gist:fea8b74fce5e1b8119e6c8dc96671f07
Created December 1, 2018 05:52 — forked from suimye/gist:d8edec74cce9ecad86cd
For retrieving fastq file from NCBI GEO using GSM number.
#!/bin/sh
if ! type "efetch" > /dev/null; then
print "Please install E-utilitie."
fi
GSM_LIST=$1
GSMs=`cat $GSM_LIST|cut -f1`
#echo $SRRs
! type "foo" > /dev/null 2>&1;
for GSM in $GSMs; do
import bpy
import mathutils
import numpy as np
def add_curve(
coords,
thickness = 0.2,
name='polymer',
resolution=4,
kind='BEZIER'):
@golobor
golobor / test_dask.ipynb
Created April 16, 2017 23:39
dask delayed+arrays for multiprocessing in hdf5
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@golobor
golobor / gridspec_inches.py
Created March 3, 2017 18:04
making grids of plots in units that make sense
def gridspec_inches(
wcols,
hrows,
fig_kwargs={}):
import matplotlib as mpl
import matplotlib.pyplot as plt
fig_height_inches = (
sum(hrows)