利用R语言计算实验数据的重复性并绘制折线图

最新推荐文章于 2025-12-05 18:19:08 发布

原创最新推荐文章于 2025-12-05 18:19:08 发布 · 925 阅读

8 ·

CC 4.0 BY-SA版权

文章标签：

#r语言 #开发语言

文章详细描述了对一组设备在四个测试点进行的三次实验，通过计算各测试点的误差值平均值和标准差，探讨设备的精度及重复性。结果以误差曲线图的形式展示并可视化.

现做了一组设备的准确度测试实验，分别在不同的4个测试点进行3次实验，每个测试点的实验最终结果采取平均值的方式体现，随后计算其重复性，重复性计算公式为：

$s = \sqrt{\frac{\sum_{i=1}^{n}{({x}_i-\bar{x})}}{n-1}}$

即第i次实验点的第j次测得误差值（%）-第i次实验j次测试测得数据平均值的集合除以测试次数-1，所得值再开根号。平均值计算公式为：

TestPoint1 <- c(0.22, 0.25, 0.31)
TestPoint2 <- c(0.78, 0.81, 0.85)
TestPoint3 <- c(0.15, 0.13, 0.22)
TestPoint4 <- c(-1.96, -1.78, -1.43)
mean_value1 <- mean(TestPoint1)
mean_value2 <- mean(TestPoint2)
mean_value3 <- mean(TestPoint3)
mean_value4 <- mean(TestPoint4)
std_dev1 <- sd(TestPoint1)
std_dev2 <- sd(TestPoint2)
std_dev3 <- sd(TestPoint3)
std_dev4 <- sd(TestPoint4)
x <- c(100, 40, 20, 10)
y <- c(mean_value1, mean_value2, mean_value3, mean_value4)
y_sd <- c(std_dev1, std_dev2, std_dev3, std_dev4)
plot(x, y, type = "o", ylim = c(-2,1), ylab = "误差值", xlab = "测试点", main = "误差曲线图")
grid(nx = NA, ny = NULL, col = "gray", lwd = 0.2)
text(x = x, y = y, labels = y, pos = 3)
for(i in 1:length(x)) {
  arrows(x[i], y[i] - y_sd[i], x[i], y[i] + y_sd[i], code=3, length = 0)
}