3 组的设计矩阵
在缺氧实验中,干细胞暴露于 3 个不同的氧气水平:1%、5% 和 21%。为了对比这 3 种环境下的细胞功能,请使用组均值参数化(group-means parametrization),指定一个对每个氧气水平各有一个系数的线性模型。
本练习是课程的一部分
在 R 中使用 limma 进行差异表达分析
练习说明
包含缺氧数据的 ExpressionSet 对象 eset 已加载到您的工作区。
- 使用
model.matrix创建一个无截距的设计矩阵。请回忆,本研究关注的变量(3 个氧气水平)位于表型数据框的oxygen列中。
交互式实操练习
通过完成这段示例代码来试试这个练习。
# Create design matrix with no intercept
design <- ___(~___ + ___, data = ___(eset))
# Count the number of samples modeled by each coefficient
colSums(design)