Answers to: indelを配列に反映させる

Comment by yag_ays on kenken's 回答

yag_ays — Sat, 26 Feb 2011 05:22:51 +0900

BioRubyを使わなくても良いように，スクリプトを書き換えました．単独で動くので，上記の使い方の項目1.は必要ありません．

使用例　：　https://skitch.com/yagays/rt2ma/indel

(ruby 1.8.7とruby 1.9.2p136で動作確認しています)

Answer by kenken

kenken — Thu, 24 Feb 2011 19:58:41 +0900

わざわざありがとうございます。感謝いたします。プログラムは初心者なのでよくわからないのですが、 biorubyは以下のホームページの通りにダウンロードしたらいいんでしょうか。過去に何回かcygwinで解析したことがあるのですが、今回は、いまいちよくわかりません。

http://dev.bioruby.org/ja/?install

Answer by yag_ays

yag_ays — Thu, 24 Feb 2011 09:52:43 +0900

Biostarに似たような質問があったので，リンクを貼っておきます．

http://biostar.stackexchange.com/questions/2529/how-to-insert-or-delete-nucleoties-at-various-locations-in-large-dna-file

Comment by yag_ays on kenken's 回答

yag_ays — Thu, 24 Feb 2011 04:11:40 +0900

取り敢えずは上のような感じで良さそうですね．スクリプトを少し書き換えたので，使い方を含めて補足します．

1．BioRubyというパッケージをインストールします

以下のコマンドを実行してください．

  $ gem install bio

(うまく行かない場合は行頭にsudoを付けて管理者権限になってください)

2．以下の3つのファイルを同じフォルダに入れてください

indel.rb

これは gist: 840281 をコピペしたファイル

indels.txt

indelの情報が入ったファイルです．タブ区切りのファイルで保存してください（もし違う形式だった場合はExcelの「別名で保存」からtsv形式で保存）

fasta.fasta

対象となるFASTAファイルです．

3．スクリプトを実行します

先ほどファイルを置いたフォルダに移動し，以下のコマンドを実行してください．

 $ ruby indel.rb indels.txt fasta.fasta

もし成功すれば，フォルダ内に「output.fasta」というファルが生成されているハズです．中身を見て，書き換えがうまくいっているか確かめてください．

以上です．

Answer by kenken

kenken — Thu, 24 Feb 2011 00:22:08 +0900

参考になるアドバイスありがとうございます。使い方と結果と全く同じことを cygwin上で行ってみました。エラーが出てしまうのですが、私のやり方の問題でしょうか。。。

$ ruby indel.rb indels.txt fasta.fasta
indel.rb:3:in `require': no such file to load -- bio (LoadError)
        from indel.rb:3

とのことです。

Answer by yag_ays

yag_ays — Wed, 23 Feb 2011 20:04:59 +0900

塩基配列を弄る系のソフトを簡単に見てみた限りだとkenkenさんの要求する操作ができるものがなかったので，適当にスクリプト書いてみました．こんなものが欲しいという理解でよろしいでしょうか．

ソースコード　：　https://gist.github.com/840281

使い方と結果　：　http://gyazo.com/d055680402948cd245eeaba2b60a373e.png

久しぶりにruby書いたので，ちょっとバグがあるかもしれません．~~動作保証もできないですが，要望があれば直します．~~

追記(2011/03/04)：取り敢えず解決された(?)のでしょうか．特に要望も無さそうなので，これで一端終わりにします．ソースコード等は自由に使ってもらって構いませんが，スクリプトの保守管理やbug fixは今後やりません．

Comment by mn3 on kenken's 質問

mn3 — Wed, 23 Feb 2011 18:55:02 +0900

真面目な話として、数百個分の操作なら、手でやるのも検討するといいかもしれません。

その操作が今回一回きりで、プログラミングの手間やソフトウェアの購入が難しいなら、決して悪くない方法だと思います。今後も同様の操作が繰り返される場合は、プログラミングすることを検討するべきだと思います。

もしくは、質問として求めている処理をより正確に記述して、これを見た方にプログラミングしてもらうという手があります。その場合は、正確に入力ファイルやパラメータや条件を定義し、出力として期待する型式を定義すると回答が得られる確率が向上します。たとえば、上記の質問では、1239000baseがなにを意味しているのか不明なので、期待されている回答が得られにくいかなとおもいます。