utils/saarclust_output_accuracy_plot.snakemake.R

log <- file(snakemake@log[[1]], open='wt')
sink(file=log, type='message')
sink(file=log, type='output')

source("postProcessing.R")

clusters <- snakemake@input[["clusters"]]
raw_data <- snakemake@input[["raw_data"]]
v <- ClustersAccuracyPerChrPerDir(Clusters2process=clusters, Quals2process=raw_data)

write.table(v, file = snakemake@output[[1]], quote = F, sep = "\t", row.names = F)
daewoooo/SaaRclust documentation built on May 28, 2019, 7:50 p.m.