[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