Mercurial > repos > iuc > deseq2
annotate deseq2.R @ 41:0a0a3388e3f2 draft
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
author | iuc |
---|---|
date | Mon, 29 Nov 2021 18:16:10 +0000 |
parents | b30a443484ab |
children | 6ef2cba4e35a |
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 |
37
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
34 options(show.error.messages = F, error = function() { |
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
35 cat(geterrmessage(), file = stderr()); q("no", 1, F) |
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
36 }) |
17
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
37 |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
38 # we need that to not crash galaxy with an UTF8 error on German LC settings. |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
39 loc <- Sys.setlocale("LC_MESSAGES", "en_US.UTF-8") |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
40 |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
41 library("getopt") |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
42 library("tools") |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
43 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
|
44 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
|
45 |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
46 # 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
|
47 # 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
|
48 spec <- matrix(c( |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
49 "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
|
50 "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
|
51 "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
|
52 "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
|
53 "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
|
54 "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
|
55 "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
|
56 "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
|
57 "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
|
58 "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
|
59 "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
|
60 "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
|
61 "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
|
62 "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
|
63 "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
|
64 "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
|
65 "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
|
66 "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
|
67 "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
|
68 "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
|
69 "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
|
70 "auto_mean_filter_off", "c", 0, "logical", |
39
b30a443484ab
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit c8e7020d39d581d7374a13cf94c64998a6481e05"
iuc
parents:
37
diff
changeset
|
71 "beta_prior_off", "d", 0, "logical", |
b30a443484ab
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit c8e7020d39d581d7374a13cf94c64998a6481e05"
iuc
parents:
37
diff
changeset
|
72 "alpha_ma", "A", 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) |
30
3b9af9a27b5b
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents:
29
diff
changeset
|
217 # estimate size factors for 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
|
218 if (!is.null(opt$esf)) { |
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
219 dds <- estimateSizeFactors(dds, type = opt$esf) |
30
3b9af9a27b5b
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents:
29
diff
changeset
|
220 } |
41
0a0a3388e3f2
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents:
39
diff
changeset
|
221 |
0a0a3388e3f2
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents:
39
diff
changeset
|
222 # 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
|
223 # - 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
|
224 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
|
225 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
|
226 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
|
227 ## 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
|
228 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
|
229 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
|
230 } else { |
0a0a3388e3f2
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents:
39
diff
changeset
|
231 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
|
232 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
|
233 ## 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
|
234 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
|
235 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
|
236 } else { |
0a0a3388e3f2
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents:
39
diff
changeset
|
237 ## 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
|
238 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
|
239 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
|
240 } |
0a0a3388e3f2
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents:
39
diff
changeset
|
241 } |
0a0a3388e3f2
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents:
39
diff
changeset
|
242 write.table(size_factors, file = opt$sizefactorsfile, sep = "\t", col.names = F, quote = FALSE) |
0a0a3388e3f2
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents:
39
diff
changeset
|
243 } |
0a0a3388e3f2
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit cbeb1c4c436be04323bd9a809a6393d00b168d07"
iuc
parents:
39
diff
changeset
|
244 |
37
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
245 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
|
246 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
|
247 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
|
248 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
|
249 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
|
250 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
|
251 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
|
252 } |
66203279224f
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents:
28
diff
changeset
|
253 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
|
254 # 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
|
255 # Set sort to False, which maintains the order in dds_data |
37
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
256 reordered_batch <- merge(dds_data, batch_factors, by.x = 0, by.y = 0, sort = F) |
29
66203279224f
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents:
28
diff
changeset
|
257 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
|
258 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
|
259 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
|
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 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
|
262 return(dds) |
66203279224f
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents:
28
diff
changeset
|
263 } |
27
dc6bc19f05ab
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents:
26
diff
changeset
|
264 |
29
66203279224f
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents:
28
diff
changeset
|
265 if (!is.null(opt$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
|
266 batch_factors <- read.table(opt$batch_factors, sep = "\t", header = T) |
29
66203279224f
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit b95582cea8320d5488056a9576474f79cec53be8
iuc
parents:
28
diff
changeset
|
267 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
|
268 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
|
269 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
|
270 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
|
271 } |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
272 |
27
dc6bc19f05ab
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents:
26
diff
changeset
|
273 if (verbose) { |
dc6bc19f05ab
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents:
26
diff
changeset
|
274 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
|
275 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
|
276 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
|
277 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
|
278 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
|
279 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
|
280 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
|
281 } |
dc6bc19f05ab
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 83eb5b2665d87c02b270596f8175499e69061032
iuc
parents:
26
diff
changeset
|
282 |
17
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
283 # optional outlier behavior |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
284 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
|
285 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
|
286 } else { |
37
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
287 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
|
288 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
|
289 } |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
290 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
|
291 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
|
292 } else { |
37
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
293 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
|
294 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
|
295 } |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
296 |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
297 # 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
|
298 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
|
299 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
|
300 } else { |
37
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
301 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
|
302 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
|
303 } |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
304 |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
305 # shrinkage of LFCs |
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$beta_prior_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 beta_prior <- TRUE |
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 beta_prior <- FALSE |
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("beta prior 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 |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
313 # dispersion fit type |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
314 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
|
315 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
|
316 } else { |
37
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
317 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
|
318 } |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
319 |
37
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
320 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
|
321 |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
322 # 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
|
323 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
|
324 |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
325 # 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
|
326 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
|
327 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
|
328 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
|
329 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
|
330 } |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
331 |
37
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
332 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
|
333 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
|
334 |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
335 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
|
336 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
|
337 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
|
338 } |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
339 |
30
3b9af9a27b5b
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents:
29
diff
changeset
|
340 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
|
341 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
|
342 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
|
343 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
|
344 } |
3b9af9a27b5b
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents:
29
diff
changeset
|
345 |
3b9af9a27b5b
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents:
29
diff
changeset
|
346 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
|
347 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
|
348 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
|
349 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
|
350 } |
3b9af9a27b5b
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents:
29
diff
changeset
|
351 |
3b9af9a27b5b
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 448dccb0c02aba00d8301247b0f0f406ab3d4fa2
iuc
parents:
29
diff
changeset
|
352 |
17
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
353 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
|
354 # 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
|
355 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
|
356 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
|
357 res <- results( |
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
358 dds, |
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
359 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
|
360 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
|
361 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
|
362 ) |
17
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
363 if (verbose) { |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
364 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
|
365 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
|
366 print(summary(res)) |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
367 } |
37
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
368 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
|
369 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
|
370 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
|
371 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
|
372 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
|
373 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
|
374 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
|
375 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
|
376 } else { |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
377 threshold <- 0 |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
378 } |
37
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
379 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
|
380 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
|
381 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
|
382 } |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
383 } else { |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
384 # 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
|
385 # 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
|
386 # 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
|
387 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
|
388 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
|
389 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
|
390 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
|
391 res <- results( |
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
392 dds, |
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
393 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
|
394 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
|
395 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
|
396 ) |
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
397 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
|
398 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
|
399 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
|
400 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
|
401 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
|
402 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
|
403 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
|
404 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
|
405 } else { |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
406 threshold <- 0 |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
407 } |
37
951876579237
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 9eb6d07510ccf27d6499172d62c81661078ec57b"
iuc
parents:
36
diff
changeset
|
408 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
|
409 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
|
410 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
|
411 } |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
412 } |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
413 } |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
414 } |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
415 |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
416 # 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
|
417 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
|
418 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
|
419 dev.off() |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
420 } |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
421 |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
422 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
|
423 |
24930eeff731
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/deseq2 commit 23559e873291b04840f186edb0b4b85e76359e05
iuc
parents:
diff
changeset
|
424 sessionInfo() |