Пытаюсь получить уравнение множественной регрессии. Использую книгу-учебник Р.И.Кабакова "R в действии. Анализ и визуализация данных в программе R", последовательно выполняю примеры из главы 8. Дошел до построения диагностических диаграмм (глава 8.3.1), применяю функцию plot()
Спойлер
Код: Выделить всё
> library(car)
> setwd("D:/data/R Project/M")
> x <- read.table('_M_.dat', sep=' ', header= TRUE)
> summary(x)
long lat M X2_2
Min. :136.8 Min. :36.67 Min. :-38.24 Min. :-13.4225
1st Qu.:147.3 1st Qu.:44.89 1st Qu.:-26.09 1st Qu.: -6.4841
Median :149.5 Median :47.67 Median :-22.77 Median : -4.3408
Mean :150.2 Mean :48.05 Mean :-21.41 Mean : -4.5497
3rd Qu.:152.9 3rd Qu.:51.01 3rd Qu.:-15.81 3rd Qu.: -2.6442
Max. :164.9 Max. :59.49 Max. : 15.00 Max. : 0.2172
X3_1 N X1 X2
Min. :-9.3983 Min. : 1.0 Min. :-282.387 Min. :-1.6837
1st Qu.:-3.8409 1st Qu.: 80.5 1st Qu.: -6.557 1st Qu.: 0.7391
Median :-1.6686 Median :160.0 Median : 12.905 Median : 1.6464
Mean :-2.5937 Mean :160.0 Mean : 13.621 Mean : 1.9559
3rd Qu.:-0.8717 3rd Qu.:239.5 3rd Qu.: 37.842 3rd Qu.: 2.7114
Max. : 0.5194 Max. :319.0 Max. : 240.795 Max. : 8.2283
X3 X2_1 X4
Min. :-9115.2 Min. :-0.5961 Min. : 12.16
1st Qu.:-4655.8 1st Qu.: 0.5199 1st Qu.: 50.02
Median :-1287.5 Median : 1.3203 Median : 65.23
Mean :-2382.7 Mean : 1.9593 Mean : 73.52
3rd Qu.: -458.2 3rd Qu.: 3.0517 3rd Qu.: 80.89
Max. : 493.9 Max. : 7.3753 Max. :654.58
> fit <- lm(I(-1*M) ~ X1 + X2 + I(X3/-1000) + X4, data=x)
> summary(fit)
Call:
lm(formula = I(-1 * M) ~ X1 + X2 + I(X3/-1000) + X4, data = x)
Residuals:
Min 1Q Median 3Q Max
-25.8916 -2.9421 -0.2681 2.2475 13.8906
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 31.346898 0.799152 39.225 < 2e-16 ***
X1 -0.041058 0.004324 -9.495 < 2e-16 ***
X2 -0.817607 0.170710 -4.789 2.58e-06 ***
I(X3/-1000) -3.040876 0.168223 -18.076 < 2e-16 ***
X4 -0.007242 0.005491 -1.319 0.188
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 5 on 314 degrees of freedom
Multiple R-squared: 0.5185, Adjusted R-squared: 0.5123
F-statistic: 84.52 on 4 and 314 DF, p-value: < 2.2e-16
> confint(fit)
2.5 % 97.5 %
(Intercept) 29.77452766 32.919268162
X1 -0.04956610 -0.032549826
X2 -1.15348683 -0.481726516
I(X3/-1000) -3.37186279 -2.709889603
X4 -0.01804681 0.003562752
> plot(fit)
Ожидаю подтверждения смены страницы...
Ожидаю подтверждения смены страницы...
Ошибка: $ operator is invalid for atomic vectors
Хотелось бы узнать в чем причина и как исправить.