Skip to content

GitLab

Commit 26436db7 authored by Sascha Herzinger's avatar Sascha Herzinger
Browse files

added outrider processing

parent ddea3d4e
Show whitespace changes
Inline Side-by-side
process.R
View file @ 26436db7
...@@ -44,6 +44,13 @@ write.table(fraserTsv, sep = "\t", file = "fraser_complete_output.tsv", row.name ...@@ -44,6 +44,13 @@ write.table(fraserTsv, sep = "\t", file = "fraser_complete_output.tsv", row.name
outriderTsvPath <- "./Output/processed_results/aberrant_expression/v29/outrider/outrider/OUTRIDER_results.tsv" outriderTsvPath <- "./Output/processed_results/aberrant_expression/v29/outrider/outrider/OUTRIDER_results.tsv"
outriderRdsPath <- "./Output/processed_results/aberrant_expression/v29/outrider/outrider/ods.Rds" outriderRdsPath <- "./Output/processed_results/aberrant_expression/v29/outrider/outrider/ods.Rds"
outriderTsv <- read.table(outriderTsvPath, sep="\t", header=TRUE) outriderTsv <- read.table(outriderTsvPath, sep="\t", header=TRUE)
ods <- readRDS(outriderRdsPath) ods <- readRDS(outriderRdsPath)
ods
gr <- as.data.frame(granges(ods)) gr <- as.data.frame(granges(ods))
sampleRank <- rank(outriderTsv["normcounts"])
outriderTsv <- cbind(outriderTsv, sampleRank)
write.table(outriderTsv, sep = "\t", file = "outrider_complete_output.tsv", row.names=FALSE)
\ No newline at end of file
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