Mercurial > repos > davidvanzessen > mutation_analysis
view pattern_plots.r @ 123:0453ea4d9f14 draft
Uploaded
author | davidvanzessen |
---|---|
date | Mon, 22 Aug 2016 07:00:23 -0400 |
parents | |
children | 4a93146f87aa |
line wrap: on
line source
library(ggplot2) library(reshape2) args <- commandArgs(trailingOnly = TRUE) input.file = args[1] #the data that's get turned into the "SHM overview" table in the html report "data_sum.txt" plot1.file = args[2] plot2.file = args[3] plot3.file = args[4] dat = read.table(input.file, header=F, sep=",", quote="", stringsAsFactors=F, fill=T) classes = c("ca", "ca1", "ca2", "cg", "cg1", "cg2", "cg3", "cg4", "cm") xyz = c("x", "y", "z") names(dat) = c("info", paste(rep(classes, each=3), xyz, sep="."), paste("un", xyz, sep="."), paste("all", xyz, sep="."))