R 使い方 | volcano plot 描き方（ボルケーノプロット）

volcano plotは、発現解析の結果を表現でよく使われるグラフです。

横軸に発現比(logFC)、縦軸にpvalue(-log10した値)をとります。

Rを使って、volcano plot （ボルケーノ プロット）を作成します。

valcano plot 描き方

単純な volcano-plot

> plot(-log10(comp_pValue) ~ comp_logFC)

または

> plot(comp_logFC, -log10(comp_pValue))

データフレームから必要な列を抽出して描く場合は、tidyverseのtransmutateを使うと便利です。

library(tidyverse)

> de_result <- read_tsv("de_result.tsv")

> vdata <- de_result %>% transmute(GeneSymbol, x=comp_logFC, y=-log10(comp_pValue))
> plot(vdata$y ~ vdata$x)

volcano-plot に色を付ける

条件（p値 0.05以下、|logFC|>1）でボルケーノプロットに色をつけました。タイトル・軸ラベルも追加しました。

> plot(y ~ x, 
       main="volcano-plot",
       xlab="log2 fold change",
       ylab="-log10 p-value",
       col=ifelse(y>-log10(0.05), 
           ifelse(abs(x)>1, ifelse(x>0, "red", "blue"), 
           "black"), "black"),
       data=vdata
      )

領域を区切る補助線をひきたい場合は、ablineで追加します。

abline(v=1, lty=2)
abline(v=-1, lty=2)

ちなみに、画像で保存する場合は、png(“volocano-plot.png”) 〜 dev.off()で囲みます。

> png("volocano-plot.png")

> plot(y ~ x, 
       main="volcano-plot",
       xlab="log2 fold change",
       ylab="-log10 p-value",
       col=ifelse(y>-log10(0.05), 
           ifelse(abs(x)>1, ifelse(x>0, "red", "blue"), 
           "black"), "black"),
       data=vdata
      )

> abline(v=1, lty=2)
> abline(v=-1, lty=2)
> abline(h=-log10(0.05), lty=2)

> dev.off()

論文・データ解析の支援

個別に直接相談できるのが良いと、たくさんの方にご利用頂いています。

関連記事１

R 使い方 ボックスプロット（box plot） グラフの描き方

データを理解する方法の１つに「グラフ」があります。詳細なデータをまとめた表に比べると「詳細さ」は劣りますが、データの全体像を理解しやすくなります。そして、目的に応じたグラフを使うことが大事です。Rを使って、色んなグラフを描画してみましょう。...

bioinfo-dojo.net

2016.05.23

R 使い方 折れ線グラフ グラフの描き方

データを理解する方法の１つに「グラフ」があります。詳細なデータをまとめた表に比べると「詳細さ」は劣りますが、データの全体像を理解しやすくなります。そして、目的に応じたグラフを使うことが大事です。Rを使って、色んなグラフを描画してみましょう。...

bioinfo-dojo.net

2016.05.27

R 使い方 パイチャート（円グラフ） グラフの描き方

データを理解する方法の１つに「グラフ」があります。詳細なデータをまとめた表に比べると「詳細さ」は劣りますが、データの全体像を理解しやすくなります。そして、目的に応じたグラフを使うことが大事です。Rを使って、色んなグラフを描画してみましょう。...

bioinfo-dojo.net

2016.05.24

R 使い方 棒グラフ グラフの描き方

データを理解する方法の１つに「グラフ」があります。詳細なデータをまとめた表に比べると「詳細さ」は劣りますが、データの全体像を理解しやすくなります。そして、目的に応じたグラフを使うことが大事です。今回は、「棒グラフ」です。とにかく早く問題解決...

bioinfo-dojo.net

2016.05.22

R 使い方 散布図 グラフの描き方

データを理解する方法の１つに「グラフ」があります。詳細なデータをまとめた表に比べると「詳細さ」は劣りますが、データの全体像を理解しやすくなります。そして、目的に応じたグラフを使うことが大事です。Rを使って、色んなグラフを描画してみましょう。...

bioinfo-dojo.net

2016.05.24

R 使い方 ヒストグラム グラフの描き方

データを理解する方法の１つに「グラフ」があります。詳細なデータをまとめた表に比べると「詳細さ」は劣りますが、データの全体像を理解しやすくなります。そして、目的に応じたグラフを使うことが大事です。Rを使って、色んなグラフを描画してみましょう。...

bioinfo-dojo.net

2016.05.23

関連記事２

404 NOT FOUND | バイオインフォ 道場 [bioinfo-Dojo]

バイオインフォ 道場

bioinfo-dojo.net

R 使い方 目盛りのべき乗回避 グラフの描き方

Rでグラフを描画する場合、位置やサイズを自動調整してくれるので大変便利です。しかし、カスタマイズしようと思えば、各種パラメータを調整する必要があります。より綺麗なグラフに仕上げるために、軸ラベルや目盛り、余白の調整についてまとめます。目盛り...

bioinfo-dojo.net

2017.04.11

R 使い方 任意の軸ラベル グラフの描き方

Rでグラフを描画する場合、位置やサイズを自動調整してくれるので大変便利です。しかし、カスタマイズしようと思えば、各種パラメータを調整する必要があります。より綺麗なグラフに仕上げるために、軸ラベルや目盛り、余白の調整についてまとめます。任意の...

bioinfo-dojo.net

2017.04.11

R 使い方 グラフにテキストを書き込む グラフの描き方

Rでグラフを描画する場合、位置やサイズを自動調整してくれるので大変便利です。しかし、カスタマイズしようと思えば、各種パラメータを調整する必要があります。より綺麗なグラフに仕上げるために、軸ラベルや目盛り、余白の調整についてまとめます。データ...

bioinfo-dojo.net

2017.04.11

