【绘图进阶】之带连线的箱体图绘制(8)

发稿时间：2020-08-21来源：天昊生物

上回我们使用ggpaired包绘制成对样品（2组）的箱体图，但ggpaired只能绘制2个箱体图，有老师希望能绘制连续时间上物种的变化。这里，我们展示4组之间的箱体图绘制结果。

In[1]:

library(ggplot2)

 

read.table读取文件，head查看文件头部

In[2]:

data = read.table('alpha.diversity.index.xls',header = T,sep = 't',row.names = 1)
head(data)

Out[2]:

查看文件列名

In[3]:

rownames(data)

Out[3]:

'A1' 'A2' 'A3' 'A4' 'A5' 'A6' 'B1' 'B2' 'B3' 'B4' 'B5' 'B6' 'C1' 'C2' 'C3' 'C4' 'C5' 'C6' 'D1' 'D2' 'D3' 'D4' 'D5' 'D6'

In [4]:

apply(data,2,sum)

Out[4]:

In [5]:

apply(data,2,function(x){ df = data.frame(x,group=rep(c('A','B','C','D'),each = 6),sub_group=as.factor(rep(1:6,times  = 4)))
                        ggplot(data = df,aes(x=group,y=x,fill=group)) + geom_boxplot() +
                        geom_point(aes(colour=sub_group),size=4) +    
                        geom_line(aes(group=sub_group,colour=sub_group) ,size=0.8,) 
                        })

Out[5]:

往期相关链接：

1、R基础篇

excel不熟练怎么办，R来帮您（一）数据分类汇总；

如何使用Rstudio练习R基础教程；

R相关软件及R包安装；

【零基础学绘图】之气泡图绘制（六）；

【零基础学绘图】之绘制venn图（五）；
【零基础学绘图】之绘制barplot柱状图图（四）；
【零基础学绘图】之绘制heatmap图（三）；

【零基础学绘图】之绘制PCA图（二）；

【零基础学绘图】之alpha指数箱体图绘制（一）；

2、R进阶

【绘图进阶】之配对箱体图绘制(七)；

【绘图进阶】之通路与菌的相关性分析热图(六)；

【绘图进阶】之lefse定制化绘图(五)；

【绘图进阶】之六种带中心点的PCA 图和三维PCA图绘制（四）；

【绘图进阶】之交互式可删减分组和显示样品名的PCA 图（三）；

【绘图进阶】之绘制PCA biplot图(二)；

【进阶篇绘图】之带P值的箱体图、小提琴图绘制（一）；

3、数据提交

3分钟学会微生物多样性云平台数据分析;

3分钟学会CHIP-seq类实验测序数据可视化 —IGV的使用手册；

10分钟搞定多样性数据提交，最快半天内获取登录号，史上最全的多样性原始数据提交教程；

20分钟搞定GEO上传，史上最简单、最详细的GEO数据上传攻略；

4、表达谱分析

表达谱分析（二）通路富集分析和基因互作网络图绘制；

如何对GEO数据进行差异分析；

miRNA靶基因预测软件__miRWalk 3.0；

5、医学数据分析

KING: 样本亲缘关系鉴定工具；

【WGS服务升级】人工智能软件SpliceAI助力解读罕见和未确诊疾病中的非编码突变；

隐性疾病trio家系别忽视单亲二倍体现象——天昊数据分析助力临床疾病诊断新添UPD（单亲二倍体）可视化分析工具；

【昊工具】Oh My God! 太好用了吧！疾病或表型的关键基因查询数据库，我不允许你不知道Phenolyzer；

天昊客户服务中心

手机/微信号：18964693703

【本群将为大家提供】

分享生信分析方案

提供数据素材及分析软件支持

定期开展生信分析线上讲座

QQ号：1040471849