Skip to content

Instantly share code, notes, and snippets.

@golobor
golobor / bottle-brush-scaling-theory.ipynb
Created Jun 22, 2020
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.
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
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.
@golobor
golobor / w303.centromeres.tsv
Last active Dec 14, 2018
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
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 / gist:fea8b74fce5e1b8119e6c8dc96671f07
Created Dec 1, 2018 — 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
@golobor
golobor / centromere_downloader.ipynb
Last active Oct 29, 2020
centromere_downloaded
View centromere_downloader.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 / test_dask.ipynb
Created Apr 16, 2017
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.