[BioRuby-ja] bioruby の.naseq メソッド

KAWASHIMA Shuichi shuichi @ kuicr.kyoto-u.ac.jp
2004年 2月 3日 (火) 05:31:14 EST


川島@京大化研です。

At Tue, 03 Feb 2004 17:49:24 +0900,
takeshi kawashima wrote:
> しかし、マルチファスタを読み込んで同じことをしようと思い、
> 
> m = Bio::FlatFile.new(Bio::FastaFormat, ARGF)
> m.each_entry do |x|
>   i = 1
>   remainder = x.naseq.window_search(60, 60) do |subseq|
>     puts subseq.to_fasta("segment #{i}", 60)
>     i += 1
>   end
>   puts remainder.to_fasta("segment #{i}", 60)      #####<=ここ
> end
> 
> とすると、ループの最後の putsの行でエラーがでてしまいます。

試しましたが、動いているようですが?

> remainderがStringクラスだからのようですが、
> うまい回避方法はないでしょうか?
> 今は、
> tmpseq = Bio::Sequence.new("#{remainder}")
> puts tmpseq.to_fasta("segment #{i}", 60)     
> としてます。

NA クラスは、String クラスを継承しているので String クラスで実装されて
いる to_fasta は、NA クラスでも使えるはずですが。

何か別の問題では?

川島



BioRuby-ja メーリングリストの案内