abline() 不适用于 weighted.hist()

Question

所以我使用 plotrix 库使用一些权重绘制直方图，直方图按预期显示，但是当我尝试绘制均值作为垂直线时，它根本不会显示
这是我的代码片段：

library("plotrix")
library("zoom")


vals = seq.int(from = 52.5 , to = 97.5 , by = 5)
weights <- c(18.01,18.26,16.42,14.07,11.67,9.19,6.46,3.85,1.71,0.34)/100

mean <- sum(vals*weights)
wh <- weighted.hist(x = vals , w = weights , freq = FALSE)
abline(v = mean)

abline() 似乎只适用于正常的 hist() 函数

如果这个问题听起来很愚蠢，我很抱歉，我是 R 新手，但是我进行了研究但找不到任何有用的信息。

提前致谢。

Answer 1

您应该提供您的数据样本。仅当您的权重之和为 1 时，您对加权平均值的计算才正确。否则，您应该使用 weighted.mean(vals, weights) 或 sum(vals * weights/sum(weights))。以下示例对 weighted.hist 手册页 (help(weighted.hist)) 中的示例稍作修改：

vals <- sample(1:10, 300, TRUE)
weights <- (101:400)/100
weighted.hist(vals, weights, breaks=1:10, main="Test weighted histogram")
(mean <- weighted.mean(vals, weights))
# [1] 5.246374

直方图从 1 开始，但 x 轴坐标为 0，因此我们需要减去 1 才能在正确的位置得到直线：

abline(v=mean-1, col="red")

使用您的数据，我们需要确定第一个边界以调整平均值，以便它绘制在正确的位置

wh$breaks[1]
# [1] 52.5
abline(v=mean - wh$breaks[1], col="red")

abline() 不适用于 weighted.hist()

abline() is not working with weighted.hist()

r

r-environment