Top2b の箱ひげ図
次に妥当性チェックとして、Top2b の発現量の箱ひげ図を作成し、null マウスで Top2b の発現が低いことを確認します。
第1章の演習4で、これと同様の箱ひげ図を作成したことを思い出してください。
この演習はコースの一部です
R での limma を用いた Differential Expression 解析
演習の手順
ドキソルビシンのデータを含む ExpressionSet オブジェクト eset がワークスペースに読み込まれています。
特徴量データの列
"symbol"を"Top2b"で検索し、Top2b のデータを含む行を見つけてください。公式記法を使って箱ひげ図を作成します。式では、
~の左側に発現量、右側に遺伝子型（genotype）変数を置きます。プロットのタイトルには、対応する特徴量データフレームの行を
mainに渡してください。
実践的なインタラクティブ演習
このサンプルコードを完成させて、この演習に挑戦してみましょう。
# Find the row which contains Top2b expression data
top2b <- which(___(eset)[___] == "Top2b")
# Plot Top2b expression versus genotype
boxplot(___(eset)[top2b, ] ~ ___(eset)[___],
main = ___(eset)[top2b, ])