首页
学习
活动
专区
圈层
工具
发布
社区首页 >问答首页 >Gggcorr图:无关紧要的值

Gggcorr图:无关紧要的值
EN

Stack Overflow用户
提问于 2017-09-06 08:06:32
回答 1查看 1.4K关注 0票数 3

我用下面的代码使用ggcorrplot创建了一个关联图:

代码语言:javascript
复制
library(ggcorrplot)

corCN <- subset(csfBL, BL_Diaggroups==1, 
                select=c("log_ab38", "log_ab40", "log_ab42", "ABratio", "log_YKL", "logNgcomplete", "log_NFL"))

colnames(corCN) <- c("Aß 1-38","Aß 1-40", "Aß 1-42", "Aß 42/40", "YKL-40", "Ng", "NFL")

corrCN <- cor(corCN, method=c("spearman"), use="complete.obs")
p.matCN <- cor_pmat(corCN, method=c("spearman"), use="complete.obs")

CNcorrplot <- ggcorrplot(corrCN, type="lower", lab=TRUE, ggtheme = ggplot2::theme_classic,
                          p.mat=p.matCN, insig ="blank")

CNcorrplot2 <- CNcorrplot + ggtitle("CN") + theme(plot.title = element_text(hjust = 0.5)) 

情节是这样的。

我唯一想要更改的是,不重要的值中有"NS“,或者完全是空的(所以没有显示0)。实际上,我需要一些代码来说明:

代码语言:javascript
复制
if insig=TRUE, lab=FALSE
EN

回答 1

Stack Overflow用户

回答已采纳

发布于 2017-09-06 15:43:15

我用新的ggcorrplot选项创建了一个名为myggcorrplot的修改版lab.notsig (表示不显著关联的标签)。

下载文件这里并将其保存在工作目录中作为myggcorrplot.r

然后,运行以下代码:

代码语言:javascript
复制
library(ggcorrplot)  
library(mvtnorm)
# A toy dataset
set.seed(1)
n <- 100
p <- 7
A <- matrix(runif(p^2)*2-1, ncol=p) 
Sigma <- cov2cor(t(A) %*% A)
corCN <- data.frame(rmvnorm(n, mean=rep(0,p), sigma=Sigma))
colnames(corCN) <- c("Aß 1-38","Aß 1-40", "Aß 1-42", "Aß 42/40", "YKL-40", "Ng", "NFL")

corrCN <- cor(corCN, method=c("spearman"), use="complete.obs")
p.matCN <- cor_pmat(corCN, method=c("spearman"), use="complete.obs")

# Load the myggcorrplot function
source("myggcorrplot.r")
# Set in the lab.notsig option the label for not significant correlations
CNcorrplot <- myggcorrplot(corrCN, type="lower", lab=TRUE, 
              ggtheme = ggplot2::theme_classic,
              p.mat=p.matCN, insig ="blank", lab.notsig="NS")

CNcorrplot2 <- CNcorrplot + ggtitle("CN") + 
                   theme(plot.title = element_text(hjust = 0.5)) 
CNcorrplot2

票数 2
EN
页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接:

https://stackoverflow.com/questions/46069907

复制
相关文章

相似问题

领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档