R 如何正确绘制长轴回归的置信区间(使用lmodel2)?
如何正确绘制缩短长轴回归的置信区间 我在中找到了一些建议,但我对CIs有点困惑。对截距使用较高的CI限值,对坡度使用较低的限值,反之亦然,这难道不是更有意义吗 基于joran在上述帖子中的示例,以下代码(注意CI的索引):R 如何正确绘制长轴回归的置信区间(使用lmodel2)?,r,R,如何正确绘制缩短长轴回归的置信区间 我在中找到了一些建议,但我对CIs有点困惑。对截距使用较高的CI限值,对坡度使用较低的限值,反之亦然,这难道不是更有意义吗 基于joran在上述帖子中的示例,以下代码(注意CI的索引): 库(lmodel2) 图书馆(GG2) 种子(43) dat library(lmodel2) library(ggplot2) set.seed(43) dat <- data.frame(a=log10(rnorm(50, 30, 10)), b=log10(rn
库(lmodel2)
图书馆(GG2)
种子(43)
dat
library(lmodel2)
library(ggplot2)
set.seed(43)
dat <- data.frame(a=log10(rnorm(50, 30, 10)), b=log10(rnorm(50, 20, 2)))
mod <- lmodel2(a ~ b, data=dat,"interval", "interval", 99)
ggplot(dat,aes(x=b,y=a)) + geom_point() +
geom_abline(intercept=mod$regression.results[4,2],
slope=mod$regression.results[4,3],colour="blue") +
geom_abline(intercept=mod$confidence.intervals[4,2],
slope=mod$confidence.intervals[4,5],colour="red") +
geom_abline(intercept=mod$confidence.intervals[4,3],
slope=mod$confidence.intervals[4,4],colour="red")