Commit e550dccd authored by Brandi Cantarel's avatar Brandi Cantarel

process scripts

parents 9d8cac39 d97dfb21
Pipeline #2335 failed with stage
in 41 minutes and 8 seconds
process.executor='slurm'
process.queue='128GB,256GB,256GBv1'
process.clusterOptions = '--hold'
trace.enabled = true
trace.enabled = false
trace.file = 'pipeline_trace.txt'
trace.field = 'task_id,native_id,process,name,status,exit,submit,start,complete,duration,realtime'
......@@ -94,6 +94,7 @@ for (i in 1:a) {
filt.out <- na.omit(output[output$fdr < 0.05,])
if (nrow(filt.out) > 2) {
subset <- logcpm[row.names(logcpm) %in% filt.out$symbol,]
subset <- subset[!apply(subset, 1, function(x) {any(x == 0)}),]
gnames <- filt.out[c('ensembl','symbol')]
s <- merge(gnames,subset,by.x="ensembl",by.y="row.names",all.x=FALSE,all.y=TRUE,sort=FALSE)
STREE <- hclust(dist(t(subset)))
......@@ -148,6 +149,7 @@ for (i in 1:a) {
filt.out <- na.omit(output[output$fdr < 0.05,])
if (nrow(filt.out) > 2) {
subset <- logcpm[row.names(logcpm) %in% filt.out$symbol,]
subset <- subset[!apply(subset, 1, function(x) {any(x == 0)}),]
gnames <- filt.out[c('ensembl','symbol')]
s <- merge(gnames,subset,by.x="ensembl",by.y="row.names",all.x=FALSE,all.y=TRUE,sort=FALSE)
STREE <- hclust(dist(t(subset)))
......
Markdown is supported
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