如何基于标签向量自动生成指定0-1对比矩阵？（R语言）

给定向量：

A <- c("CTRL", "A", "B", "C")

希望生成如下对比矩阵（注：原示例列名存在笔误，应为[CTRL]而非[CTRIL]）：

[A]   [B]  [C] [CTRL]
[A-CTRL]      1     0.   0.   -1
[B-CTRL]      0.    1.   0.   -1
[C-CTRL]      0.    0.   1.   -1

当前手动实现的代码如下：

levels(global.proteinSummarization$ProteinLevelData$GROUP)

comparison1 <- matrix(c(1,-1,0,0,0,0,0,0), nrow=1)
comparison2 <- matrix(c(0,-1,1,0,0,0,0,0), nrow=1)
comparison3 <- matrix(c(0,-1,0,1,0,0,0,0), nrow=1)

comparison <- rbind(comparison1, comparison2, comparison3)
row.names(comparison) <- c("A-CTRL", "B-CTRL", "C-CTRL")
colnames(comparison) <- levels(global.proteinSummarization$ProteinLevelData$GROUP) 

请问是否有更自动化的实现方式？

方法一：基础R原生实现 #

无需额外包，通过原生矩阵操作即可批量生成目标矩阵，适配任意数量的处理组：

# 定义分组向量（可替换为你的levels结果）
groups <- c("CTRL", "A", "B", "C")
# 提取所有非CTRL的处理组
treatments <- setdiff(groups, "CTRL")
# 创建处理组对应的单位矩阵
mat <- diag(nrow = length(treatments), ncol = length(treatments))
# 添加CTRL列，全部填充-1
mat <- cbind(mat, rep(-1, length(treatments)))
# 设置列名和行名
colnames(mat) <- c(treatments, "CTRL")
rownames(mat) <- paste(treatments, "CTRL", sep = "-")
# 调整列顺序与原groups一致（可选，若不需要可省略）
mat <- mat[, groups]

方法二：利用limma包快速生成 #

如果你的工作流涉及差异分析，limma包的makeContrasts函数可以更简洁地生成对比矩阵：

library(limma)

# 定义分组向量
groups <- c("CTRL", "A", "B", "C")
treatments <- setdiff(groups, "CTRL")
# 构造对比字符串
contrast_strings <- paste(treatments, "CTRL", sep = "-")
# 生成对比矩阵
contrast_mat <- makeContrasts(contrasts = contrast_strings, levels = groups)

该方法生成的矩阵可直接用于limma的差异分析流程，无需额外调整格式。

内容的提问来源于stack exchange，提问作者Paolo Lorenzini