Mercurial > repos > yufei-luo > s_mart

comparison SMART/DiffExpAnal/DESeqTools/HTseqClean.R @ 31:0ab839023fe4

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author m-zytnicki
date Tue, 30 Apr 2013 14:33:21 -0400
parents 94ab73e8a190
children
comparison
equal deleted inserted replaced
30:5677346472b5 31:0ab839023fe4
1 # HTseqClean
2 # remove extra counts out of genes
3 # for HTseq output
4
5 # input : rawCounts
6 # output : cleaned rawCounts
7
8 # created Feb 6th, 2012
9 # Modified Feb 16th, 2012
10 # Marie-Agnes Dillies
11
12
13 HTseqClean <- function( rawCounts ){
14
15 row2remove <- c("alignment_not_unique", "ambiguous", "no_feature", "not_aligned", "too_low_aQual")
16 rawCounts <- rawCounts[!rawCounts$Id %in% row2remove,]
17 rawCounts[is.na(rawCounts)] <- 0
18 return(rawCounts)
19 }