我想知道是否可以像在eviews中一样在R中绘制实际的,拟合值和残差。我已经尝试过xyplot,但是看起来并不一样。 我用这个代码: xyplot(实际+残差+拟合值〜时间)
答案 0 :(得分:0)
我不知道eviews图表是什么样子,但是当我想查看实际,拟合和残差时,我通常使用这样的工作流程:
library(tidyverse)
data <- tibble(y = rnorm(1000),
x1 = y + rnorm(1000, 1, 10),
x2 = y + rnorm(1000, 10, 1))
model <- lm(y ~ ., data)
model %>%
predict() %>%
as_tibble() %>%
rename(.pred = value) %>%
bind_cols(data %>%
select(y)) %>%
mutate(.resid = y - .pred) %>%
ggplot(aes(x = y, y = .pred, color = .resid)) +
geom_point() +
geom_abline(slope = 1, intercept = 0) +
labs(x = "Actual",
y = "Fitted",
color = "Residuals")