基因表达矩阵保存在txt中，行为基因，列为样本，画出单样本基因表达水平箱线图，并提取每个样本的中位数，R语言代码

假设你的基因表达矩阵保存在名为 `gene_expression.txt` 的文本文件中，其中第一行为样本名称，第一列为基因名称，剩余元素为基因表达水平。你可以使用以下 R 语言代码来绘制单个样本的基因表达水平箱线图，并提取每个样本的中位数：

# 读取基因表达矩阵
gene_exp <- read.table("gene_expression.txt", header=TRUE, row.names=1)

# 选择要绘制箱线图的样本
sample_name <- "Sample1"
sample_exp <- gene_exp[, sample_name]

# 绘制箱线图
boxplot(sample_exp, main=paste("Gene expression of", sample_name), ylab="Expression level")

# 提取中位数
sample_median <- median(sample_exp)
cat("Median expression level of", sample_name, "is", sample_median)

其中，`read.table()` 函数用于读取基因表达矩阵文件，`header=TRUE` 表示第一行是列名，`row.names=1` 表示第一列是行名。`gene_exp[, sample_name]` 用于选取要绘制箱线图的样本数据，`sample_name` 为样本名称。`boxplot()` 函数用于绘制箱线图，`main` 参数为图表标题，`ylab` 参数为 Y 轴标签。`median()` 函数用于计算样本基因表达水平的中位数，`cat()` 函数用于打印输出结果。你可以将代码中的 `sample_name` 变量替换为你要绘制箱线图和计算中位数的样本名称。如果你想提取所有样本的中位数，可以使用一个循环来遍历每个样本的列，并计算中位数。