R: Идентичные значения SE для всех склонов в lm

Question

Я относительно новичок в R, и я просто столкнулся с проблемой, которой я не был раньше. Я использовал lm с видовым богатством в качестве переменной ответа и уровня высоты (у меня пять из них в диапазоне от 3000 до 5000 м, это категориальная переменная) в качестве предиктора. Выход модели дает мне точные значения для стандартных ошибок для всех склонов (0.7709 для перехвата и 1.0903 для остальных 4-х уровней). Кажется, все в порядке ... Я не могу сказать, что может быть проблемой. Кто-нибудь из вас столкнулся с подобной проблемой? Я был бы очень признателен за предложения о том, где я должен искать потенциальные проблемы.

Answer 1

Весьма легко получить данные, обладающие этим свойством. Например:

set.seed(1) 
y <- rnorm(1e3) 
x <- rep(1:4, length.out = 1e3) 
y <- y + x 
summary(lm(y ~ factor(x))) 

Call: 
lm(formula = y ~ factor(x)) 

Residuals: 
    Min  1Q Median  3Q  Max 
-3.0078 -0.6817 -0.0251 0.6980 3.8188 

Coefficients: 
      Estimate Std. Error t value Pr(>|t|)  
(Intercept) 0.95308 0.06554 14.54 <2e-16 *** 
factor(x)2 1.04662 0.09269 11.29 <2e-16 *** 
factor(x)3 2.03841 0.09269 21.99 <2e-16 *** 
factor(x)4 3.05606 0.09269 32.97 <2e-16 *** 
--- 
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

Residual standard error: 1.036 on 996 degrees of freedom 
Multiple R-squared: 0.5468, Adjusted R-squared: 0.5454 
F-statistic: 400.6 on 3 and 996 DF, p-value: < 2.2e-16 

Возможно, что-то само по себе неправильно с вашим кодом или обязательно вашими данными. Вероятно, это, скорее всего, статистический вопрос, а не вопрос R.