像这样的东西?

或者你是在中位数的某个间隔之后,就像你得到带有缺口的箱线图(但自然适合一个样本比较)?
这是一个例子:

这使用了 McGill 等人建议的区间(参考文献中的区间?boxplot.stats
)。实际上可以使用缺口,但这可能会增加将其解释为普通缺口箱线图的机会。
当然,如果您需要一些东西来更直接地复制有符号秩检验,可以构建各种东西来做到这一点,甚至可以包括伪中位数的间隔(即单样本 Hodges-Lehmann 位置估计,中位数成对平均值)。
确实,wilcox.test
可以为我们生成必要的信息,所以这很简单:
> wilcox.test(pd,mu=1.1,conf.int=TRUE)
Wilcoxon signed rank test
data: pd
V = 72, p-value = 0.5245
alternative hypothesis: true location is not equal to 1.1
95 percent confidence interval:
0.94 1.42
sample estimates:
(pseudo)median
1.1775
这也可以绘制:

[箱线图间隔较宽的原因是,当数据相当正态时,正态中位数的标准误差(这是基于 IQR 计算的假设)往往大于伪中位数的标准误差。]
当然,人们可能希望将实际数据添加到绘图中:

Z值
R 使用正等级的总和作为其测试统计量(这与测试的维基百科页面上讨论的统计量不同)。
Hollander 和 Wolfe 给出统计量的平均值为,方差为。n(n+1)/4n(n+1)(2n+1)/24
因此,对于您的数据,平均值为 60,标准差为 17.61,z 值为 0.682(忽略连续性校正)
我用来生成第四个图的代码(也可以通过省略不需要的部分来完成前面的图)有点粗糙(它主要针对问题,而不是一般的绘图函数),但我想有人可能想要它:
notch1len <- function(x) {
stats <- stats::fivenum(x, na.rm = TRUE)
iqr <- diff(stats[c(2, 4)])
(1.96*1.253/1.35)*(iqr/sqrt(sum(!is.na(x))))
}
w <- notch1len(pd)
m <- median(pd)
boxplot(pd,horizontal=TRUE,boxwex=.4)
abline(v=1.1,col=8)
points(c(m-w,m+w),c(1,1),col=2,lwd=6,pch="|")
ci=wilcox.test(pd,mu=1.1,conf.int=TRUE)$conf.int #$
est=wilcox.test(pd,mu=1.1,conf.int=TRUE)$estimate
stripchart(pd,pch=16,add=TRUE,at=0.7,cex=.7,method="jitter",col=8)
points(c(ci,est),c(0.7,0.7,0.7),pch="|",col=4,cex=c(.9,.9,1.5))
lines(ci,c(0.7,0.7),col=4)
稍后我可能会回来发布更多功能代码。