小提琴图:单个

这篇博客展示了如何利用R语言的ggplot2和reshape2库对基因表达数据进行处理和可视化。首先,通过读取`File_tpm.csv`文件获取数据,然后筛选出表达水平介于10到500之间的基因。接着,使用`melt`函数将数据重塑为长格式,以便于绘图。最后,通过`ggplot2`创建了两种不同的图形:一种包含小提琴图和点状图,另一种则展示了数据的中位数。这些图形有助于理解不同样本间的基因表达差异。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

library(gcookbook)
library(ggplot2)
library(reshape2) 
data <- read.csv("File_tpm.csv",row.names =1)
data1 <- data[rowSums(data) > 10,]
data1 <- data1[rowSums(data1) < 500,]

data_melt <- melt(data1)
##(1)
ggplot(data_melt, aes(x = variable , y = value , fill = variable)) +
    geom_violin(alpha = 0.5,aes(linetype=NA)) + 
    geom_jitter(shape=21,aes(fill=variable),position = position_jitter(width = 0.2))+
    xlab("sample")+ylab("tpm")+
    theme_bw()+
    theme(legend.position = "none")

在这里插入图片描述

ggplot(data_melt, aes(x = variable , y = value , fill = variable)) +
  geom_violin(alpha = 0.5)+
  stat_summary(fun.y=median, geom="point", fill="white",shape=21,size=2.5)

在这里插入图片描述

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值