在R中找到y-hat

时间:2015-02-24 12:43:34

标签: r

如果给出如下虚拟数据,如何为每次观察创建和打印带y,X,y-hat和e的data.frame:

y x
17 1
22 2
29 3
29 4
38 5
39 6
45 7

1 个答案:

答案 0 :(得分:1)

假设您需要拟合值和简单线性回归模型的残差,您可以按如下方式获取:

mod <- lm(y~x, data = df)
data.frame(df, y_hat = fitted(mod), e = residuals(mod))
   y x    y_hat          e
1 17 1 17.67857 -0.6785714
2 22 2 22.21429 -0.2142857
3 29 3 26.75000  2.2500000
4 29 4 31.28571 -2.2857143
5 38 5 35.82143  2.1785714
6 39 6 40.35714 -1.3571429
7 45 7 44.89286  0.1071429