Recentemente, eu estava comparando dois exercícios de estatística e descobri que saídas diferentes para a mesma entrada em R talvez seja um comportamento não intencional do R, certo?
model1 <- lm(rent ~ area + bath, data = rent99)
coefficients1 <- coef(model1)
# Using a matrix without an intercept column
X <- cbind(rent99$area, rent99$bath)
model2 <- lm(rent99$rent ~ X[, 1] + X[, 2])
coefficients2 <- coef(model2)
# Both coefficients1 and coefficients2 should be identical
coefficients1
coefficients2
Saída:
(Intercept) area bath1
144.149195 4.587025 100.661413
(Intercept) X[, 1] X[, 2]
43.487782 4.587025 100.661413
Eu diria que os coeficientes são idênticos, porque os dados de entrada são idênticos