Mercurial > repos > yufei-luo > differential_expression_analysis_pipeline_for_rnaseq_data
view DiffExpAnal/DESeqTools/HTseqClean.R @ 0:63799b789162 draft
Uploaded
author | yufei-luo |
---|---|
date | Tue, 22 Jan 2013 10:07:03 -0500 |
parents | |
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) }