WGCNA软阈值绘图

kylin王国

已于 2022-10-08 13:02:57 修改

阅读量319

点赞数

分类专栏：数据可视化数据分析 R语言文章标签：开发语言 r语言开源软件数据分析

于 2022-10-08 09:38:55 首次发布

本文链接：https://blog.csdn.net/weixin_47634487/article/details/127203089

版权

R语言同时被 3 个专栏收录

32 篇文章 5 订阅

订阅专栏

数据可视化

18 篇文章 0 订阅

订阅专栏

数据分析

17 篇文章 0 订阅

订阅专栏

df <- data.frame(x = c(seq(1, 10, by = 1), seq(12, 30, by = 2)),
                 y = c(0,0.1,0.19,0.3,0.4,0.6,0.68,
                       0.78,0.8,0.815,0.88,0.9,0.91,
                       0.92,0.91,0.93,0.93,0.92,0.915,0.93))
plot(df$x,df$y,pch = 16,cex=0,
     xlab="Soft Power Threshold",
     ylab="Scale Free Topology Model Fit,signed R^2",type="n",
     main = paste("Scale independence"),yaxt='n')
abline(h=0.8,col="black",lty = 2)
abline(v=9,col="black",lty = 2)
points(x = 9, y = 0.8,pch = 16,cex=3)
text(df$x,df$y,
     labels=df$x,cex=1,col="black")
text(x = 9, y = 0.8,labels=9,cex=1,col="white")
axis(side = 2,at = c(0,0.2,0.4,0.6,0.8,1),labels =  c(0,0.2,0.4,0.6,0.8,1))

在这里插入图片描述

df <- data.frame(x = c(seq(1, 10, by = 1), seq(12, 30, by = 2)),
                 y = c(0,0.1,0.19,0.3,0.4,0.6,0.68,
                       0.78,0.8,0.815,0.88,0.9,0.91,
                       0.92,0.91,0.93,0.93,0.92,0.915,0.93))
plot(df$x,df$y,pch = 16,cex=0,
     xlab="Soft Power Threshold",
     ylab="Scale Free Topology Model Fit,signed R^2",type="n",
     main = paste("Scale independence"),yaxt='n')
abline(h=0.8,col="black",lty = 2,lwd = 2)
abline(v=9,col="black",lty = 2,lwd = 2)
axis(side = 2,lwd = 2,at = c(-1,3))
axis(side = 1,lwd = 2,at = c(-1,300))
axis(side = 4,lwd = 2,at = c(-1,3))
axis(side = 3,lwd = 2,at = c(-1,300))

points(x = 9, y = 0.8,pch = 16,cex=3)
text(df$x,df$y,
     labels=df$x,cex=1,col="black")
text(x = 9, y = 0.8,labels=9,cex=1,col="white")
axis(side = 2,at = c(0,0.2,0.4,0.6,0.8,1),
     labels =  c(0,0.2,0.4,0.6,0.8,1))

在这里插入图片描述

df <- data.frame(x = c(seq(1, 10, by = 1), seq(12, 30, by = 2)),
                 y = c(0,0.1,0.19,0.3,0.4,0.6,0.68,
                       0.78,0.8,0.815,0.88,0.9,0.91,
                       0.92,0.91,0.93,0.93,0.92,0.915,0.93))
plot(df$x,df$y,pch = 16,cex=0,
     xlab="Soft Power Threshold",
     ylab="Scale Free Topology Model Fit,signed R^2",
     type="n",
     main = paste("Scale independence"),
     yaxt='n',
     ylim = c(0,1))
abline(h=0.8,col="black",lty = 2,lwd = 2)
abline(v=9,col="black",lty = 2,lwd = 2)
axis(side = 2,lwd = 2,at = c(-1,3))
axis(side = 1,lwd = 2,at = c(-1,300))
axis(side = 4,lwd = 2,at = c(-1,3))
axis(side = 3,lwd = 2,at = c(-1,300))

points(x = 9, y = 0.8,pch = 16,cex=3,col = 'black')
text(df$x,df$y,
     labels=df$x,cex=1,col="black")
text(x = 9, y = 0.8,labels=9,cex=1,col="white")
axis(side = 2,at = c(0,0.2,0.4,0.6,0.8,1),
     labels =  c(0,0.2,0.4,0.6,0.8,1))