R 如何正确绘制长轴回归的置信区间(使用lmodel2)?

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

如何正确绘制缩短长轴回归的置信区间

我在中找到了一些建议,但我对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(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")