[英]put correlation coefficient on ggplot scatter plot after faceting
我有问题在facet_wrap之后用另一个变量将相关系数放在我的散点图上。 下面是我使用mtcars数据集进行说明的示例。 当我绘制出来时,两个图都具有相同的相关数。 似乎没有为每个方面计算相关系数。 我无法找到实现这一目标的方法。 真的很感激,如果有人可以帮助那...
library(ggplot2)
library(dplyr)
corr_eqn <- function(x,y, method='pearson', digits = 2) {
corr_coef <- round(cor.test(x, y, method=method)$estimate, digits = digits)
corr_pval <- tryCatch(format(cor.test(x,y, method=method)$p.value,
scientific=TRUE),
error=function(e) NA)
paste(method, 'r = ', corr_coef, ',', 'pval =', corr_pval)
}
sca.plot <- function (cor.coef=TRUE) {
df<- mtcars %>% filter(vs==1)
p<- df %>%
ggplot(aes(x=hp, y=mpg))+
geom_point()+
geom_smooth()+
facet_wrap(~cyl, ncol=3)
if (cor.coef) {
p<- p+geom_text(x=0.9*max(df$hp, na.rm=TRUE),
y=0.9*max(df$mpg, na.rm=TRUE),
label = corr_eqn(df[['hp']],df[['mpg']],
method='pearson'))
}
return (p)
}
sca.plot(cor.coef=TRUE)
通过变量inputFacet
调用facet,遍历此变量以计算corr_enq
并使用带有get
变量名称绘制facet。
在闪亮的你可能有用户输入作为input$facet
在这里它被称为inputFacet
。 我们绘制了在facet_wrap(~ get(inputFacet), ncol = 3)
获取此变量的主要情节。 接下来,我们使用for(i in seq_along(resCor$facets))
循环遍历所有facet选项for(i in seq_along(resCor$facets))
并将结果存储在rescore
。
这应解决“不计算每个方面的相关系数”问题。
library(dplyr)
library(ggplot2)
inputFacet <- "cyl"
cor.coef = TRUE
df <- mtcars
p <- df %>%
ggplot(aes(hp, mpg))+
geom_point()+
geom_smooth()+
facet_wrap(~ get(inputFacet), ncol = 3)
if (cor.coef) {
resCor <- data.frame(facets = unique(mtcars[, inputFacet]))
for(i in seq_along(resCor$facets)) {
foo <- mtcars[mtcars[, inputFacet] == resCor$facets[i], ]
resCor$text[i] <- corr_eqn(foo$hp, foo$mpg)
}
colnames(resCor)[1] <- inputFacet
p <- p + geom_text(data = resCor,
aes(0.9 * max(df$hp, na.rm = TRUE),
0.9 * max(df$mpg, na.rm = TRUE),
label = text))
}
p
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.