Два способа справиться с этим. Первый - использовать грубую силу и вычислять прогнозы с newdata=<subset of the two obs>
в качестве аргумента predict.lm
, а второй - проверять значения $residuals
в Reg
-объекте. Если первые два наблюдения NLSY_education[1:2,]
, то это может быть так просто:
#See ?lm, Values section
Reg$residuals[1:2]