comparison experimental_design/experimental_design.r @ 40:106275b54470 draft

Uploaded
author davidvanzessen
date Tue, 30 May 2017 07:26:33 -0400
parents afe85eb6572e
children 124b7fd92a3e
comparison
equal deleted inserted replaced
39:bad6a9a53ae7 40:106275b54470
12 12
13 for(current in inputs){ 13 for(current in inputs){
14 if(grepl("/", current)){ #its a path to a file 14 if(grepl("/", current)){ #its a path to a file
15 print(paste("Adding file", counter, "to", current.id)) 15 print(paste("Adding file", counter, "to", current.id))
16 dat = read.table(current, sep="\t", header=T, quote="", fill=T) 16 dat = read.table(current, sep="\t", header=T, quote="", fill=T)
17
18 if(nrow(dat) == 0){
19 print(paste(counter, "of", current.id, "has no sequences, skipping"))
20 next
21 }
17 22
18 #IMGT check 23 #IMGT check
19 24
20 dat$Sample = current.id 25 dat$Sample = current.id
21 dat$Replicate = counter 26 dat$Replicate = counter