Skip to content

Instantly share code, notes, and snippets.

View ottov's full-sized avatar

Otto Valladares ottov

  • University of Pennsylvania
  • Philadelphia, PA
View GitHub Profile
@ottov
ottov / make_intervals.sh
Last active June 8, 2020 16:13
exclude SNPs, identify INDELs they have to be skipped
#!/usr/bin/env bash
#
# make_intervals.sh
# provides a list of intervals for a given chromosome
# -c contig (of the form: chrN)
# -b VCF file in BGZIP format, must have an index
# -p number of paritions
# -d path to genome DICT file (optional)
# -v verbose (optional)
# -h show usage

Keybase proof

I hereby claim:

  • I am ottov on github.
  • I am ottov (https://keybase.io/ottov) on keybase.
  • I have a public key ASAp_1bhPl6IVHVf1Kc3yzbiu962y1DMeXl9BtPSlzoP7Qo

To claim this, I am signing this object:

@ottov
ottov / infoStat.R
Created October 27, 2015 14:30
Summary stats on .info files from Michigan HRC imputation files
library(foreach)
library(doMC)
registerDoMC(8)
std_err <- function(x) sd(x)/sqrt(length(x))
freqbin=c(0,0.0025,0.005,0.0075,0.01,0.02,0.03,0.04,0.05,0.06,0.07,0.08,0.09,0.10,0.5)
freqbin1 <- c(0, 0.01, 0.10)
freqbin2 <- c(0, 0.8, 0.9, 0.95, 0.98, 0.99)
# Multiple plot function
#
# ggplot objects can be passed in ..., or to plotlist (as a list of ggplot objects)
# - cols: Number of columns in layout
# - layout: A matrix specifying the layout. If present, 'cols' is ignored.
#
# If the layout is something like matrix(c(1,2,3,3), nrow=2, byrow=TRUE),
# then plot 1 will go in the upper left, 2 will go in the upper right, and
# 3 will go all the way across the bottom.
#