始める無料で始める

2×2 要因計画のデザイン行列

干ばつが樹木の成長に与える影響を調べるこの 2×2 要因実験では、2 種の Populus(ヤマナラシ属)を、2 通りの水量で育てました。各グループに 1 つずつモデル係数を割り当てるには、まず 2 つの変数を結合する必要があります。

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

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

コースを見る

演習の手順

Populus のデータを含む ExpressionSet オブジェクト eset がワークスペースに読み込まれています。

  • 変数 type(樹種)と water(通常 vs. 乾燥)を 1 つの因子変数に結合してください。

  • model.matrix を使って、切片なしのデザイン行列を作成してください。

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

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

# Create single variable
group <- with(___(eset), paste(___, ___, sep = "."))
group <- factor(group)

# Create design matrix with no intercept
design <- model.matrix(~___ + ___)
colnames(design) <- levels(group)

# Count the number of samples modeled by each coefficient
colSums(design)
コードを編集して実行