提前道歉 - 我是 R 的新手,但我已经用谷歌搜索了好几个星期,但我还没有找到一个直接的答案。
问题
基本上,我希望使用 R 环境在盒子和胡须图上绘制比例及其置信区间。
即绘制 binom.test(100,1000) 和 binom.test(125,1000) 的结果
什么是最好的包(我没有任何偏好)和这样做的方法是什么?
提前道歉 - 我是 R 的新手,但我已经用谷歌搜索了好几个星期,但我还没有找到一个直接的答案。
问题
基本上,我希望使用 R 环境在盒子和胡须图上绘制比例及其置信区间。
即绘制 binom.test(100,1000) 和 binom.test(125,1000) 的结果
什么是最好的包(我没有任何偏好)和这样做的方法是什么?
这确实看起来更适合 SO 中的 [r] 标记问题,但似乎对 CV 中的此类问题有惊人的广泛容忍度,所以就这样吧。答案还假设您并不想使用现有boxplot()函数,该函数将用于构建真正的盒须图(并且它并不真正使用二项式 CI,而是使用 Tukey 定义的fivenum()函数)。plotrix 和 gplots 包都具有相同的 plotCI 函数;
est1 <- binom.test(100,1000) ; est2 <- binom.test(125,1000)
estvals <- rbind( c(est1$estimate, est1$conf.int[1], est1$conf.int[2]),
c(est2$estimate, est2$conf.int[1], est2$conf.int[2]))
plotCI(x = c(1,2), estvals[,1], ui=estvals[,3], li=estvals[,2],
xlim=c(0.5,2.5), ylim=c(0, 0.15) , xaxt="n", xlab="Groups")
基本函数是带有参数的arrows()绘图函数,该参数将箭头的点在两个方向上展开 90 度。