Mercurial > repos > yufei-luo > s_mart

view SMART/DiffExpAnal/DESeqTools/HTseqClean.R @ 32:3441fe98a2ba

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Deleted selected files
author m-zytnicki
date Tue, 30 Apr 2013 14:34:10 -0400
parents 94ab73e8a190
children
line wrap: on
line source

# HTseqClean
# remove extra counts out of genes
# for HTseq output

# input : rawCounts
# output : cleaned rawCounts

# created Feb 6th, 2012
# Modified Feb 16th, 2012
# Marie-Agnes Dillies


HTseqClean <- function( rawCounts ){

  row2remove <- c("alignment_not_unique", "ambiguous", "no_feature", "not_aligned", "too_low_aQual")
  rawCounts <- rawCounts[!rawCounts$Id %in% row2remove,]
  rawCounts[is.na(rawCounts)] <- 0
  return(rawCounts)
}