プログラミングPythonでファイルの読み込みや書き込みをする際にはopen関数を使うのが基本ですが、下に挙げるような少し複雑な操作をしたい場合には頭を抱えることになります。標準入力からも読み込みたい、ただし指定したときにはファイルからも読み込んでほしい。Pythonオーソログ推定プログラムのOrthoFinderですが、便利なことにオーソログ推定の結果予測されたシングルコピー遺伝子の情報を個別に抽出してディレクトリにまとめてくれます。プログラミングPythonスクリプト内で外部プログラムを呼び出して解析を実行することが多々あるのですが、プログラムの中には出力先を指定するオプションがなくカレントディレクトリに結果を出力することを想定しているものがあります。LTSにしてから、exitで端末を閉じようとするとフリーズするようになってしまいました。プログラミングpythonでディレクトリ(フォルダ)下のファイルを一括で変更して通し番号を振るスクリプトを作ったので備忘録を兼ねてまとめます。バイオインフォマティクス複数のNewick形式の系統樹から樹形が同じものを抽出する必要に駆られたので備忘録を兼ねたまとめです。デスクトップ環境この度VimからNeovimに乗り換えたので、備忘録を兼ねてそのインストール手順をまとめます。雑記どうにもパソコンの動作がモッサリとしていたので軽くする方法を調べてみました。した試みは徒労に終わる傾向が強いような印象があったのですが、珍しく満足するレベルで高速化できましたので備忘録を兼ねてまとめます。デスクトップ環境rangerはターミナルから起動できるvimライクなキーバインドで操作できるCUIのファイラーです。ターミナル内で使えるファイラーは他にも試したのですが操作が良く分からず、デフォルトの状態でも十分使いやすいrangerに落ち着きました。バイオインフォマティクスBiopythonを使って、与えられたFASTAファイルの配列データをクエリーとしてオンラインのBLASTを実行するスクリプトを作成したので、備忘録を兼ねてまとめます。公式のチュートリアルを参考に作りましたので詳しくはそちらをご覧ください。バイオインフォマティクス解析に使うプログラムの中には、FASTA形式の配列ファイルではなくPhylip形式のものを要求するものがあります。今回はそのような事情でFASTA形式のデータをPhylip形式に変換する、ちょっとしたスクリプトを書いたので備忘録を兼ねてまとめます。Python入門書で読んだきり忘れていた問題に、この度出くわしたので備忘録を兼ねてまとめておこうと思います。以下のようなPythonのコードを書いていました。バイオインフォマティクス与えられた配列ファイルをプログラムで処理する際に、ファイルのフォーマットを調べておきたいという場合がたまにあります。プログラミング新しい環境に変わるたびに必要なパッケージなどをインストール作業が段々と苦痛になってきたので、シェルスクリプト(bash)で自動化することにしました。デスクトップ環境タイル型ウィンドウマネージャのi3-wmはデフォルトの設定でもある程度快適に使えるのが魅力ですが、カスタマイズするとより使いやすくなるのも事実です。Python依存関係を手作業で解消するのが面倒だったのでPythonにやってもらうことしました。「退屈なことはPythonにやらせよう」ですね。環境はUbuntuを想定しています。普段はVirtualbox上で仮想マシンのUbuntuを使用しているのですが、何かともたつくことも多いため、できればWindows上で操作を完結させたいのが本音です。Windowsを完全に捨ててしまうのも1つの手なのでしょうが、なかなか踏ん切りがつきません。バイオインフォマティクス他人に説明するときなどに、系統樹のノードに番号が付いていると便利な場合がたまにあります。AとBの共通祖先等の表現でも良いんですが長くなりやすいので。Toolkitを使って系統樹のノードに通し番号を付けたのでそのまとめです。Linux)は入門レベルの人間が、WSLのUbuntu(bash)上で可能な限り操作を完結させようともがいた結果をまとめたメモです。初心者なので色々大目に見ていただけるとありがたいです。バイオインフォマティクスプログラム中で、NCBIの管理するデータベースに登録された配列ファイルをダウンロードしたいことがたまにあります。手作業は何かと煩雑なので。PythonPythonのmultiprocessingモジュールを試そうとしたらエラーが発生してしまったという話。問題のコードは下の通りです。PCPythonのスクリプト中で、ディレクトリ(フォルダ)が空かどうか確認してから処理を行いたいことがたまにあります。そこで、簡単に調べたところosモジュールのlistdir関数を使う方法と、globモジュールのglob関数を利用する方法があるようでした。バイオインフォマティクスFastTreeは近似最尤法による系統樹推定プログラムです。大規模データに対応しており、他のプログラムに比べ実行時間が非常にはやいです。公式の説明によれば、PhyMLやRAxMLの100から1000倍高速な場合もあるようです。PCオーソログ推定プログラム「OrthoFinder」は、推定したオーソロググループ配列を使用してそのまま系統樹作成まで行うことができます。詳しくはGithubのページに詳しい説明がありますのでご覧ください。バイオインフォマティクスIQ-Treeは最尤法による系統樹作成プログラムです。コマンドラインから簡単に実行することができ、計算時間もはやいため便利です。配列をパーティションに区切った解析を行うことができる他、塩基・アミノ酸配列以外にも、幾つかの形式の配列データに...IQ-Treeの使い方まとめ2018-10-19JupyterNotebookにIRubyをインストールしたハナシPCJupyterNotebook上でRubyを利用するために、IRubyをインストールしたまとめです。一点だけつまづいたので記事に起こしておきます。16.04以前の方や他のOSの方は手順が若干異なるためご注意ください。プログラミングプログラムで自動で配列を集めてきたり、トリミングなどの処理を行ったりした場合に、予期せず全く同じ配列ファイルが違う名前で紛れてしまうことがあります。私のやり方がまずいのかもしれませんが。プログラミングFASTAファイルの配列が途中で改行され、複数行に渡っている場合があります。通常は問題ないのですがたまに不便なときがあるので、Biopythonを使って不要な改行を削除し配列行を一行にします。完成形のコードは以下のとおりです。デスクトップ環境タイル型ウィンドウマネージャの「i3」と一緒に使っているアプリのまとめです。バイオインフォマティクスGenbankファイル中のCDS情報から遺伝子のアミノ酸配列を抽出し、遺伝子ごとのアミノ酸配列が記載されたFASTAファイルを作ったまとめです。Genbankファイルの扱いはBiopythonを利用すると簡単です。公式のチュートリアルに詳しい説明がありますが、英語だっ...BiopythonでGenbankファイルからアミノ酸配列を抽出する次のページプロフィール取り留めのないブログですがお役に立てれば幸いです。(7)月別アーカイブスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。