我无法为您提供有效的 R 代码,因为您没有提供原始数据(箱线图需要这些数据),并且不清楚您想要显示什么,因为没有任何迹象表明您的黄金标准在给定的聚合中发挥作用数据(这些是用不同仪器重复测量吗?),除非报告的平均值代表第 i 个方法和参考方法之间的差异(在这种情况下,我看不出你如何使用箱线图)。您的数据的基本图可能看起来像
dfrm <- data.frame(method=LETTERS[1:3], lcl=c(-5,-9,-8),
mean=c(4,2,4), ucl=c(15,13,16), var=c(27,33,36))
# I use stripchart to avoid axis relabeling and casting of factor to numeric
# with default plot function
stripchart(mean ~ seq(1,3), data=dfrm, vertical=TRUE, ylim=c(-10,20),
group.names=levels(dfrm$method), pch=19)
with(dfrm, arrows(1:3, mean-lcl, 1:3, mean+lcl, angle=90, code=3, length=.1))
abline(h=0, lty=2)
但是,我可以建议您查看MethComp包,它可以帮助您将几种方法与黄金标准进行比较,无论是否重复,以及显示结果。配套教材是
Carstensen, B.比较临床测量方法。约翰威利父子公司 2010
