检查R中线性模型中的观察子集

Question

在r中，我有一个名为lmod452的线性模型。从这个线性模型，我想检查案例编号25-50的观察，以观察哪些是一个异常。如何在选定数量的观察中使用绘图功能？

谢谢！

Answer 1

您可以使用[]来分组数据。所以对于假数据

x <- rnorm(100)
y <- x + 0.5*rnorm(100)
plot(x,y)

...你可以绘制观察的子集......

plot(x[25:50], y[25:50])

...或使用不同的绘图角色或颜色覆盖它们......

plot(x,y)
points(x[25:50], y[25:50], pch=16, col="red")

很多方法，希望这能让你开始。

哦，既然你提到你的lm模型而不是你的数据，那么值得注意的是你可以从lm对象中提取数据。只需致电names(lmod452$model)即可获取姓名。

lmod452 <- lm(y~x)

# extracting stuff
xxx <- lmod452$model$x
yyy <- lmod452$model$y