1 Kallisto and salmon: how similar are they?

1.1 version: 20180306

2 Create experiments from the two sets of counts

## Reading the sample metadata.
## The sample definitions comprises: 9, 12 rows, columns.
## Reading count tables.
## Reading salmon data with tximport.
## Finished reading count tables.
## Matched 81155 annotations and counts.
## Bringing together the count matrix and gene information.
## Some annotations were lost in merging, setting them to 'undefined'.
## Reading the sample metadata.
## The sample definitions comprises: 9, 12 rows, columns.
## Reading count tables.
## Reading kallisto inputs with tximport.
## Finished reading count tables.
## Matched 93422 annotations and counts.
## Bringing together the count matrix and gene information.
## Some annotations were lost in merging, setting them to 'undefined'.
## Using limma's removeBatchEffect to visualize before/after batch inclusion.
## Finished running DE analyses, collecting outputs.
## Comparing analyses 1/1: Rv_vs_Control
## Using limma's removeBatchEffect to visualize before/after batch inclusion.
## Finished running DE analyses, collecting outputs.
## Comparing analyses 1/1: Rv_vs_Control
## Writing a legend of columns.
## Working on table 1/1: Rv_vs_Control
## Writing a legend of columns.
## Working on table 1/1: Rv_vs_Control

3 TODO

  • 2018-03-06:
    • Load the data from them and see how they look.

index.html

## If you wish to reproduce this exact build of hpgltools, invoke the following:
## > git clone http://github.com/abelew/hpgltools.git
## > git reset 1b009834267dea125ee94934203413fbd606e783
## R> packrat::restore()
## This is hpgltools commit: Mon Apr 23 14:59:56 2018 -0400: 1b009834267dea125ee94934203413fbd606e783
## Saving to 02_kallisto_vs_salmon-v20180306.rda.xz

R version 3.4.4 (2018-03-15)

Platform: x86_64-pc-linux-gnu (64-bit)

locale: LC_CTYPE=en_US.utf8, LC_NUMERIC=C, LC_TIME=en_US.utf8, LC_COLLATE=en_US.utf8, LC_MONETARY=en_US.utf8, LC_MESSAGES=en_US.utf8, LC_PAPER=en_US.utf8, LC_NAME=C, LC_ADDRESS=C, LC_TELEPHONE=C, LC_MEASUREMENT=en_US.utf8 and LC_IDENTIFICATION=C

attached base packages: stats, graphics, grDevices, utils, datasets, methods and base

other attached packages: ruv(v.0.9.7) and hpgltools(v.2018.03)

loaded via a namespace (and not attached): nlme(v.3.1-137), bitops(v.1.0-6), matrixStats(v.0.53.1), devtools(v.1.13.5), bit64(v.0.9-7), doParallel(v.1.0.11), RColorBrewer(v.1.1-2), rprojroot(v.1.3-2), GenomeInfoDb(v.1.14.0), tools(v.3.4.4), backports(v.1.1.2), R6(v.2.2.2), rpart(v.4.1-13), Hmisc(v.4.1-1), DBI(v.0.8), lazyeval(v.0.2.1), BiocGenerics(v.0.24.0), mgcv(v.1.8-23), colorspace(v.1.3-2), nnet(v.7.3-12), withr(v.2.1.2), gridExtra(v.2.3), DESeq2(v.1.18.1), bit(v.1.1-12), compiler(v.3.4.4), preprocessCore(v.1.40.0), Biobase(v.2.38.0), htmlTable(v.1.11.2), xml2(v.1.2.0), DelayedArray(v.0.4.1), checkmate(v.1.8.5), scales(v.0.5.0.9000), DEoptimR(v.1.0-8), robustbase(v.0.92-8), readr(v.1.1.1), genefilter(v.1.60.0), commonmark(v.1.4), stringr(v.1.3.0), digest(v.0.6.15), foreign(v.0.8-70), rmarkdown(v.1.9), XVector(v.0.18.0), base64enc(v.0.1-3), pkgconfig(v.2.0.1), htmltools(v.0.3.6), limma(v.3.34.9), htmlwidgets(v.1.2), rlang(v.0.2.0.9001), rstudioapi(v.0.7), RSQLite(v.2.1.0), BiocParallel(v.1.12.0), acepack(v.1.4.1), RCurl(v.1.95-4.10), magrittr(v.1.5), GenomeInfoDbData(v.1.0.0), Formula(v.1.2-2), Matrix(v.1.2-14), Rcpp(v.0.12.16), munsell(v.0.4.3), S4Vectors(v.0.16.0), stringi(v.1.1.7), yaml(v.2.1.18), edgeR(v.3.20.9), SummarizedExperiment(v.1.8.1), zlibbioc(v.1.24.0), rhdf5(v.2.22.0), plyr(v.1.8.4), grid(v.3.4.4), blob(v.1.1.1), parallel(v.3.4.4), ggrepel(v.0.7.0), lattice(v.0.20-35), splines(v.3.4.4), pander(v.0.6.1), annotate(v.1.56.2), hms(v.0.4.2), locfit(v.1.5-9.1), knitr(v.1.20), pillar(v.1.2.1), GenomicRanges(v.1.30.3), rjson(v.0.2.15), corpcor(v.1.6.9), geneplotter(v.1.56.0), codetools(v.0.2-15), stats4(v.3.4.4), XML(v.3.98-1.11), evaluate(v.0.10.1), latticeExtra(v.0.6-28), data.table(v.1.10.4-3), foreach(v.1.4.4), gtable(v.0.2.0), ggplot2(v.2.2.1), openxlsx(v.4.0.17), xtable(v.1.8-2), roxygen2(v.6.0.1), survival(v.2.42-3), tibble(v.1.4.2), iterators(v.1.0.9), AnnotationDbi(v.1.40.0), memoise(v.1.1.0), IRanges(v.2.12.0), tximport(v.1.6.0), cluster(v.2.0.7-1) and sva(v.3.26.0)

