ゴミのやま

完全自分向け雑記

NCBIからFASTAを取ってくる

以前UniprotからFastaを取ってくる方法の記事を書いた。UniprotではURLを適宜いじることでFastaを直接ダウンロードできたが, NCBIで同様の方法は使えない(ようだ)。 具体的にはBiopythonのEntrezを使う。詳細は以下参照。 https://www.biostars.org/p/1771…

Uniprotから任意のFASTAをとってくる

やりたいこと Uniprotから任意のアミノ酸配列をFASTA形式でPythonなどのプログラミング言語を使うことにより自動取得したい。 参考 以下のサイトでまずUniprotとは何かを理解した。アミノ酸配列とその機能情報を掲載しているらしい。 magattaca.hatenablog.c…

HHsuite覚書

目的 仕事でホモロジーモデリングをする必要が出てきた。それにはhhsuiteなるものが必要らしい。よくわからないので入力/出力および特性を理解する。アルゴの詳細には深入りしない。 同じようなソフトウェアは他にも有るはずなので, それらのソフトウェアと…

Linux mint 導入時にやったことメモ

Ubuntuくんが激重なのでLinux mint (Cinnamon)に乗り換えたのメモ。完全自分用 とりあえず結論 めっちゃ快適なのでさっさとやればよかった とりあえずアップデート インストール後表示されるFirst Stepsに従ってドライバをインストールしたりアップデートし…

LGAのコンパイルでSegmentation faultが出る

タンパク質のアラインメント計算ソフトウェアのLGAでコンパイルはできるのだが, Segmentation faultで実行できない問題の対応法。 LGA - Protein Structure Comparison Facility [解決策]ulimit -s unlimitedする どうもコアダンプの容量が制限されている…

Pymol 見た目を調整する

気付き次第, 順次追加予定。。。 DSSコマンド 水素結合および主鎖の構造を用いて二次構造情報を自動的に追加する。これにより二次構造情報に基づいた色分けなどが可能になる。 https://pymolwiki.org/index.php/Dss デフォルトの設定を変更する(背景など) …

Poisson-Boltzmann方程式を導出する

目的 分子シミュレーションでよく使われるPossion-Boltzmann方程式を導出する。間違ってたら教えて下さい。 出発 まずPossion equationから出発する :点rにおける誘電率; :点rにおける静電ポテンシャル; :点rにおける電荷密度 この式はある点における電荷分…

PythonのPandasで一列ずつ足す

やりたいこと a = [] for i in range(0,n): a.append("hogehoge") みたいなことをpandasでやりたい 解決策 一旦python標準配列で作ってからpd.DataFrameに変換する方法が早いらしい qiita.com

SMARTSってなんやねん(化学)

参考: magattaca.hatenablog.com www.molsis.co.jp 要は SMILES 分子の二次元構造を文字列で記述 ワイルドカード表記はできない(多分…?) 一つの化合物構造に対して複数のsmilesがありえる(一対一対応するものはcanonical SMILESと呼ばれる) SMARTS SMI…

sshfsでマウントしたフォルダが"You do not have the permissions necessary to view the contents of mounted directory"と出て開けない

Mac OSX MojaveでSshfsでマウントしたフォルダがパーミッションエラーで開けない You do not have the permissions necessary to view the contents of mounted directory. というメッセージによりsshfsが開けなくなったのでその対処法 対策 -o defer_permis…

Jupyter Notebookの導入メモ

pipでインストール qiita.com 拡張機能を入れる www.pynote.inf とりあえず上のサイトから適当につっこむ。vimmerなので"jupyter-vim-binding"は必須。 リモートのスクリプトをいじる sshfsでマウントしていじる方法。この場合, スクリプトはローカルで実行…

t検定ってなんやねん

まず仮設をたてる 参考: best-biostatistics.com 帰無仮説(Null hypothesis):差がない 対立仮説(Alternative hypothesis):差がある 帰無仮説の否定:差がないわけではない⇒差がある t検定する 帰無仮説を立てる(差がないのではという仮設) p値を算出し…

Jupyter Notebookってなんやねん

目的 同僚がJupyter Notebookを使っていたのだが何が便利なのかわからないので全体像を掴む 参考: ここを参考にしました blog.tsurubee.tech 要は Pythonスクリプト内部にLatexやmatplotlibを用いた色々な形式でメモができたり それをチームとシェアできた…

Pymol 便利なコマンド集

目的 自分用によく使うコマンドを整理しておく ファイルのリストを得る cmd.get_names() ファイルの名前を書きかえる cmd.set_name("fileA", "fileB")) [応用]"filename A"を含むファイル名のみ抽出し, 一部を書き換える for a in [file for file in cmd.get…

線形代数はなぜ便利なのか?

ここを参考にしました qiita.com なぜ便利か? 時間発展の運動方程式が解けるnstepにおける任意のベクトルに時間発展行列をかけてやることでn+1stepの結果が得られる特にベクトルを時間発展行列の固有ベクトルの線形結合で示してやることにより, 時間発展行…

線形代数を再学習する

背景 応用数学分野において行列計算のために線形代数の知識は必須。しかし著者は大学の基礎教養でやって以来ご無沙汰なので基礎計算はできるものの概念は未だによくわかっていない。 目的 線形代数の大枠の理解(どのように便利なのか) 線形空間・固有値分…

このブログについて

目的 アウトプットすることで理解を深める アウトプットすることで間違いを(望ましくはやさしく)正してもらう 基本的にゴミの山と思っているが, アウトプットしておくことで皆様に有用なものを適宜見つけて貰えるかもしれない お題目(頻度順…のはず) 応…