broom
パッケージを使用して信頼区間を計算する方法は不思議です。私がやろうとしています何R - tidy augment confidence interval
は、シンプルで標準です:
set.seed(1)
x <- runif(50)
y <- 2.5 + (3 * x) + rnorm(50, mean = 2.5, sd = 2)
dat <- data.frame(x = x, y = y)
mod <- lm(y ~ x, data = dat)
visreg
を使用して私は非常に簡単CI
と回帰モデルをプロットすることができます
library(visreg)
visreg(mod, 'x', overlay=TRUE)
私は面白いですよbroom
とggplot2
を使ってこれを再現すると、これまでのところ私は達成していますEDこれは:
library(broom)
dt = lm(y ~ x, data = dat) %>% augment(conf.int = TRUE)
ggplot(data = dt, aes(x, y, colour = y)) +
geom_point() + geom_line(data = dt, aes(x, .fitted, colour = .fitted))
augment
funcitonはconf.int
を計算しません。どのように私はいくつかのsmooth
信頼性の反転を追加することができますか? broom
出力を使用して
geom_smooth(data=dt, aes(x, y, ymin=lcl, ymax=ucl), size = 1.5,
colour = "red", se = TRUE, stat = "smooth")
CIで回帰直線を追加するデータフレーム? –
なぜなら、私はたくさんの異なる行を追加し、より複雑なものを作成する必要があるからです。だから、「ほうき」とプロットする簡単な方法があるのだろうかと思っていました。 'geom_smooth(method =" lm ")'は、変数が多く、どの行をプロットするかを制御したいときに複雑になりますか? – giacomo