R回归三倍残差以外观测值提取的问题

举个例子：resg <- lm(y ~ x1+x2, data=demo)，用这个回归模型可以得到，
residual，standard，student，hat_matrix，DFFITS，cooks_distance，COVRATIO等其它统计量，可以显示统计量
的同时把对应的x1，x2值显示出来吗？？想把三倍残差以外的观测值提取后，再拿数据来做回归，这样R平方判定系数要大些。
R的观测值提取怎么就那么难呢？？，SAS都可以把统计量输出的同时也把原始变量也显示出来，原始变量显示不出来，怎么进一步提取数据啊。求高手解答，谢谢！

对了需要说明的是，不用merge的方法，可以回归出统计量的同时显示出原始的那些自变量（预测变量）吗？因为数据量大，又要分组回归，所以merge的方法数据量过大会出错。