lapply为R中列表比较

Posted 2023-01-09

我的前3行中低于工作微小的R代码。但是，当我做这3个行成一个功能，我得到的Error in m[[i + 1]] : subscript out of bounds。

这个问题能解决吗？

m1 <- lm(hp ~ vs, mtcars)
m2 <- lm(hp ~ vs*wt, mtcars)

pchisq(2 * (logLik(m2) - logLik(m1)), df = abs(m1$df.residual - m2$df.residual), lower = F)

###### Make a function of above 3 lines: ######

 compare <- function(...)

   m <- list(...)
   L <- length(m)

lapply(1:L, function(i) pchisq(2 * (logLik(m[[i+1]]) - logLik(m[[i]])), df = abs(m[[i]]$df.residual - m[[i+1]]$df.residual), lower = F) )
 
 ## Example of use:
  compare(m1, m2)       # `Error in m[[i + 1]] : subscript out of bounds`

答案

你lapply语句是从第一个索引（1）m的最后一个索引去（L，在这种情况下）。所以，当你尝试调用m[[i+1]]，你调用一个指标，是不是m的一部分。如果m具有长度为5，那么你的呼叫m[[5+1]]，这是出界。