Skip to content

Instantly share code, notes, and snippets.

View jwintersinger's full-sized avatar

Jeff Wintersinger jwintersinger

View GitHub Profile
Initial:
[[0 1 1 0 1]
[0 0 1 0 1]
[0 0 0 0 0]
[0 1 1 0 1]
[0 0 0 0 0]]
Top:
[[0 1 1 0 1]
[0 0 1 0 1]
@jwintersinger
jwintersinger / README
Last active October 21, 2015 19:59
Feature counts for ECE 1508 assignment
These are the proportions of variants possessing features used in each question for the various questions comprising ECE1508 assignment 1.
CnvParser
Caller: main()
Then: grouper.add_cnvs(cn_regions)
CnvFormatter
MultisampleCnvCombiner
Called by grouper.add_cnvs()
Segmenter
function fetch_packages(db, on_done) {
// Your code goes here. You must return an array of objects that looks like
// this (which you can infer from looking at what fields we reference in
// index.html):
//
// [ { sent_on: '2016-06-27 14:31:03',
// due_on: '2016-06-28 08:00:00',
// received_on: null,
// shipper_name: 'Nils Olav',
// receiver_name: 'William Windsor',
Do proportion of SVs plots
See what BAF and logR are of my BP clusters that don't have SVs -- are they subclonal? Do they have smaller deviations?
We should report these nulls back to PCAWG-6 -- now?
Can we use our results to help PCAWG-6 -- ask them to look for non-consensus SVs that one or two methods called around our wll-supported BPs
Send list of BPs to PCAWG6, along with CN states upstream and downstream
Send this list to group, which can check BAF and logR
Include TRAs, but throw out small SVs
Geoff: Anything < 500 kb is often missed by CN methods for duplications or deletions
http://prune.da-shiz.net/uploads/jeff/sbf-battlefield-vietnam-world-tournament-recordings/20IDvs.sbf_-_Jun.10.mp3
http://prune.da-shiz.net/uploads/jeff/sbf-battlefield-vietnam-world-tournament-recordings/Pinc.vs.pfg_-_Apr.27.mp3
from __future__ import print_function
import sys
from collections import defaultdict
def parse(fn):
loci = set()
with open(fn) as F:
for L in F:
L = L.strip()
#!/usr/bin/env python2
import subprocess
def move_window(oldpos, newpos):
return subprocess.check_output(['tmux', 'move-window', '-s', str(oldpos), '-t', str(newpos)])
def main():
output = subprocess.check_output(['tmux', 'list-windows'])
lines = [l.strip() for l in output.split('\n') if l.strip() != '']
frags = [line.strip().split(': ', 1) for line in output.split('\n') if line.strip() != '']
@jwintersinger
jwintersinger / convert_lichee_outputs.py
Created September 18, 2019 20:44
Convert LICHeE outputs
import re
import sys
import os
import argparse
import json
sys.path.append(os.path.join(os.path.dirname(__file__), '..', '..', 'lib'))
sys.path.append(os.path.join(os.path.dirname(__file__), '..'))
import common
import inputparser