accaldwell

#!/usr/bin/env python
# quick script for demultiplexing reads from illumina fastq files
# Sean Davis <seandavi@gmail.com>
# 2012-06-29
#
import argparse
import Bio.SeqIO as SeqIO
from itertools import izip
import gzip
from string import maketrans

accaldwell

heatmapdata <- read.csv("~/Subsystem_by_Phylum-absolute.csv", sep="\t", check.names=FALSE)

row.names(heatmapdata) <- heatmapdata$Phylum

hmdata <- heatmapdata[,2:29]

myData <- data.matrix(hmdata)

# For melt() to work seamlessly, myData has to be a matrix.
longData <- melt(myData)