如何在R语言的小提琴图（Violin Plot）中获取各时间点的均值与平均表达量？

针对你的需求，这里提供两种简单的方法来获取每个时间点Cx3cl1基因的表达量均值，都是基于Seurat对象操作：

方法1：直接提取数据计算 #

从Seurat对象中提取基因表达数据和分组信息，用基础R函数计算均值：

# 提取Cx3cl1的表达值与对应的time分组
expr_df <- data.frame(
  cx3cl1_expr = spinalcord_neurons@assays$RNA@data["Cx3cl1", ],
  time_point = spinalcord_neurons@meta.data$time
)

# 按time分组计算均值
grouped_mean <- aggregate(cx3cl1_expr ~ time_point, data = expr_df, FUN = mean)
print(grouped_mean)

如果需要用原始测序计数计算，把@assays$RNA@data换成@assays$RNA@counts即可。

方法2：用Seurat内置函数计算 #

Seurat的AverageExpression函数专门用于按组计算基因平均表达量，适合后续批量处理多个基因：

# 按time分组计算指定基因的平均表达量
avg_expr <- AverageExpression(
  object = spinalcord_neurons,
  features = "Cx3cl1",
  group.by = "time"
)

# 转换为易读的数据框格式
result_df <- as.data.frame(t(avg_expr$RNA))
result_df$time <- rownames(result_df)
colnames(result_df)[1] <- "mean_expression"
print(result_df)

两种方法得到的结果本质一致，第二种更贴合Seurat的分析流程，后续加多个基因只需要把features参数改成c("基因1", "基因2")即可。

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