summaryrefslogtreecommitdiff
path: root/scripts/pa-fingerprints.rb
diff options
context:
space:
mode:
authorChristoph Helma <helma@in-silico.ch>2021-02-06 20:21:58 +0100
committerChristoph Helma <helma@in-silico.ch>2021-02-06 20:21:58 +0100
commit771a2381ae0fd5e352f23d7223baeb26e8bb4e02 (patch)
tree1a7182a83ec15b9d7a69e5508c8d3b2bfa5983df /scripts/pa-fingerprints.rb
parent241f997c3a1a6a38fa47070f5efbd23852cc432b (diff)
svm validation, tensorflow mp2d pa predictions, r results removed, cleanup
Diffstat (limited to 'scripts/pa-fingerprints.rb')
-rwxr-xr-xscripts/pa-fingerprints.rb20
1 files changed, 20 insertions, 0 deletions
diff --git a/scripts/pa-fingerprints.rb b/scripts/pa-fingerprints.rb
new file mode 100755
index 0000000..344ba86
--- /dev/null
+++ b/scripts/pa-fingerprints.rb
@@ -0,0 +1,20 @@
+#!/usr/bin/env ruby
+require_relative "../../lazar/lib/lazar.rb"
+training_fingerprints = `sed -n '1p' data/mutagenicity-fingerprints.csv`.chomp.split(",")
+training_fingerprints.pop
+puts training_fingerprints.join(",")
+training_fingerprints.shift
+
+File.readlines(File.join("pyrrolizidine-alkaloids","lazar","pa-smiles.csv")).each_with_index do |line,i|
+ if i > 0
+ (id,smiles) = line.chomp.split(",")
+ c = Compound.new(smiles)
+ out = [c.smiles]
+ fp = c.fingerprint
+ training_fingerprints.each do |frag|
+ fp.include?(frag) ? out << 1 : out << 0
+ end
+ puts out.join(",")
+ end
+end
+