[BioRuby] Bio.pdb doubt ...
K. Shameer
shameer at ncbs.res.in
Wed Feb 20 18:01:31 UTC 2008
Hi Alex,
I used this code to check chain identifiers. But I noticed the code is
giving true for multichain and single chain entries.
Can you please have a look at the code,
require 'bio'
module Bio
class PDB
def multichain?
self.chains.size > 1
end
end
end
inputpdb = ARGV[0]
if ARGV[0] == nil
puts "please give an argument"
end
test = Bio::PDB.new(IO.read(inputpdb)).multichain?
print test,"\n"
I run this code using chain.rb test/1tca.pdb -> out put was 'true'(this
shouldbe false - there is only one chain in 1TCA) and out put was 'true'
for 1a7b as well.
Thanks in advance,
Shameer
PS. Now I am going to try Goto's code.
More information about the BioRuby
mailing list