annotate deseq/differential_expression_analysis_pipeline_for_rnaseq_data-a03838a6eb54/DiffExpAnal/DESeqTools/raw2counts.R @ 10:6e573fd3c41b draft

Uploaded
author yufei-luo
date Mon, 13 May 2013 10:06:30 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
10
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
1 # raw2counts
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
2 # extract counts only from rawCounts
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
3 # and add rownames to counts
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
4
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
5 # input : rawCounts
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
6 # output : counts
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
7
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
8 # created Feb 6th, 2012
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
9 # modified April 12, 2012
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
10 # Marie-Agnes Dillies
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
11
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
12
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
13 raw2counts <- function( rawCounts, annot=1 ){
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
14
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
15 ex <- 1:annot
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
16 counts <- as.matrix( rawCounts[,-ex] )
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
17 rownames(counts) <- rawCounts[,1]
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
18 infoCounts <- rawCounts[,ex]
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
19 return( list("counts"=counts, "infoCounts"= infoCounts) )
6e573fd3c41b Uploaded
yufei-luo
parents:
diff changeset
20 }