Как я могу сделать тест Левена, используя функцию lapply? - PullRequest
1 голос
/ 15 апреля 2020

Как я могу пройти тест Левена (автомобильный пакет) моих регрессий? Для регрессий я использовал функцию lapply.

Вот мой код

library(data.table)

Regressions<- 
data.table(February)[, 
                  .(Lm = lapply(.SD, function(x) summary(lm(February$PPNA ~ February$Acum1 + x)))),
                  .SDcols = 80:157]

1 Ответ

1 голос
/ 15 апреля 2020

Если у нас есть столбец группировки ('group'), то мы l oop по интересующим столбцам из .SDcols, применяем levene.test с этим значением столбца и группой

library(data.table)
library(car)
as.data.table(February)[, 
              .(Lm = lapply(.SD, function(x) list(levene.test(x, group)))),
              .SDcols = 80:157]
...