Как определить, где ширина лесса в 95% CI достигает определенной величины? - PullRequest
Новая
       50

Как определить, где ширина лесса в 95% CI достигает определенной величины?

0 голосов
/ 08 января 2020

Я установил кривую лёсса с 95% -ным доверительным интервалом, и теперь я хотел бы иметь возможность определить, где КИ достигает определенной ширины.

Например, используя набор данных "cars" : 

plot <- ggplot (cars, aes (x=speed, y=dist)) +
geom_point() +
stat_smooth (method= "loess", se=TRUE) +
xlab("Speed")+
ylab("Distance")+
theme_bw()
plot

plot

Мне бы хотелось узнать, при каких значениях "скорости" CI равны 20 единицам расстояние. Глядя на сюжет, это может быть примерно 7 и 24.

Спасибо!

1 Ответ

1 голос
/ 08 января 2020

Вы можете использовать ggplot_build(plot) для извлечения релевантных данных о слоях, встроенных в ggplot2, и другой неверной c информации.

В этом случае пределы доверительных интервалов находятся в ymin и ymax столбцы и могут быть достигнуты с помощью: 

foo <- ggplot_build(plot)
foo[["data"]][[2]]

Затем вы можете сделать простой мутирование, чтобы изучить различия между ymax и ymin и с какой «скоростью» разрыв CI достигает 20 через столбец x. 

mutate_info <- foo[["data"]][[2]] %>% dplyr::mutate(ci_gap = ymax-ymin)
...