R中散点图标准差的计算

Question

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

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)))

除非你有充分的理由不这样做，否则你应该考虑x上y的线性模型得到任何有意义的结果。如果x和y之间的关系不是1对1，那么假设正确的模型是没有意义的。如果y上的x关系不是1对1，则y-x将不会正常分布，因此sd将难以以有意义的方式解释。