R 使い方 ラベル向き調整 グラフの描き方

Rでグラフを描画する場合、位置やサイズを自動調整してくれるので大変便利です。しかし、カスタマイズしようと思えば、各種パラメータを調整する必要があります。より綺麗なグラフに仕上げるために、軸ラベルや目盛り、余白の調整についてまとめます。ラベル...

bioinfo-dojo.net

2017.04.11

R 使い方 余白サイズ グラフの描き方

Rでグラフを描画する場合、位置やサイズを自動調整してくれるので大変便利です。しかし、カスタマイズしようと思えば、各種パラメータを調整する必要があります。より綺麗なグラフに仕上げるために、軸ラベルや目盛り、余白の調整についてまとめます。余白の...

bioinfo-dojo.net

2017.04.11

R 使い方 軸ラベルのサイズ グラフの描き方

Rでグラフを描画する場合、位置やサイズを自動調整してくれるので大変便利です。しかし、カスタマイズしようと思えば、各種パラメータを調整する必要があります。より綺麗なグラフに仕上げるために、軸ラベルや目盛り、余白の調整についてまとめます。軸ラベ...

bioinfo-dojo.net

2017.04.11

R 使い方 チャートやグラフの組合せ グラフの描き方

チャートやグラフは、データの理解を手助けしてくれます。「統計ソフト R」では様々なグラフが手軽に描けるので便利です。単体のチャートやグラフでも役に立ちますが、表現には限界があり、一長一短となるものも少なくありません。異なる種類のチャートやグ...

bioinfo-dojo.net

2016.11.21

R 使い方 軸と目盛りの調整まとめ グラフの描き方

Rでグラフを描画する場合、位置やサイズを自動で調整してくれるので大変便利です。しかし、少しカスタマイズしようと思えば、色々なパラメータを調整する必要があります。より綺麗なグラフに仕上げるために、軸の目盛りの調整についてまとめます。とにかく早...

bioinfo-dojo.net

2016.08.23

tidyverse インストールとロード tidyverseの主なパッケージ

tidyverseはRのパッケージの１つで、データ操作（抽出・連結など）を統合的に扱うためのツール群です。Rを起動後に、tidyverseのインストールを行います。インストール完了後は、使うときにパッケージをロードして使います。 とにかく早...

bioinfo-dojo.net

2019.03.01

R 使い方 カラーパレット RColorBrewer 透明な色とグラフ

RColorBrewerは便利でよく使われるパッケージです。色々な色の組み合わせを簡単に使うことができます。しかし、複数のグラフを重ねたときに、透過度を簡単に設定できないのは使いにくく感じる場合があります。色の透過度を設定する方法を調べまし...

bioinfo-dojo.net

2017.08.18

R 使い方 グラフの色 自作のカラーパレットとデフォルトのカラーパレット

グラフを描画する際、色の組み合わせを１つ１つ考えるのは面倒ですよね。Rには、あらかじめ色の組み合わせを用意したカラーパレットが準備されています。カラーパレットを使えば、洗練された色の組み合わせで綺麗なグラフを描くことができます。自分で色の組...

bioinfo-dojo.net

2017.08.14

R グラフ 凡例の調整 legendオプション

Rでグラフを描いたとき、凡例が気になる場合がありませんか？工夫できるポイントを紹介します。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談R グラフ 凡例の細かい調整(legendオプション)グラフでよく使うヒストグラムを例にし...

bioinfo-dojo.net

2017.08.17

R 使い方 複数のグラフを並べる グラフの描き方

Rでグラフを並べてみましょう。似たようなデータや条件を変えたデータを比較する場合に、わかりやすくなります。注意点や工夫すべきポイントをまとめます。並べるという観点で別の記事を整理しました。とにかく早く問題解決したい人はこちら＞＞直接、データ...

bioinfo-dojo.net

2017.04.13

R 使い方 凡例・エラーバー・グラフ保存などの図の仕上げ グラフの描き方

「統計ソフト R」によるチャートやグラフの描画は、手軽で大変便利です。ただ、仕上げには細かい調整が必要です。今回は、そのようなときに役立つ「凡例」や「エラーバー」の描き方に関するテクニックをまとめます。とにかく早く問題解決したい人はこちら＞...

bioinfo-dojo.net

2016.11.25

R 使い方 軸・ラベルの調整（向き・サイズ・色など） グラフの描き方

Rによるplot(グラフ)の描画は、手軽で大変便利です。ただ、細かい調整に関しては、全てをデフォルトに任せるわけにはいかないようです。今回は軸について、色々と制御してみます。軸の表示を制御する軸ラベルを変更する・消すx軸・y軸のラベルをそれ...

bioinfo-dojo.net

2016.11.24

R 使い方 散布図 工夫（透明度・グラフ並べるなど） グラフの描き方

統計処理などで２変数の関係を調べるときに、良く使われるのが散布図です。少し工夫するとより使いやすくなります。散布図の工夫についてまとめます。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談散布図を描く散布図はplot関数を使って...

bioinfo-dojo.net

2016.11.17

R 使い方 グラフ重ねる・グラフ分割「上下・左右」 グラフの描き方

R で複数のグラフを描画する際、いろいろな描き方を知っていると役に立つことがあります。よく用いられる「重ねる」「画面を分割する」について、まとめます。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談「R」で複数のグラフを描画する...

bioinfo-dojo.net

2016.03.17

Rの使い方 エラーバーを描く方法

Rのグラフには、エラーバーや凡例、そして文字による説明があるとわかりやすくなります。グラフを仕上げるためのテクニックをまとめます。エラーバー 編。とにかく早く問題解決したい人はこちら＞＞直接、データ解析相談エラーバーを描くRのエラーバーは、...

bioinfo-dojo.net

2017.04.15