Question

大家好，请在这里帮助我。.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))

计算a和b的平均值和95％的自举置信区间（数据不是正态分布的）。

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 =平均值和y =平均值。

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。

这是我的问题 ，有没有办法在我已经拥有的绘图上添加boostrap CI（数据点和曲线）？我想看到类似数据点周围的阴影圆区域和沿曲线的阴影上/下区域，以表示置信区间。

谢谢。

用R度的置信区间绘制三角函数

0 个答案: