R Durbin Watson测试lm对象列表

Question

我有一个包含两个（或更多）lm个对象的列表。现在，我想同时在dwtest个durbinWatsonTest对象上分别lmtest或car执行Durbin-Watson测试，即。我想做

lm

当我这样做时，我会得到以下信息：使用dwtest(reg_objects) durbinWatsonTest(reg_objects) ：值为dwtest。使用NA：

durbinWatsonTest

为了促进这个想法，我已经生成了一个5 lm对象的列表对象：

Error in durbinWatsonTest.default(reg_1y$call) : 
requires vector of residuals 

你能帮我解决这个问题吗？ 安德烈亚斯

修改：# non-sensical lm! set.seed(123) y <- 1 + rnorm(10) # this makes a list of 5 lm's using the response data 'y' data <- replicate( 5 , lm( y ~ 1 + rnorm(10) ) , simpl = FALSE ) ，bgtest和bptest的同一问题。注意：sctest来自sctest。

Answer 1

试试这个：

my_dw_stat <- list()

for(dw in seq(data)) {

  my_dw_stat[[dw]] <- dwtest(data[[dw]]$model[, 1] ~ data[[dw]]$model[, 2])

}