Last updated: 2018-07-30
workflowr checks: (Click a bullet for more information) ✔ R Markdown file: up-to-date
Great! Since the R Markdown file has been committed to the Git repository, you know the exact version of the code that produced these results.
✔ Environment: empty
Great job! The global environment was empty. Objects defined in the global environment can affect the analysis in your R Markdown file in unknown ways. For reproduciblity it’s best to always run the code in an empty environment.
✔ Seed:
set.seed(12345)
The command set.seed(12345)
was run prior to running the code in the R Markdown file. Setting a seed ensures that any results that rely on randomness, e.g. subsampling or permutations, are reproducible.
✔ Session information: recorded
Great job! Recording the operating system, R version, and package versions is critical for reproducibility.
✔ Repository version: 422a428
wflow_publish
or wflow_git_commit
). workflowr only checks the R Markdown file, but you know if there are other scripts or data files that it depends on. Below is the status of the Git repository when the results were generated:
Ignored files:
Ignored: .DS_Store
Ignored: .Rhistory
Ignored: .Rproj.user/
Ignored: output/.DS_Store
Untracked files:
Untracked: data/18486.genecov.txt
Untracked: data/APApeaksYL.total.inbrain.bed
Untracked: data/Totalpeaks_filtered_clean.bed
Untracked: data/YL-SP-18486-T_S9_R1_001-genecov.txt
Untracked: data/bedgraph_peaks/
Untracked: data/bin200.5.T.nuccov.bed
Untracked: data/bin200.Anuccov.bed
Untracked: data/bin200.nuccov.bed
Untracked: data/clean_peaks/
Untracked: data/combined_reads_mapped_three_prime_seq.csv
Untracked: data/gene_cov/
Untracked: data/leafcutter/
Untracked: data/nuc6up/
Untracked: data/reads_mapped_three_prime_seq.csv
Untracked: data/smash.cov.results.bed
Untracked: data/smash.cov.results.csv
Untracked: data/smash.cov.results.txt
Untracked: data/smash_testregion/
Untracked: data/ssFC200.cov.bed
Untracked: output/picard/
Untracked: output/plots/
Untracked: output/qual.fig2.pdf
Unstaged changes:
Modified: analysis/cleanupdtseq.internalpriming.Rmd
Modified: analysis/dif.iso.usage.leafcutter.Rmd
Modified: analysis/explore.filters.Rmd
Modified: analysis/test.max2.Rmd
Modified: code/Snakefile
Note that any generated files, e.g. HTML, png, CSS, etc., are not included in this status report because it is ok for generated content to have uncommitted changes.
File | Version | Author | Date | Message |
---|---|---|---|---|
Rmd | 422a428 | Briana Mittleman | 2018-07-30 | add peak cove pipeline and combined lane qc |
I want to use this analysis to run simple QC on the first 32 libraries now that we have 2 lanes per library.
First, I will look at the new map stats to see how many more reads/mapped reads the socond lane provided.
library(tidyr)
library(reshape2)
Attaching package: 'reshape2'
The following object is masked from 'package:tidyr':
smiths
library(ggplot2)
library(dplyr)
Attaching package: 'dplyr'
The following objects are masked from 'package:stats':
filter, lag
The following objects are masked from 'package:base':
intersect, setdiff, setequal, union
library(workflowr)
This is workflowr version 1.1.1
Run ?workflowr for help getting started
comb_map=read.csv("../data/combined_reads_mapped_three_prime_seq.csv", header = T, stringsAsFactors = T)
comb_map$line=as.factor(comb_map$line)
mapped_melt=melt(comb_map, id.vars=c("line", "fraction"), measure.vars = c( "lane1_mapped", "comb_mapped"))
mapped_melt$line=as.factor(mapped_melt$line)
ggplot(mapped_melt, aes(y=value, x=line, by=fraction,fill=fraction)) + geom_bar(stat="identity", position = "dodge") + facet_grid(.~ variable)+ labs(y="Reads Mapped")
Next I want to look at the x more mapped reads we got by line and fraction
ggplot(comb_map, aes(x=line,y=combed_xrmappedmore, fill=fraction)) + geom_bar(stat="identity", position = "dodge") + labs(title="X more mapped reads in adding second lane", y="X more mapped reads")
sessionInfo()
R version 3.5.1 (2018-07-02)
Platform: x86_64-apple-darwin15.6.0 (64-bit)
Running under: macOS Sierra 10.12.6
Matrix products: default
BLAS: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRblas.0.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRlapack.dylib
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] workflowr_1.1.1 dplyr_0.7.6 ggplot2_3.0.0 reshape2_1.4.3
[5] tidyr_0.8.1
loaded via a namespace (and not attached):
[1] Rcpp_0.12.18 pillar_1.3.0 compiler_3.5.1
[4] git2r_0.23.0 plyr_1.8.4 bindr_0.1.1
[7] R.methodsS3_1.7.1 R.utils_2.6.0 tools_3.5.1
[10] digest_0.6.15 evaluate_0.11 tibble_1.4.2
[13] gtable_0.2.0 pkgconfig_2.0.1 rlang_0.2.1
[16] rstudioapi_0.7 yaml_2.1.19 bindrcpp_0.2.2
[19] withr_2.1.2 stringr_1.3.1 knitr_1.20
[22] tidyselect_0.2.4 rprojroot_1.3-2 grid_3.5.1
[25] glue_1.3.0 R6_2.2.2 rmarkdown_1.10
[28] purrr_0.2.5 magrittr_1.5 whisker_0.3-2
[31] backports_1.1.2 scales_0.5.0 htmltools_0.3.6
[34] assertthat_0.2.0 colorspace_1.3-2 labeling_0.3
[37] stringi_1.2.4 lazyeval_0.2.1 munsell_0.5.0
[40] crayon_1.3.4 R.oo_1.22.0
This reproducible R Markdown analysis was created with workflowr 1.1.1