前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >样本相关性分析

样本相关性分析

原创
作者头像
爱学习的小明明
修改2020-09-22 17:55:54
1.4K0
修改2020-09-22 17:55:54
举报
文章被收录于专栏:R语言学习

使用corrplot包分析,使用RPKM值

用于分析的数据,行为基因名字,列为样本名字
用于分析的数据,行为基因名字,列为样本名字
代码语言:javascript
复制
> # <样本相关性>
> ## 1.如果不存在corrplot就安装这个包
> if (!requireNamespace('corrplot', quietly = TRUE))
+   install.packages('corrplot')
> library('corrplot')	#加载corrplot包用于绘制相关性矩阵热图
> 
> 
> p<-file.choose()
> 
> 
> fpkm <- read.table(p, header = T, row.names = 1)	#header=T,第一行指定为列名,row.names=1指定第一列为行名
> 
> 
> corr <- cor(fpkm, method = 'spearman')  #cor函数计算两两样本(列与列)之间的相关系数
> 
> pdf('./sample_correlation.pdf', width = 8, height = 8)	#打开绘图设备,保存为pdf文件
> 
> corrplot(corr, method="shade",type = 'upper', tl.col = 'black', hclust.method="median",order = 'hclust', cl.pos="r",tl.srt = 45, addCoef.col = 'white',diag=F)	
> # type='upper':只显示右上角相关系数矩阵
> # cl.pos=r:图例的颜色条码在右边显示
> # tl.col='black':字体颜色黑色
> # order='hclust':使用层次聚类算法
> # tl.srt = 45:x轴标签倾斜45度
> # addCoef.col='white':添加相关系数数值,颜色白色
> #method:相关图显示方式吗,默认为circle
> dev.off() #配合pdf()使用,关闭绘图设备

method=shade
method=shade
method=pie
method=pie

原创声明:本文系作者授权腾讯云开发者社区发表,未经许可,不得转载。

如有侵权,请联系 cloudcommunity@tencent.com 删除。

原创声明:本文系作者授权腾讯云开发者社区发表,未经许可,不得转载。

如有侵权,请联系 cloudcommunity@tencent.com 删除。

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 使用corrplot包分析,使用RPKM值
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档