6

我需要使用 ggplot 创建一个多变量回归线。

我的数据：

dput(head(x2,15))
structure(list(Date = structure(c(15608, 15609, 15610, 15611, 
15612, 15613, 15614, 15615, 15616, 15617, 15618, 15619, 15620, 
15621, 15622), class = "Date"), Cpu = c(77.0763, 51.8909, 59.3229, 
89.5822, 87.7448, 80.4413, 57.5009, 99.8185, 99.9969, 91.5528, 
50.0793, 56.4049, 57.808, 51.0453, 56.0505), Memory = c(369.667979452055, 
341.572253381722, 345.013066490241, 334.520135424091, 374.107056613899, 
1592.38342810723, 470.204599904169, 393.802909594735, 540.817571059432, 
425.49563812601, 438.326775174387, 614.417456359102, 1255.63550519358, 
466.993243243243, 358.445879354291), Response = c(52.25, 48.36, 
49.23, 50.99, 48.63, 46.11, 43.03, 45.35, 50.03, 46.18, 47.39, 
43.28, 55.36, 50.59, 50.44)), .Names = c("Date", "Cpu", "Memory", 
"Response"), row.names = c(1L, 4L, 6L, 7L, 9L, 10L, 13L, 16L, 
19L, 25L, 29L, 32L, 35L, 39L, 42L), class = "data.frame")

我可以在Response和之间做到这一点Cpu：

ggplot(x2, aes(Response)) + 
   geom_point(aes(y = Memory), size = 2, colour = "blue") +
   geom_point(aes(y = Cpu), size = 2, colour = "orange") + 
   geom_smooth(method = "lm", formula = "Response ~ CPU+Memory", 
               size = 1.5, colour = "red", se = T)

我收到此错误：

Error: stat_smooth requires the following missing aesthetics: y

有任何想法吗？

4

1 回答 1

9

首先，stat_smooth取公式；我不认为geom_smooth需要一个。其次，我不认为你可以输入一个超过一个predictorin的公式stat_smooth。如果我错了，请纠正我。另一种解决方案是自己拟合模型并计算预测值，并以这种方式自己计算和绘制 SE：（取自http://docs.ggplot2.org/0.9.3/geom_smooth.html）

model <- lm(data = df, Response ~ Memory + Cpu)
df$model <- stats::predict(model, newdata=df)
err <- stats::predict(model, newdata=df, se = TRUE)
df$ucl <- err$fit + 1.96 * err$se.fit
df$lcl <- err$fit - 1.96 * err$se.fit

g <- ggplot(df)
g <- g + geom_point(aes(x=Response, y = model), size = 2, colour = "blue")
g <- g + geom_smooth(data=df, aes(x=Response, y=model, ymin=lcl, ymax=ucl), size = 1.5, 
            colour = "red", se = TRUE, stat = "smooth")

这给出了如下所示的结果（不确定它是您期望的）： 在此处输入图像描述

于 2013-01-03T20:46:47.790 回答