Mercurial > repos > davidvanzessen > argalaxy_tools
view experimental_design/experimental_design.r @ 3:beaa487ecf43 draft
Uploaded
author | davidvanzessen |
---|---|
date | Wed, 07 Dec 2016 08:36:23 -0500 |
parents | afe85eb6572e |
children | 106275b54470 |
line wrap: on
line source
args <- commandArgs(trailingOnly = TRUE) print(args) inputs = args[1:(length(args) - 1)] output = args[length(args)] current.id = "" counter = 1 result = NULL for(current in inputs){ if(grepl("/", current)){ #its a path to a file print(paste("Adding file", counter, "to", current.id)) dat = read.table(current, sep="\t", header=T, quote="", fill=T) #IMGT check dat$Sample = current.id dat$Replicate = counter if(is.null(result)){ result = dat[NULL,] } result = rbind(result, dat) counter = counter + 1 } else { #its an ID of a patient print(paste("New patient", current)) current.id = current counter = 1 } } write.table(result, output, sep="\t", quote=F, row.names=F, col.names=T)