2×2 因子设计的设计矩阵
在这个 2×2 因子实验中，我们研究干旱对树木生长的影响。两种不同类型的杨树（Populus）在两种不同的供水条件下生长。若要为每个分组设置一个模型系数，您需要先将这两个变量组合起来。
本练习是课程的一部分
在 R 中使用 limma 进行差异表达分析
练习说明
工作区中已加载包含 Populus 数据的 ExpressionSet 对象 eset。
将变量
type（树的类型）和water（正常供水 vs. 干旱）合并为一个单一的因子变量。使用
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)