提取expresssion_data矩阵中的前两列为gene_mapping矩阵，回复R代码

如果你有两个矩阵，一个叫`expression_data`，它有多个列，而另一个叫`gene_mapping`，并且你想从`expression_data`中提取出前两列，同时将这两列的内容与`gene_mapping`的前两列进行关联，你可以使用`cbind()`函数结合子集操作来完成。如果`expression_data`的第一、二列恰好与`gene_mapping`的相应列相匹配，那么代码可以像这样：

# 假设expression_data的第一、二列分别是你要提取的基因ID
# gene_mapping的第一、二列分别是对应基因名称

# 提取expression_data的前两列
selected_expression_data <- expression_data[, 1:2]

# 提取gene_mapping的前两列并重命名，方便后续操作
selected_gene_mapping <- gene_mapping[, c("Gene_ID_1", "Gene_ID_2")] %>%
  setNames(c("Gene_ID", "Gene_Name"))

# 如果expression_data的第一列与gene_mapping的"Gene_ID_1"列匹配...
if (all(selected_expression_data[, 1] == selected_gene_mapping$Gene_ID)) {
  # ...则可以使用cbind将两个部分合并
  merged_data <- cbind(selected_expression_data, selected_gene_mapping[, -c(1)])  
} else {
  warning("基因ID在两个数据集中不匹配")
}

# 结果`merged_data`现在包含expression_data的前两列加上gene_mapping的相关信息

注意这里的`Gene_ID_1`和`Gene_ID_2`是我假设的`gene_mapping`中的列名，实际使用时应替换为你实际的列名。