Mercurial > repos > iuc > deseq2

annotate deseq2.R @ 43:621ddf5967d2 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
author iuc
date Tue, 18 Jul 2023 14:58:28 +0000
parents 6ef2cba4e35a
children 1cb33de18af5
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
1 #!/usr/bin/env Rscript
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
2
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
3 # A command-line interface to DESeq2 for use with Galaxy
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
4 # written by Bjoern Gruening and modified by Michael Love 2016.03.30
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
5 #
26
47267a044ef1 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit e811a7887db870f4f94f620f52bce656c8d5ba23
iuc
parents: 24
diff changeset
6 # This argument is required:
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
7 #
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
8 # 'factors' a JSON list object from Galaxy
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
9 #
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
10 # the output file has columns:
27
dc6bc19f05ab planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents: 26
diff changeset
11 #
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
12 # baseMean (mean normalized count)
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
13 # log2FoldChange (by default a moderated LFC estimate)
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
14 # lfcSE (the standard error)
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
15 # stat (the Wald statistic)
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
16 # pvalue (p-value from comparison of Wald statistic to a standard Normal)
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
17 # padj (adjusted p-value, Benjamini Hochberg correction on genes which pass the mean count filter)
27
dc6bc19f05ab planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents: 26
diff changeset
18 #
26
47267a044ef1 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit e811a7887db870f4f94f620f52bce656c8d5ba23
iuc
parents: 24
diff changeset
19 # the first variable in 'factors' will be the primary factor.
47267a044ef1 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit e811a7887db870f4f94f620f52bce656c8d5ba23
iuc
parents: 24
diff changeset
20 # the levels of the primary factor are used in the order of appearance in factors.
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
21 #
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
22 # by default, levels in the order A,B,C produces a single comparison of B vs A, to a single file 'outfile'
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
23 #
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
24 # for the 'many_contrasts' flag, levels in the order A,B,C produces comparisons C vs A, B vs A, C vs B,
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
25 # to a number of files using the 'outfile' prefix: 'outfile.condition_C_vs_A' etc.
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
26 # all plots will still be sent to a single PDF, named by the arg 'plots', with extra pages.
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
27 #
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
28 # fit_type is an integer valued argument, with the options from ?estimateDisperions
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
29 # 1 "parametric"
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
30 # 2 "local"
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
31 # 3 "mean"
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
32
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
33 # setup R error handling to go to stderr
42
6ef2cba4e35a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 6868b66f73ddbe947986d1a20b546873cbd515a9
iuc
parents: 41
diff changeset
34 options(show.error.messages = FALSE, error = function() {
6ef2cba4e35a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 6868b66f73ddbe947986d1a20b546873cbd515a9
iuc
parents: 41
diff changeset
35 cat(geterrmessage(), file = stderr())
6ef2cba4e35a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 6868b66f73ddbe947986d1a20b546873cbd515a9
iuc
parents: 41
diff changeset
36 q("no", 1, FALSE)
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
37 })
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
38
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
39 library("getopt")
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
40 library("tools")
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
41 options(stringAsFactors = FALSE, useFancyQuotes = FALSE)
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
42 args <- commandArgs(trailingOnly = TRUE)
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
43
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
44 # get options, using the spec as defined by the enclosed list.
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
45 # we read the options from the default: commandArgs(TRUE).
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
46 spec <- matrix(c(
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
47 "quiet", "q", 0, "logical",
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
48 "help", "h", 0, "logical",
33
4f497a8a72e8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 05bc13637dc9f8e523bc72844ff6eb0794f98ad3
iuc
parents: 31
diff changeset
49 "cores", "s", 0, "integer",
4f497a8a72e8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 05bc13637dc9f8e523bc72844ff6eb0794f98ad3
iuc
parents: 31
diff changeset
50 "batch_factors", "w", 1, "character",
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
51 "outfile", "o", 1, "character",
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
52 "countsfile", "n", 1, "character",
41
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
53 "sizefactorsfile", "F", 1, "character",
30
3b9af9a27b5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents: 29
diff changeset
54 "rlogfile", "r", 1, "character",
3b9af9a27b5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents: 29
diff changeset
55 "vstfile", "v", 1, "character",
27
dc6bc19f05ab planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents: 26
diff changeset
56 "header", "H", 0, "logical",
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
57 "factors", "f", 1, "character",
24
1a34aec44d60 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 62e9101c1e7b8467e395f31ccbd9321de01a6418
iuc
parents: 23
diff changeset
58 "files_to_labels", "l", 1, "character",
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
59 "plots", "p", 1, "character",
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
60 "tximport", "i", 0, "logical",
23
d1241c2b81a7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 0a56599c36b4968095ec5a3cb589f94fb139466c
iuc
parents: 21
diff changeset
61 "txtype", "y", 1, "character",
31
7a189a7634c7 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 5b6dc96c6e14582d5bb1dc213ac8d26dc7b2829e
iuc
parents: 30
diff changeset
62 "tx2gene", "x", 1, "character", # a space-sep tx-to-gene map or GTF/GFF3 file
30
3b9af9a27b5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents: 29
diff changeset
63 "esf", "e", 1, "character",
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
64 "fit_type", "t", 1, "integer",
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
65 "many_contrasts", "m", 0, "logical",
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
66 "outlier_replace_off", "a", 0, "logical",
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
67 "outlier_filter_off", "b", 0, "logical",
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
68 "auto_mean_filter_off", "c", 0, "logical",
43
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
69 "use_beta_priors", "d", 0, "logical",
42
6ef2cba4e35a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 6868b66f73ddbe947986d1a20b546873cbd515a9
iuc
parents: 41
diff changeset
70 "alpha_ma", "A", 1, "numeric",
6ef2cba4e35a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 6868b66f73ddbe947986d1a20b546873cbd515a9
iuc
parents: 41
diff changeset
71 "prefilter", "P", 0, "logical",
6ef2cba4e35a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 6868b66f73ddbe947986d1a20b546873cbd515a9
iuc
parents: 41
diff changeset
72 "prefilter_value", "V", 1, "numeric"
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
73 ), byrow = TRUE, ncol = 4)
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
74 opt <- getopt(spec)
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
75
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
76 # if help was asked for print a friendly message
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
77 # and exit with a non-zero error code
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
78 if (!is.null(opt$help)) {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
79 cat(getopt(spec, usage = TRUE))
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
80 q(status = 1)
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
81 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
82
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
83 # enforce the following required arguments
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
84 if (is.null(opt$outfile)) {
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
85 cat("'outfile' is required\n")
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
86 q(status = 1)
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
87 }
26
47267a044ef1 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit e811a7887db870f4f94f620f52bce656c8d5ba23
iuc
parents: 24
diff changeset
88 if (is.null(opt$factors)) {
47267a044ef1 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit e811a7887db870f4f94f620f52bce656c8d5ba23
iuc
parents: 24
diff changeset
89 cat("'factors' is required\n")
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
90 q(status = 1)
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
91 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
92
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
93 verbose <- is.null(opt$quiet)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
94
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
95 source_local <- function(fname) {
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
96 argv <- commandArgs(trailingOnly = FALSE)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
97 base_dir <- dirname(substring(argv[grep("--file=", argv)], 8))
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
98 source(paste(base_dir, fname, sep = "/"))
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
99 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
100
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
101 source_local("get_deseq_dataset.R")
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
102
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
103 suppressPackageStartupMessages({
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
104 library("DESeq2")
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
105 library("RColorBrewer")
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
106 library("gplots")
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
107 })
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
108
33
4f497a8a72e8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 05bc13637dc9f8e523bc72844ff6eb0794f98ad3
iuc
parents: 31
diff changeset
109 if (opt$cores > 1) {
4f497a8a72e8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 05bc13637dc9f8e523bc72844ff6eb0794f98ad3
iuc
parents: 31
diff changeset
110 library("BiocParallel")
4f497a8a72e8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 05bc13637dc9f8e523bc72844ff6eb0794f98ad3
iuc
parents: 31
diff changeset
111 register(MulticoreParam(opt$cores))
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
112 parallel <- TRUE
33
4f497a8a72e8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 05bc13637dc9f8e523bc72844ff6eb0794f98ad3
iuc
parents: 31
diff changeset
113 } else {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
114 parallel <- FALSE
33
4f497a8a72e8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 05bc13637dc9f8e523bc72844ff6eb0794f98ad3
iuc
parents: 31
diff changeset
115 }
4f497a8a72e8 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 05bc13637dc9f8e523bc72844ff6eb0794f98ad3
iuc
parents: 31
diff changeset
116
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
117 # build or read sample table
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
118
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
119 trim <- function(x) gsub("^\\s+|\\s+$", "", x)
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
120
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
121 # switch on if 'factors' was provided:
26
47267a044ef1 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit e811a7887db870f4f94f620f52bce656c8d5ba23
iuc
parents: 24
diff changeset
122 library("rjson")
47267a044ef1 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit e811a7887db870f4f94f620f52bce656c8d5ba23
iuc
parents: 24
diff changeset
123 parser <- newJSONParser()
47267a044ef1 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit e811a7887db870f4f94f620f52bce656c8d5ba23
iuc
parents: 24
diff changeset
124 parser$addData(opt$factors)
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
125 factor_list <- parser$getObject()
26
47267a044ef1 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit e811a7887db870f4f94f620f52bce656c8d5ba23
iuc
parents: 24
diff changeset
126 filenames_to_labels <- fromJSON(opt$files_to_labels)
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
127 factors <- sapply(factor_list, function(x) x[[1]])
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
128 primary_factor <- factors[1]
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
129 filenames_in <- unname(unlist(factor_list[[1]][[2]]))
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
130 labs <- unname(unlist(filenames_to_labels[basename(filenames_in)]))
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
131 sample_table <- data.frame(
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
132 sample = basename(filenames_in),
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
133 filename = filenames_in,
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
134 row.names = filenames_in,
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
135 stringsAsFactors = FALSE
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
136 )
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
137 for (factor in factor_list) {
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
138 factor_name <- trim(factor[[1]])
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
139 sample_table[[factor_name]] <- character(nrow(sample_table))
26
47267a044ef1 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit e811a7887db870f4f94f620f52bce656c8d5ba23
iuc
parents: 24
diff changeset
140 lvls <- sapply(factor[[2]], function(x) names(x))
47267a044ef1 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit e811a7887db870f4f94f620f52bce656c8d5ba23
iuc
parents: 24
diff changeset
141 for (i in seq_along(factor[[2]])) {
47267a044ef1 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit e811a7887db870f4f94f620f52bce656c8d5ba23
iuc
parents: 24
diff changeset
142 files <- factor[[2]][[i]][[1]]
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
143 sample_table[files, factor_name] <- trim(lvls[i])
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
144 }
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
145 sample_table[[factor_name]] <- factor(sample_table[[factor_name]], levels = lvls)
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
146 }
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
147 rownames(sample_table) <- labs
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
148
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
149 design_formula <- as.formula(paste("~", paste(rev(factors), collapse = " + ")))
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
150
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
151 # these are plots which are made once for each analysis
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
152 generate_generic_plots <- function(dds, factors) {
24
1a34aec44d60 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 62e9101c1e7b8467e395f31ccbd9321de01a6418
iuc
parents: 23
diff changeset
153 library("ggplot2")
1a34aec44d60 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 62e9101c1e7b8467e395f31ccbd9321de01a6418
iuc
parents: 23
diff changeset
154 library("ggrepel")
1a34aec44d60 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 62e9101c1e7b8467e395f31ccbd9321de01a6418
iuc
parents: 23
diff changeset
155 library("pheatmap")
1a34aec44d60 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 62e9101c1e7b8467e395f31ccbd9321de01a6418
iuc
parents: 23
diff changeset
156
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
157 rld <- rlog(dds)
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
158 p <- plotPCA(rld, intgroup = rev(factors))
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
159 print(p + geom_text_repel(aes_string(x = "PC1", y = "PC2", label = factor(colnames(dds))), size = 3) + geom_point())
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
160 dat <- assay(rld)
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
161 dists_rl <- dist(t(dat))
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
162 mat <- as.matrix(dists_rl)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
163 colors <- colorRampPalette(rev(brewer.pal(9, "Blues")))(255)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
164 pheatmap(
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
165 mat,
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
166 clustering_distance_rows = dists_rl,
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
167 clustering_distance_cols = dists_rl,
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
168 col = colors,
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
169 main = "Sample-to-sample distances"
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
170 )
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
171 plotDispEsts(dds, main = "Dispersion estimates")
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
172 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
173
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
174 # these are plots which can be made for each comparison, e.g.
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
175 # once for C vs A and once for B vs A
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
176 generate_specific_plots <- function(res, threshold, title_suffix) {
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
177 use <- res$baseMean > threshold
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
178 if (sum(!use) == 0) {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
179 h <- hist(res$pvalue, breaks = 0:50 / 50, plot = FALSE)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
180 barplot(
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
181 height = h$counts,
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
182 col = "powderblue",
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
183 space = 0,
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
184 xlab = "p-values",
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
185 ylab = "frequency",
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
186 main = paste("Histogram of p-values for", title_suffix)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
187 )
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
188 text(x = c(0, length(h$counts)), y = 0, label = paste(c(0, 1)), adj = c(0.5, 1.7), xpd = NA)
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
189 } else {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
190 h1 <- hist(res$pvalue[!use], breaks = 0:50 / 50, plot = FALSE)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
191 h2 <- hist(res$pvalue[use], breaks = 0:50 / 50, plot = FALSE)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
192 colori <- c("filtered (low count)" = "khaki", "not filtered" = "powderblue")
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
193 barplot(
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
194 height = rbind(h1$counts, h2$counts),
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
195 beside = FALSE,
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
196 col = colori,
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
197 space = 0,
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
198 xlab = "p-values",
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
199 ylab = "frequency",
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
200 main = paste("Histogram of p-values for", title_suffix)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
201 )
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
202 text(x = c(0, length(h1$counts)), y = 0, label = paste(c(0, 1)), adj = c(0.5, 1.7), xpd = NA)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
203 legend("topright", fill = rev(colori), legend = rev(names(colori)), bg = "white")
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
204 }
39
b30a443484ab "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit c8e7020d39d581d7374a13cf94c64998a6481e05"
iuc
parents: 37
diff changeset
205 plotMA(res, main = paste("MA-plot for", title_suffix), ylim = range(res$log2FoldChange, na.rm = TRUE), alpha = opt$alpha_ma)
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
206 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
207
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
208 if (verbose) {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
209 cat(paste("primary factor:", primary_factor, "\n"))
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
210 if (length(factors) > 1) {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
211 cat(paste("other factors in design:", paste(factors[-length(factors)], collapse = ","), "\n"))
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
212 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
213 cat("\n---------------------\n")
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
214 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
215
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
216 dds <- get_deseq_dataset(sample_table, header = opt$header, design_formula = design_formula, tximport = opt$tximport, txtype = opt$txtype, tx2gene = opt$tx2gene)
43
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
217
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
218 # use/estimate size factors with the chosen method
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
219 if (!is.null(opt$esf)) {
43
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
220 if (opt$esf %in% list("ratio", "poscounts", "iterate")) {
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
221 cat("Calculating size factors de novo\n")
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
222 dds <- estimateSizeFactors(dds, type = opt$esf)
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
223 } else {
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
224 sf_table <- read.table(opt$esf)
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
225 # Sort the provided size factors just in case the order differs from the input file order.
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
226 merged_table <- merge(sample_table, sf_table, by.x = 0, by.y = 1, sort = FALSE)
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
227 sf_values <- as.numeric(unlist(merged_table[5]))
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
228 "sizeFactors"(dds) <- sf_values
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
229
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
230 cat("Using user-provided size factors:\n")
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
231 print(sf_values)
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
232 }
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
233 } else {
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
234 cat("No size factor was used\n")
30
3b9af9a27b5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents: 29
diff changeset
235 }
41
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
236
43
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
237
41
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
238 # estimate size factors for each sample
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
239 # - https://support.bioconductor.org/p/97676/
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
240 if (!is.null(opt$sizefactorsfile)) {
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
241 nm <- assays(dds)[["avgTxLength"]]
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
242 if (!is.null(nm)) {
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
243 ## Recommended: takes into account tximport data
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
244 cat("\nsize factors for samples: taking tximport data into account\n")
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
245 size_factors <- estimateSizeFactorsForMatrix(counts(dds) / nm)
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
246 } else {
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
247 norm_factors <- normalizationFactors(dds)
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
248 if (!is.null(norm_factors)) {
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
249 ## In practice, gives same results as above.
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
250 cat("\nsize factors for samples: no tximport data, using derived normalization factors\n")
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
251 size_factors <- estimateSizeFactorsForMatrix(norm_factors)
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
252 } else {
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
253 ## If we have no other information, estimate from raw.
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
254 cat("\nsize factors for samples: no tximport data, no normalization factors, estimating from raw data\n")
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
255 size_factors <- estimateSizeFactorsForMatrix(counts(dds))
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
256 }
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
257 }
42
6ef2cba4e35a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 6868b66f73ddbe947986d1a20b546873cbd515a9
iuc
parents: 41
diff changeset
258 write.table(size_factors, file = opt$sizefactorsfile, sep = "\t", col.names = FALSE, quote = FALSE)
41
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
259 }
0a0a3388e3f2 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents: 39
diff changeset
260
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
261 apply_batch_factors <- function(dds, batch_factors) {
29
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
262 rownames(batch_factors) <- batch_factors$identifier
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
263 batch_factors <- subset(batch_factors, select = -c(identifier, condition))
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
264 dds_samples <- colnames(dds)
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
265 batch_samples <- rownames(batch_factors)
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
266 if (!setequal(batch_samples, dds_samples)) {
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
267 stop("Batch factor names don't correspond to input sample names, check input files")
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
268 }
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
269 dds_data <- colData(dds)
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
270 # Merge dds_data with batch_factors using indexes, which are sample names
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
271 # Set sort to False, which maintains the order in dds_data
42
6ef2cba4e35a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 6868b66f73ddbe947986d1a20b546873cbd515a9
iuc
parents: 41
diff changeset
272 reordered_batch <- merge(dds_data, batch_factors, by.x = 0, by.y = 0, sort = FALSE)
29
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
273 batch_factors <- reordered_batch[, ncol(dds_data):ncol(reordered_batch)]
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
274 for (factor in colnames(batch_factors)) {
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
275 dds[[factor]] <- batch_factors[[factor]]
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
276 }
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
277 colnames(dds) <- reordered_batch[, 1]
29
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
278 return(dds)
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
279 }
27
dc6bc19f05ab planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents: 26
diff changeset
280
29
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
281 if (!is.null(opt$batch_factors)) {
42
6ef2cba4e35a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 6868b66f73ddbe947986d1a20b546873cbd515a9
iuc
parents: 41
diff changeset
282 batch_factors <- read.table(opt$batch_factors, sep = "\t", header = TRUE)
29
66203279224f planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents: 28
diff changeset
283 dds <- apply_batch_factors(dds = dds, batch_factors = batch_factors)
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
284 batch_design <- colnames(batch_factors)[-c(1, 2)]
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
285 design_formula <- as.formula(paste("~", paste(c(batch_design, rev(factors)), collapse = " + ")))
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
286 design(dds) <- design_formula
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
287 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
288
27
dc6bc19f05ab planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents: 26
diff changeset
289 if (verbose) {
dc6bc19f05ab planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents: 26
diff changeset
290 cat("DESeq2 run information\n\n")
dc6bc19f05ab planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents: 26
diff changeset
291 cat("sample table:\n")
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
292 print(sample_table[, -c(1:2), drop = FALSE])
27
dc6bc19f05ab planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents: 26
diff changeset
293 cat("\ndesign formula:\n")
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
294 print(design_formula)
27
dc6bc19f05ab planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents: 26
diff changeset
295 cat("\n\n")
dc6bc19f05ab planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents: 26
diff changeset
296 cat(paste(ncol(dds), "samples with counts over", nrow(dds), "genes\n"))
dc6bc19f05ab planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents: 26
diff changeset
297 }
dc6bc19f05ab planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents: 26
diff changeset
298
42
6ef2cba4e35a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 6868b66f73ddbe947986d1a20b546873cbd515a9
iuc
parents: 41
diff changeset
299 # minimal pre-filtering
6ef2cba4e35a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 6868b66f73ddbe947986d1a20b546873cbd515a9
iuc
parents: 41
diff changeset
300 if (!is.null(opt$prefilter)) {
6ef2cba4e35a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 6868b66f73ddbe947986d1a20b546873cbd515a9
iuc
parents: 41
diff changeset
301 keep <- rowSums(counts(dds)) >= opt$prefilter_value
6ef2cba4e35a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 6868b66f73ddbe947986d1a20b546873cbd515a9
iuc
parents: 41
diff changeset
302 dds <- dds[keep, ]
6ef2cba4e35a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 6868b66f73ddbe947986d1a20b546873cbd515a9
iuc
parents: 41
diff changeset
303 }
6ef2cba4e35a planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 6868b66f73ddbe947986d1a20b546873cbd515a9
iuc
parents: 41
diff changeset
304
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
305 # optional outlier behavior
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
306 if (is.null(opt$outlier_replace_off)) {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
307 min_rep <- 7
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
308 } else {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
309 min_rep <- Inf
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
310 if (verbose) cat("outlier replacement off\n")
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
311 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
312 if (is.null(opt$outlier_filter_off)) {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
313 cooks_cutoff <- TRUE
27
dc6bc19f05ab planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents: 26
diff changeset
314 } else {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
315 cooks_cutoff <- FALSE
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
316 if (verbose) cat("outlier filtering off\n")
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
317 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
318
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
319 # optional automatic mean filtering
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
320 if (is.null(opt$auto_mean_filter_off)) {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
321 independent_filtering <- TRUE
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
322 } else {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
323 independent_filtering <- FALSE
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
324 if (verbose) cat("automatic filtering on the mean off\n")
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
325 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
326
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
327 # shrinkage of LFCs
43
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
328 if (is.null(opt$use_beta_priors)) {
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
329 beta_prior <- FALSE
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
330 if (verbose)
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
331 cat("Applied default - beta prior off\n")
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
332 } else {
43
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
333 beta_prior <- opt$use_beta_priors
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
334 }
43
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
335 sprintf("use_beta_prior is set to %s", beta_prior)
621ddf5967d2 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 469558ddf5bc6249874fe5826637fd6ee81588cf
iuc
parents: 42
diff changeset
336
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
337
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
338 # dispersion fit type
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
339 if (is.null(opt$fit_type)) {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
340 fit_type <- "parametric"
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
341 } else {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
342 fit_type <- c("parametric", "local", "mean")[opt$fit_type]
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
343 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
344
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
345 if (verbose) cat(paste("using disperion fit type:", fit_type, "\n"))
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
346
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
347 # run the analysis
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
348 dds <- DESeq(dds, fitType = fit_type, betaPrior = beta_prior, minReplicatesForReplace = min_rep, parallel = parallel)
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
349
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
350 # create the generic plots and leave the device open
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
351 if (!is.null(opt$plots)) {
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
352 if (verbose) cat("creating plots\n")
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
353 pdf(opt$plots)
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
354 generate_generic_plots(dds, factors)
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
355 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
356
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
357 n <- nlevels(colData(dds)[[primary_factor]])
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
358 all_levels <- levels(colData(dds)[[primary_factor]])
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
359
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
360 if (!is.null(opt$countsfile)) {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
361 normalized_counts <- counts(dds, normalized = TRUE)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
362 write.table(normalized_counts, file = opt$countsfile, sep = "\t", col.names = NA, quote = FALSE)
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
363 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
364
30
3b9af9a27b5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents: 29
diff changeset
365 if (!is.null(opt$rlogfile)) {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
366 rlog_normalized <- rlogTransformation(dds)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
367 rlog_normalized_mat <- assay(rlog_normalized)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
368 write.table(rlog_normalized_mat, file = opt$rlogfile, sep = "\t", col.names = NA, quote = FALSE)
30
3b9af9a27b5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents: 29
diff changeset
369 }
3b9af9a27b5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents: 29
diff changeset
370
3b9af9a27b5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents: 29
diff changeset
371 if (!is.null(opt$vstfile)) {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
372 vst_normalized <- varianceStabilizingTransformation(dds)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
373 vst_normalized_mat <- assay(vst_normalized)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
374 write.table(vst_normalized_mat, file = opt$vstfile, sep = "\t", col.names = NA, quote = FALSE)
30
3b9af9a27b5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents: 29
diff changeset
375 }
3b9af9a27b5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents: 29
diff changeset
376
3b9af9a27b5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents: 29
diff changeset
377
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
378 if (is.null(opt$many_contrasts)) {
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
379 # only contrast the first and second level of the primary factor
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
380 ref <- all_levels[1]
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
381 lvl <- all_levels[2]
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
382 res <- results(
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
383 dds,
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
384 contrast = c(primary_factor, lvl, ref),
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
385 cooksCutoff = cooks_cutoff,
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
386 independentFiltering = independent_filtering
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
387 )
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
388 if (verbose) {
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
389 cat("summary of results\n")
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
390 cat(paste0(primary_factor, ": ", lvl, " vs ", ref, "\n"))
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
391 print(summary(res))
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
392 }
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
393 res_sorted <- res[order(res$padj), ]
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
394 out_df <- as.data.frame(res_sorted)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
395 out_df$geneID <- rownames(out_df) # nolint
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
396 out_df <- out_df[, c("geneID", "baseMean", "log2FoldChange", "lfcSE", "stat", "pvalue", "padj")]
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
397 filename <- opt$outfile
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
398 write.table(out_df, file = filename, sep = "\t", quote = FALSE, row.names = FALSE, col.names = FALSE)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
399 if (independent_filtering) {
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
400 threshold <- unname(attr(res, "filterThreshold"))
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
401 } else {
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
402 threshold <- 0
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
403 }
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
404 title_suffix <- paste0(primary_factor, ": ", lvl, " vs ", ref)
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
405 if (!is.null(opt$plots)) {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
406 generate_specific_plots(res, threshold, title_suffix)
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
407 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
408 } else {
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
409 # rotate through the possible contrasts of the primary factor
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
410 # write out a sorted table of results with the contrast as a suffix
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
411 # add contrast specific plots to the device
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
412 for (i in seq_len(n - 1)) {
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
413 ref <- all_levels[i]
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
414 contrast_levels <- all_levels[(i + 1):n]
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
415 for (lvl in contrast_levels) {
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
416 res <- results(
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
417 dds,
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
418 contrast = c(primary_factor, lvl, ref),
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
419 cooksCutoff = cooks_cutoff,
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
420 independentFiltering = independent_filtering
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
421 )
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
422 res_sorted <- res[order(res$padj), ]
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
423 out_df <- as.data.frame(res_sorted)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
424 out_df$geneID <- rownames(out_df) # nolint
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
425 out_df <- out_df[, c("geneID", "baseMean", "log2FoldChange", "lfcSE", "stat", "pvalue", "padj")]
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
426 filename <- paste0(primary_factor, "_", lvl, "_vs_", ref)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
427 write.table(out_df, file = filename, sep = "\t", quote = FALSE, row.names = FALSE, col.names = FALSE)
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
428 if (independent_filtering) {
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
429 threshold <- unname(attr(res, "filterThreshold"))
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
430 } else {
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
431 threshold <- 0
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
432 }
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
433 title_suffix <- paste0(primary_factor, ": ", lvl, " vs ", ref)
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
434 if (!is.null(opt$plots)) {
37
951876579237 "planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents: 36
diff changeset
435 generate_specific_plots(res, threshold, title_suffix)
17
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
436 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
437 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
438 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
439 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
440
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
441 # close the plot device
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
442 if (!is.null(opt$plots)) {
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
443 cat("closing plot device\n")
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
444 dev.off()
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
445 }
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
446
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
447 cat("Session information:\n\n")
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
448
24930eeff731 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff changeset
449 sessionInfo()