cooper56 发表于 2017-8-13 15:25 
能解释一下吗,按我的理解您这边的残差E和F都是和B不相关的,这两个残差进行回归会是完全线性关系,没有残 ...
多元回归对系数的理解是ceteris paribus,在控制B效应之后的回归。由单变量回归建立多变量回归的时候方程两边都要partial out掉B的效应。可以用数据验证我的这两种方法得到的C的系数和残差都是一样的
- dat <- as.data.frame(a = rnorm(100),
- b = rnorm(100),
- c = rnorm(100))
- reg11 <- lm(a~b, dat)
- dat$r1 <- reg11$residuals
- reg12 <- lm(c~b, dat)
- dat$r2 <- reg12$residuals
- reg13 <- lm(r1~r2, dat)
- reg2 <- lm(a~b+c, dat)
- ## check coef
- reg13$coefficients['r2']
- reg2$coefficients['c']
- ## check residual
- resi <-data.frame(r1 = reg13$residuals,
- r2 = reg2$residuals)
- View(resi)
复制代码