Для простого примера ниже вы можете видеть, что есть определенные точки, которые идентифицируются на последующих графиках. Как я могу извлечь номера строк, указанные на этих графиках, особенно график Normal Q-Q?
set.seed(2016)
maya <- data.frame(rnorm(100))
names(maya)[1] <- "a"
maya$b <- rnorm(100)
mara <- lm(b~a, data=maya)
plot(mara)
Я попытался использовать str (мара), чтобы увидеть, могу ли я найти список там, но я не вижу ни одного из чисел из графика Normal Q-Q. Мысли?
Я редактировал ваш вопрос, используя set.seed(2016) для воспроизводимости. Чтобы ответить на ваш вопрос, мне нужно объяснить, как создать график Q-Q, который вы видите.
se <- sqrt(sum(mara$residuals^2) / mara$df.residual) ## Pearson residual standard error
hii <- lm.influence(mara, do.coef = FALSE)$hat ## leverage
std.resi <- mara$residuals / (se * sqrt(1 - hii)) ## standardized residuals
## these three lines can be replaced by: std.resi <- rstandard(mara)
Теперь сравним график Q-Q, который мы сгенерируем и сгенерируем с помощью plot.lm:
par(mfrow = c(1,2))
qqnorm(std.resi, main = "my Q-Q"); qqline(std.resi, lty = 2)
plot(mara, which = 2) ## only display Q-Q plot
То же самое, правильно?
Теперь единственная проблема заключается в том, как обозначаются цифры. Эти помеченные точки отмечают наибольшие 3 абсолютные стандартизированные остатки. Рассмотрим:
x <- sort(abs(std.resi), decreasing = TRUE)
id <- as.integer(names(x))
id[1:3]
# [1] 23 8 12
Теперь, если вы внимательно посмотрите на график, вы увидите, что эти три числа - это именно то, что показано. Зная это, вы также можете проверить, например, id[1:5].