一
加载R包
In [1]:
library(ggplot2)
二
数据读取
data = read.table('alpha.diversity.index.xls',header = T,sep = ' ',row.names = 1)
head(data)
三
绘图
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;
【本群将为大家提供】
分享生信分析方案
提供数据素材及分析软件支持
定期开展生信分析线上讲座
QQ号:1040471849
作者:大熊
审核:有才
来源:天昊生信团
微信扫一扫
关注该公众号
前往“发现”-“看一看”浏览“朋友在看”