Эй, ребята, пожалуйста, помогите мне здесь .. совершенно новый для R и программного обеспечения для кодирования, поэтому, пожалуйста, потерпите меня
Допустим, у меня есть 2 набора переменных (в градусах) ...
data = data.frame(a = c(-75, 7, -78, -56, -67, -82, -2, -34, -49, 83, -38),
b = c(78, -75, 69, 58, 62, 58, -71, 69, 61, -56, 70))
Вычислены средние значения и 95% -ные доверительные интервалы как для a, так и для b (данные ненормально распределен).
meana = mean(data$a)
meanb = mean(data$b)
mean.fun <- function (d, i)
{ m <- mean(d[i])
n <- length(i)
v <- (n-1)*var(d[i])/n^2
c(m, v)
}
set.seed(1111)
boota <- boot(data = data$a, mean.fun, R = 1000)
bootcia = boot.ci(boota, type = "all")
set.seed(1112)
bootb <- boot(data = data$b, mean.fun, R = 1000)
bootcib = boot.ci(bootb, type = "all")
И затем у меня есть эта тригонометрическая функция, которую я хочу построить с R, ось X будет в градусах от -90 до 90, ось Y будет ответом от функции также в -От 90 до 90 градусов.Я установил пакет REdaS, чтобы иметь функцию rad2deg для преобразования выходных данных в радианы в градусы.Я также отображаю точку данных = x = meana и y = meanb на графике.
library(REdaS)
ggplot(data.frame(x = c(-90,90)), aes(x)) +
stat_function(fun = function(x){
rad2deg(-atan(cospi((x-meana)/180)/tanpi(meanb/180)))
}) + ylim(c(-90,90)) + geom_point(x = meana, y = meanb)
С помощью кода выше я получил этот график: [1]: https://i.stack.imgur.com/kplhJ.png.
Итак, вот мой вопрос , есть ли способ добавить CI Boostrap на этот график, который у меня уже есть (как для точки данных, так и для кривой)?Я хочу видеть что-то вроде заштрихованной области круга вокруг точки данных и заштрихованной верхней / нижней области вдоль кривой как представление доверительного интервала.
Заранее спасибо.