Skip to content

Instantly share code, notes, and snippets.

View mpkocher's full-sized avatar

M. Kocher mpkocher

View GitHub Profile
@mpkocher
mpkocher / abstract_graph_example.xml
Last active August 29, 2015 13:56
Tasks XML Manifest
<?xml version="1.0" encoding="UTF-8"?>
<workflow>
<!-- Not sure this make sense to be here -->
<settings>
<param id="pbsmrtpipe.task_modules.filter.my_option">
<value>12</value>
</param>
</settings>
<abstract-graph id="pbsmrtpipe.abstract_graphs.rs_resequencing">
<file-node id="2">
@mpkocher
mpkocher / test_filter.py
Last active August 29, 2015 13:56
New Pysiv API to support file Validation
"""Unittests for P_Filter module"""
import logging
from pysiv.core.base import TestBase, TestModuleBase
from pysiv.core.utils import monkey_patch_task, monkey_patch_module
from pysiv.validators import FastaValidator, FastqValidator, CsvValidator
log = logging.getLogger(__name__)
@mpkocher
mpkocher / example_validation_subparsers.py
Last active August 29, 2015 13:56
Commandline Subparser with Validation
import os
import sys
import argparse
import logging
__version__ = '0.1'
log = logging.getLogger(__name__)
@mpkocher
mpkocher / summarizeCompareByMovie.py
Created February 25, 2014 18:51
SMRTAnalysis 2.2 summarizeCompareByMovie that supports expired SMRT Cell lots (bug 24721).
#!/usr/bin/env python
from pprint import pformat
__doc__ = """Script to break down post-mapping statistics by movie.
OLD Pre-bax files
m130715_185638_SMRT1_c000000062559900001500000112311501_s1_p0.bas.h5
m{DATE_DATE}_{INSTRUMENT_NAME}_c{CHIP_STRIP_BARCODE}{CELL_NUMBER}_s{
@mpkocher
mpkocher / Notes.md
Last active August 29, 2015 13:56
P_Control.py Disable task (which calls summarizeCompareByMovie.py) to handle expired SMRT Cell lots.

Updating 2.1

Replace $SEYMOUR_HOME/analysis/lib/python2.7/pbpy-0.1-py2.7.egg/pbpy/smrtpipe/modules/P_Control.py in SMRTAnalysis 2.1.

@mpkocher
mpkocher / bug_status_by_state.py
Created March 17, 2014 23:15
Generate bugs to verify from bugzilla CSV by state
import os
import pandas as pd
import pysiv.toolbook as tb
_STATES = ("RESOLVED", "ASSIGNED", "FIXED")
def make_bug_summary(csv_file_name, output_image, state="RESOLVED"):
if state not in _STATES:
raise ValueError("Unsupported state {s}".format(s=state))
@mpkocher
mpkocher / tool_runner.py
Created March 29, 2014 23:39
General Pacbio Tool test runner
import os
import logging
import abc
import functools
import networkx as nx
from fabric.api import lcd, local, shell_env, prefix, settings
log = logging.getLogger(__name__)
@mpkocher
mpkocher / compose_functions.py
Created March 31, 2014 16:18
compose functions in python
def compose(*funcs):
"""Functional composition
[f, g, h] will be f(g(h(x)))
"""
def compose_two(f, g):
def c(x):
return f(g(x))
return c
return functools.reduce(compose_two, funcs)
@mpkocher
mpkocher / mhouse_runner.py
Created May 2, 2014 17:31
Submission to Milhouse via commandline tool.
"""
Usage:
python mhouse_runner.py submit my_awesome.csv --debug -m "http://milhouse:9999" internal_test.csv -a '.m2__m_loading_stats' '.m2__gAccuracy' '.m2__gAccuracyByReadLengthBoxplot' -n 'mk_test_01' --description "services-submitted-job"
This is included in the milhouse code (milhouse/src/milhouse/scripts/mhouse_runner.py
"""
import os
from pprint import pformat
import sys
@mpkocher
mpkocher / add_to_bashrc.sh
Last active August 29, 2015 14:02
Get current changelist from SEYMOUR_HOME
# Add to .bashrc
function seamour() {
#source $SEYMOUR_HOME/etc/setup.sh
x=`env | grep SEYMOUR_HOME | awk -F "=" '{print $2}'`
source $x/etc/setup.sh
current_version=`get_current_version.sh`
export PS1="[${current_version}]$PS1"
}