annotate rmarkdown_deseq2_count_matrix.Rmd @ 1:4a7131658ca6 draft

update tool
author mingchen0919
date Sat, 30 Dec 2017 12:35:31 -0500
parents 867a77a7a380
children 18e3fc69da76
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
1 ---
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
2 title: 'DESeq2 analysis'
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
3 output:
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
4 html_document:
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
5 number_sections: true
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
6 toc: true
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
7 theme: cosmo
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
8 highlight: tango
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
9 ---
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
10
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
11 ```{r setup, include=FALSE, warning=FALSE, message=FALSE}
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
12 knitr::opts_chunk$set(
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
13 echo = opt$echo,
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
14 error = TRUE
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
15 )
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
16 ```
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
17
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
18
1
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
19 # User input
0
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
20
1
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
21 ```{r 'user input'}
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
22 df = data.frame(name = names(opt)[-1],
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
23 value = unlist(opt))
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
24 datatable(df, rownames = FALSE)
0
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
25 ```
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
26
867a77a7a380 planemo upload
mingchen0919
parents:
diff changeset
27
1
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
28 # Count Matrix
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
29
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
30 Display the first 100 rows of count data matrix.
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
31
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
32 ```{r 'count matrix'}
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
33 count_data = read.table(opt$count_data)
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
34 col_names = trimws(strsplit(opt$count_matrix_column_names, ',')[[1]])[1:ncol(count_data)]
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
35 col_names = col_names[!is.na(col_names)]
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
36 colnames(count_data)[1:length(col_names)] = col_names
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
37 datatable(head(count_data, 100))
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
38 ```
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
39
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
40 # Column Data
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
41
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
42 ```{r 'column data'}
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
43 col_data = read.table(opt$col_data,
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
44 stringsAsFactors = FALSE, sep=',', header = TRUE, row.names = 1)
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
45 datatable(col_data)
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
46 ```
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
47
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
48 # Match sample names
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
49
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
50 The goal of this step is to rearrange the rows of the column data matrix so that the samples rows in the count data matrix and the sample columns in the count data matrix are in the same order.
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
51
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
52 ```{r 'match sample names'}
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
53 col_data = col_data[col_names, ]
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
54 datatable(col_data)
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
55 ```
4a7131658ca6 update tool
mingchen0919
parents: 0
diff changeset
56