Replace DEanalysis.R

scripts/DEanalysis.R

@@ -4,6 +4,8 @@ padj <- read.csv("padjwoquote.csv", sep="\t", header=TRUE, row.names=1)
 
 # prep data
 padj[is.na(padj)] = 1
+logf[is.na(logf)] = 0
+
 
 # print # of DEG values for samples
 for (i in 1:ncol(padj)) {
@@ -18,4 +20,5 @@ DEGs <-as.matrix(DEGs)
 DEGs <- DEGs[,order(colnames(DEGs))]
 # plot data
 barplot(DEGs, main="Differentially expressed genes (DEGs)")
+boxplot(logf)