R中散点图标准差的计算

Question

我使用 R 与代表 x=y 对角线的线（使用 abline）结合创建了两个向量的散点图。我希望计算点与对角线的标准偏差，并对第一和第三分位数之间找到的区域进行着色。
我不知道该怎么做，非常感谢所有帮助！ 提前致谢。 朝觐。

Answer 1

好吧，你想要做的是这样的：

# sample data
x <- rnorm(50,0,2)
y <- x+rnorm(50,0,2)

# construct polygons
div <- quantile(y-x,c(0.25,0.75))
x1 <- min(c(x,y))
x2 <- max(c(x,y))


plot(x,y,type="n")
polygon(x=c(x1,x1,x2,x2),y=c(x1+div,(x2+div)[c(2,1)]),col="grey")
abline(0,1)
points(x,y)

我会做的是这样的：

qplot(x,y,geom="point") + stat_smooth(method="lm")

你想要计算的标准偏差是

sd(y-x)

你可能正在寻找的正确度量是：

sd(residuals(lm(y~x)))

你应该根据 y 的线性模型来思考x 以获得任何有意义的结果，除非您有充分的理由不这样做。如果 x 和 y 之间的关系不是 1 on 1，那么假设正确的模型是没有意义的。如果 x 与 y 的关系不是 1 on 1，则 yx 将不是正态分布，因此 sd 将难以以有意义的方式解释。