The next Helmholtz Hacky Hour will take place on Wednesday, December 9, 2020 from 2PM to 3PM! Topic: Sustainable Programming! more...

Commit 64b2097b authored by Florian Centler's avatar Florian Centler

Update AnnotationCombine.r

parent 36b31899
......@@ -59,3 +59,15 @@ myTable <- table(data$Species, data$Level1, useNA="always")
#colnames(myTable)[ncol(myTable)] <- "unclassified"
write.table(myTable, "myTable_sampleID_taxa-vs-fct.txt", sep = "\t")
# Visualization starts here ...
plotData <- as.data.frame.matrix(myTable)
rownames(plotData)[nrow(plotData)] <- "unclassified"
colnames(plotData)[ncol(plotData)] <- "unclassified"
colFunc <- colorRampPalette(c("white", "steelblue"))
#heatmap(as.matrix(plotData), Rowv=NA, Colv=NA, col=colFunc(256), scale="none")
#heatmap(as.matrix(plotData), col=colFunc(256), scale="row")
heatmap(as.matrix(plotData), Rowv=NA, Colv=NA, col=colFunc(256), scale="row")
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment