始める無料で始める

ボルケーノプロット

各コントラストにおける遺伝子発現変動の大きさを、ボルケーノプロットで可視化します。ボルケーノプロットでは、y軸に差次的発現である対数オッズ、x軸に対数フォールドチェンジを表示します。

この演習はコースの一部です

R での limma を用いた Differential Expression 解析

コースを見る

演習の手順

当てはめ済みのモデルオブジェクト fit2 がワークスペースに読み込まれています。limma パッケージもすでにロード済みです。

  • データフレーム fit2$genes から遺伝子シンボル(symbol)にアクセスしてください。

  • 3つのコントラストそれぞれについて、上位5遺伝子をハイライトしたボルケーノプロットを作成してください。

実践的なインタラクティブ演習

このサンプルコードを完成させて、この演習に挑戦してみましょう。

# Extract the gene symbols
gene_symbols <- fit2$genes[___]

# Create a volcano plot for the contrast dox_wt
___(fit2, coef = ___, highlight = 5, names = gene_symbols)

# Create a volcano plot for the contrast dox_top2b
___(fit2, coef = ___, highlight = ___, names = gene_symbols)

# Create a volcano plot for the contrast interaction
___(fit2, coef = ___, highlight = ___, names = ___)
コードを編集して実行