pruebas que se basan en residuos

Question

Esta podría ser una pregunta r, que creo que no está permitida aquí, pero pensé que lo intentaría en caso de que esté bien. Hay una serie de pruebas que, según tengo entendido, se basan en residuos. Un ejemplo de esto es la prueba de Breusch-Pagan. Esto genera resultados aparentes a partir de eso (y la documentación parece requerir un modelo de regresión en lugar de residuos).

# load the dataset
data(mtcars)
# fit a regression model
model <- lm(mpg~disp+hp, data=mtcars)
library(lmtest)
# perform Breusch-Pagan Test
bptest(model)

¿El código está extrayendo los residuos o me equivoco al suponer que necesita residuos para esta prueba (los comentarios que leí parecen bastante claros de que se ejecutan en los residuos).

Answer 1

Para agregar al comentario de @ChristophHanck, puede verificar el código fuente de bptest y en él verá que se puede aplicar una fórmula y un marco de datos:

bptest(mpg~disp+hp, data=mtcars)

studentized Breusch-Pagan test

data:  mpg ~ disp + hp
BP = 4.0861, df = 2, p-value = 0.1296

Debajo del código, hay dos partes. Si proporciona un objeto lm, extraerá x e y, reajustará la fórmula usando lm.fit. Después de lo cual, extrae los residuos del objeto lm.fit y calcula las estadísticas relevantes. Estas son las líneas correspondientes:

resi <- lm.fit(X,y)$residuals
  sigma2 <- sum(resi^2)/n

  if(studentize)
  {
    w <- resi^2 - sigma2
    aux <- lm.fit(Z, w)    
    bp <- n * sum(aux$fitted.values^2)/sum(w^2)
    method <- "studentized Breusch-Pagan test"
  }
  else
  {
    f <- resi^2/sigma2 -1
    aux <- lm.fit(Z, f)
    bp <- 0.5 * sum(aux$fitted.values^2)
    method <- "Breusch-Pagan test"
  }

  names(bp) <- "BP"
  df <- c("df" = aux$rank - 1)
  RVAL <- list(statistic = bp,
      parameter = df,
      method = method,
      p.value= pchisq(bp, df, lower.tail = FALSE),
      data.name = dname)