LS0tCnRpdGxlOiAiTS50dWJlcmN1bG9zaXMgMjAxODogQ29tcGFyZSBrYWxsaXN0byBhbmQgc2FsbW9uIgphdXRob3I6ICJhdGIgYWJlbGV3QGdtYWlsLmNvbSIKZGF0ZTogImByIFN5cy5EYXRlKClgIgpvdXRwdXQ6CiBodG1sX2RvY3VtZW50OgogIGNvZGVfZG93bmxvYWQ6IHRydWUKICBjb2RlX2ZvbGRpbmc6IHNob3cKICBmaWdfY2FwdGlvbjogdHJ1ZQogIGZpZ19oZWlnaHQ6IDcKICBmaWdfd2lkdGg6IDcKICBoaWdobGlnaHQ6IHRhbmdvCiAga2VlcF9tZDogZmFsc2UKICBtb2RlOiBzZWxmY29udGFpbmVkCiAgbnVtYmVyX3NlY3Rpb25zOiB0cnVlCiAgc2VsZl9jb250YWluZWQ6IHRydWUKICB0aGVtZTogY29zbW8KICB0b2M6IHRydWUKICB0b2NfZmxvYXQ6CiAgICBjb2xsYXBzZWQ6IGZhbHNlCiAgICBzbW9vdGhfc2Nyb2xsOiBmYWxzZQotLS0KCjxzdHlsZT4KICBib2R5IC5tYWluLWNvbnRhaW5lciB7CiAgICBtYXgtd2lkdGg6IDE2MDBweDsKfQo8L3N0eWxlPgoKYGBge3Igb3B0aW9ucywgaW5jbHVkZT1GQUxTRX0KbGlicmFyeSgiaHBnbHRvb2xzIikKdHQgPC0gZGV2dG9vbHM6OmxvYWRfYWxsKCJ+L2hwZ2x0b29scyIpCmtuaXRyOjpvcHRzX2tuaXQkc2V0KHByb2dyZXNzPVRSVUUsCiAgICAgICAgICAgICAgICAgICAgIHZlcmJvc2U9VFJVRSwKICAgICAgICAgICAgICAgICAgICAgd2lkdGg9OTAsCiAgICAgICAgICAgICAgICAgICAgIGVjaG89VFJVRSkKa25pdHI6Om9wdHNfY2h1bmskc2V0KGVycm9yPVRSVUUsCiAgICAgICAgICAgICAgICAgICAgICBmaWcud2lkdGg9OCwKICAgICAgICAgICAgICAgICAgICAgIGZpZy5oZWlnaHQ9OCwKICAgICAgICAgICAgICAgICAgICAgIGRwaT05NikKb3B0aW9ucyhkaWdpdHM9NCwKICAgICAgICBzdHJpbmdzQXNGYWN0b3JzPUZBTFNFLAogICAgICAgIGtuaXRyLmR1cGxpY2F0ZS5sYWJlbD0iYWxsb3ciKQpnZ3Bsb3QyOjp0aGVtZV9zZXQoZ2dwbG90Mjo6dGhlbWVfYncoYmFzZV9zaXplPTEwKSkKc2V0LnNlZWQoMSkKdmVyIDwtICIyMDE4MDMwNiIKcHJldmlvdXNfZmlsZSA8LSAiMDFfYW5ub3RhdGlvbi5SbWQiCgp0bXAgPC0gdHJ5KHNtKGxvYWRtZShmaWxlbmFtZT1wYXN0ZTAoZ3N1YihwYXR0ZXJuPSJcXC5SbWQiLCByZXBsYWNlPSIiLCB4PXByZXZpb3VzX2ZpbGUpLCAiLXYiLCB2ZXIsICIucmRhLnh6IikpKSkKcm1kX2ZpbGUgPC0gIjAyX2thbGxpc3RvX3ZzX3NhbG1vbi5SbWQiCmBgYAoKS2FsbGlzdG8gYW5kIHNhbG1vbjogaG93IHNpbWlsYXIgYXJlIHRoZXk/Cj09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PQoKIyMgdmVyc2lvbjogYHIgdmVyYAoKIyBDcmVhdGUgZXhwZXJpbWVudHMgZnJvbSB0aGUgdHdvIHNldHMgb2YgY291bnRzCgpgYGB7ciBleHB0c30Kcm93bmFtZXMoaHNfYW5ub3QpIDwtIG1ha2UubmFtZXMocGFzdGUwKGhzX2Fubm90W1siZW5zZW1ibF90cmFuc2NyaXB0X2lkIl1dLCAiLiIsIGhzX2Fubm90W1sidHJhbnNjcmlwdF92ZXJzaW9uIl1dKSwgdW5pcXVlPVRSVUUpCnNhbG1vbl9leHB0IDwtIGNyZWF0ZV9leHB0KG1ldGFkYXRhPSJzYW1wbGVfc2hlZXRzL3JuYXNlcV9zYW1wbGVzLnhsc3giLAogICAgICAgICAgICAgICAgICAgICAgICAgICBnZW5lX2luZm89aHNfYW5ub3QsIGZpbGVfY29sdW1uPSJzYWxtb25maWxlIikKa2FsbGlzdG9fZXhwdCA8LSBjcmVhdGVfZXhwdChtZXRhZGF0YT0ic2FtcGxlX3NoZWV0cy9ybmFzZXFfc2FtcGxlcy54bHN4IiwKICAgICAgICAgICAgICAgICAgICAgICAgICAgICBnZW5lX2luZm89aHNfYW5ub3QsIGZpbGVfY29sdW1uPSJrYWxsaXN0b2ZpbGUiKQoKc2FsbW9uX2RlIDwtIGFsbF9wYWlyd2lzZShzYWxtb25fZXhwdCkKa2FsbGlzdG9fZGUgPC0gYWxsX3BhaXJ3aXNlKGthbGxpc3RvX2V4cHQpCnNhbG1vbl90YWIgPC0gY29tYmluZV9kZV90YWJsZXMoc2FsbW9uX2RlKQprYWxsaXN0b190YWIgPC0gY29tYmluZV9kZV90YWJsZXMoa2FsbGlzdG9fZGUpCgpzaW1pbGFycCA8LSBjb21wYXJlX2RlX3Jlc3VsdHMoc2FsbW9uX3RhYiwga2FsbGlzdG9fdGFiLCBjb3JfbWV0aG9kPSJzcGVhcm1hbiIpCgpgYGAKCiMgVE9ETwoKKiAyMDE4LTAzLTA2OgogIC0gIExvYWQgdGhlIGRhdGEgZnJvbSB0aGVtIGFuZCBzZWUgaG93IHRoZXkgbG9vay4KCltpbmRleC5odG1sXShpbmRleC5odG1sKQoKYGBge3Igc2F2ZW1lfQptZXNzYWdlKHBhc3RlMCgiVGhpcyBpcyBocGdsdG9vbHMgY29tbWl0OiAiLCBnZXRfZ2l0X2NvbW1pdCgpKSkKdGhpc19zYXZlIDwtIHBhc3RlMChnc3ViKHBhdHRlcm49IlxcLlJtZCIsIHJlcGxhY2U9IiIsIHg9cm1kX2ZpbGUpLCAiLXYiLCB2ZXIsICIucmRhLnh6IikKbWVzc2FnZShwYXN0ZTAoIlNhdmluZyB0byAiLCB0aGlzX3NhdmUpKQp0bXAgPC0gc20oc2F2ZW1lKGZpbGVuYW1lPXRoaXNfc2F2ZSkpCnBhbmRlcjo6cGFuZGVyKHNlc3Npb25JbmZvKCkpCmBgYAo=