Commit fc3930a5 authored by Schoorlemmer, Joran's avatar Schoorlemmer, Joran
Browse files

Upload New File

parent 02c1cb67
# import data
logf <- read.csv("log2fwoquote.csv", sep="\t", header=TRUE, row.names=1)
padj <- read.csv("padjwoquote.csv", sep="\t", header=TRUE, row.names=1)
# prep data
padj[is.na(padj)] = 1
# print # of DEG values for samples
for (i in 1:ncol(padj)) {
print(colnames(padj)[i])
print(length(padj[,i][padj[,i]<0.05]))
}
# manually write data to csv file in notepad
# read DEG data
DEGs <-read.csv("nDEGs.csv", sep="\t", header=TRUE, row.names=1)
DEGs <-as.matrix(DEGs)
DEGs <- DEGs[,order(colnames(DEGs))]
# plot data
barplot(DEGs, main="Differentially expressed genes (DEGs)")
Supports Markdown
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