我有一个线性模型,我使用函数plot(model.lm$residuals)
绘制残差。但是,我的数据集由两大组观察(领土)组成,我想按组识别残差。我的数据集中的二进制变量已经通过向一个组提供“1”而向另一个组提供“0”来识别两个组。我怎么能得到R来绘制我的残差并给与我的1和0相关的残差提供不同的颜色?
谢谢!
答案 0 :(得分:4)
data(iris)
iris$isVersicolor <- factor(as.numeric(iris$Species == "versicolor"))
fit <- lm(Sepal.Length ~ Sepal.Width + Petal.Length, data = iris)
plot(residuals(fit), col = iris$isVersicolor)