wflynny

import scanpy.api as sc
import matplotlib.pyplot as plt
import seaborn as sns

def get_cluster_proportions(adata,
                            cluster_key="cluster_final",
                            sample_key="replicate",
                            drop_values=None):
    """
    Input

wflynny

# Find myself referring to this thread a lot:
# https://www.biostars.org/p/198143/
# However updating codes with what I see at JAX
@Mxxxx - MiSeq
@Dxxxx - HiSeq 2500
@Kxxxx - HiSeq 4000
@NSxxx - NextSeq 500/550
@Axxxxx - NovaSeq

wflynny

Introduction

This process has gotten simpler with 10X Genomics Loupe Browser 5.0.0+. I try to keep a cache of all 10X downloads internal to JAX here (must be on the VPN):

http://singlecell-software.s3-far.jax.org/list.html

where you can download the latest version of Loupe Browser for your OS. Otherwise, head over to the 10X site:

https://support.10xgenomics.com/single-cell-gene-expression/software/downloads/latest#loupe

wflynny

import scanpy as sc
import pandas as pd
import GEOparse

def gse_to_adata(gse):
    X = gse.pivot_samples("VALUE")
    key = list(gse.gpls.keys())[0]
    annotated = X.index.to_frame(index=False).merge(
        gse.gpls[key].table[["ID", "Gene Symbol"]], 
        left_on="ID_REF", right_on="ID"

wflynny

# I usually put this in my ~/.bash_aliases
# A portion of our GPFS storage removes files after 21 days of creation.
# `stat` does not show creation time, so we have to resort to parsing the
# output of `mmlsattr`

ftime() {
    # Usage: 
    #   ftime path/to/file
    # 
    # Outputs:

wflynny

###################
# chord diagram
import matplotlib.pyplot as plt
from matplotlib.path import Path as mplPath
import matplotlib.patches as patches

import numpy as np
import pandas as pd

LW = 0.3

wflynny

from matplotlib.text import Annotation
from matplotlib.legend_handler import HandlerBase
from matplotlib import patheffects
from matplotlib.colors import CSS4_COLORS, hex2color

class AnnotationHandler(HandlerBase):
    def create_artists(self, legend, artist, xdescent, ydescent,
                       width, height, fontsize, trans):
        a = Annotation(
            artist.get_text(), [width/2, height/2],

wflynny

import re
from pathlib import Path
from xml.etree import ElementTree as ET

import numpy as np
from skimage import exposure
from skimage.io import imread, imsave, imshow
from skimage.filters import threshold_otsu, gaussian

wflynny

#!/usr/bin/env bash

TEMP=$(getopt -o hsg --long help,snapshot,gpu -n 'susuage' -- "$@")

if [ $? != 0 ] ; then echo "Terminating..." >&2 ; exit 1 ; fi

# Note the quotes around `$TEMP': they are essential!
eval set -- "$TEMP"

SNAPSHOT=false

wflynny

import re
import argparse

parser = argparse.ArgumentParser()
parser.add_argument("-i", "--infile", required=True)
parser.add_argument("-o", "--outfile", required=True)
args = parser.parse_args()

gene_matcher = re.compile('\tgene\t.*gene_id (".*?");.*Name (".*?");')
parent_matcher = re.compile('gene_id (".*?");.*Parent (".*?");')