# HG changeset patch # User yhoogstrate # Date 1389268528 18000 # Node ID c672e76503b2733726cd2d58f31dc7d9ea4e2f77 # Parent b8b3d63b60abbccad6d8cd5beac0ba76c49ddc42 Uploaded diff -r b8b3d63b60ab -r c672e76503b2 edgeR_DGE.xml --- a/edgeR_DGE.xml Thu Jan 09 06:48:17 2014 -0500 +++ b/edgeR_DGE.xml Thu Jan 09 06:55:28 2014 -0500 @@ -68,7 +68,12 @@ raw_data <- read.delim(designmatrix,header=T,stringsAsFactors=T) ## Obtain read-counts -read_counts = read.delim(as.character(raw_data[1,1]),header=F,stringsAsFactors=F,row.names=1) + +header = read.delim(as.character(raw_data[1,1]),header=F,stringsAsFactors=F,row.names=1,nrows=1) +has_header = (class(header[1,1]) == "character") + +read_counts = read.delim(as.character(raw_data[1,1]),header=has_header,stringsAsFactors=F,row.names=1) + for(i in 2:length(raw_data[,1])) { print("parsing counts from:") print(raw_data[i,1])