rbbt-marq 1.0.9 → 1.1.0
Sign up to get free protection for your applications and to get access to all the features.
- data/R/GEO.R +68 -124
- data/R/GEOquery_patch.R +44 -0
- data/R/MARQ.R +76 -0
- data/install_scripts/GEO/Rakefile +111 -74
- data/lib/MARQ/GEO.rb +459 -451
- data/lib/MARQ/ID.rb +1 -0
- data/lib/MARQ/main.rb +82 -2
- metadata +13 -2
data/lib/MARQ/ID.rb
CHANGED
data/lib/MARQ/main.rb
CHANGED
@@ -22,9 +22,9 @@ module MARQ
|
|
22
22
|
def self.platform_organism(platform)
|
23
23
|
if platform_type(platform) == :GEO
|
24
24
|
if platform.match(/^GPL/)
|
25
|
-
GEO::
|
25
|
+
GEO::SOFT::GPL(platform)[:organism]
|
26
26
|
else
|
27
|
-
GEO::
|
27
|
+
GEO::SOFT::GPL(GEO::dataset_platform(platform))[:organism]
|
28
28
|
end
|
29
29
|
else
|
30
30
|
CustomDS::organism(platform)
|
@@ -50,6 +50,7 @@ module MARQ
|
|
50
50
|
CustomDS::is_cross_platform?(platform)
|
51
51
|
end
|
52
52
|
end
|
53
|
+
|
53
54
|
def self.complete_positions(positions, matched, genes)
|
54
55
|
|
55
56
|
pos = Hash[*matched.zip(positions).flatten]
|
@@ -62,6 +63,83 @@ module MARQ
|
|
62
63
|
end
|
63
64
|
}
|
64
65
|
|
66
|
+
end
|
67
|
+
|
68
|
+
module Dataset
|
69
|
+
|
70
|
+
def self.exists?(dataset)
|
71
|
+
MARQ::dataset_path(dataset) != nil
|
72
|
+
end
|
73
|
+
|
74
|
+
def self.read_file(dataset, ext)
|
75
|
+
Open.read(MARQ::dataset_path(dataset) + '.' + ext)
|
76
|
+
end
|
77
|
+
|
78
|
+
def self.read_values(dataset, file)
|
79
|
+
result = {}
|
80
|
+
|
81
|
+
experiments = experiments(dataset)
|
82
|
+
experiments.each{|experiment| result[experiment] = [] }
|
83
|
+
read_file(dataset, file).split(/\n/).each do |line|
|
84
|
+
values = line.chomp.split(/\t/)
|
85
|
+
values.each_with_index{|value, i| result[experiments[i]] << (value == 'NA' ? nil : value.to_f) }
|
86
|
+
end
|
87
|
+
|
88
|
+
result
|
89
|
+
end
|
90
|
+
|
91
|
+
def self.read_values_t(dataset, file)
|
92
|
+
result = {}
|
93
|
+
|
94
|
+
experiments = experiments(dataset).select{|experiment| experiment !~ /\[ratio\]$/ }
|
95
|
+
return {} if experiments.empty?
|
96
|
+
experiments.each{|experiment| result[experiment] = [] }
|
97
|
+
|
98
|
+
read_file(dataset, file).split(/\n/).each do |line|
|
99
|
+
values = line.chomp.split(/\t/)
|
100
|
+
values.each_with_index{|value, i| result[experiments[i]] << (value == 'NA' ? nil : value.to_f) }
|
101
|
+
end
|
102
|
+
|
103
|
+
result
|
104
|
+
end
|
105
|
+
|
106
|
+
|
107
|
+
def self.platform_codes(platform)
|
108
|
+
if MARQ::is_cross_platform? platform
|
109
|
+
file = 'cross_platform'
|
110
|
+
else
|
111
|
+
file = 'codes'
|
112
|
+
end
|
113
|
+
|
114
|
+
Open.read(File.join(MARQ::platform_path(platform), file))
|
115
|
+
end
|
116
|
+
|
117
|
+
def self.experiments(dataset)
|
118
|
+
read_file(dataset, 'experiments').split(/\n/)
|
119
|
+
end
|
120
|
+
|
121
|
+
def self.codes(dataset)
|
122
|
+
read_file(dataset, 'codes').split(/\n/)
|
123
|
+
end
|
124
|
+
|
125
|
+
def self.orders(dataset)
|
126
|
+
read_values(dataset, 'orders')
|
127
|
+
end
|
128
|
+
|
129
|
+
def self.logratios(dataset)
|
130
|
+
read_values(dataset, 'logratios')
|
131
|
+
end
|
132
|
+
|
133
|
+
def self.pvalues(dataset)
|
134
|
+
read_values_t(dataset, 'pvalues')
|
135
|
+
end
|
136
|
+
|
137
|
+
def self.t(dataset)
|
138
|
+
read_values_t(dataset, 't')
|
139
|
+
end
|
140
|
+
|
141
|
+
|
142
|
+
|
65
143
|
end
|
66
144
|
|
67
145
|
def self.platform_scores_up_down(platform, up, down)
|
@@ -189,6 +267,8 @@ module MARQ
|
|
189
267
|
end
|
190
268
|
|
191
269
|
if __FILE__ == $0
|
270
|
+
p MARQ::Dataset.orders('GDS2419')
|
271
|
+
exit
|
192
272
|
#puts MARQ::organism_platforms('human')
|
193
273
|
#puts MARQ.platform_organism("HaploidData")
|
194
274
|
#puts MARQ::platform_scores_up_down("HaploidData",%w( YMR261c YDL140c YIL122w YPL093w YHR211w YDL142c YHR106w YOR103c YDR233c YLR181c),%w()).keys
|
metadata
CHANGED
@@ -1,7 +1,7 @@
|
|
1
1
|
--- !ruby/object:Gem::Specification
|
2
2
|
name: rbbt-marq
|
3
3
|
version: !ruby/object:Gem::Version
|
4
|
-
version: 1.0
|
4
|
+
version: 1.1.0
|
5
5
|
platform: ruby
|
6
6
|
authors:
|
7
7
|
- Miguel Vazquez
|
@@ -9,7 +9,7 @@ autorequire:
|
|
9
9
|
bindir: bin
|
10
10
|
cert_chain: []
|
11
11
|
|
12
|
-
date: 2009-12-
|
12
|
+
date: 2009-12-16 00:00:00 +01:00
|
13
13
|
default_executable: marq_config
|
14
14
|
dependencies:
|
15
15
|
- !ruby/object:Gem::Dependency
|
@@ -62,6 +62,16 @@ dependencies:
|
|
62
62
|
- !ruby/object:Gem::Version
|
63
63
|
version: "0"
|
64
64
|
version:
|
65
|
+
- !ruby/object:Gem::Dependency
|
66
|
+
name: png
|
67
|
+
type: :runtime
|
68
|
+
version_requirement:
|
69
|
+
version_requirements: !ruby/object:Gem::Requirement
|
70
|
+
requirements:
|
71
|
+
- - ">="
|
72
|
+
- !ruby/object:Gem::Version
|
73
|
+
version: "0"
|
74
|
+
version:
|
65
75
|
description: Find microarray experiments with similar or opposite signature to a given query. A SOAP interface and a merb portal can be found in rbbt-marq-www.
|
66
76
|
email: miguel.vazquez@fdi.ucm.es
|
67
77
|
executables:
|
@@ -74,6 +84,7 @@ extra_rdoc_files:
|
|
74
84
|
files:
|
75
85
|
- R/CustomDS.R
|
76
86
|
- R/GEO.R
|
87
|
+
- R/GEOquery_patch.R
|
77
88
|
- R/MA.R
|
78
89
|
- R/MARQ.R
|
79
90
|
- bin/marq_config
|