x <- matrix(c(rep(1,11),rep(0,8),rep(1,3),rep(0,8),rep(1,2)),ncol = 4)
y <- matrix(c( 10,11.4,8.5,11.2,11.9,10.7,7.5,11.2),ncol = 1)
model <- lm(y~x[,c(2,3,4)])
summary(model)
上面的代码给出了第三系数的NA,是否是因为在模型中也考虑了截距,然后将第三系数写成截距和前两个矩阵的线性组合?
x <- matrix(c(rep(1,11),rep(0,8),rep(1,3),rep(0,8),rep(1,2)),ncol = 4)
y <- matrix(c( 10,11.4,8.5,11.2,11.9,10.7,7.5,11.2),ncol = 1)
model <- lm(y~x[,c(2,3,4)])
summary(model)
上面的代码给出了第三系数的NA,是否是因为在模型中也考虑了截距,然后将第三系数写成截距和前两个矩阵的线性组合?
是的,如果要获取所有三个变量的估计值,则需要省略截距。这可以通过以下方式实现: