如何显示线性回归的去趋势数据?

机器算法验证 数据可视化
2022-03-27 15:03:25

我只是在看一下哈德利关于他的 ggplot2 R 包的非常优秀的书。

他有一些代码可以删除 diamonds 数据集中的线性趋势,如下所示:

d <- subset(diamonds, carat < 2.5 & rbinom(nrow(diamonds), 1, 0.2) == 1)
d$lcarat <- log10(d$carat)
d$lprice <- log10(d$price)

detrend <- lm(lprice ~ lcarat, data = d)
d$lprice2 <- resid(detrend)

qplot(lcarat, lprice, data = d)
qplot(lcarat, lprice2, data = d)

生成这些图表

未调整...

在此处输入图像描述

去势...

在此处输入图像描述

我想看看没有 lcarat 影响的 lprice 的实际值是多少。绘制残差与 lcarat 显示正确的形状,但点向 y = 0 移动(查看 y 轴单位的范围)。

为了得到我想要的,简单地绘制残差+均值(lprice)是否有意义?即按均值(lprice)向上移动上一个图。

qplot(lcarat, lprice2 + mean(lprice), data = d)

在此处输入图像描述

这样做有意义吗?我想做的事情有名字吗?

1个回答

对我来说,这非常令人困惑,尤其是当你可以做更简单的事情时——计算price/carat得到一克拉的价格,这将更容易解释。