我是ggplot2的新手。是否有可能只为数据添加更平滑的数据,在0.05的水平上存在显着的相关性?
我尝试使用以下代码提取重要关联数据,然后为子集添加更平滑的数据。但它没有用。
my_custom_df <- function(p) {
function(df) {
fit <- lm(y ~ x, data = df)
fs <- summary(fit)$fstatistic
pValue <- pf(fs[1], fs[2], fs[3], lower.tail = FALSE)
if (pValue < p) {
df
}
}
}
dat.sig <- ddply(dat, 'g', my_custom_df(0.05))
gp <- ggplot(data = dat, aes(x = c, y = v)) + geom_point() + geom_smooth(data = dat.sig, method = 'lm') + facet_wrap( ~ g, scales = 'free', ncol = 2)
预期结果
我的数据:
c = rep(c(3,2,1.5,1.25,1,0.65),4)
g = c(rep('narG',6), rep('napA',6), rep('nirK',6), rep('nirS',6))
v = c(0.5, 8.0, 19.1, 34.8, 60.6, 63.5, 1.5, 0, 3.7, 22.4, 4.7, 7.9, 2780.7, 179.2, 327.6, 2213, 370.6, 151.4, 13.6, 0.7, 0.1, 1846.9, 132.4, 79)*100000