基本 R 图中的 lines() 函数给出了几条线而不是平滑线

Question

我正在尝试为我拥有的回归模型图添加平滑器。 我只是使用基数 R 来绘制我的 X 和 Y 向量，然后使用 plot() 和lines() 添加一个平滑器。 我以前做过这个，它奏效了，但今天我得到了一个多条线连接点的图，而不是一条穿过所有数据的平滑线。 我无法弄清楚我编写的这段代码有什么不同，所以我希望这里有人可以帮助我确定问题。

这是我的代码。 我正在使用我随机生成的数据来练习其他东西：

X 和随机变量向量来创建“Y”：

X <- rnorm(100, mean = 10, sd = 1)
epsilon <- rnorm(100, 0, 1)

是：

b0 <- 0.27
b1 <- 0.49
b2 <- 0.62
b3 <- 0.8

Y <- b0 + b1*X + b2*2^2 + b3*X^3 + epsilon

创建 df 和 reg 模型/Yhat：

df = data.frame(Y,X,epsilon)
reg <- lm(Y ~ I(X^3), data = df)
Yhat <- fitted.values(reg)
cbind(df, Yhat) -> df

阴谋：

plot(X, Y)
lines(X, Yhat, col = "blue", lwd = 0.5)

Answer 1

为此，必须对X值进行排序，并根据对应的X值对Y值进行排序：

X <- rnorm(100, mean = 10, sd = 1)
epsilon <- rnorm(100, 0, 1)

b0 <- 0.27
b1 <- 0.49
b2 <- 0.62
b3 <- 0.8

Y <- b0 + b1*X + b2*2^2 + b3*X^3 + epsilon

df = data.frame(Y,X,epsilon)
reg <- lm(Y ~ I(X^3), data = df)
Yhat <- fitted.values(reg)
cbind(df, Yhat) -> df

plot(X, Y)
lines(X[order(X)], Yhat[order(X)], col = "blue", lwd = 0.5)