Mercurial > repos > artbio > gsc_gene_expression_correlations
annotate correlation_with_signature.R @ 2:b49295546f29 draft default tip
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
author | artbio |
---|---|
date | Wed, 18 Oct 2023 10:00:34 +0000 |
parents | 8ad272e0b640 |
children |
rev | line source |
---|---|
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
1 # Performs multi-correlation analysis between the vectors of gene expressions |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
2 # in single cell RNAseq libraries and the vectors of signature scores in these |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
3 # same single cell RNAseq libraries. |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
4 # Example of command |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
5 # Rscript correlations_with_signature.R --expression_file <expression_data.tsv> |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
6 # --signatures_file <signature_scores.tsv> |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
7 # --sep "\t" |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
8 # --colnames "T" |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
9 # --gene_corr <gene-gene corr file> |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
10 # --gene_corr_pval <gene-gene corr pvalues file> |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
11 # --sig_corr <genes correlation to signature file> |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
12 |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
13 options(show.error.messages = FALSE, |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
14 error = function() { |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
15 cat(geterrmessage(), file = stderr()) |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
16 q("no", 1, FALSE) |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
17 } |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
18 ) |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
19 loc <- Sys.setlocale("LC_MESSAGES", "en_US.UTF-8") |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
20 |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
21 library(optparse) |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
22 library(Hmisc) |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
23 |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
24 # Arguments |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
25 option_list <- list( |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
26 make_option( |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
27 "--sep", |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
28 default = "\t", |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
29 type = "character", |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
30 help = "File separator, must be the same for all input files [default : '%default' ]" |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
31 ), |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
32 make_option( |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
33 "--colnames", |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
34 default = TRUE, |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
35 type = "logical", |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
36 help = "Consider first lines as header (must stand for all input files) [default : '%default' ]" |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
37 ), |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
38 make_option( |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
39 "--expression_file", |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
40 default = NA, |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
41 type = "character", |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
42 help = "Input file that contains log2(CPM +1) expression values" |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
43 ), |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
44 make_option( |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
45 "--signatures_file", |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
46 default = NA, |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
47 type = "character", |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
48 help = "Input file that contains cell signature" |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
49 ), |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
50 make_option( |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
51 "--sig_corr", |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
52 default = "sig_corr.tsv", |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
53 type = "character", |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
54 help = "signature correlations output [default : '%default' ]" |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
55 ), |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
56 make_option( |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
57 "--gene_corr", |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
58 default = "gene_corr.tsv", |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
59 type = "character", |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
60 help = "genes-genes correlations output [default : '%default' ]" |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
61 ), |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
62 make_option( |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
63 "--gene_corr_pval", |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
64 default = "gene_corr_pval.tsv", |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
65 type = "character", |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
66 help = "genes-genes correlations pvalues output [default : '%default' ]" |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
67 ) |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
68 ) |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
69 |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
70 opt <- parse_args(OptionParser(option_list = option_list), |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
71 args = commandArgs(trailingOnly = TRUE)) |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
72 |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
73 if (opt$sep == "tab") { |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
74 opt$sep <- "\t" |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
75 } |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
76 if (opt$sep == "comma") { |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
77 opt$sep <- "," |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
78 } |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
79 |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
80 # Open files |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
81 data <- read.delim( |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
82 opt$expression_file, |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
83 header = opt$colnames, |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
84 row.names = 1, |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
85 sep = opt$sep, |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
86 check.names = FALSE |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
87 ) |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
88 signature <- read.delim( |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
89 opt$signatures_file, |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
90 header = TRUE, |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
91 stringsAsFactors = FALSE, |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
92 row.names = 1, |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
93 sep = opt$sep, |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
94 check.names = FALSE |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
95 ) |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
96 |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
97 |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
98 # keep only signatures that are in the expression dataframe |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
99 signature <- subset(signature, rownames(signature) %in% colnames(data)) |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
100 |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
101 # Add signature score to expression matrix |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
102 data <- rbind(t(signature), data) |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
103 |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
104 # Gene correlation |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
105 gene_corr <- rcorr(t(data), type = "pearson") # transpose because we correlate genes, not cells |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
106 |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
107 # Gene correlation with signature score |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
108 gene_signature_corr <- cbind.data.frame(gene = colnames(gene_corr$r), |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
109 Pearson_correlation = gene_corr$r[, 1], |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
110 p_value = gene_corr$P[, 1]) |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
111 gene_signature_corr <- gene_signature_corr[order(gene_signature_corr[, 2], decreasing = TRUE), ] |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
112 |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
113 |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
114 ### Save files ### |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
115 |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
116 write.table( |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
117 format(gene_signature_corr, digits = 2), |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
118 file = opt$sig_corr, |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
119 sep = "\t", |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
120 quote = FALSE, |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
121 col.names = TRUE, |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
122 row.names = FALSE |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
123 ) |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
124 |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
125 r_genes <- data.frame(gene = rownames(gene_corr$r), gene_corr$r) # add rownames as a variable for output |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
126 write.table( |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
127 format(r_genes[-1, -2], digits = 2), |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
128 file = opt$gene_corr, |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
129 sep = "\t", |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
130 quote = FALSE, |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
131 col.names = TRUE, |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
132 row.names = FALSE |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
133 ) |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
134 |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
135 p_genes <- data.frame(gene = rownames(gene_corr$P), gene_corr$P) # add rownames as a variable for output |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
136 write.table( |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
137 format(p_genes[-1, -2], digits = 2), |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
138 file = opt$gene_corr_pval, |
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
139 sep = "\t", |
2
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
140 quote = FALSE, |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
141 col.names = TRUE, |
b49295546f29
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 91d59a3a90a9bdb64ec70000b69a864285411d9a
artbio
parents:
1
diff
changeset
|
142 row.names = FALSE |
0
734ab9c5595a
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/gsc_gene_expression_correlations commit 09dcd74dbc01f448518cf3db3e646afb0675a6fe
artbio
parents:
diff
changeset
|
143 ) |