Mercurial > repos > mingchen0919 > aurora_htseq_site

annotate 01_htseq_count_analysis.Rmd @ 5:cf3fde5d0834 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
author mingchen0919
date Tue, 01 May 2018 11:27:32 -0400
parents
children a40d5db8687e
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
5
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
1 ---
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
2 title: 'HTSeq-count Analysis'
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
3 output:
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
4 html_document:
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
5 highlight: pygments
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
6 ---
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
7
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
8 ## Job script
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
9
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
10 ```{bash, echo=FALSE}
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
11 sh ${TOOL_INSTALL_DIR}/build-and-run-job-scripts.sh
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
12 ```
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
13
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
14 ```{r echo=FALSE, comment='', results='asis'}
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
15 cat('```bash\n')
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
16 cat(readLines(paste0(Sys.getenv('REPORT_FILES_PATH'), '/htseq-count.sh')), sep = '\n')
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
17 cat('\n```')
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
18 ```
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
19
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
20 ## Counts
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
21
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
22 Only the first few rows are displayed.
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
23
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
24 ```{r, echo=FALSE}
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
25 count_data = read.table(paste0(opt$X_d, '/counts.txt'), row.names = 1)
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
26 sample_names = trimws(strsplit(opt$X_B, ',')[[1]])
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
27 colnames(count_data) = rep(sample_names, length = ncol(count_data))
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
28 kable(head(count_data))
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
29
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
30 # modify column names
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
31 count_data = data.frame(feature_id = rownames(count_data), count_data)
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
32 write.csv(count_data,
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
33 file = paste0(Sys.getenv('REPORT_FILES_PATH'), '/count_data.csv'),
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
34 quote = FALSE, row.names = FALSE)
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
35 ```
cf3fde5d0834 planemo upload commit 104ae24ee30761a0099eeb91362ed1e3e13aba4b-dirty
mingchen0919
parents:
diff changeset
36