Hi-Roto's Note

博士(理学)筑波大学。職業:バイオ系研究者(糖鎖関連疾患の研究に従事)。 研究や実験に…

Hi-Roto's Note

博士(理学)筑波大学。職業:バイオ系研究者(糖鎖関連疾患の研究に従事)。 研究や実験に関するメモ,pythonによるデータ処理や解析についてのメモを記録する

記事一覧

Pythonを使ったプロテオミクスデータのVolcano plot例

はじめにコントロール細胞(WT細胞)と,ある遺伝子をノックアウトした細胞(KO細胞)に発現するタンパク質群について,増減を比較し,その増減が統計的に有意なのかを,ボルケ…

2

Scipyのcurve_fitを使用した非線形曲線Y=ax/x+bへの近似方法

はじめに 生化学でよく出てくる酵素反応の経時変化を,下記の式への近似する。 まずはモジュールのimport # import moduleimport numpy as npimport matplotlib.pyplot …

ウエスタンブロットの画像トリミング方法

はじめにウエスタンブロットを含めた画像データは画像編集ソフトでトリミングを行なった後に,Adbe IllustratorやInkscapeなどのベクター系の描画ツールで作図を行うことが…

外部ファイルのクラスを読み込む方法

ipython上などで下記プログラムを実行する場合 class Hello(): def __init__(self, name): self.name = name def sayhello(self): print ('Hello {} !!!'…

1

Numpy memo (随時更新)

Numpyで今まで使ったものを記録(随時更新) 一次元配列のarrayを二次元配列にして指定の数に分割する:np.array_splita = np.array([1,2,3,4,5,6,7,8,9])を二次元に変換後,…

1

Pandasを使って複数の独立した実験データから平均値と標準偏差を算出する方法

三回の独立した実験(exp_number, 1-3) で,サンプルa-cの酵素活性(relateive_activity) を算出したデータ例を使用し,平均値(mean)と標準偏差(std)を計算する。  自作関…

6

Pandas でよく忘れてしまう事柄(随時追加)

Pythonを使ったデータ解析をする際によく忘れてgoogle検索しがちなものを下記に書き留めておく(随時追加予定) Excelファイルのimport(read_excel)'xlrd'をインストール…

細胞内に発現する3つのタンパク質の相対量変化を例にした折れ線グラフの描画 (Matplotlib, Seaborn)

3種類のタンパク質 Protein1, Protein2, Protein3について,あるタイムポイント(0分)から始めて90分まで,30分おきに細胞内発現量を経時的に定量・算出したデータを想定し…

酵素活性データを例にした棒グラフの描画その2(Matplotlib, Seaborn)

前回に引き続き今回は酵素1-3の酵素活性について,阻害剤のある場合ない場合に分け,3回繰り返し測定したデータ例を用意。 平均値,標準偏差を計算後,matplotlib, seaborn…

酵素活性データを例にした棒グラフの描画(Matplotlib, Seaborn)

バイオ実験で出てきそうな酵素活性を測定した実験データをmatplotlibとseabornで棒グラフにする例を書き留める。 (データ例の概要)酵素1-3の酵素活性を3回繰り返し測定…

4

matplotlibでsubplot

subplotの方法について覚えたことを記録。 suplotは複数のグラフを一つの図として表示する方法。ネットで検索しながら何となくコピペして作っていていたのだが,どうやらや…

1

Matplotlib(Seaborn)の基本的な使い方

最近,研究活動(おもに実験)で得たデータを計算したり,可視化するのに使っているPython。 可視化(グラフ作成)に使うMatplotlib(Seaborn)の使い方をすぐ忘れてしまうの…

21
Pythonを使ったプロテオミクスデータのVolcano plot例

Pythonを使ったプロテオミクスデータのVolcano plot例

はじめにコントロール細胞(WT細胞)と,ある遺伝子をノックアウトした細胞(KO細胞)に発現するタンパク質群について,増減を比較し,その増減が統計的に有意なのかを,ボルケーノプロットにより描画・検討する。

ボルケーノプロット(Volcano plot):
縦軸が群間比 (WTとKOのlog2(fold change)) ,横軸がp値 (-log10(p値) )の散布図。

データ例のDLこちらのリ

もっとみる
Scipyのcurve_fitを使用した非線形曲線Y=ax/x+bへの近似方法

Scipyのcurve_fitを使用した非線形曲線Y=ax/x+bへの近似方法

はじめに

生化学でよく出てくる酵素反応の経時変化を,下記の式への近似する。

まずはモジュールのimport

# import moduleimport numpy as npimport matplotlib.pyplot as plt

数値例を作成
(X値,0から360までの範囲で,上記のような曲線に近いY値が得られたという想定)。

X = [0,10,20,30,60,120,180

もっとみる
ウエスタンブロットの画像トリミング方法

ウエスタンブロットの画像トリミング方法

はじめにウエスタンブロットを含めた画像データは画像編集ソフトでトリミングを行なった後に,Adbe IllustratorやInkscapeなどのベクター系の描画ツールで作図を行うことが理想である。
Power Pointを用いることで画像データを作成することも出来るが,本来Power Pointはプレゼンテーションツールであり,画像編集のためのツールではない。また,pptファイルにコピーされた画像

もっとみる
外部ファイルのクラスを読み込む方法

外部ファイルのクラスを読み込む方法

ipython上などで下記プログラムを実行する場合

class Hello(): def __init__(self, name): self.name = name def sayhello(self): print ('Hello {} !!!'.format(self.name))

いつも通り

greeting = Hello('Python')gree

もっとみる
Numpy memo (随時更新)

Numpy memo (随時更新)

Numpyで今まで使ったものを記録(随時更新)

一次元配列のarrayを二次元配列にして指定の数に分割する:np.array_splita = np.array([1,2,3,4,5,6,7,8,9])を二次元に変換後,3分割して
np.array([[1,2,3],[4,5,6],[7,8,9]])にする

aをn分割する際は np.array_aplit(a, n)
三分割なら np.arr

もっとみる
Pandasを使って複数の独立した実験データから平均値と標準偏差を算出する方法

Pandasを使って複数の独立した実験データから平均値と標準偏差を算出する方法

三回の独立した実験(exp_number, 1-3) で,サンプルa-cの酵素活性(relateive_activity) を算出したデータ例を使用し,平均値(mean)と標準偏差(std)を計算する。 
自作関数を作ってそれぞれを計算する方法も記録しておく。

データの準備# import modulesimport numpy as npimport pandas as pdimport ma

もっとみる
Pandas でよく忘れてしまう事柄(随時追加)

Pandas でよく忘れてしまう事柄(随時追加)

Pythonを使ったデータ解析をする際によく忘れてgoogle検索しがちなものを下記に書き留めておく(随時追加予定)

Excelファイルのimport(read_excel)'xlrd'をインストールすると下記方法でxlsxファイルをデータフレームとして読み込み可能

pd.read_excel('filename.xlsx')

Excelファイルのexport(to_excel)下記でデータ

もっとみる
細胞内に発現する3つのタンパク質の相対量変化を例にした折れ線グラフの描画 (Matplotlib, Seaborn)

細胞内に発現する3つのタンパク質の相対量変化を例にした折れ線グラフの描画 (Matplotlib, Seaborn)

3種類のタンパク質 Protein1, Protein2, Protein3について,あるタイムポイント(0分)から始めて90分まで,30分おきに細胞内発現量を経時的に定量・算出したデータを想定して,折れ線グラフを作成する方法を書き留める。

全てのデータは独立した三回の実験(exp_number:1-3)によって算出されたものと想定し,平均値をプロットし,エラーバーは標準偏差値を基に描く。
pl

もっとみる
酵素活性データを例にした棒グラフの描画その2(Matplotlib, Seaborn)

酵素活性データを例にした棒グラフの描画その2(Matplotlib, Seaborn)

前回に引き続き今回は酵素1-3の酵素活性について,阻害剤のある場合ない場合に分け,3回繰り返し測定したデータ例を用意。
平均値,標準偏差を計算後,matplotlib, seabornの両方で棒グラフを描いてみる。
なお,このデータではEnzyme2にのみ阻害剤が効果的に効き酵素活性が低下しているという結果を想定している。

前回

仮想データの準備

# 仮想データの準備a,b,c,d,e,f

もっとみる
酵素活性データを例にした棒グラフの描画(Matplotlib, Seaborn)

酵素活性データを例にした棒グラフの描画(Matplotlib, Seaborn)

バイオ実験で出てきそうな酵素活性を測定した実験データをmatplotlibとseabornで棒グラフにする例を書き留める。

(データ例の概要)酵素1-3の酵素活性を3回繰り返し測定したデータ(下表)。 

平均値,標準偏差を計算後, matplotlib, seabornの両方で棒グラフを描く

データの準備# import modulesimport numpy as npimport pa

もっとみる
matplotlibでsubplot

matplotlibでsubplot

subplotの方法について覚えたことを記録。
suplotは複数のグラフを一つの図として表示する方法。ネットで検索しながら何となくコピペして作っていていたのだが,どうやらやり方が一つではないらしく,よく混乱する。。。 今回subplotの方法についてまとめることで,理解を深めようと思う。

Subplotの方法についての概要自分が調べた限りsubplotする方法は大きく分けて以下の二つがあるよう

もっとみる
Matplotlib(Seaborn)の基本的な使い方

Matplotlib(Seaborn)の基本的な使い方

最近,研究活動(おもに実験)で得たデータを計算したり,可視化するのに使っているPython。 可視化(グラフ作成)に使うMatplotlib(Seaborn)の使い方をすぐ忘れてしまうので,自分のために書き留めておこうと思う。

ごく基本的な使い方モジュールのインポート

import numpy as npimport matplotlib.pyplot as pltimport seaborn

もっとみる