【R】数据相关性的可视化

一千零一技|相关性分析及其可视化:copy&paste,搞定

r 复制代码
.libPaths(c("/bioinfo/home/software/miniconda3/envs/R4.0/lib/R/library"))
#data("mtcars")
library("PerformanceAnalytics")
# pdf("test.pdf")
# my_data <- mtcars[, c(1,3,4,5,6,7)]
# print (head(my_data))
# chart.Correlation(my_data, histogram=TRUE, pch=19)
# dev.off()

args    <- commandArgs(trailingOnly = TRUE)
infile  <- args[1]
outdir <- args[2]
names<-basename(infile)
df <- read.delim(infile, header = T, stringsAsFactors = F,row.names = NULL)
#df <- df[,-ncol(df)]
print (df)
index <- row.names(df)
print(index)
# q()
# data_T=as.data.frame((data))
# df1 = anno_col[,"types",drop=FALSE]
pdf(paste0(outdir,"/",names,".correlation_chart.pdf"))
chart.Correlation(df, histogram=TRUE, pch=19,method='spearman')
dev.off()

输入:

输出:

相关推荐
ZZZKKKRTSAE10 分钟前
RHCE认证题解
linux·运维·服务器·ansible·rhce·rhel9·红帽
Lbwnb丶14 分钟前
openssl创建公钥私钥 pem+ der
linux·运维·服务器
❀͜͡傀儡师22 分钟前
达梦数据库日常运维命令
运维·数据库
喵王叭34 分钟前
【大模型实战】向量数据库实战 - Chroma & Milvus
数据库·人工智能·langchain
2501_920047031 小时前
linux-系统性能监控
linux·运维·服务器
云台0951 小时前
HIVE 窗口函数处理重复数据
数据库
星石传说2 小时前
使用R将nc文件转换为asc文件或者tif文件
r语言·生信
得物技术2 小时前
Valkey 单点性能比肩 Redis 集群了?Valkey8.0 新特性分析|得物技术
大数据·数据库·redis
白白白鲤鱼2 小时前
Vue2项目—基于路由守卫实现钉钉小程序动态更新标题
服务器·前端·spring boot·后端·职场和发展·小程序·钉钉
倾听醉梦语2 小时前
Redis作为MySQL缓存的完整指南:从原理到实战
数据库·redis·mysql·缓存