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)