Skip to content

Instantly share code, notes, and snippets.

@golobor
golobor / gridspec_inches.py
Created March 3, 2017 18:04
making grids of plots in units that make sense
View gridspec_inches.py
def gridspec_inches(
wcols,
hrows,
fig_kwargs={}):
import matplotlib as mpl
import matplotlib.pyplot as plt
fig_height_inches = (
sum(hrows)
@golobor
golobor / test_dask.ipynb
Created April 16, 2017 23:39
dask delayed+arrays for multiprocessing in hdf5
View test_dask.ipynb
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
View polender.py
import bpy
import mathutils
import numpy as np
def add_curve(
coords,
thickness = 0.2,
name='polymer',
resolution=4,
kind='BEZIER'):
@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.
View gist:fea8b74fce5e1b8119e6c8dc96671f07
#!/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
View sra_downloader.py
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 / w303.centromeres.tsv
Last active December 14, 2018 05:50
Yeast cohesin supplemental files (w303 genome)
View w303.centromeres.tsv
chrom start end
chrVII 496916 497034
chrVI 148509 148626
chrV 151987 152104
chrIV 449712 449822
chrIII 114385 114501
chrII 238206 238322
chrI 151465 151582
chrIX 355629 355745
chrVIII 105591 105708
@golobor
golobor / yeast_cohesin_project_fig6.ipynb
Created December 14, 2018 06:09
Yeast dCohesin Hi-C, fig6
View yeast_cohesin_project_fig6.ipynb
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
View my_savez.py
# from https://stackoverflow.com/questions/22712292/how-to-use-numpy-savez-in-a-loop-for-save-more-than-one-array/22716159#22716159
import numpy as np
import tempfile
class my_savez(object):
def __init__(self, file, compression=None, compresslevel=None):
# Import is postponed to here since zipfile depends on gzip, an optional
# component of the so-called standard library.
import zipfile
@golobor
golobor / bottle-brush-scaling-theory.ipynb
Created June 22, 2020 20:26
Bottle Brush P(s) theory for Gibcus et at 2018
View bottle-brush-scaling-theory.ipynb
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.