Mercurial > repos > mingchen0919 > aurora_fastqc_site

comparison 08_sequence_length_distribution.Rmd @ 2:c64267b9f754 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload commit 841d8b22bf9f1aaed6bfe8344b60617f45b275b2-dirty
author mingchen0919
date Sun, 30 Dec 2018 12:48:14 -0500
parents b7c115edd970
children
comparison
equal deleted inserted replaced
1:92e3de11b9f8 2:c64267b9f754
1 --- 1 ---
2 output: html_document 2 output: html_document
3 --- 3 ---
4 4
5 ```{r setup, include=FALSE, warning=FALSE, message=FALSE} 5 ```{r setup, include=FALSE, warning=FALSE, message=FALSE}
6 knitr::opts_chunk$set( 6 knitr::opts_knit$set(progress = FALSE)
7 echo = as.logical(opt$X_e), 7 knitr::opts_chunk$set(error = TRUE, echo = FALSE)
8 error = TRUE,
9 eval = TRUE
10 )
11 ``` 8 ```
12 9
13 # Sequence Length Distribution 10 ### Sequence Length Distribution
14 11
15 ```{r 'Sequence Length Distribution', fig.width=10} 12 ```{r 'Sequence Length Distribution', fig.width=10}
16 ## reads 1 13 ## reads 1
17 sld_1 = extract_data_module(paste0(opt$X_d, '/read_1_fastqc/fastqc_data.txt'), 'Sequence Length Distribution') 14 sld_1 = extract_data_module(paste0(opt$X_d, '/read_1_fastqc/fastqc_data.txt'), 'Sequence Length Distribution')
18 sld_1$id = 1:length(sld_1$X.Length) 15 sld_1$id = 1:length(sld_1$X.Length)
30 geom_line(color = 'red') + 27 geom_line(color = 'red') +
31 scale_x_continuous(breaks = sld_2$id, labels = sld_2$X.Length) + 28 scale_x_continuous(breaks = sld_2$id, labels = sld_2$X.Length) +
32 facet_grid(. ~ trim) + 29 facet_grid(. ~ trim) +
33 xlab('Sequence Length (bp)') + 30 xlab('Sequence Length (bp)') +
34 ylab('') + 31 ylab('') +
35 theme(axis.text.x = element_text(angle=45)) 32 theme(axis.text.x = element_text(size = 5),
33 panel.background = element_rect(fill = NA),
34 axis.line = element_line(),
35 plot.margin = margin(2,2,2,10) )
36 ggplotly(p) 36 ggplotly(p)
37 ``` 37 ```