这是我为线性建模函数构建设计矩阵的代码:
f <- factor(targets$Sample.Name, levels = unique(targets$Sample.Name))
design <- model.matrix(~0 + f)
colnames(design) <- levels(f)
我不确定如何解释公式“~0”。我查找了 ?lm() 并发现如果一个公式具有隐含的截距项,则可以使用 y ~ x - 1 或 y ~ 0 + x 来删除它,但我不确定这是否是相同的情况。